Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,070,410	C→A	9.1%	intergenic (‑9/+105)	yihR ← / ← yihS	putative sulphoquinovose mutarotase/sulphoquinovose isomerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,070,410	0	C	A	9.1%	297.0 / 17.5	88	intergenic (‑9/+105)	yihR/yihS	putative sulphoquinovose mutarotase/sulphoquinovose isomerase
Reads supporting (aligned to +/- strand):  ref base C (42/38);  new base A (4/4);  total (46/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.40e-01

GATATTCCTGCCCCTGATAGCGGTAACAGCCGTTGGCGATCCGGTTCGGCCACGGAATAAGCACTTTACCCAGATGCGCCAGCGGCATCTCTTCGGGTTTGTGCGGGATGACTAAATGACAACCCTGAAAGGTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACG‑ACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACGGGATACGCGGGATCACCAGGCAATGCAGCAGGTGATAAATATCCTGTTTGCCGTCCCAGACTTTGGTGGTGACCTTATTGTCCGCATCCAGCTCCTGCC  >  NC_000913/4070143‑4070675                                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                                           gATATTCCGGCCCCTGATAGCGGTAACAGCCGTTGGCGATCCGGTTCGGCCACGGAATAAGCACTTTACCCAGATGCGCCAGCGGCATCTCTTCGGGTTTGTGCGGGATGACTAAATGACAACCCTGAAAGGTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACG‑ACATACAGCCTCCCGTGAgg                                                                                                                                                                                                                                                                     <  4:334757/275‑1 (MQ=255)  atatTCCTGCCCCTGATAGCGGTAACAGCCGTTGGCGATCCGGTTTGGCCACGGAATCAGCACTTTCCCCAGATGCGCCAGCCGCATCTCTTCGGGTTTGTGCGGCATGTCTACATGACAACCCTGAATGGTTAATTCAGCCCGACCACCACCGACGGTCACAATGGTCGCGTGATACTCCCCCGCGGCAACTGAAACGGTCTGTCCACGTCAACTCATATTGGTGATTTGCAGATAGGTCACCTTAATTAACC‑ACCTACCCCCTcccagcagca                                                                                                                                                                                                                                                                    >  4:110238/1‑268 (MQ=255)          gCCCCTGATAGCGGTAACAGCCGTTGGCGATCCGGTTCGGCCACGGAATAAGCACTTTACCCAGATGCGCCAGCGGCATCTCTTCGGGTTTGTGCGGGATGACTAAATGACAACCCTGAAAGGTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCGTAATTCACG‑ACATAAAGCCTCCCGTGAGGCGG‑ATGTc                                                                                                                                                                                                                                                            >  2:421451/1‑274 (MQ=255)                                                ggCCACGGAATAAGCACTTTACCCAGATGCGCCAGCGGCATCTCTTCGGGTTTGTGCGGGATGACTAAATGACAACCCTGAAAGGTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACG‑ACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCtt                                                                                                                                                                                                                      <  2:151523/274‑1 (MQ=255)                                                    aCGGAATAAGCACTTTACCCAGATGCGCCAGCGGCATCTCTTCGGGTTTGTGCGGGATGACTAAATGACAACCCTGAAAGGTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACG‑ACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTggg                                                                                                                                                                                                                  <  4:53692/274‑1 (MQ=255)                                                     cGGAATAAGCACTTTACCCAGATCCCCCAGCGGCATCTCTTCGGGTTTGTCCGGGATTACTAAATGACCCCCCTGAAAGGTTAATTCAGCCAGCCCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGAACCTTAATTAACG‑ACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGga                                                                                                                                                                                                                <  1:75077/275‑1 (MQ=255)                                                      ggAATAAGCACTTTACCCAGATGCGCCAGCGGCATCTCTTCGGGTTTGTGCGGGATGACTAAATGACAACCCTGAAAGGTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACG‑ACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGgag                                                                                                                                                                                                               <  2:131988/275‑1 (MQ=255)                                                                  ttACCCAGATGCGCCAGCGGCATCTCTTCGGGTTTGTGCGGGATGACTAAATGACAACCCTGAAAGGTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGCTTTGCATACAAGGCACCTTAATTAACG‑ACATACAGCCTCCCGTGAGGCTT‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGCGAGAGGGTTAGGga                                                                                                                                                                                                    >  2:72682/1‑274 (MQ=255)                                                                             cgcCAGCGGCATCTCTTCGGGTTTGTGCGGGATGACTAAATGACAACCCTGAAAGGTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACG‑ACATACAGCCTCCCGTGAGGCTG‑ATGTCaa                                                                                                                                                                                                                                                          >  1:154838/1‑209 (MQ=255)                                                                             cgcCAGCGGCATCTCTTCGGGTTTGTGCGGGATGACTAAATGACAACCCTGAAAGGTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACG‑ACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTg                                                                                                                                                                                                                    >  4:237616/1‑247 (MQ=255)                                                                                 aGCGGCATCTCTTCGGGTTTGTGCGGGATGACTAAATGACAACCCTGAAAGGTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACG‑ACATACAGCCTCCCGTGAGGCTG‑ATGTCAAAccc                                                                                                                                                                                                                                                      <  3:154838/209‑1 (MQ=255)                                                                                 aGCGGCATCTCTTCGGGTTTGTGCGGGATGACTAAATGACAACCCTGAAAGGTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACG‑ACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGga                                                                                                                                                                                                                <  2:237617/247‑1 (MQ=255)                                                                                                    tgtgCGGGATGACTAAATGACAACCCTGAAAGGTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACG‑ACATACAGCCTCCCGTGAGGCTG‑At                                                                                                                                                                                                                                                               >  4:538582/1‑181 (MQ=255)                                                                                                      tgCGGGATGACTAAATGACAACCCTGAAAGGTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACG‑ACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAACACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTa                                                                                                                                                                >  2:524348/1‑274 (MQ=255)                                                                                                        cGGGATGACTAAATGACAACCCTGAAAGGTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACG‑ACATACAGCCTCCCGTGAGGCTG‑ATGTCa                                                                                                                                                                                                                                                           <  2:538583/181‑1 (MQ=255)                                                                                                           gATGACTAAATGACAACCCTGAAAGGTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACG‑ACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGc                                                                                                                                                          >  2:87625/1‑275 (MQ=255)                                                                                                                         aaCCCTGAAAGGTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACG‑ACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTg                                                                                                                                                                                                                    >  4:429662/1‑203 (MQ=255)                                                                                                                           cccTGAAAGGTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACG‑ACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCACTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGATAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCcagcag                                                                                                                                           >  2:445036/1‑274 (MQ=255)                                                                                                                             cTGAAAGGTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACG‑ACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGga                                                                                                                                                                                                                <  2:429663/203‑1 (MQ=255)                                                                                                                              tGAAAGGTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACG‑ACATACAGCCTCCCGTGa                                                                                                                                                                                                                                                                       <  5:1076039/148‑1 (MQ=255)                                                                                                                                     ttAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACG‑ACATACAGCCTCCCGTGAGGc                                                                                                                                                                                                                                                                    <  8:110245/144‑1 (MQ=255)                                                                                                                                     ttAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACG‑ACATACAGCCTCCCGTGAGGc                                                                                                                                                                                                                                                                    >  7:110245/1‑144 (MQ=255)                                                                                                                                      tAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACG‑ACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAAc                                                                                                                                                                                         >  2:314276/1‑217 (MQ=255)                                                                                                                                        aTTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACG‑ACATACAGCCTCCCGTGAg                                                                                                                                                                                                                                                                      >  6:996749/1‑139 (MQ=255)                                                                                                                                        aTTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACG‑ACATACAGCCTCCCGTGAg                                                                                                                                                                                                                                                                      <  5:996749/139‑1 (MQ=255)                                                                                                                                              ccAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCCTATTGGTGATTTGCATACAAGGCACCTTAATTAACG‑ACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGTTCCCCTCGCCCCTTTGGGGACAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCCGATCAGTCATTTCACATTAATATCCACCAGACCCGCTGCAACCCCTGGAg                                                                                                                       >  1:147008/1‑275 (MQ=255)                                                                                                                                                    cAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACG‑ACAT‑CAGCCTCACGGGAGGCTgtatg                                                                                                                                                                                                                                                              >  5:490588/1‑130 (MQ=255)                                                                                                                                                    cAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACG‑ACAT‑CAGCCTCACGGGAGGCTgtatg                                                                                                                                                                                                                                                              <  6:490588/134‑5 (MQ=255)                                                                                                                                                     aGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACG‑ACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACt                                                                                                                                                                                                                                               >  8:44140/1‑148 (MQ=255)                                                                                                                                                            aCGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACG‑ACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGCGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATTCAGAAGCACCGATGCACCCGCTGGCAGCATTCCTGgcgcc                                                                                                         >  2:537731/1‑274 (MQ=255)                                                                                                                                                                  acaATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACG‑ACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTc                                                                                                                                                                                                                                  >  8:96848/1‑148 (MQ=255)                                                                                                                                                                   caATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACG‑ACATACAGCCTCCCGTGAGGCTG‑ATGTCaaa                                                                                                                                                                                                                                                         <  5:784275/124‑1 (MQ=255)                                                                                                                                                                   caATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACG‑ACATACAGCCTCCCGTGAGGCTG‑ATGTCaaa                                                                                                                                                                                                                                                         >  6:784275/1‑124 (MQ=255)                                                                                                                                                                     aTGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACG‑ACATACAGCCTccc                                                                                                                                                                                                                                                                           <  8:894148/105‑1 (MQ=255)                                                                                                                                                                     aTGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACG‑ACATACAGCCTccc                                                                                                                                                                                                                                                                           >  7:894148/1‑105 (MQ=255)                                                                                                                                                                            cgTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACG‑ACATACAGCCTCCCGTGAGGCTG‑ATGTCaa                                                                                                                                                                                                                                                          >  6:1126903/1‑114 (MQ=255)                                                                                                                                                                            cgTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACG‑ACATACAGCCTCCCGTGAGGCTG‑ATGTCaa                                                                                                                                                                                                                                                          <  5:1126903/114‑1 (MQ=255)                                                                                                                                                                             gTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACG‑ACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCt                                                                                                                                                                                                                       <  7:44140/148‑1 (MQ=255)                                                                                                                                                                              ggATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACG‑ACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCtt                                                                                                                                                                                                                      <  7:96848/147‑1 (MQ=255)                                                                                                                                                                                 tAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACG‑ACATACAGCCTCCCGTGAGGCTG‑At                                                                                                                                                                                                                                                               >  5:465230/1‑104 (MQ=255)                                                                                                                                                                                 tAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACG‑ACATACAGCCTCCCGTGAGGCTG‑At                                                                                                                                                                                                                                                               <  6:465230/104‑1 (MQ=255)                                                                                                                                                                                 tAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACG‑ACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTg                                                                                                                                                                                                                    >  5:923484/1‑147 (MQ=255)                                                                                                                                                                                 tAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACG‑ACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTg                                                                                                                                                                                                                    <  6:923484/147‑1 (MQ=255)                                                                                                                                                                                  aaTCCCCCGGGGCGATAGAACCGTACTGTCCAATACAATGCAAATGGGTGATATGCATCAGACCCACCTTAATAAACA‑ACAGACATCCTCCCGTGAGGCGT‑GTGTAAAACCCAAAAACTGCACCGACCGGTCCCCTCCCCCTTTTGGGGAGAGGGTTGGGGAGGGGGAAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACGGGATACGCGGGATCACCa                                                                                   <  4:445035/275‑1 (MQ=255)                                                                                                                                                                                     ccccctggcaatgTGAAACCGTTCGTCCACTACAGTGACTTTTGGTGTTTTGCATACAATGCCCCTTAATTAACG‑CCTTACAGCCCCCCGTGAGGCTG‑TTGTCAAACCCAAAACCTGCCCCGGACGGCCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGGGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGATCCATGCCGGGCGCTAACGGGATACGCGGGATCACCAgg                                                                                 <  4:87624/262‑1 (MQ=255)                                                                                                                                                                                     cccccGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACG‑ACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTg                                                                                                                                                                                                                    <  8:876529/143‑1 (MQ=255)                                                                                                                                                                                     cccccGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACG‑ACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTg                                                                                                                                                                                                                    >  7:876529/1‑143 (MQ=255)                                                                                                                                                                                        ccGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACG‑ACATACAGCCTCACGGGAGGCt                                                                                                                                                                                                                                                                   >  7:435844/1‑94 (MQ=255)                                                                                                                                                                                        ccGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACG‑ACATACAGCCTCACGGGAGGCt                                                                                                                                                                                                                                                                   <  8:435844/94‑1 (MQ=255)                                                                                                                                                                                        ccGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACG‑ACATACAGCCTCACGGGAGGCTgta                                                                                                                                                                                                                                                                <  6:239365/97‑3 (MQ=255)                                                                                                                                                                                        ccGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACG‑ACATACAGCCTCACGGGAGGCTgta                                                                                                                                                                                                                                                                >  5:239365/1‑95 (MQ=255)                                                                                                                                                                                         cgcgGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACG‑ACATACAGCCTCCCGTGAGGCTG‑ATGTCAAAcc                                                                                                                                                                                                                                                       <  5:914830/104‑1 (MQ=255)                                                                                                                                                                                         cgcgGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACG‑ACATACAGCCTCCCGTGAGGCTG‑ATGTCAAAcc                                                                                                                                                                                                                                                       >  6:914830/1‑104 (MQ=255)                                                                                                                                                                                             tgagtagAAACGGTGTGCCCACAGAAAAGAAGAGTGGTTAGTTGGTTGCAAAGCACCTGAATTCAGGCACA‑GCAGCACCCCGTGGGCCGG‑AGGTGAAACCAAGAACCGGCGCCGGCCGGGCCCCGGGCCCCTTGAGGGGGAGGGTTAGGGGGACGGGAAGACCCGCACCCTATTTTCAGGTCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCATCCGCTGGAGCCATGCCGGGCGCTAACGGGATACGCGGGATCACCAGGCAATgcag                                                                         <  4:281740/268‑1 (MQ=255)                                                                                                                                                                                               aaGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACG‑ACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGCGAGAGGGTTAGGCAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCATCCCCCCTTCAACCGCTGGACCCATGCCGTGCGCTAACGGGATACCCGTGCTCCCCCGGCAAtcccg                                                                         >  4:250076/1‑268 (MQ=255)                                                                                                                                                                                                            tGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACG‑ACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTggg                                                                                                                                                                                                                  >  4:88582/1‑122 (MQ=255)                                                                                                                                                                                                                cACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACG‑ACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGgaga                                                                                                                                                                                                              <  2:88583/122‑1 (MQ=255)                                                                                                                                                                                                                  cTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACG‑ACATACAGCCTCCCGTGAgg                                                                                                                                                                                                                                                                     >  6:563514/1‑66 (MQ=255)                                                                                                                                                                                                                  cTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACG‑ACATACAGCCTCCCGTGAgg                                                                                                                                                                                                                                                                     <  5:563514/66‑1 (MQ=255)                                                                                                                                                                                                                  cTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACG‑ACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAggg                                                                                                                                                                                                     <  7:191727/129‑1 (MQ=255)                                                                                                                                                                                                                  cTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACG‑ACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAggg                                                                                                                                                                                                     >  8:191727/1‑129 (MQ=255)                                                                                                                                                                                                                        tGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACG‑ACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACt                                                                                                                                                                                                                                               >  1:27585/1‑81 (MQ=255)                                                                                                                                                                                                                        tGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACG‑ACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTg                                                                                                                                                                                                                    <  7:444499/108‑1 (MQ=255)                                                                                                                                                                                                                        tGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACG‑ACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTg                                                                                                                                                                                                                    >  8:444499/1‑108 (MQ=255)                                                                                                                                                                                                                         gCATATTGGTGATTTGCATACAAGGCACCTTAATTAACG‑ACATACAGCCTCCCGTGAGGCTG‑ATg                                                                                                                                                                                                                                                              >  5:1086383/1‑65 (MQ=255)                                                                                                                                                                                                                         gCATATTGGTGATTTGCATACAAGGCACCTTAATTAACG‑ACATACAGCCTCCCGTGAGGCTG‑ATg                                                                                                                                                                                                                                                              <  6:1086383/65‑1 (MQ=255)                                                                                                                                                                                                                         gCATATTGGTGATTTGCATACAAGGCACCTTAATTAACG‑ACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAATCTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAggg                                                                                                                                                                                                           <  5:893759/116‑1 (MQ=255)                                                                                                                                                                                                                         acATATTGGTGATTTGCATACAAGGCACCTTAATTAACG‑ACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAggg                                                                                                                                                                                                           >  6:893759/2‑116 (MQ=255)                                                                                                                                                                                                                          cATATTGGTGATTTGCATACAAGGCACCTTAATTAACG‑ACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTg                                                                                                                                                                                                                    <  7:200494/106‑1 (MQ=255)                                                                                                                                                                                                                          cATATTGGTGATTTGCATACAAGGCACCTTAATTAACG‑ACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTg                                                                                                                                                                                                                    >  8:200494/1‑106 (MQ=255)                                                                                                                                                                                                                            tatTGGTGATTTGCATACAAGGCACCTTAATTAACG‑ACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCAc                                                                                                                                                                                                                                           <  3:27585/81‑1 (MQ=255)                                                                                                                                                                                                                                   gATTTGCATACAAGGCACCTTAATTAACG‑ACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTgg                                                                                                                                                                                                                   <  7:658136/98‑1 (MQ=255)                                                                                                                                                                                                                                   gATTTGCATACAAGGCACCTTAATTAACG‑ACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTgg                                                                                                                                                                                                                   >  8:658136/1‑98 (MQ=255)                                                                                                                                                                                                                                    aTTTGCATACAAGGCACCTTAATTAACG‑ACATACAGCCTCCCGTGAGGCTG‑a                                                                                                                                                                                                                                                                >  6:885603/1‑52 (MQ=255)                                                                                                                                                                                                                                    aTTTGCATACAAGGCACCTTAATTAACG‑ACATACAGCCTCCCGTGAGGCTG‑a                                                                                                                                                                                                                                                                <  5:885603/52‑1 (MQ=255)                                                                                                                                                                                                                                         cATACAAGGCACCTTAATTAACG‑ACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGACGGGAAAACCCGCACCCTATTTTCAGATCNGTGATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCGAACGGGAGACGCGGGNACACCAGGCAATGCACCAGGGGATAAAGATNNTGTTTCGTGTCCCACACTTTGgtggtg                              >  3:64913/1‑273 (MQ=255)                                                                                                                                                                                                                                               aGGCACCTTAATTAACG‑ACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGgaga                                                                                                                                                                                                              >  7:140545/1‑91 (MQ=255)                                                                                                                                                                                                                                               aGGCACCTTAATTAACG‑ACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGgaga                                                                                                                                                                                                              <  8:140545/91‑1 (MQ=255)                                                                                                                                                                                                                                               aGGCACCTTAATTAACG‑ACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATNTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACGGGATACGCGGGATACCCNGGCAATGCAGACGGGGATAAATATCCCGTTNNCCATCCCAGACTTTGGCGGTTActgtc                        >  3:58543/1‑269 (MQ=255)                                                                                                                                                                                                                                                    ccTTAATTAACG‑ACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGgaga                                                                                                                                                                                                              <  7:602876/86‑1 (MQ=255)                                                                                                                                                                                                                                                    ccTTAATTAACG‑ACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGgaga                                                                                                                                                                                                              >  8:602876/1‑86 (MQ=255)                                                                                                                                                                                                                                                       taattaaCG‑ACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCa                                                                                                                                                                                                                                            >  5:258963/1‑53 (MQ=255)                                                                                                                                                                                                                                                       taattaaCG‑ACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCa                                                                                                                                                                                                                                            <  6:258963/53‑1 (MQ=255)                                                                                                                                                                                                                                                       taattaaCG‑ACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGccc                                                                                                                                                                                                                         >  7:202905/1‑72 (MQ=255)                                                                                                                                                                                                                                                       taattaaCG‑ACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGccc                                                                                                                                                                                                                         <  8:202905/72‑1 (MQ=255)                                                                                                                                                                                                                                                               g‑ACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGAGCAGTTATTTCGCATTAATATCCAGCAGATCCGCTGCACCCGCTGGAGCCATGCCGGGCGCTAACGGGATACGCGGGATCACCAGGCAATGCAGCAGGTGATAAATATACTGTTTTCCGTCCCAGACTTTGGTGGTGACCTTATTGTCCgatttcatctc       >  1:22945/1‑264 (MQ=255)                                                                                                                                                                                                                                                                    tACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACGGGATACGCGGGATCACCAGGCAATGCAGCAGGTGATAAATATCCTGTTTGCCGTCCCAGACTTTGGTGGTGACCTTATTGTCCGCATCCAGCTCCTGcc  >  2:523476/1‑275 (MQ=255)                                                                                                                                                                                                                                                                      cAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAgg                                                                                                                                                                                                <  6:859602/83‑1 (MQ=255)                                                                                                                                                                                                                                                                      cAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAgg                                                                                                                                                                                                >  5:859602/1‑83 (MQ=255)                                                                                                                                                                                                                                                                        gCCTCACGGGAGGCTGTATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCtttt                                                                                                                                                                                                                    <  6:892725/62‑2 (MQ=25)                                                                                                                                                                                                                                                                        gCCTCACGGGAGGCTGTATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTg                                                                                                                                                                                                                    >  5:892725/1‑62 (MQ=37)                                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                                           GATATTCCTGCCCCTGATAGCGGTAACAGCCGTTGGCGATCCGGTTCGGCCACGGAATAAGCACTTTACCCAGATGCGCCAGCGGCATCTCTTCGGGTTTGTGCGGGATGACTAAATGACAACCCTGAAAGGTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACG‑ACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACGGGATACGCGGGATCACCAGGCAATGCAGCAGGTGATAAATATCCTGTTTGCCGTCCCAGACTTTGGTGGTGACCTTATTGTCCGCATCCAGCTCCTGCC  >  NC_000913/4070143‑4070675

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑