Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,264,460	2 bp→AT	coding (598‑599/1182 nt)	setB →	sugar exporter SetB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,264,460	0	G	A	100.0%	134.8 / NA	40	E200M (GAG→ATG)	setB	sugar exporter SetB
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (23/17);  total (23/17)
*	NC_000913	2,264,461	0	A	T	95.0%	129.9 / ‑3.7	40	E200M (GAG→ATG)	setB	sugar exporter SetB
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base T (22/16);  minor base G (1/1);  total (23/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.36e-01

CGACCGCTAACCCGCAAATGTTTGCCCTTGCCCGTGAACATGCCGACAAAACCGGACGTGAGGCGGTGATGTTCAGCTCTTTTTTACGCGCTCAGGTTTCACTGGCATGGGTCATTGGCCCACCGCTGGCTTATGCCTTAGCGATGGGTTTCAGCTTTACGGTAATGTATCTGAGCGCAGCGGTAGCGTTTATTGTTTGCGGTGTGATGGTGTGGCTGTTTTTACCGTCGATGCGAAAAGAGCTTCCGCTGGCGACCGGCACGATCGAAGCGCCGCGCCGTAACCGTCGCGATACGCTGCTGCTGTTTGTCATTTGTACATTGATGTGGGGCTCGAACAGCCTGTACATCATCAACATGCCGCTATTTATTATCAACGAACTGCATCTTCCCGAGAAACTGGCCGGTGTGATGATGGGGACCGCCGCCGGGCTGGAAATCCCGACGATGTTGATTGCCGGATATTTCGCC  >  NC_000913/2264221‑2264690                                                                                                                                                                                                                                                ||                                                                                                                                                                                                                                      cGACCGCTAACCCGCAAATGTTTGCCCTTGCCCGTGAACATGCCGACAAAACCGGACGTGAGGCGGTGATGTTCAGCTCTTTTTTACGCGCTCAGGTTTCACTGGCATGGGTCATTGGCCCACCGCTGGCTTATGCCTTAGCGATGGGTTTCAGCTTTACGGTAATGTATCTGAGCGCAGCGGTAGCGTTTATTGTTTGCGGTGTGATGGTGTGGCTGTTTTTACCGTCGATGCGAAAAATGCTTCCGCTg                                                                                                                                                                                                                             >  1:274353/1‑251 (MQ=255)      gCTAACCCGCAAATGTTTGCCCTTGCCCGTGAACATGCCGACAAAACCGGACGTGAGGCGGTGATGTTCAGCTCTTTTTTACGCGCTCAGGTTTCACTGGCATGGGTCATTGGCCCACCGCTGGCTTATGCCTTAGCGATGGGTTTCAGCTTTACGGTAATGTATCTGAGCGCAGCGGTAGCGTTTATTGTTTGCGGTGTGATGGTGTGGCTGTTTTTACCGTCGATGCGAAAAATGCTTCCGCTGGCGAc                                                                                                                                                                                                                        <  1:255000/251‑1 (MQ=255)         aaCCCGCAAATGTTTGCCCTTGCCCGTGAACATGCCGACAAAACCGGACGTGAGGCGGTGATGTTCAGCTCTTTTTTACGCGCTCAGGTTTCACTGGCATGGGTCATTGGCCCACCGCTGGCTTATGCCTTAGCGATGGGTTTCAGCTGTACGGTAATGTATCTGAGCGCAGCGGTAGCGTTTATTGTTTGCGGTGTGATGGTGTGGCTGTTTTTACCGTCGATGCGAAAAATGCTTCCGCTGGCGACCg                                                                                                                                                                                                                      <  2:228963/250‑1 (MQ=255)                aaaTGTTTGCCCTTGCCCGTGAACATGCCGACAAAACCGGACGTGAGGCGGTGATGTTCAGCTCTTTTTTACGCGCTCAGGTTTCACTGGCATGGGTCATTGGCCCACCGCTGGCTTATGCCTTAGCGATGGGTTTCAGCTTTACGGTAATGTATCTGAGCGCAGCGGTAGCGTTTATTGTTTGCGGTGTGATGGTGTGGCTGTTTTTACCGTCGATGCGAAAAATGCTTCCGCTGGCGACCGGCACGATc                                                                                                                                                                                                              <  1:373442/251‑1 (MQ=255)                        gcccttgcccGTGAACATGCCGACAAAACCGGACGTGAGGCGGTGATGTTCAGCTCTTTTTTACGCGCTCAGGTTTCACTGGCATGGGTCATTGGCCCACCGCTGGCTTATGCCTTAGCGCTGGGTTTCAGCTTTACGGTAATGTATCTGAGCGCAGCGGTAGCGTTTATTGTTTGCGGTGTGATGGTGTGGCTGTTTTTACCGTCGATGCGAAAAATGCTTCCGCTGGCGACCGGCACGATCGAAgcgcc                                                                                                                                                                                                      >  1:62677/1‑251 (MQ=255)                             tgcccGTGAACATGCCGACAAAACCGGACGTGAGGCGGTGATGTTCAGCTCTTTTTTACGCGCTCAGGTTTCACTGGCATGGGTCATTGGCCCACCGCTGGCTTATGCCTTAGCGATGGGTTTCAGCTTTACGGTAATGTATCTGAGCGCAGCGGTAGCGTTTATTGTTTGCGGTGTGATGGTGTGGCTGTTTTTACCGTCGATGCGAAAAATGCTTCCGCTGGCGACCGGCACGATCGAAgcgccgcg                                                                                                                                                                                                   <  2:101926/249‑1 (MQ=255)                                       cATGCCGACAAAACCGGACGTGAGGCGGTGATGTTCAGCTCTTTTTTACGCGCTCAGGTTTCACTGGCATGGGTCATTGGCCCACCGCTGGCTTATGCCTTAGCGATGGGTTTCAGCTTTACGGTAATGTATCTGAGCGCAGCGGTAGCGTTTATTGTTTGCGGTGTGATGGTGTGGCTGTTTTTACCGTCGATGCGAAAAATGCTTCCGCTGGCGACCGGCACGATCGAAGCGCCGCGCCGTAACCggcc                                                                                                                                                                                       >  1:127621/1‑248 (MQ=255)                                       cATGCCGACAAAACCGGACGTGAGGCGGTGATGTTCAGCTCTTTTTTACGCGCTCAGGTTTCACTGGCATGGGTCATTGGCCCACCGCTGGCTTATGCCTTAGCGATGGGTTTCAGCTTTACGGTAATGTATCTGAGCGCAGCGGTAGCGTTTATTGTTTGCGGTGTGATGGTGTGGCTGTTTTTACCGTCGATGCGAAAAATGCTTCCGCTGGCGACCGGCACGATCGAAGCGCCGCGCCGTAACCGTcg                                                                                                                                                                                       <  1:352690/251‑1 (MQ=255)                                                  aaCCGGACGTGAGGCGGTGATGTTCAGCTCTTTTTTACGCGCCCAGGTTTCACTGGCATGGGTAATTGGCCCACCGCTGGATTATGCCTTAGCGATGGGGTTAAGCTTTACGGTAATGTATCTGAGCGCAGCGGCAGCGTTTATTGTTTGCGGTGGGATGGTGTGGGTGTTTTTACAGGCGATGCGAAAAAGGCTTACGCGGGCGACCGGCACGAACGAAGCGCCGCGCCGAAACCGTCGCGATACgatgg                                                                                                                                                                            >  2:158644/1‑247 (MQ=255)                                                        aCGTGAGGCGGTGATGTTCAGCTCTTTTTTACGCGCTCAGGTTTCACTGGCATGGGTCATTGGCCCACCGCTTGCTTATGCCTTAGCGATGGGTTTCAGCTTTACGGTAATGTATCTGAGCGCAGCGGTAGCGTTTATTGTTTGCGGTGTGATGGTGTGGCTGTTTTTACCGTCGATGCGAAAAATGCTTCCGCTGGCGACCGGCACGATCGACGCGCCGCGCCGTAACCGTCGCGATACCCGGCTGCTGt                                                                                                                                                                      >  2:374605/1‑251 (MQ=255)                                                         cGTGAGGCGGTGATGTTCAGCTCTTTTTTACGCGCTCAGGTTTCACTGGCATGGGTCATTGGCCCACCGCTGGCTTATGCCTTAGCGATGGGTTTCAGCTTTACGGTAATGTATCTGAGCGCAGCGGTAGCGTTTATTGTTTGCGGTGTGATGGTGTGGCTGTTTTTACCGTCGATGCGAAAAATGCTTCCGCTGGCGACCGGCACGATCGAAGCGCCGCGCCGTAACCGTCGCGATACGCTGCTGCTGtt                                                                                                                                                                     <  2:282014/251‑1 (MQ=255)                                                         cGTGAGGCGGTGATGGTCAGCTCTTTTTTACGCGCTCAGGTTTCACTGGCATGGGTCATTGGCCCACCCCTGGCTTATGCCTTAGCGATGGGTTTCAGCTTTACGGTAATGTATCTGAGCGCAGCGGTAGCGTTTCTTGTTTGCGGTGTGATGGTGTGGCTGGTTTTACCGTCGATGCGAAAAATGCTTCCGCTGGCGACCGGCACGATCGAAGCGCCGCGCCGTAAACGTCGCGATACGCTGCTGCTGtt                                                                                                                                                                     >  2:412210/1‑251 (MQ=255)                                                                cGGTGATGTTCAGCTCTTTTTTACGCGCTCAGGTTTCACTGGCATGGGTCATTGGCCCACCGCTGGCTTATGCCTTAGCGATGGGTTTCAGCTTTACGGTAATGTATCTGAGCGCAGCGGTAGCGTTTATTGTTTGCGGTGTGATGGTGTGGCTGTTTTTACCGTCGATGCGAAAAATGCTTCCGCTGGCGACCGGCACGATCGAAGCGCCGCGCCGTAACCGTCGCGATACGCTGCTGCTGTTTGTCAtt                                                                                                                                                              <  2:25471/251‑1 (MQ=255)                                                                       gTTCAGCTCTTTTTTACGCGCTCAGGTTTCACTGGCATGGGTCATTGGCCCACCGCTGGCTTATGCCTTAGCGATGGGTTTCAGCTTTACGGTAATGTATCTGAGCGCAGCGGTAGCGTTTATTGTTTGCGGTGTGATGGTGTGGCTGTTTTTACCGTCGATGCGAAAAATGCTTCCGCTGGCGACCGGCACGATCGAAGCGCCGCGCCGTAACCGTCGCGATACGCTGCTGCTGTTTGTCATTTGTACAt                                                                                                                                                       <  2:274353/251‑1 (MQ=255)                                                                           aGCTCTTTTTTACGCGCTCAGGTTTCACTGGCATGGGTCATTGGCCCACCGCTGGCTTATGCCTTAGCGATGGGTTTCAGCTTTACGGTAATGTATCTGAGCGCAGCGGTAGCGTTTATTGTTTGCGGTGTGATGGTGTGGCTGTTTTTACCGTCGATGCGAAAAATGCTTCCGCTGGCGACCGGCACGATCGAAGCGCCGCGCCGTAACCGTCGCGGTACGCTGCTGCTGGTTGTCATTTGTTCATTGAt                                                                                                                                                   <  1:115940/251‑1 (MQ=255)                                                                               ctTTTTTACGCGCTCAGGTTTCACTGGCATGGGTCATTGGCCCACCGCTGGCTTATGCCTTAGCGATGGGTTTCAGCTTTACGGTAATGTATCTGAGCGCAGCGGTAGCGTTTATTGTTTGCGGTGTGATGGTGTGGCTGTTTTTACCGTCGATGCGAAAAATGCTTCCGCTGGCGACCGGCACGATCGAAGCGCCGCGCCGTAACCGTCGCGATACGCTGCTGCTGTTTGTCATTTGTACATTGATGTgg                                                                                                                                               >  2:267353/1‑251 (MQ=255)                                                                                      aCGCGCTCAGGTTTCACTGGCATGGGTCATTGGCCCACCGCTGGCTTATGCCTTAGCGATGGGTTTCAGCTTTACGGTAATGTATCTGAGCGCAGCGGTAGCGTTTATTGTTTGCGGTGTGATGGTGTGGCTGTTTTTACCGTCGATGCGAAAAATGCTTCCGCTGGCGACCGGCACGATCGAAGCGCCGCGCCGTAACCGTCGCGATACGCTGCTGCTGTTTGTCATTTGTACATTGATGTGGGGCTCg                                                                                                                                         <  2:62677/250‑1 (MQ=255)                                                                                       cgcgCTCAGGTTTCACTGGCATGGGTCATTGGCCCACCGCTGGCTTATGCCTTAGCGATGGGTTTCAGCTTTACGGTAATGTATCTGAGCGCAGCGGTAGCGTTTATTGTTTGCGGTGTGATGGTGTGGCTGTTTTTACCGTCGATGCGAAAAATGCTTCCGCTGGCGACCGGCACGATCGAAGCGCCGCGCCGTAACCGTCGCGATACGCTGCTGCTGTTTGTCATTTGTACATTGATGTGGGGCTCGaa                                                                                                                                       <  1:180703/251‑1 (MQ=255)                                                                                              aGGTTTCACTGGCATGGGTCATTGGCCCACCGCTGGCTTATGCCTTAGCGATGGGTTTCAGCTTTACGGTAATGTATCTGAGCGCAGCGGTAGCGTTTATTGTTTGCGGTGTGATGGTGTGGCTGTTTTTACCGTCGATGCGAAAAATGCTTCCGCTGGCGACCGGCACGATCGAAGCGCCGCGCCGTACCCGTCGCGCTACGCTGCTGCTGTTTGTCATTTGTACATTGCTGTGGGGCTCGAACAGCCTg                                                                                                                                >  2:112407/1‑251 (MQ=255)                                                                                                gTTTCACTGGCATGGGTCATTGGCCCACCGCTGGCTTATGCCTTAGCGATGGGTTTCAGCTTTACGGTAATGTATCTGAGCGCAGCGGTAGCGTTTATTGTTTGCGGTGTGATGGTGTGGCTGTTTTTACCGTCGATGCGAAAAATGCTTCCGCTGGCGACCGGCACGATCGAAGCGCCGCGCCGTAACCGTCGCGATACGCTGCTGCTGTTTGTCATTTGTACATTGATGTGGGGCTCGAACAGCCTGt                                                                                                                               >  1:205244/1‑250 (MQ=255)                                                                                                gTTTCACTGGCATGGGTCATTGGCCCACCGCTGGCTTATGCCTTAGCGATGGGTTTCAGCTTTACGGTAATGTATCTGAGCGCAGCGGTAGCGTTTATTGTTTGCGGTGTGATGGTGTGGCTGTTTTTACCGTCGATGCGAAAAATGCTTCCGCTGGCGACCGGCACGATCGAAGCGCCGCGCCGTAACCGTCGCGATACGCTGCTGCTGTTTGTCATTTGTACATTGATGTGGGGCTCGAACAGCCTGt                                                                                                                               >  1:396084/1‑250 (MQ=255)                                                                                                           aTGGGTCATTGGCCCACCGCTGGCTTATGCCTTAGCGATGGGTTTCAGCTTTACGGTAATGTATCTGAGCGCAGCGGTAGCGTTTATTGTTTGCGGTGTGATGGTGTGGCTGTTTTTACCGTCGAGGCGAAAAATGCTTCCGCTGGCGACCGGCACGATCGAAGCGCCGCGCCGTAACCGTCGCGATACGCTGCTGCTGTTTGTCATTTGTACATTGATGTGGGGCTCGAACAGCCTGTACATCATCAACa                                                                                                                   <  1:398248/251‑1 (MQ=255)                                                                                                            tGGGTCATTGGCCCACCGCTGGCTTATGCCTTAGCGATGGGTTTCAGCTTTACGGTAATGTATCTGAGCGCAGCGGTAGCGTTTATTGTTTGCGGTGTGATGGTGTGGCTGTTTTTACCGTCGATGCGAAAAATGCTTCCGCTGGCGACCGGCACGATCGAAGCGCCGCGCCGTAACCGTCGCGATACGCTGCTGCTGTTTGTCATTTGTACATTGATGTGGGGCTCGAACAGCCTGTACATCATCAACAt                                                                                                                  >  1:230868/1‑251 (MQ=255)                                                                                                                             gCTGGCTTATGCCTTAGCGATGGGTTTCAGCTTTACGGTAATGTATCTGAGCGCAGCGGTAGCGTTTATTGTTTGCGGTGTGATGGTGTGGCTGTTTTTACCGTCGATGCGAAAAATGCTTCCGCTGGCGACCGGCACGATCGAAGCGCCGCGCCGTAACCGTCGCGATACGCTGCTGCTGTTTGTCATTTGTACATTGATGTGGGGCTCGAACAGCCTGTACATCATCAACATGCCGCTATTTATTATCa                                                                                                 <  2:328743/251‑1 (MQ=255)                                                                                                                                                aTGGGTTTCAGCTTTACGGTAATGTATTGGAGCGCAGCGGTAGCGTTTATGGGTTGCTGTGTGATGGTGTGGCTGTTTGTACCGTCGATGCGAAAAAGGCTTCCGCTGGCGACCGGCACGATCGAAGCGTCGCTCCGTACCCTGCGCGATACGCTGCTGCGGTTTGTCAGTTGTACAGTGATGTGGTGCTCGAACAGCCTGTACATCATCAACATGCCGCTATTTATTATCAACGAACTGCATCTTccc                                                                                <  1:397424/249‑1 (MQ=255)                                                                                                                                                 tGGGTTTCAGCTTTACGATAATGTATCTGAGCGCAGCGGTAGCGTTTATTGTTTGCGGTGTGATGGTGTGGCTGTTTTTACCGTCGATGCGAAAAATGCTTCCGCTGGCGACCGGCACGATCGAAGCGCCGCGCCGTCACCGCCGCGATACGCTGCTGCTGTTTGTCATTTGTCCATTGATGTGGGGCTCGAACAGCCTGTACATCATCCACATGCCGCTCTTTATTATCAACGAACTGCATCTTCCCg                                                                               >  2:196985/1‑249 (MQ=255)                                                                                                                                                         aGCTTTACGGTAATGTATCTGAGCGCAGCGGTAGCGTTTATTGTTTGCGGTGTGATGGTGTGGCTGTTTTTACCGTCGATGCGAAAAATGCTTCCGCTGGCGACCGGCACGATCGAAGCGCCGCGCCGTAACCGTCGCGATACGCTGCTGCTGTTTGGCATTTGTACATTGATGTGGGGCTCGAACAGCCTGTACATCATCAACATGCCGCTATTTATTATCAACGAACTGCATCTTCCCGAGAAACTGGc                                                                     <  1:290778/251‑1 (MQ=255)                                                                                                                                                                  gTAATGTATCTGAGCGCAGCGGTAGCGTTTATTGTTTGCGGTGTGATGGTGTGGCTGTTTTTACCGTCGATGCGAAAAATGCTTCCGCTGGCGACCGGCACGATCGAAGCGCCGCGCCGTAACCGTCGCGATACgct                                                                                                                                                                              >  2:107674/1‑137 (MQ=255)                                                                                                                                                                  gTAATGTATCTGAGCGCAGCGGTAGCGTTTATTGTTTGCGGTGTGATGGTGTGGCTGTTTTTACCGTCGATGCGAAAAATGCTTCCGCTGGCGACCGGCACGATCGAAGCGCCGCGCCGTAACCGTCGCGATACgct                                                                                                                                                                              <  1:107674/137‑1 (MQ=255)                                                                                                                                                                   tAATGTATCTGAGCGCAGCGGTAGCGTTTATTGTTTGCGGTGTGATGGTGTGGCTGTTTTTACCGTCGATGCGAAAAATGCTTCCGCTGGCGACCGGCACGATCGAAGCGCCGCGCCGTAACCGTCGCGATACGCTGCTGCTGTTTGTCATTTGTACATTGATGTGGGGCTCGAACAGCCTGTACATCATCAACATGCCGCTATTTATTATCAACGAACTGCATCTTCCCGAGAAACTGGCCGGTGtgatg                                                           >  2:251010/1‑251 (MQ=255)                                                                                                                                                                   tAATGTATCTGAGCGCAGCGGTAGCGTTTATTGTTTGCGGTGTGATGGTGTGGCTGTTTTTACCGTCGATGCGAAAAATGCTTCCGCTGGCGACCGGCACGATCGAAGCGCCGCGCCGTAACCGTCGCGATAAGCTGCTGCTGTTTGTCATTTGTACATTGATGTGGGGCTCGAACAGCCTGTACATCATCAACATGCCGCTATTTATTATCAACGAACTGCATCTTCCCGAGAAACTGGCCGGTGtgatg                                                           >  2:379131/1‑251 (MQ=255)                                                                                                                                                                                                        cGGTGTGATGGTGTGGCTGTTTTTACCGTCGATGCGAAAAATGCTTCCGCTGGCGACCGGCACGATCGAAGCGCCGCGCCGTAACCGTCGCGATACGCTGCTGCTGTTTGTCATTTGTACATTGATGTGGGGCTCGAACAGCCTGTACATCATCAACATGCCGCTATTtattat                                                                                                   >  2:126389/1‑174 (MQ=255)                                                                                                                                                                                                        cGGTGTGATGGTGTGGCTGTTTTTACCGTCGATGCGAAAAATGCTTCCGCTGGCGACCGGCACGATCGAAGCGCCGCGCCGTAACCGTCGCGATACGCTGCTGCTGTTTGTCATTTGTACATTGATGTGGGGCTCGAACAGCCTGTACATCATCAACAGGCCGCTAGTtattat                                                                                                   <  1:126389/174‑1 (MQ=255)                                                                                                                                                                                                            gtgATGGTGTGGCTGTTTTTACCGTCGATGCGAAAAATGCTTCCGCTGGCGACCGGCACGATCGAAGCGCCGCGCCGTAACCGTCGCGATACGCTGCTGCTGTTTGTCATTTGTACATTGATGTGGGGCTCGAACAGCCTGTACATCATCAACATGCCGCTATTTATTATCAACGAACTGCATCTTCCCGAGAAACTGGCCGGTGTGATGATGGGGACCGCCGCCGGGCTGGAAATCCCGACGATGTTGAt                  >  1:107878/1‑251 (MQ=255)                                                                                                                                                                                                                       gCTGTTTTTACCGTCGATGCGAAAAATGCTTCCGCTGGCGACCGGCACGATCGAAGCGCCGCGCCGTAACCGTCGCGATACGCTGCTGCTGTTTGTCATTTGTACATTGATGTGGGGCTCGAACAGCCTGTACATCATCAACATGCCGCTATTTATTATCAACGAACTGCATCTTCCCGAGAAACTGGCCGGTGTGATGATGGGGACCGCCGCCGGGCTGGAAATCCCGACGATGTTGATTGCCGGATAtt       >  1:206870/1‑251 (MQ=255)                                                                                                                                                                                                                        cTGTTTTTACCGTCTAAGCGAAAAATGCTTCCGCTGGCGACCGGCACGATCGAAGCGCCGCGCCGTCACCGTCGCGATACGCTGCTGCTGTTTGTCATTTGTACATTGATGTGGGGCTCGAACAGCCTGTCCATCATCAACATGCCGCTATTTATTATCAACGACCTGCATCTTCCCGAGAAACTGGCCGGGGTGATGATGGGGACCGCCGCCGGGCTGGCAATCCCGACGATGTTGTTTGCCGGATAttt      >  2:323401/1‑251 (MQ=255)                                                                                                                                                                                                                        cTGTTTTTACCGTCGATGCGAAAAATGCTTCCGCTGGCGACCGGCACGATCGAAGCGCCGCGCCGTAACCGTCGCGATACGCTGCTGCTGTTTGTCATTTGTACATTGATGTGGGGCTCGAACAGCCTGTACATCATCAACATGCCGCTATTTATTATCAACGAACTGAATCTTCCCGAGAAACTGGCCGGTGTGATGATGGGGACCGCCGCCGGGCTGGAAATCCCGACGATGTTGATTGCCGGATAttt      >  1:224568/1‑251 (MQ=255)                                                                                                                                                                                                                          gTTTTTACCGTCGATGCGAAAAATGCTTCCGCTGGCGACCGGCACGATCGAAGCGCCGCGCCGTAACCGTCGCGATACGCTGCTGCTGTTTGTCATTTGTACATTGATGTGGGGCTCGAACAGCCTGTACATCATCAACATGCCGCTATTTATTATCAACGAACTGCATCTTCCCGAGAAACTGGCCGGTGTGATGATGGGGACCGCCGCCGGGCTGGAAATCCCGACGATGTTGATTGCCGGATATTTCg    >  1:152026/1‑251 (MQ=255)                                                                                                                                                                                                                           tttttACCGTCGATGCGAAAAATGCTTCCGCTGGCGACCGGCACGATCGAAGCGCCGCGCCGTAACCGTCGCGATACGCTGCTGCTGTTTGTCATTTGTACATTGATGTGGGGCTCGAACAGCCTGTACATCATCAACATGCCGCTATTTATTATCAACGAACTGCATCTTCCCGAGAAACTGGCCGGTGTGATGATGGGGACCGCCGCCGGGCTGGAAATCCCGACGATGTTGATTGCCGGATATTTCGc   >  1:376700/1‑251 (MQ=255)                                                                                                                                                                                                                             tttACCGTCGATGCGAAAAATGCTTCCGCTGGCGACCGGCACGATCGAAGCGCCGCGCCGTAACCGTCGCGATACGCTGCTGCTGTTTGTCATTTGTACATTGATGTGGGGCTCGAACAGCCTGTACATCATCAACATGCCGCTATTTATTATCAACGAACTGCATCTTCCCGAGAAACTGGCCGGTGTGATGATGGGGACCGCCGCCGGGCTGGAAATCCCGACGATGTTGATTGCCGGATATTTCGcc  >  1:167949/1‑250 (MQ=255)                                                                                                                                                                                                                                                ||                                                                                                                                                                                                                                      CGACCGCTAACCCGCAAATGTTTGCCCTTGCCCGTGAACATGCCGACAAAACCGGACGTGAGGCGGTGATGTTCAGCTCTTTTTTACGCGCTCAGGTTTCACTGGCATGGGTCATTGGCCCACCGCTGGCTTATGCCTTAGCGATGGGTTTCAGCTTTACGGTAATGTATCTGAGCGCAGCGGTAGCGTTTATTGTTTGCGGTGTGATGGTGTGGCTGTTTTTACCGTCGATGCGAAAAGAGCTTCCGCTGGCGACCGGCACGATCGAAGCGCCGCGCCGTAACCGTCGCGATACGCTGCTGCTGTTTGTCATTTGTACATTGATGTGGGGCTCGAACAGCCTGTACATCATCAACATGCCGCTATTTATTATCAACGAACTGCATCTTCCCGAGAAACTGGCCGGTGTGATGATGGGGACCGCCGCCGGGCTGGAAATCCCGACGATGTTGATTGCCGGATATTTCGCC  >  NC_000913/2264221‑2264690

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