Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	89,089	C→A	intergenic (+57/‑545)	cra → / → mraZ	DNA‑binding transcriptional dual regulator Cra/DNA‑binding transcriptional repressor MraZ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	89,089	0	C	A	95.4%	167.7 / ‑0.7	65	intergenic (+57/‑545)	cra/mraZ	DNA‑binding transcriptional dual regulator Cra/DNA‑binding transcriptional repressor MraZ
Reads supporting (aligned to +/- strand):  ref base C (2/1);  new base A (21/41);  total (23/42)
Fisher's exact test for biased strand distribution p-value = 2.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.78e-01

TTGGCGATAACGAACTGCTCGACTTCTTACAGTGTCCGGTGCTG‑GCAGTGGCTCAACGTCACCGCGATGTCGCAGAGCGTGTGCTGGAGATTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTG‑AATTTCCCTCTGGCGCGTAGAGTACGGGACTGGACATCAATATGCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAATGCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCGGCGATGAATTGCCGATTTAACAAACACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCACACAGACGTAACAAGGACT  >  NC_000913/88845‑89326                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                              ttGGCGATAACGAACTGTGCGACTTCTTACAGGGTCCGGTGATG‑GGAGTGGCTCAACGTCGCCGCGAAGTCGCAGAGCGTGAGCTGTAGATTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCAAAATCCTGGTTTAACGCGCATTAAACGTAATCTCTACCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAAGACGCGCCAGGTG‑AATGTCCCTCTGGCGCGTAGAGTACGGGAATgggc                                                                                                                                                                                                                                          <  2:116099/250‑3 (MQ=255) ttGGCGATAACGAACTGCTCGACTTCTTTCAGTGTCCGGTGCTG‑GCAGTGGCTCAACGTCACCGCGATGTCGCAGAGCGTGTGCTGGAGATTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTG‑AATTTCCCTCTGGCGCGTAGAGTACGGGAATGGAc                                                                                                                                                                                                                                          <  2:342934/250‑1 (MQ=255) tgggctaTAACGAACTGCTCGACTTCTTACAGGGTCCGGGGCTG‑GCAGTGGCTCAACGTCACCGCGATGTCGCAGAGCGTGTGCTGGAGATTGCCCTGGCTAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACTCGCCAGGTG‑AATTTCCCTCTGTCTCGTAGCGTACGGGAATGGa                                                                                                                                                                                                                                           <  2:574478/243‑1 (MQ=255)   ggCGATAACGAACTGCTCGACTTCTTACCGTGTCCGGTGCTG‑GCAGGGGCTCAACGTCACCGCGATGTCGCAGAGCGTGTGCTGGAGATTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGTCAGGTG‑AATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACAt                                                                                                                                                                                                                                        <  2:768886/250‑1 (MQ=255)   ggCGATAACGAACTGCTCGACTTCTTACAGTGTCCGGTGCTG‑GCAGTGGCTCAACGTCACCGCGATGTCGCAGAGCGTGTGCTGGAGATTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTG‑AATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACAt                                                                                                                                                                                                                                        <  2:863822/250‑1 (MQ=255)    tcGATAACGAACTGCTCGACTTCTTACAGTGTCCGGTGCTG‑GCAGTGGCTCAACGTCACCGCGATGTCGCAGAGCGTGTGCTGGAGATTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTG‑ACTTTCCCTCTGGCGCGTAGAGTACGGGAATGGATATc                                                                                                                                                                                                                                       >  1:780517/2‑250 (MQ=255)      gATAACGAACTGCTCGACTTCTTACAGTGTCCGGTGCTG‑GCAGTGGCTCAACGTCACCGCGATGTCGCAGAGCGTGTGCTGGAGATTGTCCTGGCAAGCCTGGACGACCCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACCAAAATACGCGCCAGGTG‑AATTTCCCTCTAGCGCGTAGAGTACGGGAAGGGACATCaa                                                                                                                                                                                                                                     >  2:453537/1‑250 (MQ=255)       aTAACGACCTGGTCGACTTTTTACATTGGCCGGTG‑TGTGCGGTGTCTCAACGTCACCGCGATGTCGCAGAGCGGGTGCTGGAGATTGTCCTGGCAAGCCTGGACTAACCGCGTAAGCCAAACCCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGACGTAGCTAAGCCGTGAACAAAAATACGCGCCAGGTG‑AATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAata                                                                                                                                                                                                                                   <  2:706685/251‑1 (MQ=255)       aTAACGAACTGCTCGACTTCTTACAGTGTCCGGTGCTG‑GCAGTGGCTCAACGTCACCGCGATGTCGCAGAGCGTGTGCTGGAGATTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTG‑AATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAata                                                                                                                                                                                                                                   <  1:212122/251‑1 (MQ=255)       aTAACGAACTGCTCGACTTCTTACAGTGTCCGGTGCGG‑GCAGTGGCTAAACGGCACCGCTATGTCGCAGAGCGTGTGCTGGAGATTGTCCTTGCAAGCCGGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGGAATCTCTATCGCCGCGGCGTTCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTG‑AATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAata                                                                                                                                                                                                                                   <  2:663902/251‑1 (MQ=255)         aaCGAACTGCTCGACTTCTTACAGTGTCCGGTGCTG‑GAAGTGGCTCAACGTCACCGCGATGTCGCAGAGCGTGTGCTGGAGATTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAATCCGCGAACAAAAATACGCGCCAGGTG‑AATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAatat                                                                                                                                                                                                                                  <  1:545707/250‑1 (MQ=255)           cGAACTGCTCGACTTCTTACAGTGTCAGGTGCTG‑GCAGTGGCTCAACGTCACCGCGATGACGCAGAGCGTGTGCTGGAGATTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTG‑ACTTTCCCTCTGGCGCGTCGCGTACGGGAATGGACCTCAatatcct                                                                                                                                                                                                                               >  2:603214/1‑248 (MQ=255)            gAACTGCTCGACTTCTTACAGTGTCCGGTGCTG‑GCAGTGGCTCAACGTCACCGCGATGTCGCAGAGCGTGTGCTGGAGATTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTG‑AATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCtt                                                                                                                                                                                                                              <  1:459199/251‑1 (MQ=255)              acggCTCGACTTCTTACAGTGTCCGGTGCTG‑GAAGTGGCTAAACGTCACCGCGATGGCGCAGAGCGTGTGCTGGAGATTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCGATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTG‑AATTTCCCGCTGGCGCGTAGAGTACGGGAATGGACAGCAATATGCTTaa                                                                                                                                                                                                                            <  2:724662/248‑1 (MQ=255)                tGCTCGACTTCTTACAGTGTCCGGTGCTG‑GCAGTGGCTCAACGTCACCGCGATGTCGCAGAGCGTGTGCTGGAGATTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTG‑AATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTaa                                                                                                                                                                                                                            >  1:774132/1‑249 (MQ=255)                tGCTCGACTTCTTACAGTGTCCGGAGCTG‑GCAGTGGCTCAACGTCACCGCGATGTCGCAGAGCGTGTGCTGGAGATTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTG‑AATTTCCCTCTGGCGCGTAGAGTCCGGGAATGGACCTCAATATCCTTaa                                                                                                                                                                                                                            >  2:436915/1‑249 (MQ=255)                     gACTTCTTACAGGGTCCGGTGCTG‑GCAGTGGCTCAACGTCACCGCGATGTCGCAGAGCGTGTGCTGGAGATTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTG‑AATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAAt                                                                                                                                                                                                                     <  1:96537/251‑1 (MQ=255)                       cttcttACAGTGTCCGGTGCTG‑GCAGTGGCTCAACGTCACCGCGATGTCGCAGAGCGTGTGCTGGAGATTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTG‑AATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATa                                                                                                                                                                                                                    >  1:766889/1‑250 (MQ=255)                            tACAGTGTCCGGTGCTT‑GCAGTGGTTCAACGTCCCCGCGATTGCGCATGGCGTGTGCTGGAGATTGTCCTGGCAAGCCTGGACGAACCGCTTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCATCCGTAGCTAAGCTGCGAACAAAAATACGCGCCAG‑TGTAATTTCACTCTGGCTCGTAGAGTACGGGAATTGACATCAATATCCTTAAAGTAAAGAAGACt                                                                                                                                                                                                               <  2:524778/250‑1 (MQ=255)                            tACAGTGTCCGGTGCTG‑GCAGTGGCTCAACGTCACCGCGATGTCGCAGAGCGTGTGCTGGAGATTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTG‑AATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGa                                                                                                                                                                                                                 <  1:594090/248‑1 (MQ=255)                                gtgtCCGGTGCTG‑GCAGTGGCTCACCGTCACCGCGATGTCGCAGAGCGTGTGCTGGAGATTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAATCCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTG‑AATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACAGCAATATCCTTAAAGTAAATAAGACTATTcc                                                                                                                                                                                                          <  2:476617/251‑1 (MQ=255)                                      ggTGCTG‑GCAGTGGCTCAACGTCACCGCGATGTCGCAGAGCGTGTGCTGGAGATTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTG‑AATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGCCTATTCCTgcctt                                                                                                                                                                                                    >  1:346025/1‑247 (MQ=255)                                            g‑gCAGTGGCTCAACGTCACCGCGATGTCGCAGAGCGTGTGCTGGAGATTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTG‑AATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTg                                                                                                                                                                                              <  1:436915/251‑1 (MQ=255)                                                aGTGGCTCAACGTCACCGCGATGTCGCAGAGCGTGTGCTGGAGATTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTG‑AATTTCCCTCTGGCGCGTAGAGTACGGGACTGGACATCAATATGCTTAAAGTAAATAAGACTATTCCTGACTATTATTGAt                                                                                                                                                                                            >  1:397238/1‑250 (MQ=255)                                                  tGGCTCAACGTCACCGCGATGTCGCAGAGCGTGTGCTGGAGATTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTG‑AATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATa                                                                                                                                                                                           >  1:359102/1‑249 (MQ=255)                                                   ggCTCAACGTCACCGCGATGTCGCAGAGCGTGTGCTGGAGATTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTG‑AATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATaagt                                                                                                                                                                                        <  1:561376/251‑3 (MQ=255)                                                    gCTCAACGTCACCGCGATGTCGCAGAGCGTGTGCTGGAGATTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTG‑AATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATaagtt                                                                                                                                                                                       >  1:465036/1‑248 (MQ=255)                                                        aaCGTCACCGCGATGTCGCAGAGCGTGTGCTGGAGATTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTG‑AATTTCCCTCTGGCGCGTCGAGTACGGGAATGGCCATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATacgttcttt                                                                                                                                                                                   >  2:233832/1‑243 (MQ=255)                                                        aaCGTCACCGCGATGTCGCAGAGCGTGTGCTGGAGATTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTG‑AATTTCCCTCTGGCGCGTAGAGTACGGGACTGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATaagttcttt                                                                                                                                                                                   >  2:818029/1‑244 (MQ=255)                                                             cACCGCGATGTCGCAGAGCTGGTGCTGGAGATTGTCCTGGTAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTG‑AATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAAc                                                                                                                                                                              <  2:343393/251‑1 (MQ=255)                                                                cgcgATGTCGCAGGGCGTGTGCTGGAGATTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTG‑AATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAAcccg                                                                                                                                                                           <  2:81791/251‑1 (MQ=255)                                                                   gATGTCGCAGAGTGTGGGCTTGAGATTGTCCTGGCAAGCATGGACGACCCGCGTAAGCCAAAACCTGGTTTAACTCGCATTAAACGTAAGCTCGATCGCGGCGGCGTGCTCAGCCGGAGCTAAGCCGCGAACAAAAATACGCGCCAGGTG‑AATGTCCCTATGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTAGTGATAAGTTCTTTTAAAcccgccc                                                                                                                                                                        <  2:827764/251‑1 (MQ=255)                                                                     tGTCGCAGAGCGTGTGCTGGAGATTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAACACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTG‑AATTTCCCTCTGGCGCGGAGAGTACGGGAATGGACATCAAGATCCTTAAAGTAAATAAGACTATTACTGACAAATATTGaaaaggtctttaaaaccagacagg                                                                                                                                                                      >  2:661965/1‑228 (MQ=255)                                                                                  gtgCTGGAGATTGTCCTGGCAATCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTG‑AATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCAcc                                                                                                                                                         <  2:263692/251‑1 (MQ=255)                                                                                         agaTTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAATCGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTG‑AATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGaa                                                                                                                                                  <  1:16643/251‑1 (MQ=255)                                                                                                ccTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTG‑AATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTcac                                                                                                                                            <  1:699800/250‑1 (MQ=255)                                                                                                   ggCAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTG‑AATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACaat                                                                                                                                         <  1:862926/250‑1 (MQ=255)                                                                                                    gCAAGCCAGGACACACCGCGTAAGCAAAAACCTTGTTTAACGCGCATTAAACGTAATCTCTAACGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCAAACAAAAATACGCGCCATGTG‑AATTTCCCTCTTGCGCGTAGAGTACGGGACTGGACATCAATATGCTTAAAGTAAATAAGACTATTCCTGACTATTATTTATAAATGCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTTAAATTCACaataa                                                                                                                                       <  2:375728/251‑1 (MQ=255)                                                                                                     cAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTG‑AATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATACGACTATTCCTGACTATTATTGCTAAGTTCTTTTAAACCCGCCCGTTAATTAACTCCCCAGCTGAAATTCACattaat                                                                                                                                      >  2:183955/1‑251 (MQ=255)                                                                                                      aaGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTG‑AATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAAtt                                                                                                                                     <  1:29695/251‑1 (MQ=255)                                                                                                              aCGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTG‑AATTTCCCTCTGGCGCTTAGAGTTCGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGata                                                                                                                             <  2:848709/251‑1 (MQ=255)                                                                                                                 aaCCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTG‑AATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAAGTCACAATAATTAAGTGATATCg                                                                                                                          <  2:289382/251‑1 (MQ=255)                                                                                                                   ccGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTG‑AATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGAc                                                                                                                        <  1:315910/251‑1 (MQ=255)                                                                                                                    cgcgTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTG‑ACTTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTTAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTAACCAGCTGAAATTCCCAATAATTAAGTGATATCGaa                                                                                                                        >  2:22558/1‑249 (MQ=255)                                                                                                                        tAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTG‑AATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAgcgc                                                                                                                   >  1:84650/1‑251 (MQ=255)                                                                                                                               aaaCCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTG‑AATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTg                                                                                                            <  1:850668/251‑1 (MQ=255)                                                                                                                               aaaCCTGGTTTAACACGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTG‑AATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTg                                                                                                            <  2:766889/251‑1 (MQ=255)                                                                                                                                  ccTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCTCGGCGTGCTCAGCCGTAGCTAAGCCGAGAACAAAAATACGCGCCAGGTG‑AATTTCCCTCTGGCGCGTAGAGTACGGGAATGGATATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATAGCGACATCGCGTTTTTGCat                                                                                                         <  2:780517/251‑1 (MQ=255)                                                                                                                                     ggTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTG‑AATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCatta                                                                                                       <  1:332828/250‑1 (MQ=255)                                                                                                                                          aaCGGGCATTTAACGTAAACTCCATCGCCGCGGCGTGCTTAGCCGTAGCTAAGCCGCGAAAAAAAATACGCGCCAGGTG‑AATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTGTTTGCATTATTTTGt                                                                                                 <  2:346025/251‑1 (MQ=255)                                                                                                                                                       cGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTG‑AATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCCTTATTTTGTTACATgcag                                                                                        >  2:853206/1‑245 (MQ=255)                                                                                                                                                       cGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTG‑AATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATgcagagat                                                                                    <  1:453537/251‑7 (MQ=255)                                                                                                                                                             cgctATCGCCGCGGCGTGCTCAGCCGTAGCGAAGCCGCGAACAAAAAGACGCGCCAGGTT‑AATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAAtt                                                                               <  2:667460/248‑1 (MQ=255)                                                                                                                                                               ctATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTG‑AATTTCCCTCTGGCGCGTAGCGTACGGGAATGGACATCAATATACTTAAAGGAAATAAGCCTATTCCTGACTATTATTGTTACGTTCTTTTACACCCGCCCGTTAATTAACTCACCAGCTGAAATTTACAATACTTATGTGATATCGCCAGCGCGTTTTTGCATTATTTTGTTACATGCCAAGATGAATTGcc                                                                            >  2:865544/1‑251 (MQ=255)                                                                                                                                                                       gcggcgTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTG‑AATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATc                                                                                                                           >  2:246892/1‑196 (MQ=255)                                                                                                                                                                       gcggcgTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTG‑AATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATc                                                                                                                           <  1:246892/196‑1 (MQ=255)                                                                                                                                                                         ggcgTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTG‑AATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTaaca                                                                   >  1:169844/1‑250 (MQ=255)                                                                                                                                                                         ggcgTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGAG‑AATTTCCATCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGCTACGTTCTTTTAAACCCGCCCGTTAATTAAATCACCAGCTGTAATTCACAATAATTAAGTGATATCGAAAGCTCGTTTTTGCATTATTTTGTTACATACAGAGATGAATTgacgcttaaact                                                                   >  2:118913/1‑239 (MQ=255)                                                                                                                                                                                 cAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTG‑AATGTCCCTCTGGCGCGTAGAGTACGTGAATGGTCATCAATATCCTTAAAGTAAATAAGACTATTCCTGTCTATTATGGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGGGATATCGGCAGCGCGTTTTTGCAGTATGTTGTTACGTGCAGAGATGAATTGCCGATTGAACAAGCACtt                                                            <  2:846675/249‑1 (MQ=255)                                                                                                                                                                                                 gcgAACAAAAATACGCGCCAGGTG‑AATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCa                                           >  1:64366/1‑250 (MQ=255)                                                                                                                                                                                                   gAACAAAAATACGCGCCAGGTG‑AATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAAc                                        >  2:193337/1‑251 (MQ=255)                                                                                                                                                                                                             aCGCGCCAGGTG‑AATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCa                               <  1:865540/250‑1 (MQ=255)                                                                                                                                                                                                                     ggTG‑AATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGcc                      <  1:786963/251‑1 (MQ=255)                                                                                                                                                                                                                               cccTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCACACAGAc              <  1:603214/250‑1 (MQ=255)                                                                                                                                                                                                                                          gTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCACACAGACGCAACAGGGACt  <  1:818029/251‑1 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                              TTGGCGATAACGAACTGCTCGACTTCTTACAGTGTCCGGTGCTG‑GCAGTGGCTCAACGTCACCGCGATGTCGCAGAGCGTGTGCTGGAGATTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTG‑AATTTCCCTCTGGCGCGTAGAGTACGGGACTGGACATCAATATGCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAATGCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCGGCGATGAATTGCCGATTTAACAAACACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCACACAGACGTAACAAGGACT  >  NC_000913/88845‑89326

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