Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	70,740	T→G	G118G (GGT→GGG)	araC →	DNA‑binding transcriptional dual regulator AraC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	70,740	0	T	G	91.2%	175.4 / 5.8	68	G118G (GGT→GGG)	araC	DNA‑binding transcriptional dual regulator AraC
Reads supporting (aligned to +/- strand):  ref base T (6/0);  new base G (36/26);  total (42/26)
Fisher's exact test for biased strand distribution p-value = 7.53e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

GACCGCTGGGAATGAAAGGTTATATTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCT‑GTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGGTATCACCAGTGGGTTTACTTTC‑GTCCGCGCGCCTACTGGCATGAATGGCTTAAC‑TGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCAT‑TAACGCCGGGCAAGGGGAAGGGCGCTATTCGGAGCTGCTGGCGATAAATCTGCTTGAGCAATTGTTACTGCGGCG‑CATGGAAGCGATTAACGAGTCGCTCCATCCACCGATGGATAATCGGGTACGCGAGGCTTGTCAGTACATCAGCGATCACCTGGCAGACAGCAATTTTGATATC  >  NC_000913/70499‑70977                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                gaccgCTGGGAATGAAAGGTTATATTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCT‑GTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGGTATCACCAGTGGGTTTACTTTC‑GTCCGCGCGCCTACTGGCATGAATGGCTTAAC‑TGGCCGTCAATATTTGCCAATACGGGGTTCTTTCGc                                                                                                                                                                                                                                        >  1:17939/1‑251 (MQ=255)  cccgCTGGGTAGGAAAGGTTATTTTCTCAATCTCACCTTGCCGGGTCAGGGGGGGGTAAAAAATCAGGGACGATAATTTGTTTACCGACCGGGTGATATTTTGCT‑TTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGGTATCACCAGTGGGTTTACTTTC‑GTCCGCGCGCCTACTGGCATGAATGGCTTAAC‑TGGCCGTCAATATTTGCCAATACGGGGTTCTTTCGcc                                                                                                                                                                                                                                       <  2:476797/250‑1 (MQ=255)     gCTGGGAATGAAAGGTTATATTCTCAATCTCACCATTCGCGGTTAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCT‑GTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGGTATCACAAGTGGGTTTACTTTC‑GTCCGCGCGCATACTGGCATGAATGGCTTAAC‑TGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGCCCgg                                                                                                                                                                                                                                    >  1:356394/1‑251 (MQ=255)       tGGGAATGAAAGGTTATATTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCT‑GTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGGTATCACCAGTGGGTTTACTTTC‑GTCCGCGCGCCTACTGGCATGAATGGCTTAAC‑TGGCCGTCAATATTTGCCAATACGGGGTTCTTTCGCCCGGAt                                                                                                                                                                                                                                  >  1:9993/1‑251 (MQ=255)           aaTGAAGGGTTATATTCTCAATTCCACCATTCGCGGTCAGGGGTTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCT‑GTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGGTATCACCAGTGGGTTTACTTTC‑GTCCGCGCGCCTACTGGCATGAATGGCTTAAC‑TGGCCGTCAATATTTGCCAATACGGGGTTCTTTCGCCCGGATGAAg                                                                                                                                                                                                                              <  1:469344/251‑1 (MQ=255)           aaTGAAAGGTTATATACACAAGCTCACCATTCGCGGTCAGGGGGTGGTAAAAAAACAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCT‑GTTCCCGCCAGGAGAGATTCATCTCTACGGTCGTCATCCGGAGGCTCGCGAATGGTATCACCAGTGGGTTTACTTTC‑GTCCGCGCGCCTACTGGCATGAATGGCTTAAC‑TGGCCGTCAATATTTGCCAATACGGGGTTCTTTCGCCCGGATGAAg                                                                                                                                                                                                                              <  1:431184/251‑1 (MQ=255)                        atTCTCAATCTCACCATTCGCGGTCAGGGGGAGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCT‑GTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGGTATCACCAGTGGGTTTACTTTC‑GTCCGCGCGCCTACTGGCATGAATGGCTTAAC‑TGGCCGTCAATATTTGCCAATACGGGGTTCTTTCGCCCGGATGAAGCGCACCAGCCGc                                                                                                                                                                                                                  >  1:295647/1‑250 (MQ=255)                         ttCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCT‑GTTACCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGGTATCACCAGTGGGTTTACTTTC‑GTCCGCGCGCCTACTGGCATGAATGGATTAAC‑TGGCAGTCAATATTTGCCAATACGGGGTTCTTTCGCCCGGATGAAgcgccaccgccgcct                                                                                                                                                                                                                >  1:50450/1‑240 (MQ=255)                          tctcAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCT‑GTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGGTATCACCAGTGGGTTTACTTTC‑GTCCGCGCGCCTACTGGCATGAATGGCTTAAC‑TGGCCGTCAATATTTGCCAATACGGGGTTCTTTCGCCCGGCTGAAGCGCACCAGCCGCAtt                                                                                                                                                                                                               >  1:182478/1‑251 (MQ=255)                           ctcAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCT‑GTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGGTATCACCAGTGGGTTTACTTTC‑GTCCGCGCGCCTACTGGCATGAATGGCTTAAC‑TGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCCCCAGCCGCAttt                                                                                                                                                                                                              >  1:202501/1‑251 (MQ=255)                                tatcACCATACGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGACTTTGTTTGCCGAACGGATGATATTTTACT‑GTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATTGTATCACCAGTGGGTTTACTTTC‑GTCCGCGCGCCATCTGGAATGAAGGTCTTAACGTGGC‑GTCAATTTTTGACAAAACGGGGCTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTAAGCg                                                                                                                                                                                                         >  2:533456/3‑251 (MQ=255)                                   cACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCT‑GTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGGTATCACCAGTGGGTTTACTTTC‑GTCCGCGCGCCTACTGGCATGAATGGCTTAAC‑TGGCCGTCAATATTTGCCAATACGGGGTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGAcc                                                                                                                                                                                                      >  1:574027/1‑251 (MQ=255)                                        ttCGCGGTCAGGGGGTGGTGAAAATTCAGGGACGAGATTTTTTTGGCCTACGGGGTGATATTTTGTT‑GTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGGTATCACCAGTGGGTTTACTTTC‑GTCCGCGCGCCTACTGGCATGAATGGCTTAAC‑TGGCCGTCAATATTTGCCAATACGGGGTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGtt                                                                                                                                                                                                  <  1:320954/250‑1 (MQ=255)                                             ggTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCT‑GTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGGTATCACCAGTGGGTTTACTTTC‑GTCCGCGCGCCTACTGGCATGAATGGCTTAAC‑TGGCCGTCAATATTTGCCAATACGGGGTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGc                                                                                                                                                                                             >  2:705362/1‑250 (MQ=255)                                             ggTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCT‑GTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGGTATCACCAGTGGGTTTACTTTC‑GTCCGCGCGCCTACTGGCATGAATGGCTTAAC‑TGGCCGTCAATATTTGCCAATACGGGGTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGc                                                                                                                                                                                             >  1:29091/1‑250 (MQ=255)                                                 aGGGGGTGGTGAACAATCAGGGACGAGAATTTTTTTGCCGACCGGGTGATATTTTGAT‑GTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGGGATCACCAGTGGGTTTACTTTC‑GACCGCGCGCCTACTGGCATGCATGGCTTAAC‑TGGCCGTCAATATTTGCCAATACGGGGTTCTTTCGCCCGGATGAAACGCCCCCACCGCATTTCACCGCCCTGTTTTGGCCAATc                                                                                                                                                                                        >  2:580660/1‑251 (MQ=255)                                                        ggtgAAAAATAAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCT‑GTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGGTATCACCAGTGGGTTTACTTTC‑GTCCGCGCGCCTACTGGCATGAATGGCCTCAC‑TGGCCGTCACTATTTGCCCATACGGGGTTTTTTCGCACGGATGAAGCGCCCCAGCCGACTTTCAGCGACCTGTTTGGGCAAATCAT‑TAACg                                                                                                                                                                                >  2:317415/1‑251 (MQ=255)                                                                   aGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCT‑GTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGGTATCACCAGTGGGTTTACTTTC‑GTCCGCGCGCCTACTGGCATGAATGGCTTAAC‑TGGCCGTCAATATTTGCCAATACGGGGTTCTTTCGCCCGGATGAAGCGCACCCGCCGCATTTCAGCGACCTGTTTGGGCAAATCAT‑TAACGCCGGGCAAggg                                                                                                                                                                     >  2:601780/1‑251 (MQ=255)                                                                    gggatggatttgtgcttGCCGACGGGGTGATATTTGGCT‑GTTTCCGCCAGGAGAGATTCAGCACTACGGTCGGCATCCGGAGGCTCGTGAATGGTATCACCAGTGGGTTTACTTTC‑GTCCGCGCGCCTACGTGCATGAATGGCTTAAC‑TGGCCGTCAATATTTGCCAATACGGGGTTCTGTCGGCCGGATGAAGCTCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAAGCAT‑TAACGCCGGGCAAgggg                                                                                                                                                                    <  2:40661/235‑1 (MQ=255)                                                                       aCGAGAATTTGTTTGCCGACCGGGTGATATTTTGCT‑GTTACCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGCGGCTCGCGAATGGTATCACCAGTGGGTTTACTTTC‑GTCCGCGCGCCTACTGGCATGAATGGCTTAAC‑TGGCCGTCAATATTTGCCAATACGGGGTTCTTTCGCCCGGATGAAGCGCACCAGCAGCATTTCAGCGACCTGTTTGGGCAAATCAT‑TAACGCCGGGCAAGGGGa                                                                                                                                                                   >  2:415032/1‑249 (MQ=255)                                                                                       cGACCGGGTGATATTTTGCT‑GTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGGTATCACCAGTGGGTTTACTTTC‑GTCCGCGCGCCTACTGGCATGAATGGCTTAAC‑TGGCCGTCAATATTTGCCATTACGGGGTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCAT‑TAACGCCGGGCAAGGGGAGGGGCGCTGTTCGGAgct                                                                                                                                                 >  2:269951/1‑251 (MQ=255)                                                                                       cGACCGGGTAATATTT‑GCTTGTCCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGTGGCTCGCGAATGGTATCACCAGTGGGTTTACTTTC‑GTCCGCGCGCCTACTGGCATGAATGGCTTAAC‑TGGCCGTCAAGATTTGCCAATACGGGGTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCAT‑TAACGCCGGGCAAGGGGAAGGGCGCTATTCGGAgct                                                                                                                                                 <  1:302573/251‑1 (MQ=255)                                                                                         aCCGGGGGATTTTTTGCT‑GTTTCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGGTATCACCAGTGGGTTTACTTTC‑GTCCGCGCGCCTACTGGCATGAATGGCTTAAC‑TGGCCGTCAATATTTGCCAATACGGGGTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCAT‑TAACGCCGGGCAAGGGGAAGGGCGCTATTCGGAgctgc                                                                                                                                               <  1:424457/251‑1 (MQ=255)                                                                                          ccGGGTGATATTTTGCT‑GTTCCCGCCAGGAGAGATTCATCACTACGATCGTCATCCGGAGGCTCGCGAATGGTATCACCAGTGGGTTTACTTTC‑GTCCGCGAGCCTACTGGCATGAATGGCTTAAC‑TGGCCGTCAATATTTACCCATACGGGGTTCTTTCGCCCGGTAGAAGGGCAACAGCCGCATTTCAGCGACCTGTTTGGGCAAATAAT‑TAACGCAGGGCAAGGGGAAGGGCGCAATTCGGAgaggcg                                                                                                                                              >  2:539842/1‑246 (MQ=255)                                                                                                      ttGCT‑GTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGGTATCACCAGTGGGTTTACTTTC‑GTCCGCGCGCCTACTGGCATGAATGGCTTAAC‑TGGCCGTCAATATTTGCCAATACGGGGTTCTTTCGCCCGGGTGAAGCGCCCCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCAT‑TAACGCCGGGCACGGGGAAGGGCGCTATTCGGGGCTGCTGGCGAAAAATCt                                                                                                                                  >  2:299749/1‑251 (MQ=255)                                                                                                       tGCT‑GTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGGTATCACCAGTGGGTTTACTTTC‑GTCCGCGCGCCTACTGGCATGAATGGCTTAAC‑TGGCCGTCAATATTTGCCAATACGGGGTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCAT‑TAACGCCGGGCAAGGGGAAGGGCGCTATTCGGAGCTGCTGGCGGTAAATCTg                                                                                                                                 >  1:493651/1‑251 (MQ=255)                                                                                                       tGCT‑GTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGGTATCACCAGTGGGTTTACTTTC‑GTCCGCGCGCCTACTGGCATGAATGGCTTAAC‑TGGCCGTCAATATTTGCCAATACGGGGTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCAT‑TAACGCCGGGCAAGGGGAAGGGCGCTATTCGGAGCTGCTGGCGATAAATCTg                                                                                                                                 >  1:562125/1‑251 (MQ=255)                                                                                                        gCT‑GTGCCCGCCAGGAGATATTCATCACTACGGTTGTCATCCGGAGGCTCGCGAATGGTATCACCAGTGGGTTTACGTTC‑GTCCGCGCGCCTACTGGCATGAATGGCTTAAC‑TGGCCGTCAATATTTGCCAATACGGGGTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCAT‑TAACGCCGGGCAAGGGGAAGGGCGCTATTCGGAGCTGCTGGCGATAAATCTGc                                                                                                                                <  1:377047/251‑1 (MQ=255)                                                                                                                    cgtgggagaTTCATACCTACGTCCGTCATCCGGAGGTTCGCTAAGGGTAGGACCAGTGTGGTTACTTTC‑GTCCGCGCGCCTACTTGCATGATGGGCTTAAC‑TGGCCTTCAATATTTGCCAATACGGGGTTCTTTCGCCCGGATGAAGGTCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAAGCAT‑TAACGCCGGGCAAGGGGAAGGTCGCTATTCGGAGCTGCTGGCGATAAACCTGCTTGAGCAATTg                                                                                                                     <  2:434449/246‑1 (MQ=255)                                                                                                                      ggAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGGTATCACCCGTGGGTTTACTTTC‑GTCCGCGCGCCTACTGGCATGAATGGCTTAAC‑TGGCCGTCAATATTTGCCAATACGGGGTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCAT‑TAACGCCGGGCAAGGGGAAGGGCGCTATTCGGAGCTGCTGGCGATAAATCTGCTTGAGCAATTg                                                                                                                     >  1:434449/1‑249 (MQ=255)                                                                                                                       gagagaTTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGGTATCACCAGTGGGTTTACTTTC‑GTCCGCGCGCCTACTGGCATGAATGGCTTAAC‑TGGCCGTCAATATTTGCCAATACGGGGTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCAT‑TAACGCCGGGCAAGGGGAAGGGCGCTATTCGGAGCTGCTGGCGATAAATCTGCTTGAGCAATTGtt                                                                                                                   >  1:672794/1‑250 (MQ=255)                                                                                                                         gagaTTCCTCCCTCCGGTCGTCATCCGGAGGCTCGCGAATGGTATCACCAGTGGGTTTCCTTTC‑GTCCGCGCGCCTACTGGCATGAATGGCTTAAC‑TGGCCGTCACTCTTTGCCAATACGGGGTTCTTTCGCCCGGCTGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCATATAA‑GCCGGGCAGGGGGGAAGGCGCTATTAGGAGCTGCCGTCGCTAAATCTGCTTGCGCCATTGTTACt                                                                                                                >  2:313346/1‑251 (MQ=255)                                                                                                                           gatgcatcaCAACGGTGGTCATCGTGAGGCTATCTAATGGTATCCCCAGGGGGTTTACTGTC‑GTCCGCGCGATTACTGGCATGAATGGCTCAAC‑TGGCCGTCAATAGTTGCCAATACGGGGTTCTTTCGCCCGGAGGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCAT‑TAACGCCGGGCAAGGGGAAGGGCGCTATGCGGAGCTGCTGGCGATAAATCTGCTTGAGCAATTGTTAcggc                                                                                                              <  2:182478/247‑4 (MQ=255)                                                                                                                                   acgaCGGTCGTCATCCGGAGGCTCGCGAATGGTATCACCAGTGGGTTTACTTTC‑GTCCGCGCGCCTACTGGCATGAATGGCTTAAC‑TGGCCGTCAATATTTGCCAATACGGGGTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCAT‑TAACGCCGGGCAAGGGGAAGGGCGCTATTCGGAGCTGCTGGCGATAAATCTGCTTGAGCAATTGTTACTGCGGCG‑CATg                                                                                                     <  1:479319/248‑1 (MQ=255)                                                                                                                                         gtcgtcATCCGGAGGCTCGCGAATGGTATCACCAGTGGGTTTACTTTC‑GTCCGCGCGCCTACTGGCATGAATGGCTTAAC‑TGGCCGTCAATATTTGCCAATACGGGGTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCAT‑TAACGCCGGGCAAGGGGAAGGGCGCTATTCGGAGCTGCTGGCGATAAATCTGCTTGAGCAATTGTTACTGCGGCG‑CATGGAAGCg                                                                                               >  2:202939/1‑251 (MQ=255)                                                                                                                                         gctgtctaccGGGGGCTCGCGAGTGGTCTCCCCAGTGGGTTTACTTTC‑GTCCGCGCGCCTACTGGCATGAATGGCTTAAC‑TGGCCGGCAAGATTTGCCAAGCTGGGGTTCTTTCGCCCGGATGAAGCGCCCCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCAT‑TAACGCCGGGCAAGGGGAAGGGCGCTATTCGGAGCTGCTGGCGATAAATCTGCTTGAGCAATTGTTACTGCGGAG‑CATGGAAGCg                                                                                               <  2:493651/243‑1 (MQ=255)                                                                                                                                           cgtcATCCGGAGGCTCGCGCATGGTATCACCAGTGGGTTTACTTTC‑GTCCGCGCGCCTACTGGCATGAATGGCTTAAC‑TGGCCGTCAATATTTGCCAATACGGGGTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCAT‑TAACGCCGGGCAAGGGGAAGGGCGCTATTCGGAGCTGCTGGCGATAAATCGGCTTGAGCAATTGTTAATGCGGCG‑CATGGAAGCGAt                                                                                             >  2:678569/1‑251 (MQ=255)                                                                                                                                           cgtcACCCGGAGGCCAGCGATTGGTAGCACCAGTGGGTTT‑CTTTCCGGCCGCGTGCCTACTGGCATGAATGGCTTAAC‑TGGCCGACAATATTTGCCAATACGGGGTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCAT‑TAACGCCGGGCAAGGGGAAGGGCGCTAGTCGGAGCTGCTGGCGATAAATCTGCTTGAGCAATTGTTACTGCGGCG‑CATGGAAGCGa                                                                                              <  2:91434/250‑1 (MQ=255)                                                                                                                                            gtcATCCGGAGGCTCGCGAATGGTATCACCAGTGGGTTTACTTTC‑ATCCGCGCGCCTACTGGCATGAATGGCTTAAC‑TGGCCGTCAATATTTGCCAATACGGGGTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCACATCAT‑TAACGCCGGGCAAGGGGCAGGGCGCTATTCGGAGCGGCAGGCCATAAATCTGCTTGGGCAATTGTTACTGCGGCG‑CCTGGacgccatt                                                                                            >  2:141403/1‑244 (MQ=255)                                                                                                                                                  cGGAGGCTCGCGAATGGTATCACCAGTGGGTTTACTTTC‑GTCCGCGCGCCTACTGGCATGAATGGCTTAAC‑TGGCCGTCAATATTTGCCCATACGGGGTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCAT‑TAACGCCGGGCAAGGGGAAGGGCGCTATTCGGAGCTGCTGGCGATAAATCTGCTTGAGACATTGTTACTGCGGCG‑CCTGGAAGCGGTTAACg                                                                                        >  2:339029/1‑249 (MQ=255)                                                                                                                                                   ggaggctgtcgaaggggatccccggTGGGTTTACTTAC‑GTCCGAGCGCCTATTGGCATGAATTGCTTAAC‑TGGCCGTCAATATTTGCCAATACGGGGTTCTTTCGCCCGGATGAGTCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCAT‑TAACGCCGGGCAAGGGGAAGGGCGCTATTCGGAGCTGCTGGCGATAAATCTGCTTGAGCAATTGTTACTGCGGCG‑CATGGAAGCGATTAACGAgg                                                                                     <  2:595927/227‑2 (MQ=255)                                                                                                                                                            cgAATGGTATCACCAGTGGGTTTACTTTC‑GTCCGCGCGCCTACTGGCATGAATGGCTTAAC‑TGGCCGTCAATATTTGCCAATGCGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCAT‑TAACGCCGGGCAAGGGGAAGGGCGCTATTCGGAGCTGCTGGCGATAAATCTGCTTGAGCAATTGTTACTGCGGCG‑CATGGAAGCGATTAACGAGTCGCtccatc                                                                            >  1:162183/1‑251 (MQ=255)                                                                                                                                                            cgAATGGTATCACCAGTGGGTTTACTTTC‑GTCCGCGCGCCTACTGGCATGAATGGCTTAAC‑TGGCCGTCAATATTTGCCAATACGGGGTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCAT‑TAACGCCGGGCAAGGGGAAGGGCGCTATTCGGAGCTGCTGGCGATAAATCTGCTTGAGCAATTGTTACTGCGGCG‑CATGGAAGCGATTAACGAGTCGCtccatc                                                                            <  1:159624/251‑1 (MQ=255)                                                                                                                                                             gAATGGTAGCCCCAGTGGGTTTACTTTC‑GTCCGCGCGCCTACTGGCATGAAGGGCTTAAC‑TGGCCGTCAATATTTGCCAATACGGGGTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCAT‑TAACGCCGGGCAAGGGGAAGGGCGCTATTCGGAGCTGCTGGCGATAAATCTGCTTGAGCAATTGTTACTGCGGCG‑CATGGAAGCGATTAACGAGTCGCtccagcc                                                                           <  2:574027/251‑4 (MQ=255)                                                                                                                                                                 ggTATCACCAGTGGGTTTACTTTC‑GTCCGCGCGCCTACTGGCATGAATGGCTTAAC‑TGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCAT‑TAACGCCGGGCAAGTGGAAGGGCGCTATTCGGAGCTGCTGGCGATAAATCTGCTTGGGCAATTGTTACTGCGGCG‑CATGGAAGCGATTAACGAGTCGCTCCATCCACCg                                                                       >  1:208687/1‑251 (MQ=255)                                                                                                                                                                        ccAGTGGGTTTACTTTC‑GTCCGCGCGCCTACTGGCATGAATGGCTTAAC‑TGGCCGTCAATATTTGCCAATACGGGgttctt                                                                                                                                                                                                                                            <  1:675886/81‑1 (MQ=255)                                                                                                                                                                        ccAGTGGGTTTACTTTC‑GTCCGCGCGCCTACTGGCATGAATGGCTTAAC‑TGGCCGTCAATATTTGCCAATACGGGgttctt                                                                                                                                                                                                                                            >  2:675886/1‑81 (MQ=255)                                                                                                                                                                         cAGTGGGTTTACTTTC‑GTCCGCGCGCCTACTGGCATGAATGGCTTAAC‑TGGCCGTCAATATTTGCCAATACGGGGTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCAT‑TAACGCCGGGCAAGGGGAAGGGCGCTATTCGGAGCTGCTGGCGATAAATCTGCTTGAGCAATTGTTACTGCGGCG‑CATGGAAGCGATTAACGAGTCGCTCCATCCACCGATGGATaa                                                               <  1:705362/251‑1 (MQ=255)                                                                                                                                                                             gggTTGACTTCC‑TTCCCCGCGCCTACGGGCAGAAATGTCTTAAC‑TGGCCGTCAAAATTTCCCATGACGGGGTTCTTGCGCCCGGATGAAGCGCCCCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCAT‑TAACGCCGGGCAAGGGGAAGGGCGCTATTCGGAGCTGCTGGCGATAAATCTGCTTGAGCAATTGTTACTGCGGCG‑CATGGAAGCGATTAACGAGTCGCTCCATCCACCGATGGATAATCg                                                            <  2:447482/250‑1 (MQ=255)                                                                                                                                                                              ggtttaaattC‑GTCCGCGCGCCTACTGGCATGAATGGCTTAAC‑TGACCGTCAATATTTGCCAATACGGGGTTATTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCAT‑TAACGCCGGGCACGGGGAAGGGCGCTCTTCGGAGCTGCTGGCGATAAATCTGCTTGAGCAATTGTTAATGCGGCGACA‑GGAAGCGAGTAACGCGTCGCTCCATCAACCGATGGAAAATCggg                                                          >  2:98586/9‑251 (MQ=255)                                                                                                                                                                                 ttACTTTC‑GTCCGCGGACCTACTGGCATGAATGGCTTAAC‑TGGCCGGCAATATTTGCCAATACGGGGTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCAT‑TAACGCCGGGCAAGGGGAAGGGCGCTATTCGGAGCTGCTGGCGATAAATCTGCTTGAGCAATTGTTACTGCGGCG‑CATGGAAGCGATTAACGAGTCGCTCCATCCACCGATGGATAATCGGGTa                                                        <  2:9993/250‑1 (MQ=255)                                                                                                                                                                                    cTTTC‑GTCCGCGCGCCTACTGGCATGAATGGCTTAAC‑TGGCCGTCAATATTTGCCAATACGGGGTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCAT‑TAACGCCGGGCAAGGGGAAGGGCGCTATTCGGAGCTGCTGGCGATAAATCTGCTTGAGCAATTGTTACTGCGGCG‑CATGGAAGCGATTAACGAGTCGCTCCATCCACCGATGGATAATCGGGTAcgcg                                                    >  1:383279/1‑251 (MQ=255)                                                                                                                                                                                          gTCCGCGCGCCTACTGGCATGAATGGCTTAAC‑TGGCCGTCAATATTTGCCAATACGGGGTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCAT‑TAACGCCGGGCAAg                                                                                                                                                                       >  2:578579/1‑132 (MQ=255)                                                                                                                                                                                          gTCCGCGCGCCTACTGGCATGAATGGCTTAAC‑TGGCCGTCAATATTTGCCAATACGGGGTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCAT‑TAACGCCGGGCAAg                                                                                                                                                                       <  1:578579/132‑1 (MQ=255)                                                                                                                                                                                             cgcgcgCCTACTGGCATGAATGGCTTAAC‑TGGCCGTCAATATTTGCCAATACGGGGTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCAT‑TAACGCCGGGCAAGGGGAAGGGCGCTATTCGGAGCTGCTGGCGATAAATCTGCTTGAGCAATTGTTACTGCGGCG‑CATGGAAGCGATTAACGAGTCGCTCCATCCACCGATGGATAATCGGGTACGCGAGGCTTGt                                            <  1:601780/251‑1 (MQ=255)                                                                                                                                                                                             cgcgcgCCTACTGGCATGAATGGCTTAAC‑TGGCCGTCAATATTTGCCAATACGGGGTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCAT‑TAACGCCGGGCAAGGGGAAGGGCGCTATTCGGAGCTGCTGGCGATAAATCTGCTTGAGCAATTGTTACTGCGGCG‑CATGGAAGCGATTAACGAGTCGCTCCAACCACCGATGGATAATCGGGGACGCGAGGCTTGt                                            >  2:667811/1‑251 (MQ=255)                                                                                                                                                                                                gcgcCTACTGGCATGAATGGCTTAAC‑TGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCAT‑TAACGCCGGGCAAGGGGAAGGGCGCTATTCGGAGCTGCTGGCGATAAATCTGCTTGAGCAATTGTTACTGCGGCG‑CATGGAAGCGATTAACGAGTCGCTCCATCCACCGATGGATAATCGGGTACGCGCGGCTTGTc                                           >  2:175268/1‑249 (MQ=255)                                                                                                                                                                                                   ccTACTGGCATGAATGGCTTAAC‑TGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCAT‑TAACGCCGGGCAAGGGGAAGGGCGCTATTCGGAGCTGCTGGCGATAAATCTGCTTGAGCAATTGTTACTGCGGCG‑CATGGAAGCGATTAACGAGTCGCTCCATCCACCGATGGATAATCGGGTAGGCGGGGCTTGTCAATAc                                      >  2:583255/1‑251 (MQ=255)                                                                                                                                                                                                                aTGGCTTAAC‑TGGCCGTCAATATTTGCCAATACGGGGTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCAT‑TAACGCCGGGCAAGGGGAAGGGCGCTATTCGGAGCTGCTGGCGATAAATCTGCTTGAGCAATTGTTACTGCGGCG‑CATGGAAGCGATTAACGAGTCGCTCCATCCACCGATGGATAATCGGGTACGCGAGGCTTGTCAGTACATCAGCGAt                             <  1:520765/247‑1 (MQ=255)                                                                                                                                                                                                                aTGGCTTAAC‑TGGCCGTCAATATTTGCCAATACGGGGTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCAT‑TAACGCCGGGCAAGGGGAAGGGCGCTATTCGGAGCTGCTGGCGATAAATCTGCTTGAGCAATTGTTACTGCGGCG‑CATGGAAGCGATTAACGAGTCGCTCCATCCACCGATGGATAATCGGGTACGCGAGGCTTGTCAGTACATCAGCGAt                             >  2:520765/1‑247 (MQ=255)                                                                                                                                                                                                                aTGGCTTAAC‑TGGCCGTCAATATTTGCCAATACGGGGTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCAT‑TAACGCCGGGCAAGGGGAAGGGCGCTATTCGGAGCTGCTGGCGATAAATCTGCTTGAGCAATTGTTACTGCGGCG‑CATGGAAGCGATTAACGAGTCGCTCCATCCACCGATGGATAATCGGGTACGCGAGGCTTGTCAGTACATCAGCGATCAcc                         >  1:32092/1‑251 (MQ=255)                                                                                                                                                                                                                    cTTAAC‑TGGCCGTCAATATTTGCCAATACGGGGTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCAT‑TAACGCCGGGCAAGGGGAAGGGCGCTATTCGGAGCTGCTGGCGATAAATCTGCTTGAGCAATTGTTACTGCGGCG‑CATGGAAGCGATTAACGAGTCGCTCCATCCACCGATGGATAATCGGGTACGCGAGGCTTGTCAGTACATCAGCGATCACCTGGc                     >  1:527625/1‑251 (MQ=255)                                                                                                                                                                                                                            ggCCGTCAATATTTGCCAATACGGGGTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCAT‑TAACGCCGGGCAAGGGGAAGGGCGCTATTCGGAGCTGCTGGCGATAAATCTGCTTGAGCAATTGTTACTGCGGCG‑CATGGAAGCGATTAACGAGTCGCTCCATCCACCGATGGATAATCGGGTACGCGAGGCTTGTCCGTACATCAGCGATCAcc                         >  2:222859/1‑240 (MQ=255)                                                                                                                                                                                                                            ggCCGTCAATATTTGCCAATACGGGGTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCAT‑TAACGCCGGGCAAGGGGAAGGGCGCTATTCGGAGCTGCTGGCGATAAATCTGCTTGAGCAATTGTTACTGCGGCG‑CATGGAAGCGATTAACGAGTCGCTCCATCCACCGATGGATAATCGGGTACGCGAGGCTTGTCAGTACATCAGCGATCAcc                         <  1:222859/240‑1 (MQ=255)                                                                                                                                                                                                                            cgCCGTCAATATTTGCCACAACGGGGCTCTTTCGCCCGGCTGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCAT‑TCACGCCGGGCAAGGGGACGGGCGCTATTCGGAGCTGCTGGCGATCAATCTGCTTGAGCAATTGTTACTGAGGCG‑CATGGAAGCGCTTAACGAGTCGCACCATCCCCCGCTGGCTAATCTGGTACGCGAGGCTTGTCAGTACATCAGAGCTCACCTGGAAGACAGc              >  2:78638/2‑251 (MQ=255)                                                                                                                                                                                                                              ccGTCAATATTTGCCAATACGGGGTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCAT‑TAACGCCGGGCAAGGGGAAGGGCGCTATTCGGAGCTGCTGGCGATAAATCTGCTTGAGCAATTGTTACTGCGGCG‑CATGGAAGCGATTAACGAGTCGCTCCATCCACCGATGGATAATCGGGTACGCGAGGCTTGTCAGTACATCAGCGATCACCTGGCAGACAGCa             >  1:423538/1‑250 (MQ=255)                                                                                                                                                                                                                                  cAATATTTGCCAATACGGGGTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCAT‑TAACGCCGGGCAAGGGGAAGGGCGCTATTCGGAGCTGCTGGCGATAAATCTGCTTGAGCAATTGTTACTGCGGCG‑CATGGAAGCGATTAACGAGTCGCTCCATCCACCGATGGATAATCGGGTACGCGAGGCTTGTCAGTACATCAGCGATCACCTGGCAGACAGCAAttt         <  1:386945/250‑1 (MQ=255)                                                                                                                                                                                                                                      atTTGCCAATACGGGGTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCAT‑TAACGCCGGGCAAGGGGAAGGGCGCTATTCGGAGCTGCTGGCGATAAATCTGCTTGAGCAATTGTTACTGCGGCG‑CATGGAAGCGATTAACGAGTCGCTCCATCCACCGATGGATAATCGGGTACGCGAGGCTTGTCAGTACATCAGCGATCACCTGGCAGACAGCAATTTTGat     <  1:339029/250‑1 (MQ=255)                                                                                                                                                                                                                                        ttGCCAATACGGGGTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCAT‑TAACGCCGGGCAAGGGGAAGGGCGCTATTCGGAGCTGCTGGCGATAAATCTGCTTGAGCAATTGTTACTGCGGCG‑CATGGAAGCGATTAACGAGTCGCTCCATCCACCGATGGATAATCGGGTACGCGAGGCTTGTCAGTACATCAGCGATCACCTGGCAGACAGCAATTTTGATATc  >  1:530398/1‑251 (MQ=255)                                                                                                                                                                                                                                         gtccattACGGGGTTATTTCGCCCGGAGGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCAT‑TAACGCCGGGCAAGGGGAAGGGCGCTATTCGGAGCTGCTGGCGATAAATCTGCTTGAGCAATTGTTACTGCGGCG‑CATGGAAGCGATTAACGAGTCGCTCCATCCACCGATGGATAATCGGGTACGCGAGGCTTGTCAGTACATCAGCGATCACCTGGCAGACAGCAATGTTGATATc  <  2:174484/244‑1 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                GACCGCTGGGAATGAAAGGTTATATTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCT‑GTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGGTATCACCAGTGGGTTTACTTTC‑GTCCGCGCGCCTACTGGCATGAATGGCTTAAC‑TGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCAT‑TAACGCCGGGCAAGGGGAAGGGCGCTATTCGGAGCTGCTGGCGATAAATCTGCTTGAGCAATTGTTACTGCGGCG‑CATGGAAGCGATTAACGAGTCGCTCCATCCACCGATGGATAATCGGGTACGCGAGGCTTGTCAGTACATCAGCGATCACCTGGCAGACAGCAATTTTGATATC  >  NC_000913/70499‑70977

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑