Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,721,136	C→T	L38F (CTT→TTT)	ydhI →	DUF1656 domain‑containing protein YdhI

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,721,136	0	C	T	92.1%	102.5 / 1.5	40	L38F (CTT→TTT)	ydhI	DUF1656 domain‑containing protein YdhI
Reads supporting (aligned to +/- strand):  ref base C (2/1);  new base T (21/14);  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 = 8.02e-01

AAATGTATTGCTGCATTTGGTTACCGTCAGACATATTTTTCAGAAATTGCGCGTAAATTTTTCGCACTTAAAGAATATTTATTAATCTAACGCAATATATTCGGTCGTAAAGGAATCTACTTTGTGAAGTTTATGCTCAATGCAACAGGATTGCCCTTACAAGATCT‑GGTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGCTTGT‑CGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCTGGCAATACTGATTGCGTGGTAACTATGTCAATTAAAACAATTAAGTATTTCTCAACAATCATTGTAGCGGTAGTTGCGGTTCTTGCCGGATGGTGGCTGTGGAATTATTACATGCAATCTCCGTGGACACGCGACGGGAAAAT  >  NC_000913/1720902‑1721382                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                        aaaTGTATTGCTGCATTTGGTTACCGTCAGACATCTTTTTCCGAAATTGCGCGTAAATTTTTCGCACTTAAAGAATATTTCTTAATCTAACGCAATATATTCGGTCGTAAAGGAATCTACTTTGTGACGTTTATGCTCAATGCAACAGGATTGCCCTTACAAGATCT‑GGTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGCTTGT‑CGTACAcctcct                                                                                                                                                                                                                                        >  2:186669/1‑251 (MQ=255) aaaTGTATTGCTGCATTTGGTTACCGTCAGACATATTTTTCAGAAATTGCGCGTAAATTTTTCGCACTTAAAGAATATTTATTAATCTAACGCAATATATTCGGTCGTAAAGGAATCTACTTTGTGAAGTTTATGCTCAATGCAACAGGATTGCCCTTACAAGATCT‑GGTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGCTTGT‑CGTACACCGCCt                                                                                                                                                                                                                                        >  2:59358/1‑251 (MQ=255)                             aGACATATTTTTCAGAAATTGCGCGTAAATTTTTCGCACTTAAAGAATATTTATTAATCTAACGCAATATATTCGGTCGTAAAGGAATCTACTTTGTGAAGTTTATGCTCAATGCAACAGGATTGCCCTTACAAGATCT‑GGTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGTTTGT‑CGTACACCGCCTGCTTCGTGGCTGGTTCTACGCCGGTGAc                                                                                                                                                                                                            >  2:228785/1‑251 (MQ=255)                              gACATATTTTTCAGAAATTGCGCGTAAATTTTTCGCACTTAAAGAATATTTATTAATCTAACGCAATATATTCGGTCGTAAAGGAATCTACTTTGTGAAGTTTATGCTCAATGCAACAGGATTGCCCTTACAAGATCT‑GGTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGTTTGT‑CGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGAc                                                                                                                                                                                                            >  2:114901/1‑250 (MQ=255)                                   atTTTTCAGAAATTGCGCGTAAATTTTTCGCACTTAAAGAATATTTATTAATCTAACGCAATATATTCGGTCGTAAAGGAATCTACTTTGTGAAGTTTATGCTCAATGCAACAGGATTGCCCTTACAAGATCT‑GGTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGTTTGT‑CGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTg                                                                                                                                                                                                       <  1:370289/250‑1 (MQ=255)                                     ttttCAGAAATTGCGCGTAAATTTTTCGCACTTAAAGAATATTTATTAATCTAACGCAATATATTCGGTCGTAAAGGAATCTACTTTGTGAAGTTTATGCTCAATGCAACAGGATTGCCCTTACAAGATCT‑GGTGTTCGGTGCGTCCGTCTTCTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGTTTGT‑CGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGc                                                                                                                                                                                                     <  1:304091/250‑1 (MQ=255)                                           gAAATTGCGCGTAAATTTTTCGCACTTAAAGAATATTTATTAATCTAACGCAATATATTCGGTCGTAAAGGAATCTACTTTGTGACGTTTATGCTCAATGCAACAGGATTGCCCTTACAAGATCT‑GGTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGTTTGT‑CGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCtt                                                                                                                                                                                              >  1:175862/1‑251 (MQ=255)                                             aaTTGCGCGTAAATTTTTCGCACTTAAAGAATATTTATTAATCTAACGCAATATATTCGGTCGTAAAGGAATCTACTTTGTGAAGTTTATGCTCAATGCAACAGGATTGCCCTTACAAGATCT‑GGTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGTTTGT‑CGTACACCGCCTGCTTCGTGGCTGGATCTACGCCTGTGACATCTGGCATCCCTTGt                                                                                                                                                                                            >  1:150107/1‑251 (MQ=255)                                                      tAAATTTTTCGCACTTAAAGAATATTTATTAATCTAACGCAATATATTCGGTCGTAAAGGAATCTACTTTGTGAAGTTTATGCTCAATGCAACAGGATTGCCCTTACAAGATCT‑GGTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGTTTGT‑CGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGAtt                                                                                                                                                                                   >  1:165567/1‑251 (MQ=255)                                                            tttCGCACTTAAAGAATATTTATTAATCTAACGCAATATATTCGGTCGTAAAGGAATCTACTTTGTGAAGTTTATGCTCAATGCAACAGGATTGCCCTTACAAGATCT‑GGTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGTTTGT‑CGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGc                                                                                                                                                                             >  1:406673/1‑251 (MQ=255)                                                               cGCACTTAAAGAATATTTATTAATCTAACGCAATATATTCGGTCGTAAAGGAATCTACTTTGTGAAGTTTATGCTCAATGCAACAGGATTGCCCTTACAAGATCT‑GGTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGTTTGT‑CGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGt                                                                                                                                                                          >  1:285588/1‑251 (MQ=255)                                                                gCACTTAAAGAATATTTATTAATCTAACGCAATATATTCGGTCGTAAAGGAATCTACTTTGTGAAGTTTATGCTCAATGCAACAGGATTGCCCTTACAAGATCT‑GGTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGTTTGT‑CGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGtt                                                                                                                                                                         >  1:356094/1‑251 (MQ=255)                                                                  aCTTAAAGAATATTTATTAATCTAAAGCAATATATTCGGTCGTAAAGGAATCTACTTTGTGAAATTTATGCACAATGCAACAGGATTGCCCTTACAAGATCT‑GGTTTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAATGCATTCGCGTTTGGATTCGTCATCTGGTTTGT‑CGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATATGGCATCCCTTTTTAATGGATTTATTGCTGTTTg                                                                                                                                                                       >  2:68045/1‑251 (MQ=255)                                                                       aaGAAAATTTATTAATCTAACGCAATATATTCGGTCGTAAAGGAATCTACTTTGTGAATTTTATGCTCAATGCAACAGGATTGCCCTTACAAGATCT‑GGTGTTCGGTGCGTCCGTCTACTTTCCTCCGTGTTTCAAAGCATTCGCGTTGGTATTCGTCATCTGGTTGGT‑CGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGAtt                                                                                                                                                                  <  2:47243/251‑1 (MQ=255)                                                                                   taATCTAACGCAGTATATTCGGTCGTAAAGGAATCTACTTTGTGAAGTTTATGCTCAATGCAACAGGATTGCCCTTACAAGATCT‑GGTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGTTTGT‑CGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCtt                                                                                                                                                      >  2:224366/1‑251 (MQ=255)                                                                                    aaTCTAACGCAATATATTCGGTCGTAAAGGAATCTACTTTGTGAAGTTTATGCTCAATGCAACAGGATTGCCCTTACAAGATCT‑GGTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGTTTGT‑CGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTaa                                                                                                                                                                                         >  2:414672/1‑215 (MQ=255)                                                                                    aaTCTAACGCAATATATTCGGTCGTAAAGGAATCTACTTTGTGAAGTTTATGCTCAATGCAACAGGATTGCCCTTACAAGATCT‑GGTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGTTTGT‑CGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTaa                                                                                                                                                                                         <  1:414672/215‑1 (MQ=255)                                                                                                             aacggAATCTACTTTGTGAAGTTTATGCTCAATGCAACAGGATTGCCCTTACAAGATCT‑GGTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGTTTGT‑CGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCGTGCTCTTGCAATACTGATTGCGTGGt                                                                                                                             >  2:278643/4‑250 (MQ=255)                                                                                                             aaaGGAATCTACTTTGTGAAGTTTATGCTCAATGCAACAGGATTGCCCTTACAAGATCT‑GGTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGTTTGT‑CGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCTGGCAATACTGATTGCGTGGt                                                                                                                             >  2:167493/1‑250 (MQ=255)                                                                                                                   atttACTTTGTGAAGTTTATGCTCAATGCAACCTGATTGCCTTTACAAGATCT‑GGTTTTCGGTGCGGCCGTCTACTTTCCTCCGTTTTTCAAAGCATGCGCGTGTGGATTCGTAACCTGGTTTGT‑CGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCTGGCAATACTGATTGCGTGGTAACTATg                                                                                                                      <  2:108697/248‑1 (MQ=255)                                                                                                                       aCTTTGTTAAGTTTATTCTCACTGCAACAGGATTGCCCTGACAAGATCTTGGT‑TTCTGTGCGTCCGTTTCCTTTCTGCCGTTTCTCAAAGCATTCGCGTTTTGATTCGTCTTCGGGTTTTT‑CGTACACCGCCTGCTTCGTGGCTGGAGCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCTATTTGCGTTTGCCTTGCACTGGCAAGACTGATTGCGTGGTTTCTATGTCa                                                                                                                   <  2:294506/250‑1 (MQ=255)                                                                                                                       aCTTTGTGAAGTTTATGCTCAATGCAACAGGATTGCCCTTACAAGATCT‑GGTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGTTTGT‑CGTACACCGACTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCTGGCAATACTGATTGCGTGGTAACTATGTCa                                                                                                                   <  1:35909/250‑1 (MQ=255)                                                                                                                       aCTTTGTGAAGTTTATGCTCAATGCAACAGGATGGCCCTTACAAGATCT‑GGTGTTCGGTGCGTCCGTCTACTTTCCTCCTTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGTTTGT‑CGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCTGGCAATACTGATTGCGTGGTAACTATg                                                                                                                      <  1:298834/247‑1 (MQ=255)                                                                                                                          ttGTGAAGTTTATGCTCAATGCAACAGGATTGCCCTTACAAGATCT‑GGTGTTCGGTGCGTCCGTCTACTTTCCTGCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGTTTGT‑CGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCTGGCAATACTGATTGCGTGGTAACTATGTCAAt                                                                                                                 <  1:167493/249‑1 (MQ=255)                                                                                                                                         tCAATGCAACAGGATTGCCCTTACAAGATCT‑GGTGTTCGGTGCGTCCGTCTTCTTTCCCCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGTTTGT‑CGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCTGGCAATACTGATTGCGTGGt                                                                                                                             <  2:97562/222‑1 (MQ=255)                                                                                                                                         tCAATGCAACAGGATTGCCCTTACAAGATCT‑GGTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGTTTGT‑CGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCTGGCAATACTGATTGCGTGGt                                                                                                                             >  1:97562/1‑222 (MQ=255)                                                                                                                                         tCAATGCAACAGGATTGCCCTTACAAGATCT‑GGTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATACGCGTTTGGATTCGTCATCTGGTTTGT‑CGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGACATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCTGGCAATACTGATTGCGAGGTAACTATGTCAATTAAAACAATTAAGTAtt                                                                                                >  1:79548/1‑251 (MQ=255)                                                                                                                                          cAATGCAACAGGATTGCCCTTACAAGATCT‑GGTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGTTTGT‑CGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCTGGCAATACTGATTGCGTGGTAACTATGTCAATTAAAACAATTAAGTAttt                                                                                               >  2:225756/1‑251 (MQ=255)                                                                                                                                                                                      cccgtcTAGTTTCCTCCGTTGTTCAAAGCATTCGCGTTTGGAGTCGTCATCTGGTTGGT‑CGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCTGGCAATACTGATTGCGTGGTAACTATGTCAATTAAAACAATTAAGTATTTCTCAACAATCATTGTAGCGGTAGTTGCGGTTCTTGCCGGAtgg                                                    <  2:165567/250‑1 (MQ=255)                                                                                                                                                                                             actgtcctccGTTTTTCAATGCATTCGCGTTTGTATTCGTCATCTGGTTTGT‑CGTCCACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCTGGCAATACTGATTGCGGGGTAACGATGTCAATTAAAACAATTAAGTATTTCTCAACAATCATTGTAGCTGTAGTGGCGGTTCTTGCTGGATGGTGGCtg                                              <  1:107569/246‑1 (MQ=255)                                                                                                                                                                                                       gTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGTTTGT‑CGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCTGGCAATACTGATTGCGTGGTAACTATGTCAATTAAAACAATTAAGTATTTCTCAACAATCATTGTAGCGGTAGTTGCGGTTCTTGCCGGATGGTGGCTGTGGAattatt                                    >  1:324505/1‑250 (MQ=255)                                                                                                                                                                                                           ttCAAAGCATTCGCGTTTGGATTCGTCATCTGGTTTGT‑CGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCTGGCAATACTGATTGCGTGGTAACTATGTCAATTAAAACAATTAAGTATTTCTCAACAATCATTGTAGCGGTAGTTGCGGTTCTTGCCGGATGGTGGCTGTGGAATTATTACATg                               >  1:353755/1‑251 (MQ=255)                                                                                                                                                                                                                     gcgcgTTTGGATTCGTTATCTGG‑TTGTTCTTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCTGGCAATACTGATTGCGTGGTAACTATGTCAATTAAAACAATTAAGTATTTCTCAACAATCATTGTAGCGGTAGTTGCGGTTCTTGCCGGATGGTGGCTGTGGAATTATTACATGCAATCTCCg                      <  2:11791/249‑1 (MQ=255)                                                                                                                                                                                                                           ttGGATTCGTCATCTGGTTTGT‑CGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTTACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCTGGCAATACTGATTGCGTGGTAACTATGTCACTTAAAACAATTAAGTATTTCTCAACAATCATTGTAGCGGAAGTTGCGGATCTTACCGGTTGGGGGATGTGGGATTTTTACATGACATCTCCGTGGaacc               >  2:7920/1‑248 (MQ=255)                                                                                                                                                                                                                            tGGATTCTTCATCTGGTTTGT‑CGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCTGGCAATACTGATTGCGTGGTAACTATGTCAATTAAAACAATTAAGTATTTCTCAACAATCATTGTAGCGGTAGTTGCGGTTCTTGCCGGATGGTGGCTGTGGAATTATTACATGCAATCTCCGTGGacac               <  1:204565/250‑1 (MQ=255)                                                                                                                                                                                                                            tGGATTCGTCATCTGGGTTGT‑CGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGAGTTGCGTTTGCCTTGCTCTGGCAATACTGATTGCGTGGTAACTATGTCAATTAAAACAATTAAGTATTTCTCAACAATCATTGTAGCGGTAGTTGCGGTTCTTGCCGGATGGTGGCTGTGGAATTATTACATGCAATCTCCGTGGacac               <  1:186669/250‑1 (MQ=255)                                                                                                                                                                                                                            tGGATTCGTCATCTGGCTTGT‑CGTCCACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCTGGCAATACTGATTGCGTGGTAACTATGTCAATTAAAACAATTAAGTATTTCTCAACAATCATTGTAGCGGTAGTTGCGGTTCTTGCCGGATGGTGGCTGTGGAATTATTACATGCAATCTCCGTGGacac               <  1:59358/250‑1 (MQ=255)                                                                                                                                                                                                                              gATTCGTCATCTGGTTTGT‑CGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGGTTGCCTTGCTCTGGCAATACTGATTGCGTGGTAACTATGTCAATTAAAACAATTAAGTATTTCTCAACAATCATTGTAGCGGTAGTTGCGGTTCTTGCCGGATGGTGGCTGTGGAATTATTACATGCAATCTCCGTGGACAcgcg            <  1:274636/251‑1 (MQ=255)                                                                                                                                                                                                                                   gTCATCTGGTTTGT‑CGTACACCGCTTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCTGGCAATACTGATTGCGTGGTAACTATGTCAATTAAAACAATTAAGTATTTCTCAACAATCATTGTAGCGGTAGTTGCGGTTCTTGCCGGATGGTGGCTGTGGAATTATTACATGCAATCTCCGTGGACACGCGACggg       >  1:37000/1‑251 (MQ=255)                                                                                                                                                                                                                                        cTGGTTTGT‑CGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCTGGCAATACTGATTGCGTGGTAACTATGTCAATTAAAACAATTAAGTATTTCTCAACAATCATTGTAGCGGTAGTTGCGGTTCTTGCCGGATGGTGGCTGTGGAATTATTACCTGCAATCTCCGTGGACACACGACGGGAAAAt  >  1:154506/1‑251 (MQ=255)                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                        AAATGTATTGCTGCATTTGGTTACCGTCAGACATATTTTTCAGAAATTGCGCGTAAATTTTTCGCACTTAAAGAATATTTATTAATCTAACGCAATATATTCGGTCGTAAAGGAATCTACTTTGTGAAGTTTATGCTCAATGCAACAGGATTGCCCTTACAAGATCT‑GGTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGCTTGT‑CGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCTGGCAATACTGATTGCGTGGTAACTATGTCAATTAAAACAATTAAGTATTTCTCAACAATCATTGTAGCGGTAGTTGCGGTTCTTGCCGGATGGTGGCTGTGGAATTATTACATGCAATCTCCGTGGACACGCGACGGGAAAAT  >  NC_000913/1720902‑1721382

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