Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,041,451	G→A	A47T (GCA→ACA)	ygfZ →	folate‑binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,041,451	0	G	A	95.3%	133.2 / ‑1.7	44	A47T (GCA→ACA)	ygfZ	folate‑binding protein
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (23/18);  total (25/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.95e-01

TCGCCTCGCAGAAGGCACAGATCGCGTAGTGAAAGCACTTGCATTGCCTCATAGCTCTTTTACCATTAGTCATTAATACGCCGTTAAGCAACTCAGGACTCTATTATGGCTTTTACACCTTTTCCTCCCCGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCTTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCAGATGGCAGAAGATCAGCACCTGCTCGCCGCCCATTGCGACGCCAAAGGTAAAATGTGG‑AGCAATTTACGTCTGTTCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGC  >  NC_000913/3041208‑3041688                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                               tCGCCTCGCAGAAGGCACAGATCGCGTAGTGAAAGCACTTGCATTGCCTCATAGCTCTTTTACCATTAGTCATTAATACGCCGTTAAGCAACTCAGGACTCTATTATGGCTTTTACACCTTTTCCTCCCCGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCTTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGTCAACAGATGt                                                                                                                                                                                                                                         >  2:12704/1‑251 (MQ=255)        gCAGAAGGCACAGATCGCGTAGTGAAAGCACTTGCATTGCCTCATAGCTCTTTTACCATTAGTCATTAATACGCCGTTAAGCAACTCAGGACTCTATTATGGCTTTTACACCTTTTCCTCCCCGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCTTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAACAGATGTCAGcc                                                                                                                                                                                                                                    >  2:328640/1‑249 (MQ=255)                 acaGATCGCGTAGTGAAAGCACTTGCATTGCTTACTAGCTCTTTTACAATTAGAAATTACTACGCCGATAAGCAACTCAGGACTCTATTATGGCTTTTACACCTTTTCCTCCCCGTCAACCTACGGCTTATGCCCGTCTACCGCTGACGCTGATGACGCTTGATGACTGGGGGCTTGCCACCATTACTGGCGCGGACAGCGAAACATATATACAGGGTAAGGTGACAACAGATGTCAGCCAGATGGCaga                                                                                                                                                                                                                          >  1:145694/1‑250 (MQ=255)                 acaGATCGCGTAGTGAAAGCACTTGCATTGCCTCATAGCTCTTTTACCATTAGTCATTAATACGCCGTTAAGCAACTCAGGACTCTATTATGGCTTTTACACCTTTTCCTCCCCGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCTTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAACAGATGTCAGCCAGATGGCagaa                                                                                                                                                                                                                         <  1:238482/251‑1 (MQ=255)                                    cACTTGCATTGCCTCATAGCTCTTTTACCATTAGTCATTAATACGCCGTTAAGCAACTCAGGACTCTATTATGGCTTTTACACCTTTTCCTCCCCGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCTTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAACAGATGTCAGCCAGATGGCAGAAGATCAGCACCTGCTcgccg                                                                                                                                                                                                      >  1:246114/1‑251 (MQ=255)                                          cATTGCCTCATCGCTCTTTTACCATTAGTCATTGATACGCCGTTAAGCAACTCAGGACTCTATTATGGCTTTTACACCTTTTCCTCCCCGTCAGACTACGGCTTCTGCCCGTTTACCGTTTCCGCTGCTGACGCTTGGTGACTGGGCGCTTGCCACCTTTACTGCAGCGGACAGCGAAAATTATATGCAGGGTCAGGTGACACCCGCTGTCACCCAGATGGCAGTATATCAGCACCTGCTCGCCGcccctt                                                                                                                                                                                                >  2:292543/1‑248 (MQ=255)                                               ccTCATTGCGCTTTTACCATTAGTCAGTAATACGCCGTTAAGCAACTCAGGACTCTATTCTGGCTTTTACACCTTTTCCGCCCCGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACGGGGCGCTTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAACAGATGTCAGCCAGATGGCAGAAGATCAGCACCTGCTCGCCGCCCATTGCGAc                                                                                                                                                                                           <  2:280322/251‑1 (MQ=255)                                                  cATAGCTCTTTTACCATTAGTCATTAATACGCCGTTAAGCAACTCAGGACTCTATTATGGCTTTTACACCTTTTCCTCCCCGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCTTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAACAGATGTCAGCCAGATGGCAGAAGATCAGCACCTGCTCGCCGCCCATTGCGACg                                                                                                                                                                                          <  1:292543/249‑1 (MQ=255)                                                         ctTTTACCATTAGTCATTAATACGCCGTTAAGCAACTCAGGCCTCTATTATGGCTTTTACACCTTTTCCTCCCCGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCTTGCCACCATTACTGGCGCGGACAGCGAAACATATATGCAGGGTCAGGTGACAACAGATGTCAGCCAGATGGCAGCAGATCAGCACCTGCTCGCCGCCCATTGCGACGCCAAAGGt                                                                                                                                                                                  >  1:127867/1‑250 (MQ=255)                                                         ctTTTACCATTAGTCATTAATACGCCGTTAAGCAACTCAGGACTCTATTATGGCTTTTACACCTTTTCCTCCCCGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCTTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAACAGATGTCAGCCAGATGGCAGAAGATCAGCACCTGCTCGCCGCCCATTGCGACGCCAAAGGt                                                                                                                                                                                  >  2:89455/1‑250 (MQ=255)                                                           tttACCATTAGTCATTAATACGCCGTTAAGCAACTCAGGCCTCTATTATGGCTTTTACACCTTTTCCTCCCCGTCAGCCTACGGCTTCTGCCCGTTTCCCGCTGACGCTGCTGACGCTTGTTGACTGGGCGCTTGCCACCATTACTGGCGCGGACAGCGAAAACTATATGCAGGGTCAGGTGACAACAGATGTCAGACAGATGGCAGAAACTCAGCACCTGCTCTCCGCACATTCGGACGCCAAAGGTaa                                                                                                                                                                                >  1:190676/1‑250 (MQ=255)                                                           tttACCATTAGTCATTAATACGCCGTTAAGCAACTCAGGACTCTATTATGGCTTTTACACCTTTTCCTCCCCGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCTTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAACAGATGTCAGCCAGATGGCAGAAGATCAGCACCTGCTCGCCGCCCATTGCGACGCCAAAGGTaa                                                                                                                                                                                >  1:62090/1‑250 (MQ=255)                                                                 aTTAGGCATTAATACGCCGTTAAGCAACTCATGACTCTATTATGGCTTTTACACCTTTTCCTCCCCGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCTTGCCTCCATTACTGGCGCTGACAGCGTAAAATATAGGCAGGGTCAGGTGACAACAGATGTCAGTGAGTTGGCAGAAGATCAGCACCTGCTCGCCGCCCATTGCGACGCCAAAGGTAAAATGTgg                                                                                                                                                                         <  2:361998/251‑1 (MQ=255)                                                                     gTCATTAATACGCCGTTAAGCAACTCAGGACTCTATTATGGCTTTTACACCTTTTCCTCCCCGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCTTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAACAGATGTCAGCCAGATGGCAGAAGATCAGCACCTGCTCGCCGCCCATTGCGACGCCAAAGGTAAAATGTGG‑AGc                                                                                                                                                                     >  1:287425/1‑250 (MQ=255)                                                                               cGCCGTTAAGCAACTCAGGACTCTATTATGGCTTTTACACCTTTTCCTCCCCGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCTTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAACAGATGTCAGCCAGATGGCAGAAGATCAGAACCTGCTCGCCGCCCATTGCGACGCCAAAGGTAAAATGTGG‑AGCAATTTACGTCt                                                                                                                                                          >  1:272639/1‑251 (MQ=255)                                                                                  cGTTAAGCAACTCAGGACTCTATTATTGCTTTTACACCTTTTCCTCCCCTTCAGCCTACGTCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCTTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAACAGATGTCAGCCAGATGGCAGAAGATCAGCACCTGCTCGCCGCCCAGTGCGTCTCCAAAGGTAAAATGTGG‑AGCAATTTACGTCTGtt                                                                                                                                                       <  1:340983/251‑1 (MQ=255)                                                                                            cTCAGGACTCTATTATGGCTTTTACACCTTTTCCTCCCCGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCTTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAACAGATGTCAGCCAGATGGCAGAAGATCAGCACCTGCTCGCCGCCCATTGCGACGCCAAAGGTAAAATGTGG‑AGCAATTTACGTCTGTTCCGCGACGGc                                                                                                                                             >  2:354605/1‑251 (MQ=255)                                                                                            cTCAGGACTCTATTATGGCTTTTACACCTTTTCCTCCCCGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCTTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAACAGATGTCAGCCAGATGGCAGAAGATCAGCACCTGCTCGCCGCCCATTGCGACGCCAAAGGTAAAATGTGG‑AGCAATTTACGTCTGTTCCGCGACGGc                                                                                                                                             >  1:38033/1‑251 (MQ=255)                                                                                                  aCTCTATTATGGCTTTTACCCCTTTTCCTCCCCGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCTTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAACAGATGTCAGCCAGATGGCAGAAGATCAGCACCTGCTCGCCGCCCATTGCGACGCCAAAGGTAAAATGTGG‑AGCAATTTACGTCTGTTCCGCGACGGCGATGGc                                                                                                                                       <  2:38033/251‑1 (MQ=255)                                                                                                             gCTTTTACACCTTTTCCTCCCCGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCTTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAACAGATGTCAGCCAGATGGCAGAAGCTCAGCACCTGCTCGCCGCCCATTGCGACGCCAAAGGGAAAATGTGG‑AGCAATTTACGTCTGTTCCGCGACGGCGATGGCTTTGCATGGAt                                                                                                                            >  2:221509/1‑251 (MQ=255)                                                                                                                tttACACCTTTTCCTCCCCGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCTTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAACAGATGTCAGCCAGATGGCAGAAGATCAGCACCTGCTCGCCGCCCATTGCGACGCCAAAGGTAAAATGTGG‑AGCAATTTACGTCTGTTCCGCGACGGCGATGGCTTTGCATGGATTg                                                                                                                          <  1:328640/250‑1 (MQ=255)                                                                                                                 ttACACCTTTTCCTCCCCGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCTTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAACAGATGTCAGCCAGATGGCAGAAGATCAGCACCTGCGCGCCGCCCATTGCGACGCCAAAGGTAAAATGTGG‑AGCAATTTACGTCTGTTCCGCGACGGCGATGGATTTGCATGGATTGa                                                                                                                         >  2:307281/1‑250 (MQ=255)                                                                                                                   accccTTTTCCTCCCCGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCTTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAACAGATGTCAGCCAGATGGCAGAAGATCAGCACCTGCTCGCCGCCCATTGCGACGCCAAAGGTAAAATGTGG‑AGCAATTTACGTCTGTTCCGCGACGGCGATGGCTTTGCATGGATTGAACg                                                                                                                      <  2:246114/248‑1 (MQ=255)                                                                                                                          ttCCTCCCCGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCTTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAACAGATGTCAGCCAGATGGCAGAAGATCAGCACCTGCTCGCCGCCCATTGCGACGCCAAAGGTAAAATGTGG‑AGCAATTTACGTCTGTTCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAgc                                                                                                               >  1:349571/1‑251 (MQ=255)                                                                                                                           tcctccCCGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCTCGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAACAGATGTCAGCCAGATGGCAGAAGATCAGCACCTGCTCGCCGCCCATTGCGACGCCAAAGGTAAAATGTGG‑AGCAATTTACGTCTGTTCCGCGACGGCGATGGCTTTGCCTGGATTGAACGGCGCAgcg                                                                                                              >  1:258421/1‑251 (MQ=255)                                                                                                                                   gTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCTTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCAGATGGCAGAAGATCAGCACCTGCTCGCCGCCCATTGCGACGCCAAAGGTAAACTGTGG‑AGCAATTTACGTCTGTTCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGa                                                                                                       >  2:110727/1‑250 (MQ=255)                                                                                                                                    tCAGCCTACGGCTTCTGCCCGTTTACAGCTGACTCTGATGACGCTTGATGACTGGGCGCTTGCCACCATTACTGGCGAGGACAGCGACCAATATATGCAGGGTCAGGAGACAACAGATGGCAGCCAGATGGCAGAAGATCAGCACCTGCTCGCCGCCCATTGCGACGCCAAAGGTAAAATGTGG‑AGCAATTTCCGTCTGTTCCGCGACGGCGATGGCTTTGCATGGATTGTACGGCGCAGCTCGCGTGAAc                                                                                                     >  2:1178/1‑251 (MQ=255)                                                                                                                                     cAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCTTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAACAGATGTCAGCCAGATGGCAGAAGATCAGCACCTGCTCGCCGCCCATTGCGACGCCAAAGGTAAAATGTGG‑AGCAATTTACGTCTGTTCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAAcc                                                                                                    <  2:62090/251‑1 (MQ=255)                                                                                                                                           aCGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCTTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAACAGATGTCAGCCAGATGGCAGAAGATCAGCACCTGCTCGCCGCCCATTGCGACGCCAAAGGTAAAATGTGG‑AGCAAtt                                                                                                                                                                 >  2:251633/1‑184 (MQ=255)                                                                                                                                           aCGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCTTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAACAGATGTCAGCCAGATGGCAGAAGATCAGCACCTGCTCGCCGCCCATTGCGACGCCAAAGGTAAAATGTGG‑AGCAAtt                                                                                                                                                                 <  1:251633/184‑1 (MQ=255)                                                                                                                                             ggCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCTTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAACAGATGTCAGCCAGATGGCAGAAGATCAGCACCTGCTCGCCGCCCATTGCGACGCCAAAGGTAAAATGTGG‑AGCAATTTACGTCTGTTCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGctg                                                                                             >  1:295076/1‑250 (MQ=255)                                                                                                                                             ggCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCTTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAACAGATGTCAGCCAGATGGCAGAAGATCAGCACCTGCTCGCCGCCCATTGCGACGCCAAAGGTAAAATGTGG‑AGCAATTTACGTCTGTTCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGctg                                                                                             >  2:258527/1‑250 (MQ=255)                                                                                                                                                    ggccGTTTGCCGCTTACGCTGATGACGCGTGATGACTGGGCGCTTGCCACCATTACTGGCGCGGTCAGCGAAAAATAGATGCCTGGTCACGTGACCACAGCGGCCAGCCAGCGGGCAGAAGTTCAGCACCTGCTCGCCGCCCACTGCGACGCCAAAGGTAAAA‑GTGGGAGCAATTTACGCCTTCTTCGCGACCGCGATGGCTTCGCATGGATAGACCGGCCCAGCGAGCGCGAACCGCAGATGaatgaact                                                                                     <  2:212495/249‑1 (MQ=255)                                                                                                                                                    gCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCTTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAACAGATGTCAGCCAGATGGCAGAAGATCAGCACCTGCTCGCCGCCCATTGCGACGCCAAAGGTAAAATGTGG‑AGCAATTTACGTCTGTTCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGATCCGCAGCTGactgaa                                                                                       <  2:36430/249‑1 (MQ=255)                                                                                                                                                               gctgacgctgaTGACGCTTGATGACTGGGCGCTTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAACAGATGTCAGCCAGATGGCAGAAGATCAGCACCTGCTCGCCGCCCATTGCGACGCCAAAGGTAAAATGTGG‑AGCAATTTACGTCTGTTCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATg                                                                          <  1:291109/251‑1 (MQ=255)                                                                                                                                                                              gCTTGATGACTGGGCGCTTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAACAGATGTCAGCCAGATGGCAGAAGATCAGCACCTGCTCGCCGCCCATTGCGACGCCAAAGGTAAAATGTGG‑AGCAATTTACGTCTGTTCCGCGACGGCGATGGCTTTGCATGGAGTGAACGGCGCAGCGTGCGTGAACGGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAg                                                           <  1:307281/251‑1 (MQ=255)                                                                                                                                                                                     gaCTGGGCGCTTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCAGATGGCAGAAGATCAGCACCTGCTCGCCGCCCATTGCGACGCCAAAGGTAAAATGTGG‑AGCAATTTACGTCTGTTCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGACCa                                                     <  1:110727/250‑1 (MQ=255)                                                                                                                                                                                         gggCGCTTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAACAGATGTCAGCCAGATGGCAGAAGATCAGCACCTGCTCGCCGCCCATTGCGACGCCAAAGGTAAAATGTGG‑AGCAATTTACGTCTGTTCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGACCATcg                                                  <  1:89455/249‑1 (MQ=255)                                                                                                                                                                                                tGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAACAGATGTCAGCCAGATGGCAGAAGATCAGCACCTGCTCGCCGCCCATTGCGACGCCAAAGGTAAAATGTGG‑AGCAATTTACGTCTGTTCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGACCATCGCGCCAgac                                          <  1:106391/250‑1 (MQ=255)                                                                                                                                                                                                    accaTTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAACAGATGTCAGCCAGATGGCAGAAGATCAGCACCTGCTCGCCGCCCATTGCGACGCCAAAGGTAAAATGTGG‑AGCAATTTACGTCTGTTCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGACCATCGCGCCAGacgacg                                      >  1:210264/1‑250 (MQ=255)                                                                                                                                                                                                       aTTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAACAGATGTCAGCCAGATGGCAGAAGATCAGCACCTGCTCGCCGCCCATTGCGACGCCAAAGGTAAAATGTGG‑AGCAATTTACGTCTGTTCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGACCATCGCGCCAGACGACGAGCg                                  >  1:262676/1‑251 (MQ=255)                                                                                                                                                                                                             ggCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAACAGATGTCAGCCAGATGGCAGAAGATCAGCACCTGCTCGCCGCCCAGTGCGACGCCAAAGGTAAAATTTGG‑AGCAATTTACGGCTGTTCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCTTGAACCGCAGCTGACTGAACTGAAAAAATATGCGTTATTCTCTAAAGTGACCATCGCGCCATACGACGAGCGTGtgct                            <  2:321448/251‑1 (MQ=255)                                                                                                                                                                                                                                      aGGGTCAGGTGACAACAGATGTCAGCCAGATGGCAGAAGATCAGCACCTGCTCGCCGCCCATTGCGACGCCAAAGGTAAAATGTGG‑AGCAATTTACGTCTGTTCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGgcg   <  2:85630/251‑1 (MQ=255)                                                                                                                                                                                                                                       gggTCAGGTGACAACAGATGTCAGCCAGATGGCAGAAGATCAGCACCTGCTCGCCGCCCATTGCGACGCCAAAGGTAAAATGTGG‑AGCAATTTACGTCTGTTCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGGGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGgcgc  >  2:360709/1‑251 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                               TCGCCTCGCAGAAGGCACAGATCGCGTAGTGAAAGCACTTGCATTGCCTCATAGCTCTTTTACCATTAGTCATTAATACGCCGTTAAGCAACTCAGGACTCTATTATGGCTTTTACACCTTTTCCTCCCCGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCTTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCAGATGGCAGAAGATCAGCACCTGCTCGCCGCCCATTGCGACGCCAAAGGTAAAATGTGG‑AGCAATTTACGTCTGTTCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGC  >  NC_000913/3041208‑3041688

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