Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,146,389	(T)6→7	coding (132/879 nt)	pflC →	putative pyruvate formate‑lyase 2 activating enzyme PflC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,146,389	1	.	T	94.2%	174.3 / 1.8	52	coding (132/879 nt)	pflC	putative pyruvate formate‑lyase 2 activating enzyme PflC
Reads supporting (aligned to +/- strand):  ref base . (1/2);  new base T (17/32);  total (18/34)
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.81e-01

TGGTGAACGCCGACACGTTGCGGGAAGCGCAACAGCGCCCACAAGATTATGCCGGGCTGGTGGTGCGCGTTGCCGGATACAGCGCCTTCTTTGTCGAACTGTCGAAGGAGATCCAGGATGACATCATCCGCCGGACAGCGCATCAGCTGTAACGTTGTGGAAACGCGCCGCGATGATGTGGCGCGCATTTTCAACATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGT‑CT‑TTT‑TT‑AAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAA‑TGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATGCCCCTCCGGGGCGTTTGAACGGATTGGTCGCGATATCAGCCTTGACGCTCTGGAACGGGAAGTGATGAAAGATGACATTTTCTTTCGCACGTCCGGCGGCGGCGTCACGCTTTCTGGCGGCGA  >  NC_000913/4146141‑4146637                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                          tggtgAACGCCGACACGTTGCGGGAAGCGCAACAGCGCCCACAAGATTATGCCGGGCTGGTGGTGCGCGTTGCCGGATACAGCGCCTTCTTTGTCGAACTGTCGAAGGAGATCCAGGATGACATCATCCGCCGGACAGCGCATCAGCTGTAACGTTGTGGAAACGCGCCGCGATGATGTGGCGCGCATTTTCAACATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGT‑Ct‑ttt‑ttt                                                                                                                                                                                                                                                           >  2:162143/1‑249 (MQ=255)       aCGCCGACACGGTGCGGGAAGCGCAACAGCGCCCACAAGATTATGCCGGGCTGGTGGGGCGCGTTGCCGGATACAGCGCCTTCTTTGTCGAACTGTCGAAGGAGATCCAGGATGACATCATCCGCCGGACAGCGCATCAGCTGTAACGTTGTGGAAACGCGCCGCGATGATGTGGCGCGCATTTTCAACATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGTTGGT‑CT‑TTT‑TTTAAAGGCt                                                                                                                                                                                                                                                    <  2:445449/251‑1 (MQ=255)          ccGACACGTTGCGGGAAGCGCAACAGCGCCCACAAGATTATGCCGGGCTGGTGGTGCGCGTTGCCGGATACAGCGCCTTCTTTGTCGAACTGTCGAAGGAGATCCAGGATGACATCATCCGCCGGACAGCGCATCAGCTGTAACGTTGTGGAAACGCGCCGCGATGATGTGGCGCGCATTTTCAACATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGT‑CT‑TTT‑TTTAAAGGCTGTc                                                                                                                                                                                                                                                 <  1:298371/251‑1 (MQ=255)          ccGACACGTTGCGGGAAGCGCAACAGCGCCCACAAGATTATGCCGGGCTGGTGGTGCGCGTTGCCGGATACAGCGCCTTCTTTGTCGAACTGTCGAAGGAGATCCAGGATGACATCATCCGCCGGACAGCGCATCAGCTGTAACGTTGTGGAAACGCGCCGCGATGATGTGGCGCGCATTTTCAACATTCAGCGTTATTCAATGAATGACGGGGAGGGCATTAGTACGGTGGT‑CT‑TTT‑TTTAAAGGCTGTc                                                                                                                                                                                                                                                 >  1:14799/1‑251 (MQ=255)           cGACACGTTGCGGGAAGCGCAACAGCGCCCACAAGATTATGCCGGGCTGGTGGTGCGCGTTGCCGGATACAGCGCCTTCTTTGTCGAACTGTCGAAGGAGATCCAGGATGACATCATCCGCCGGACAGCGCATCAGCTGTAACGTTGTGGAAACGCGCCGCGATGATGTGGCGCGCATTTTCAACATTAAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGCGGTGTT‑TTT‑TT‑AAAGGCTGTcc                                                                                                                                                                                                                                                >  2:442591/1‑251 (MQ=255)               acGTTGCGGGAAGCGCAACAGCGCCCACAAGATTATGCCGGGCTGGTGGTGCGCGTTGCCGGATACAGCGCCTTCTTTGTCGAACTGTCGAAGGAGATCCAGGATGACATCATCCGCCGGACAGCGCATCAGCTGTAACGTTGTGGAAACGCGCCGCGATGATGTGGCGCGCATTTTCAACATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGT‑CT‑TTT‑TTTAAAGGCTGTCCGCAt                                                                                                                                                                                                                                            <  2:197240/251‑1 (MQ=255)                cGTTGCGGGAAGCGCAACAGCGCCCACAAGATTATGCCGGGCTGGTGGTGCGCGTTGCCGGATACAGCGCCTTCTTTGTCGAACTGTCGAAGGAGATCCAGGATGACATCATCCGCCGGACAGCGCATCAGCTGTAACGTTGTGGAAACGCGCCGCGATGATGTGGCGCGCATTTTCAACATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGT‑CT‑TTT‑TTTAAAGGCTGTCCGCATc                                                                                                                                                                                                                                           <  1:291976/251‑1 (MQ=255)                      gggAAGCGCAACAGCGCCCACAAGATTATGCCGGGCTGGTGGTGCGCGTTGCCGGATACAGCGCCTTCTTTGTCGAACTGTCGAAGGAGATCCAGGATGACATCATCCGCCGGACAGCGCATCAGCTGTAACGTTGTGGAAACGCGCCGCGATGATGTGGCGCGCATTTTCAACATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGT‑CT‑TTT‑TTTAAAGGCTGTCCGCATCTTTGcc                                                                                                                                                                                                                                     >  1:283095/1‑251 (MQ=255)                         aaGCGCAACAGCGCCCACAAGATTATGCCGGGCTGGTGGTGCGCGTTGCCGGATACAGCGCCTTCTTTGTCGAACTGTCGAAGGAGATCCAGGATGACATCATCCGCCGGACAGCGCATCAGCTGTAACGTTGTGGAAACGCGCCGCGATGATGTGGCGCGCATTTTCAACATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGT‑CT‑TTT‑TTTAAAGGCTGTCCGCATCTTTGCCCgt                                                                                                                                                                                                                                  <  2:349139/251‑1 (MQ=255)                          aGCGCAACAGCGCCCACAAGATTATGCCGGGCTGGTGGTGCGCGTTGCCGGATACAGCGCCTTCTTTGTCGAACTGTCGAAGGAGATCCAGGATGACATCATCCGCCGGACAGCGCATCAGCTGTAACGTTGTGGAAACGCGCCGCGATGATGTGGCGCGCATTTTCAACATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGT‑CTATTT‑TT‑AAAGGCTGTCCGCATCTTTGCCCgtg                                                                                                                                                                                                                                 <  1:402513/251‑1 (MQ=255)                                      cccACAAGATTATGCCGGGCTGGTGGTGCGCGTTGCCGGATACAGCGCCTTCTTTGTCGAACTGTCGAAGGAGATCCAGGATGACATCATCCGCCGGACAGCGCATCAGCTGTAACTTTGTGGAAACGCGCCGCGATGATGTGGCGCGCATTTTCAACATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGT‑CT‑TTT‑TTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATcc                                                                                                                                                                                                                     <  2:183860/251‑1 (MQ=255)                                         acaAGATTATGCCGGGCTGGTGGTGCGCGATGACGGATACAGCGCATTCTTTGACGAACTGTCGAAGGAGATCCAGGAAGACATCATCCGCTGGACAGCGCATCAGCTGTAACGATGTGGAAACGCGCCGCGATGATGTGGCGCGCATTTTCAACATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGAGGA‑CT‑TTT‑TTTAAAGGCTGTCCGCTTCTTTGCCCGTGGTGTGCTAATCCggg                                                                                                                                                                                                                  >  1:144913/1‑250 (MQ=255)                                           aaGATTATGCCGGGCTGGTGGTGCGCGTTGCCGGATACAGCGCCTTCTTTGTCGAACTGTCGAAGGAGATCCAGGATGACATCATCCGCCGGACAGCGCATCAGCTGTAACGTTGTGGAAACGCGCCGCGATGATGTGGCGCGCATTTTCAACATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGT‑CT‑TTT‑TTTAACGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGt                                                                                                                                                                                                                >  2:86751/1‑251 (MQ=255)                                              aTTATGCCGGGCTGGTGGTGCGCGTTGCCGGATACAGCGCCTTCTTTGTCGAACTGGCGAAGGAGATCCAGGATGACATCATCCGCCGGACAGCGCATCAGCTGTAACGTTGTGGAAACGCGCCGCGATGATGTGGCGCGCATTTTCAACATTCAGCGTTATTCACTGAATGACGGCGCGGGCTTTCGTACGGTGGT‑CT‑TTT‑TTTAAAGGCTGCCCGCATCTTTGCCCGTGGTGTGTTAATCCGGAGTCgc                                                                                                                                                                                                             >  2:282421/1‑250 (MQ=255)                                                       ggCTGGTGGTGCGCGTTGCCGGATACAGCGCCTTCTTTGTCGAACTGTCGAAGGAGATCCAGGATGACATCATCCGCCGGACAGCGCATCAGCTGTAACGTTGTGGAAACGCGCCGCGATGATGTGGCGCGCATTTTCAACATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGT‑CT‑TTT‑TTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCa                                                                                                                                                                                                    <  2:359798/251‑1 (MQ=255)                                                                          cGGATACAGCGCCTTCTTTGTCGAACTGTCGAAGGAGATCCAGGATGACATCATCCGCCGGACAGCGCATCAGCTGTAACGTTGTGGAAACGCGCCGCGATGATGTGGCGCGCATTTTCAACATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGT‑CT‑TTT‑TTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCag                                                                                                                                                                                 <  2:169642/251‑1 (MQ=255)                                                                                 aGCGCCTTCTTTGTCGAACTGTCGAAGGAGATCCAGGATGACATCATCCGCCGGACAGCGCATCAGCTGTAACGTTGTGGAAACGCGCCGCGATGATGTGGCGCGCATTTTCAACATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGT‑CT‑TTT‑TTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCg                                                                                                                                                                          >  1:115917/1‑251 (MQ=255)                                                                                             gTCGAACTGTCGAAGGAGATCCAGGATGACATCATCCGCCGGACAGCGCATCAGCTGTAACGTTGTGGAAACGCGCCGCGATGATGTGGCGCGCATTTTCAACATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGT‑CT‑TTT‑TTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAA‑TGTCTGCAc                                                                                                                                                             <  2:41132/251‑1 (MQ=255)                                                                                                        gAAGGAGATCCAGGATGACATCATCCGCCGGACAGCGCATCAGCTGTAACGTTGTGGAAACGCGCCGCGATGATGTGGCGCGCATTTTCAACATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGT‑CT‑TTT‑TTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAA‑TGTCTGCACTGTGCGAAATg                                                                                                                                                  <  1:412106/251‑1 (MQ=255)                                                                                                          aGGAGATCCAGGATGACATCATCCGCCGGACAGCGCATCAGCTGTAACGGTGTGGAAACGCGCCGCGATGATGTGGCGCGCATTTTCAACATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGT‑CT‑TTT‑TTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAA‑TGTCTGCACTGTGCGAAATGtt                                                                                                                                                <  2:214295/251‑1 (MQ=255)                                                                                                                tCCAGGATGACATCATCCGCCGGACAGCGCATCAGCTGTAACGTTGTGGAAACGCGCCGCGATGATGTGGCGCGCATTTTCAACATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGT‑CT‑TTT‑TTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAA‑TGTCTGCACTGTGCGAAATGTTTGCGTg                                                                                                                                          >  1:398932/1‑251 (MQ=255)                                                                                                                        gACNTCAGCCGCCGGACAGCGCATCAGCTGTAACGTGGTGGAAACGCGCCGCGATGATGTGGCGCGCATTTTCAACATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGT‑CT‑TTT‑TTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAA‑TGTCTGCACTGTTCGAAATGTTTGCGTGATGCGGAt                                                                                                                                  <  2:406950/251‑1 (MQ=255)                                                                                                                             catcCGCCGGACAGCGCATCAGCTGTAACGTTGTGGAAACGCGCCGCGATGATGTGGCGCGCATTTTCAACATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGT‑CT‑TTT‑TTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAA‑TGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAAt                                                                                                                              <  2:300194/250‑1 (MQ=255)                                                                                                                                  gccgGACAGCGCATCAGCTGTAACGTTGTGGAAACGCGCCGCGATGATGTGGCGCGCATTTTCAACATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGT‑CT‑TTT‑TTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGATGCGAAA‑TGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATGCCCCt                                                                                                                        <  1:252221/251‑1 (MQ=255)                                                                                                                                       aCAGCGCATCAGCTGTAACGTTGTGGAAACGCGCCGCGATGATGTGGCGCGCATTTTCAACATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGT‑CT‑TTT‑TTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAA‑TGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATGCCCCTCCggg                                                                                                                   <  2:283095/251‑1 (MQ=255)                                                                                                                                             cATCAGCTGTAACGTTGTGGAAACGCGCCGCGATGATGTGGCGCGCATTTTCAACATTCAGCGTTATTCACTGAATGACGGCGCGGGCATTCGTACGGTGGT‑CT‑TTT‑TTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAA‑TGTCTGCACTGTGCGAACTGTTTGCGTGATGCGGATGAATGCCCCTCCGGGGCGttt                                                                                                             >  2:221886/1‑251 (MQ=255)                                                                                                                                               tCAGCTGTAACGTTGTGGAAACGCGCCGCGATGATGTGGCGCGCATTTTCAACATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGT‑CT‑TTT‑TTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAA‑TGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATGCCCCTCCGGGGCGTTTg                                                                                                            >  1:442712/1‑250 (MQ=255)                                                                                                                                                     gTAACGTTGTGGAAACGCGCCGCGATGATGTGGCGCGCATTTTCAACATTCAGCGTTATGCACTGAATGACGGCGAGGGCATTCGTACGGTGGT‑CT‑TTT‑TTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAA‑TGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATGCCCCTCCGGTGCGTTTGAACGGa                                                                                                      <  2:177606/250‑1 (MQ=255)                                                                                                                                                                     gcgcCGCGATGATGTGGCGCGCATTTTCAACATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGT‑CT‑TTT‑TTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAA‑TGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATGCCCCTCCGGGGCGTTTGAACGGATTGGTCGCGATATCAGc                                                                                     <  1:164648/251‑1 (MQ=255)                                                                                                                                                                        ccgcGATGATGTGGGGCGCATTTTCAACATTCAGCGTGATTCACGTAATGACGGCGATGGAATTCGTACGGTGGT‑CT‑TTTGTT‑AAAGGCGGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGTAGTCGATCTCCGGCACAATCCAGACGGTACGCAGATAGGCGAAA‑TGTCTGCACTTTGCGAAATGTTTGCGTGATGCGGATGAATGCCTGTCCGGGGTGTTTGAACGGATTGGTGGTGATATCAGCCtt                                                                                  <  2:263168/251‑1 (MQ=255)                                                                                                                                                                         cgcgATGATGTGGCGCGCATTTTCAACATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGT‑CT‑TTT‑TTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAA‑TGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATGCCCCTCCGGGGCGTTTGAACGGATTGGTCGCGATATCAGCCTTg                                                                                 <  1:162143/251‑1 (MQ=255)                                                                                                                                                                          gcgATGATGGGGCGCGCATTTTCAACATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGT‑CT‑TTT‑TTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGTCGGTACGCAGAGAGGCGAAA‑TGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATGCCCCTCCGGGGCGTTTGAACGGATTGGTCGCGATATCAGCCTTGa                                                                                <  2:227641/251‑1 (MQ=255)                                                                                                                                                                                  gtgGCGCGCATTTTCAACATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGT‑CT‑TTT‑TTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAA‑TGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATGCCCCTCCGGGGCGTTTGAACGGATTGGTCGCGATATCAGCCTTGACGCTCTg                                                                         <  1:12913/250‑1 (MQ=255)                                                                                                                                                                                      cgcgcATTTTCAACATTAAGCGTTATTCACTGAAGGACGGCGAGGGCATTCGTACGGTGGT‑CT‑TTT‑TTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAA‑TGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATGCCCCTCCGGGGCGTTTGAACGGATTGGTCGCGATATCAGCCTTGACGCGCTGGAACg                                                                    <  2:46649/251‑1 (MQ=255)                                                                                                                                                                                       gcgcATTTTCAACATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGT‑CT‑TTT‑TTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAA‑TGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATGCCCCTCCGGGGCGTTTGAACGGATTGGTCGCGATATCAGCCTTGACGCTCTGGAAc                                                                     <  2:326594/249‑1 (MQ=255)                                                                                                                                                                                       gcgcATTTTCAACATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGT‑CT‑TTT‑TTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAA‑TGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATGCCCCTCCGGGGCGTTTGAACGGATTGGTCGCGATATCAGCCTTGACGCTCTGGAAc                                                                     >  1:326594/1‑249 (MQ=255)                                                                                                                                                                                       gcgcATTTTCAACATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGT‑CT‑TTT‑TTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAA‑TGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATGCCCCTCCGGGGCGTTTGAACGGATTGGTCGCGATATCAGCCTTGACGCTCTGGAACgg                                                                   >  2:304548/1‑251 (MQ=255)                                                                                                                                                                                        cgcATTTTCAACATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGT‑CT‑TTT‑TTTAAAGGCTGTCCGCATCTTTGCACGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAA‑TGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATGCCCCTCCGGGGCGTTTGAACGGATTGGTCGCGATATCAGCCTTGACGCTCTGGAACggg                                                                  >  1:165997/1‑251 (MQ=255)                                                                                                                                                                                         gcATTTTCAACATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGT‑CT‑TTT‑TTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAA‑TGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATGCCCCTCCGGGGCGTTTGAACGGATTGGTCGCGATATCAGCCTTGACGCTCTGGAACGGGa                                                                 <  2:232054/251‑1 (MQ=255)                                                                                                                                                                                                      tCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGT‑CT‑TTT‑TTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAA‑TGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATGCCCCTCCGGGGCGTTTGAACGGATTGGTCGCGATATCAGCCTTGACGCTCTGGAACGGGAAGTGATGAAAGa                                                     <  2:306921/250‑1 (MQ=255)                                                                                                                                                                                                      tCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGT‑CT‑TTT‑TTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAA‑TGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATGCCCCTCCGGGGCGTTTGAACGGATTGGTCGCGATATCAGCCTTGACGCTCTGGAACGGGAAGTGATGAAAGa                                                     <  1:217729/250‑1 (MQ=255)                                                                                                                                                                                                        aGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGT‑CT‑TTT‑TTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAA‑TGTCTGCACTGTGCGAAATGTTTGCGTGATTCGGATGAATGCCCCTCCGGGGCGTTTTAACGGATTGGTCGCGATATCAGCCTTGACGCTCTGGAACGGGAAGTGATGAAAGATGa                                                  <  1:64138/251‑1 (MQ=255)                                                                                                                                                                                                         gCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGT‑CT‑TTT‑TTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAA‑TGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATGCCCCTCCGGGGCGTTTGAACGGATTGGTCGCGATATCAGCCTTGACGCTCTGGAACGGGAAGTGATGAAAGATGAc                                                 <  1:117081/251‑1 (MQ=255)                                                                                                                                                                                                          cGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGT‑CT‑TTT‑TTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAA‑TGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATGCCCCTCCGGGGCGTTTGAACGGATTGGTCGCGATATCAGCCTTGACGCTCTGGAACGGGAAGTGATGAAAGATGACa                                                <  1:40015/251‑1 (MQ=255)                                                                                                                                                                                                           gTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGT‑CT‑TTT‑TTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAA‑TGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATGCCCCTCCGGGGCGTTTGAACGGATTGGTCGCGATATCAg                                                                                      >  1:95996/1‑212 (MQ=255)                                                                                                                                                                                                           gTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGT‑CT‑TTT‑TTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAA‑TGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATGCCCCTCCGGGGCGTTTGAACGGATTGGTCGCGATATCAg                                                                                      <  2:95996/212‑1 (MQ=255)                                                                                                                                                                                                           gTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGT‑CT‑TTT‑TTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAA‑TGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATGCCCCTCCGGGGCGTTTGAACGGATTGGTCGCGATATCAGCCTTGACGCTCTGGAACGGGAAGTGATGAAAGATGACAt                                               <  1:86751/251‑1 (MQ=255)                                                                                                                                                                                                                            ggCGAGGGCATTCGTACGGTGGT‑CT‑TTT‑TTTAAAGGCTGTACGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAAATGTCTGCACTGTGCGAAATGTTGGCGTGATGCg                                                                                                                                     >  2:91115/1‑149 (MQ=255)                                                                                                                                                                                                                             gCGAGGGCATTCGTACGGTGGT‑CT‑TTT‑TTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACTGTACGCAGAGAGGCGAAA‑TGTCTGCACTGTGCGAAATGTTTGCGTTATGCGGATGAATGCCCCTCCGGTGCGTTTGAACGGATTGGTCGCGATATCAGCCTTGACGCTCTGGAACGGGAAGTGATGAAAGATGACATTTTCTTTCGCACGTCcgg                             <  1:442591/251‑1 (MQ=255)                                                                                                                                                                                                                               gAGGGCATTCGTACGGTGGT‑CT‑TTT‑TTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAA‑TGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATGCCCCTCCGGGGCGTTTGAACGGATTGGTCGCGATATCAGCCTTGACGCTCTGGAACGGGAAGTGATGAAAGATGACATTTTCTTTCGCACGTCcggcg                           <  2:115917/251‑1 (MQ=255)                                                                                                                                                                                                                               gAGGGCATTCGTACGGTGGT‑CT‑TTT‑TTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAA‑TGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATGCCCCTCCGGGGCGTTTGAACGGATTGGTCGCGATATCAGCCTTGACGCTCTGGAACGGGAAGTGATGAAAGATGACATTTTCTTTCGCACGTCcggcg                           >  1:438075/1‑251 (MQ=255)                                                                                                                                                                                                                                        cGTACGGTGGT‑CT‑TTT‑TTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAA‑TGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATGCCCCTCCGGGGCGTTTGAACGGATTGGTCGCGATATCAGCCTTGACGCTCTGGAACGGGAAGTGATGAAAGATGACATTTTCTTTCGCACGTCCGGCGGCGGCGTc                   >  1:248150/1‑250 (MQ=255)                                                                                                                                                                                                                                            cggtgggctt‑ttt‑tt‑AAAGGCTGTCCGCATCTTTGCCAGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAA‑TGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATGCCC‑TCCGGGGCGTTTGAACGGATTGGTCGCGATATCAGCCTTGACGCTCTGGAACGGGAAGTGATGAAAGATGACATTTTCTTTCGCACGTCCGGCGGCGGCGTCACGCtt             >  1:192768/10‑251 (MQ=255)                                                                                                                                                                                                                                               tggt‑CT‑TTT‑TTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAA‑TGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATGCCCCTCCGGGGCGTTTGAACGGATTGGTCGCGATATCAGCCTTGACGCTCTGGAACGGGAAGTGATGAAAGATGACATTTTCTTTCGCACGTCCGTCGGCGGCGTCACGCTTTc           >  2:19598/1‑251 (MQ=255)                                                                                                                                                                                                                                                  tctt‑ttt‑tt‑AAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAA‑TGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATGCCCCTCCGGGGCGTTTGAACGGATTGGTCGCGATATAAGCCTTGACGCTCTGGAACGGGAAGTGATGAAAGATGACATTTTCTTTCGCACGTCCGGCGGCGGAGTCACGCTTTCTgg        >  1:72935/4‑251 (MQ=255)                                                                                                                                                                                                                                                     t‑ttt‑tt‑AAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAA‑TGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATGCCCCTCCGGGGCGTTTGAACGGATTGGTCGCGATATCAGCCTTGACGCTCTGGAACGGGAAGTGATGAAAGATGACATTTTCTTTCGCACGTCCGGCGGCGGCGTCACGCTTTCTgg        <  2:282491/248‑1 (MQ=255)                                                                                                                                                                                                                                                           tt‑AAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAA‑TGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATGCCCCTCCGGGGCGTTTGAACGGATTGGTCGCGATATCAGCCTTGACGCTCTGGAACGGGAAGTGATGAAAGATGACATTTTCTTTCGCACGTCCGGCGGCGGCGTCACGCTTTCTGGCGGCGa  <  2:403167/250‑1 (MQ=255)                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                          TGGTGAACGCCGACACGTTGCGGGAAGCGCAACAGCGCCCACAAGATTATGCCGGGCTGGTGGTGCGCGTTGCCGGATACAGCGCCTTCTTTGTCGAACTGTCGAAGGAGATCCAGGATGACATCATCCGCCGGACAGCGCATCAGCTGTAACGTTGTGGAAACGCGCCGCGATGATGTGGCGCGCATTTTCAACATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGT‑CT‑TTT‑TT‑AAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAA‑TGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATGCCCCTCCGGGGCGTTTGAACGGATTGGTCGCGATATCAGCCTTGACGCTCTGGAACGGGAAGTGATGAAAGATGACATTTTCTTTCGCACGTCCGGCGGCGGCGTCACGCTTTCTGGCGGCGA  >  NC_000913/4146141‑4146637

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