Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	10,840	G→A	L173L (CTG→TTG)  S4N (AGT→AAT)	yaaW ← mbiA →	putative enzyme‑specific chaperone YaaW uncharacterized protein MbiA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	10,840	0	G	A	98.0%	142.4 / NA	49	L173L (CTG→TTG)  S4N (AGT→AAT)	yaaW mbiA	putative enzyme‑specific chaperone YaaW uncharacterized protein MbiA
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (24/24);  minor base C (1/0);  total (25/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTA‑GTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCAAACCATGCCCAAGTACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAGTTGGCTGGAAAGCAAATGCGACACACCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGACCGCGTGCAGAAACTCCTGCTTATGTTCCTCGTCCATTTTCTTCCAGGTATTACGCAGAAATTGTTCCAGTAACTGTTGCTCAATTTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGG  >  NC_000913/10595‑11085                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                      tCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTACGGTGTCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTA‑ATGCTGCGCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCACGCTCATTGCTGCGGGGGTGCGGAAAATGCGGGTCAATTGGc                                                                                                                                                                                                                                                  >  2:86445/1‑251 (MQ=255)                                     cGGTTTACGGTGTCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGACAGGCACTGCCGCTGGCCGCTTTAACCCCATTTA‑ATGCTGCGCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATACCCAAGCACGCCAATTGCCGCGTGGGAGCGGCAACTGCGGGTCAATTGCCTGGAAAGCAAATGCGACACCCCtttacccatttatt                                                                                                                                                                                                              >  2:580999/1‑240 (MQ=255)                                      ggTTTACGGTGTCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTA‑ATGCTGCGCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGCCACGCCTTTTGCCAATAAttt                                                                                                                                                                                                             >  1:477370/1‑251 (MQ=255)                                                          aGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTA‑ATGCTGCGCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGgcc                                                                                                                                                                                         >  1:421635/1‑250 (MQ=255)                                                                 cTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCTCATTTA‑ATGCTGCGCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCAAGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCcgcagcgc                                                                                                                                                                                  >  2:303943/1‑249 (MQ=255)                                                                  tAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCAGACCGCTTTAACCCCATTTA‑ATGCTGCGCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGAGACCCGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGctct                                                                                                                                                                                 >  2:169097/1‑251 (MQ=255)                                                                         ccGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTA‑ATGCTGCGCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCt                                                                                                                                                                          >  1:562704/1‑251 (MQ=255)                                                                           ggCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTA‑ATGCTGCGCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTc                                                                                                                                                                         >  1:242494/1‑250 (MQ=255)                                                                                         aTTTGAAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCGGACCGCTTTAACCCCATTTA‑ATGCTTCGCCCACCGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAATCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCAt                                                                                                                                                           <  2:200148/250‑1 (MQ=255)                                                                                             caaGTACGGCTGGAAGTGTCACGCGATAGGCACGGCCGCTGACCGCTTTAACCCCATCTA‑ATGCTGCGCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCTCGCTCAGTGCTGCGTGGGTGCGTAAAATGCGTGTCAAGGGGCTGGAAAGCAAATGCGACACGCCGTTTGCCAATAATTTGTCTTTCATCAGCCGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGAATCGAc                                                                                                                                                       <  2:242494/248‑1 (MQ=255)                                                                                                    ggCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTA‑ATGCTGCGCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGACCGCGTGCa                                                                                                                                               <  1:93534/251‑1 (MQ=255)                                                                                                           attgTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTA‑ATGCTGCGCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGACCGCGTGCAGAAACTc                                                                                                                                        >  1:654241/4‑251 (MQ=255)                                                                                                              tctccgcgATAGGCACTACCTCTAACCGCTTTAACCACATTTACATCTCGC‑CCTCCAGGGCCGCCCATGCCCGCGCCTCACACAAAACCATGCCCAAGCACGCTCATTGCTGCGTGGTTGCTTAAAATGCGGGTCAATTGGCTTGAAAGCAAATGCGACACGCCTTTTGCCAACTATTTGTCTTCCATCAGCAGCGGCAGCAGCTCGTCCAGCGTAGTCACCCTGGCATCGACCGCGTGCAGAAAcgcca                                                                                                                                      <  2:344785/246‑5 (MQ=255)                                                                                                                cACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTA‑ATGCTGCGCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGACCGCGTGCAGAAACTCCTGc                                                                                                                                    <  1:92522/250‑1 (MQ=255)                                                                                                                    cgATAGGCACTGCCGCTGACCGCTTTAACCCCATTTA‑ATGCTGCGCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCCGCTCTTCCAGCTCATTCACCATGGCATCGACCGCGTGCAGAAACTCCTGCTTATg                                                                                                                               >  2:9616/1‑251 (MQ=255)                                                                                                                         ggcacTGCCGCTGACCGCTTTAACCCCATTTA‑ATGCTGCGCCCCCAGCGCCTCCCAGCCCCGCGCCGAGCAGCAAACCATGCCCAAGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCACTTGGCTGGAAAGCAACTGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGACCGCGTGCAGAAACTCCTGCTTAtactccc                                                                                                                          >  2:679355/5‑245 (MQ=255)                                                                                                                             cTGCCGCTGACCGCTTTAACCCCATTTA‑ATGCTGCGCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACAATGCCCAAGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGACCGCGTGCAGAAACTACTGCTACTGTTCCTCGTc                                                                                                                      >  2:447729/1‑251 (MQ=255)                                                                                                                               gCCGCTGACCGCCTTAACCCCATTTA‑ATGCTGCGCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGACCGCGTGCAGAAACTCCTGCTTATGTTCCTCGTCCa                                                                                                                    <  1:17093/251‑1 (MQ=255)                                                                                                                                ccGCTGACCGCTTTAACGCTAATTA‑ATGCGGCGCCCTCAGGGCCTCCCAGACCCGCGCCTCGCAGCAAACCATGCCCAAGCACGCTCATTGCTTCGTGGGTGCGTAAAATGCGGGTCAATGGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGACCGCGTGCAGAAACTCCTGCTTATGTTCCTCGTCCAt                                                                                                                   <  2:489802/251‑1 (MQ=255)                                                                                                                                 cGCTGAACGCTTTAACCCCATTTA‑ATGCTGCGCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCACGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGACCGCGTGCAGAAACTCCTGCTTATGTTTCTCGTCCAtt                                                                                                                  >  2:662177/1‑251 (MQ=255)                                                                                                                                  gCTGACCGCTTTAACCCCATTTA‑ATGCTGCGCCCACAGGGCCTCCCAGACCCGCGCCGCGGAGCAAACCATGCCCAAGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGACCGCGTGCAGAAACTCCTGCTTATGTTCCTCGTCCAtat                                                                                                                 <  1:575848/251‑3 (MQ=255)                                                                                                                                           tttAACCCCATTTA‑ATGCTGCGCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGACCGCGTGCAGAAACTCCTGCTTATGTTCCTCGTCCATATTCTTCCAgg                                                                                                        >  1:278852/1‑251 (MQ=255)                                                                                                                                                ccccATTTA‑ATGCTGCGCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGACCGCGTGCAGAAACTCCTGCTTATGTTCCTCGTCCATATTCTTCCAGGTGTTa                                                                                                   <  1:594424/251‑1 (MQ=255)                                                                                                                                                 cccATTTA‑ATGCTGCGCCCACAGGGCATCCCAGACCCGCGCCGCGCAGCAAACCACGCCCAAGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCCTCGCCCGCGCGCAGACACTCCAGCTTATGTTCCTCGTCCATATTCTTCCAGGTGTTAc                                                                                                  >  2:366859/1‑251 (MQ=255)                                                                                                                                                                                cagacccGCGCCGCGCAGCAAACCATGCCCAAGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGACCGCGTGCAGAAACTCCTGCTTATGTTCCTCGTCCATATTCTTCCAGGTGTTACGCAGAAATTGTTCCAGTAACTGTTGCTCaa                                                                    <  1:589966/247‑1 (MQ=255)                                                                                                                                                                                   ccccGCGCCGCGCAGCAAACCATGCCCAAGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGACCGCGTGCAGAAACTCCTGCTTATGTTCCTCGTCCATATTCTTCCAGGTGTTACGCAGAAATTGTTCCAGTAACTGTTGCTCAAtct                                                                 <  1:366859/251‑3 (MQ=255)                                                                                                                                                                                       gcgccgcgcAGCAAACCATGCCCAAGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGACCGCGTGCAGAAACTCCTGCTTATGTTCCTCGTCCATATTCTTCCAGGTGTTACGCAGAAATTGTTCCAGTAACTGTTGCTCAATCTCaaa                                                             <  1:86445/251‑1 (MQ=255)                                                                                                                                                                                           cgcgcAGCAAACCATGCCCAAGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGACCGCGTGCAGAAACTCCTGCTTATGTTCCTCGTCCATATTCTTCCAGGTGTTACGCAGAAATTGTTCCAGTAACTGTTGCTCAATCTCAAACGt                                                          >  1:298110/1‑250 (MQ=255)                                                                                                                                                                                               cagcaAACCATGCCCAAGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGCCCGCGTGCAGAAACTCCTGCTTATGTTCCTCGTCCACATTCTTCCAGGTGTTACGCAGAAATTGTTCCAGTAACTGCTGCTCAATATCAAACGTAGAc                                                      >  1:66848/1‑250 (MQ=255)                                                                                                                                                                                               cagcaAACCATGCCCAAGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGACCGCGTGCAGAAACTCCTGCTTATGTTCCTCGTCCATATTCTTCCAGGGGTTACGCAGCACTTGTTCCAGTAACTGTTGCTCAATCTCAAACTTAGAc                                                      >  2:342096/1‑250 (MQ=255)                                                                                                                                                                                                   aaaCCATGCCCAAGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGCCCGCGTGCAGAAACTCCTGCTTATGTTCCTCGTCCATATTCTTCCAGGTGTTACGCAGAAATTGTTCCAGTAACTGTTGCTCAATCTCAAACGTAGACAtctc                                                 >  2:264673/1‑251 (MQ=255)                                                                                                                                                                                                          gCCCAAGCACGCTCATTGCTGCGTGGGTGCGTCAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGACCGCGTGCAGCAACTCCTGCTTATGTTCCTCGTCCATATTCTTCCAGGTGTTACGCAGAAATTGTTCCAGTAACTGTTGCTCAATCTCAAACGTAGACATATCTTTGTCg                                          >  2:60940/1‑251 (MQ=255)                                                                                                                                                                                                          gCCCAAGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGACCGCGTGCAGAAACTCCTGCTTATGTTCCTCGTCCATATTCTTCCAGGTGTTACGCAGAAATTGTTCCAGTAACTGTTGCTCAATCTCAAACGTAGACATCTCTTTGTCg                                          <  2:421635/251‑1 (MQ=255)                                                                                                                                                                                                               agcaCGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGACCGCGTGCAGAAACTCCTGCTTATGTTCCTCGTCCATATTCTTCCAGGtg                                                                                                      <  2:208257/183‑2 (MQ=255)                                                                                                                                                                                                               agcaCGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGACCGCGTGCAGAAACTCCTGCTTATGTTCCTCGTCCATATTCTTCCAGGtg                                                                                                      >  1:208257/4‑185 (MQ=255)                                                                                                                                                                                                                 caCGCTCATTGCTGCGTGGGGGCTTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGAAACGCCTTTTGCCAATAATTTGGCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGACCGCGTGCAGAAACTCCTGCTTATGTTCCTCGTCCATATTCTTCCAGGTGTTACGAAGAAATTGTTCCAGTAACTGTTGCTCAATCTCAAACGTAGACATCTCTTTGTCGTCTTTCa                                   <  2:477370/250‑1 (MQ=255)                                                                                                                                                                                                                    gCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGACCGCGTGCAGAAACTCCTGCTTATGTTCCTCGTCCATATTCTTCCAGGTGTTACGCAGAAATTGTTCCAGTAACTGTTGCTCAATCTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCt                                >  2:218760/1‑251 (MQ=255)                                                                                                                                                                                                                      tCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGACCGCGTGCAGAAACTCCTGCTTATGTTCCTCGTCCATATTCTTCCAGGTGTTACGCCGAAATTGTTCCAGTAACTGTTGCTCAATCTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTc                              >  2:617259/1‑251 (MQ=255)                                                                                                                                                                                                                       cATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGACCGCGTGCAGAAACTCCTGCTTATGTTCCTCGTCCATATTCTTCCAGGTGTTACGCAGAAATTGTTCCAGTAACTGTTGCTCAATCTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCtt                               <  1:202153/249‑1 (MQ=255)                                                                                                                                                                                                                        aTTGCTGCGTGGGTGCGAAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGACCGCGTGCAGAAACTCCTGCTTATGTTCCTCGTCCATATTCTTCCAGGTGTTACGCAGAAATTGTTCCAGTAACTGTTGCTCAATCTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCa                             <  2:400146/250‑1 (MQ=255)                                                                                                                                                                                                                           gctgcGTGGGTGCGTAAAATGCGGGTCAATTGTCTGGAAAGCAAAGGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGACCGCGTGCAGAATCTCCTGCTTATGTTCCTCGTCCATATTCTTCCAGGTGTTACGCAGAAATTGTTCCAGTAACGGTTGCTCAATCTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCg                         <  1:601720/251‑1 (MQ=255)                                                                                                                                                                                                                            ctgcGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGACCGCGTGCAGCAACTCCTGCTTATGTTCCTCGTCCATATTCTTCCAGGTGTTACGCAGAAATTGTTCCAGTAACTGTTGCTCAATATCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCgt                        >  1:323740/1‑250 (MQ=255)                                                                                                                                                                                                                            ctgcGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGACCGCGTGCAGAAACTCCTGCTTATGTTCCTCGTCCATATTCTTCCAGGTGTTACGCAGAAATTGTTCCAGTAACTGTTGCTCAATCTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCgt                        >  1:456523/1‑250 (MQ=255)                                                                                                                                                                                                                             gtcgggggggtgTAAAATGCGTGTCAAATGGCTGGAGAGCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTTCAACAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGACCGCGTGCAGAAACTCCTGCTTATGTTCCTCGTCCATATTCTTCCAGGTGTTACGCAGAAATTGTTCCAGTAACTGTTGCTCAATCTCAAACGGAGACATCACTTTGTCGGCTTTCAGCTTCAATCgtt                       <  2:351970/240‑3 (MQ=255)                                                                                                                                                                                                                               cGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGACCGCGTGCAGAAACTCCTGCTTATGTTCCTCGTCCATATTCTTCCAGGTGTTACGCAGAAATTGTTCCAGTAACTGTTGCTCAATCTCAAACGTAGACATCTCTGTGTCGGCTTTCAGCTTCAATCGtttt                     <  1:264673/251‑1 (MQ=255)                                                                                                                                                                                                                                  gtgTGCGTAAAATGCGGGACAAATGTCTGGAAAGCAAATGCGACACGTCTTTTGCCAATTAATTGTTTGTCATCAGCAGCGGCAGCAGCTCTTCCAGGTCATTCACCCTGGCATCGACCGCGTGCAGACACTCCGTCTTATGTTCCTCGTCCTTATTCTTCCAGGTGTTACGCAGAAATTGTTCCAGTAACTGTTGCTCAATCTCAAACTTAGACATCTCTTTTTCGTATTTCAGCTTCAATCGTTTTGaa                  <  2:579250/249‑1 (MQ=255)                                                                                                                                                                                                                                          aaaaTGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAAGTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGACCGCGTGCAGAAACTCCTGCTTATGTTCCTCGTCCATATTCTTCCAGGTGTTACGCAGAAATTGTTCCAGTAACTGTTGCTCAATCTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGTTTTGAAACATCGa           <  1:344835/250‑1 (MQ=255)                                                                                                                                                                                                                                               gCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGACCGCGTGCAGAAACTCCTGCTTATGTTCCTCGTCCATATTCTTCCAGGTGTTACGCAGAAATTGTTCCAGTAACTGTTGCTCAATCTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGTTTTGAAACATCGAGCaaaa     <  1:616718/251‑1 (MQ=255)                                                                                                                                                                                                                                                  ggTCAATTGGCTGGAAAGCAAATTCTACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGAATCTCTTCCAGCTCATGCACCCTGGCATCGACCGCGTTCAGAAACTCCTGCTTATGTTCCTCTTCCATATTCTTCCAGGTGTTACGCAGAAATTGTTCCAGTAACTGTTGCGCAATCTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGTTTTGAAACAGCGAGCAAAATgg  <  2:671098/251‑1 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                      TCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTA‑GTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCAAACCATGCCCAAGTACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAGTTGGCTGGAAAGCAAATGCGACACACCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGACCGCGTGCAGAAACTCCTGCTTATGTTCCTCGTCCATTTTCTTCCAGGTATTACGCAGAAATTGTTCCAGTAACTGTTGCTCAATTTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGG  >  NC_000913/10595‑11085

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: ‑