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	92.1%	125.8 / 3.1	38	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 (18/17);  total (19/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 = 9.95e-01

GTGAACGCCGACACGTTGCGGGAAGCGCAACAGCGCCCACAAGATTATGCCGGGCTGGTGGTGCGCGTTGCCGGATACAGCGCCTTCTTTGTCGAACTGTCGAAGGAGATCCAGGATGACATCATCCGCCGGACAGCGCATCAGCTGTAACGTTGTGGAAACGCGCCGCGATGATGTGGCGCGCATTTTCAACATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGTCT‑TTTTT‑AAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATGCCCCTCCGGGGCGTTTGAACGGATTGGTCGCGATATCAGCCTTGACGCTCTGGAACGGGAAGTGATGAAAGATGACATTTTCTTTCGCACGTCCGGCGGCGGCGTCACGCTTTCTGGC  >  NC_000913/4146143‑4146632                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                    gtgAACGCCGACACGTTGCGGGAAGCGCAACAGCGCCCACAAGATTATGCCGGGCTGGTGGTGCGCGTTGCCGGATACAGCGCCTTCTTTGTCGAACTGTCGAAGGAGATCCAGGATGACATCATCCGCCGGACAGCGCATCAGCTGTAACGTTGTGGAAACGCGCCGCGATGATGTGGCGCGCATTTTCAACATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGTCt‑ttttttaa                                                                                                                                                                                                                                                   >  1:379454/1‑247 (MQ=255)       gCCGACACGTTGCGGGAAGCGCAACAGCGCCCACAAGATTATGCCGGGCTGGTGGTGCGCGTTGCCGGATACAGCGCCTTCTTTGTCGAACTGTCGAAGGAGATCCAGGATGACATCATCCGCCGGACAGCGCATCAGCTGTAACGTTGTGGAAACGCGCCGCGATGATGTGGCGCGCATTTTCAACATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGTCT‑TTTTTTAAAGGCTGt                                                                                                                                                                                                                                            >  1:12615/1‑251 (MQ=255)       gCCGACACGTTGAGGGAAGCGCACCAGCGCCCACAAGATTATGCCGGGCTGGTGGTGCGCGTTGCCGGAGACAGCGCCTTCTTTGTCCAACTGTCGAAGGAGATCCAGGATGACATCATCCGCCGGACAGCGCATCAGCTGTAACGTTGTGGAAACGCGCCGCGATGATGTGGCGCGCCTTTTCAACATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGGGGTCT‑TTTTT‑TAAgggctgg                                                                                                                                                                                                                                           >  2:214202/1‑246 (MQ=255)            cacGTTGCGGGAAGCGCAACAGCGCCCACAAGATTATGCCGGGCTGGTGGTGCGCGTTGCCGGATACAGCGCCTTCTTTGTCGAACTGTCGAAGGAGATCCAGGATGACATCATCCGCCGGACAGCGCATCAGCTGTAACGTTGTGGAAACGCGCCGCGATGATGTGGCGCGCATTTTCAACATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGTCT‑TTTTTTAAAGGCTGTCCGCa                                                                                                                                                                                                                                       <  2:349414/251‑1 (MQ=255)                                   gcCCACAAGATTATGCCGGGCTGGTGGTGCGCGTTGCCGGATACAGCGCCTTCTTTGTCGAACTGTCGAAGGAGATCCAGGATGACATCATCCGCCGGACAGCGCATCAGCTGTAACGTTGTGGAAACGCGCCGCGATGATGTGGCGCGCATTTTCAACATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTAAGGTGGTCT‑TTTTTTAAAGGCTGTCCGCATCTGTGCCCGTGGTATGCTAATc                                                                                                                                                                                                                >  2:325414/1‑251 (MQ=255)                                                 gCCGGGCTGGTGGTGCGTGTTGCCGGATACAGCGCCTTCTTTGTCGAACTGTCGAAGGAGATCCAGGATGACATCATCCGCCGGACAGCGCATCAGTTGTAACGTTGTGGAAACGCGCCGCGATGATGTGGCGCGCATTTTCAACATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTGCGGTGGTCTCTTTTT‑AAAGGCTGTCCGCATCTTTGCACGTGTTGTGCTAATCCGGAGTCGATCTcc                                                                                                                                                                                                  <  2:351365/251‑1 (MQ=255)                                                  ccGGGCTGGTGGTGCGCGTTGCCGGATACAGCGCCTTCTTTGTCGAACTGTCGAAGGAGATCCAGGATGACATCATCCGCCGGACAGCGCATCAGCTGTAACGTTGTGGAAACGCGCCGCGATGATGTGGCGCGCATTTTCAACATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGTCT‑TTTTTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCg                                                                                                                                                                                                 <  1:142627/251‑1 (MQ=255)                                                     ggCTGGTGGTGCGCGTTGCCGGATACAGCGCCTTCTTTGTCGAACTGTCGAAGGAGATCCAGGATGACATCATCCGCCGGACAGCGCATCAGCTGTAACGTTGTGGAAACGCGCCGCGATGATGTGGCGCGCATTTTCAACATTCAGCGTTATTCACTGACTGACGGCGAGGGCATTCGTACGGTGGTCT‑TTTTTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCgac                                                                                                                                                                                               >  2:210928/1‑248 (MQ=255)                                                                 gcgNTGCCGGATACAGCGCCTTCTTTGTCGAACTGTCGAAGGAGATCCAGGATGACATCATCCGCCGGACAGCGCATCAGCTGTAACGTTGTGGAAACGCGCCGCGATGATGTGGCGCGCATTTTCAACATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGTCT‑TTTTTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACgg                                                                                                                                                                                  <  2:356579/251‑1 (MQ=255)                                                                              cAGCGCCTTCTTTGTCGAACTGTCGAAGGAGATCCAGGATGACATCATCCGCCGGACAGCGCATCAGCTGTAACGTTGTGGAAACGCGCCGCGATGATGTGGCGCGCATTTTCAACATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGTCT‑TTTTTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGc                                                                                                                                                                     <  2:292878/251‑1 (MQ=255)                                                                               aGCGCCTTCTTTGTCGAACTGTCGAAGGAGATCCAGGATGACATCATCCGCCGGACAGCGCATCAGCTGTAACGTTGTGGAAACGCGCCGCGATGATGTGGCGCGCATTTTCAACATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGTCT‑TTTTTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCg                                                                                                                                                                    >  2:127381/1‑251 (MQ=255)                                                                                  gcCTTCTTTGTCGAACTGTCGAAGGAGATCCAGGATGACATCATCCGCCGGACAGCGCATCAGCTGTAACGTTGTGGAAACGCGCCGCGATGATGTGGCGCGCATTTTCAACATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGTCT‑TTTTTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGaaa                                                                                                                                                                 <  1:399294/251‑1 (MQ=255)                                                                                                  ggTCGAAGGAGATCCAGGATGACATCATCCGCCGGACAGCGCATCAGCTGTAACGTTGTGGAAACGCGCCGCGATGATGTGGCGCGCATTTTCAACATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGTCT‑TTTTTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCGCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGa                                                                                                                                                 <  1:273932/250‑1 (MQ=255)                                                                                                   gTCGAAGGAGATCCAGGATGACATCATCCGCCGGACAGCGCATCAGCTGTAACGTTGTGGAAACGCGCCGCGATGATGTGGCGCGCATTTTCAACATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGTCT‑TTTTTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGaa                                                                                                                                                <  1:334347/251‑1 (MQ=255)                                                                                                                cAGGATGCCATCATCCGCCGGACAGCGCATCAGCTGTAACGTTGTGGAAACGCGCCGCGATGATGTGGCGCGCATTTTCAACATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGTCT‑TTTTTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGTTTGCGTGAt                                                                                                                                   <  2:103661/251‑1 (MQ=255)                                                                                                                      gACATCATCCGCCGGACAGCGCATCAGCTGTAACGTTGTGGAAACGCGCCGCGATGATGTGGCGCGCATTTTCAACATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGTCT‑TTTTTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGa                                                                                                                              <  1:334351/250‑1 (MQ=255)                                                                                                                              ccgccgGACAGCGCATCAGCTGTAACGTTGTGGAAACGCGCCGCGATGATGTGGCGCGCATTTTCAACATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGTCT‑TTTTTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATGccc                                                                                                                     >  2:378955/1‑251 (MQ=255)                                                                                                                               cgccgGACAGCGCATCAGCTGTAACGTTGTGGAAACGCGCCGCGATGATGTGGCGCGCATTTTCAACATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGTCT‑TTTTTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATGcccc                                                                                                                    >  1:300296/1‑251 (MQ=255)                                                                                                                                gccgGACAGCGCATCAGCTGTAACGTTGTGGAAACGCGCCGCGATGATGTGGCGCGCATTTTCAACATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGTCT‑TTTTTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATGCCCCt                                                                                                                   <  2:125381/251‑1 (MQ=255)                                                                                                                                 ccgGACAGCGCATCAGCTGTAACGTTGTGGAAACGCGCCGCGATGATGTGGCGCGCATTTTCAACATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGTCT‑TTTTTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTc                                                                                                                                                             <  1:203021/208‑1 (MQ=255)                                                                                                                                 ccgGACAGCGCATCAGCTGTAACGTTGTGGAAACGCGCCGCGATGATGTGGCGCGCATTTTCAACATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGTCT‑TTTTTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTc                                                                                                                                                             >  2:203021/1‑208 (MQ=255)                                                                                                                                  cgGACAGCGCATCAGCTGTAACGTTGTGGAAACGCGCCGCGATGATGTGGCGCGCATTTTCAACATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGTCT‑TTTTTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATGCCCNTcc                                                                                                                 >  2:355610/1‑251 (MQ=255)                                                                                                                                                    tAACGTTGTGGAAACGCGCCGCGATGATGTGGCGCGCATTTTCAACATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGTCT‑TTTTTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGTTTGCGGGAGGCGGATGAATGCCCCTCCGGGGCGGTTGAACGGAtt                                                                                               >  2:134805/1‑251 (MQ=255)                                                                                                                                                       cGTTGTGGAAACGCGCCGCGATGATGTGGCGCGCATTTTCAACATTCAGCGATATTCACTGAATGACGGCGAGGGCCTTCGTCCGGTGGTCT‑TTTTTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTACTCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGAAGAGAGGCGAAATGTCTGCACTGTGCCAAATGTTTGCGTGATGCGGATGAATGCCCCTCCGGGGCGTTTGAACGGATTggg                                                                                            >  1:19403/1‑250 (MQ=255)                                                                                                                                                        gTTGTGGAAACGCGCCGCGATGATGTGGCGCGCATTTTCAACATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGTCT‑TTTTTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATGCCCCTCCGGGGCGTTTGAACGGATTGGTc                                                                                           <  1:378955/251‑1 (MQ=255)                                                                                                                                                             ggAAACGCGCCGCGATGATGTGGCGCGCATTTTCAACATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGTCT‑TTTTTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATGCCCCTCCGGGGCGTTTGAACGGATTGGTCGCGat                                                                                      >  1:315453/1‑251 (MQ=255)                                                                                                                                                             ggAAACGCGCAGCGATGATGTGGCGCGCATTTTCAACATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGTCT‑TTTTTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATGCCCCTCCGGGGCGTTTGAACGGATTGGTCGCGat                                                                                      >  2:132548/1‑251 (MQ=255)                                                                                                                                                                              atgTGGTGCGCATTTTCAACATTAAGCGTTATGCACTGAGTGACGGGGAGGGCATTCGTAGGGGGGTCT‑TTTTTTAAAGGCTGTCCGCATATTTGCCAGTGGTGTGCTAATACGGAGTCGATCTCCGGCAAAATCCAGACGGTACGGAGAGAGGCGAACTGACTGCACTGTGGGAAATGGTTGCGTGAGGCGGGTGAATAGCCCTCCGGGGCGTTTGAACGGATTGGTAGAGATATCAGCCTAGACGctct                                                                     >  2:170087/1‑251 (MQ=255)                                                                                                                                                                                gtgGCGCGCATTTTCAACATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGTCT‑TTTTTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATGCCCCTCCGGGGCGTTTGAACGGATTGGTCGCGATATCAGCCTTGACGCTCTgg                                                                   <  1:274343/251‑1 (MQ=255)                                                                                                                                                                                   gcgcgcATTTTCAACATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGCTGGTCT‑TTTTTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATGCCCCTCCGGGGCGTTTGAACGGCTTGGTCGCGATATCAGCCTTGACGCTCTGGAAc                                                                >  2:251774/1‑251 (MQ=255)                                                                                                                                                                                              cAACATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGTCT‑TTTTTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATGCCCCTCCGGGGCGTTTGAACGGATTGGTCGCGATATCAGCCTTGACGCTCTGGAACGGGAAGtgatg                                                     <  2:131775/251‑1 (MQ=255)                                                                                                                                                                                              cAACATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGTCT‑TTTTTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATGCCCCTCCGGGGCGTTTGAACGGATTGGTCGCGATATCAGCCTTGACGCTCTGGAACGGGAAGtgat                                                      <  1:244320/250‑1 (MQ=255)                                                                                                                                                                                                        cGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGTCT‑TTTTTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATGCCCCTCCGGGGCGTTTGAACGGATTGGTCGCGATATCAGCCTTGACGCTCTGGAACGGGAAGTGATGAAAGATGAc                                            >  2:396368/1‑250 (MQ=255)                                                                                                                                                                                                         gTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGTCT‑TTTTTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCCAAATGTTTGCGTGATGCGGATGAATGCCCCTCCGGGGCGTTTGAACGGATTGGTCGCGATATCAGCCTTGACGCTCTGGAACGGGAAGTGATGAAAGATGACAt                                          >  1:144703/1‑251 (MQ=255)                                                                                                                                                                                                           tagtCACTGAATGACGGCGAGGGCATGCGTACGGTGGTCTGTTTTT‑AAATGCTTTCCGCATCTTTGCCCGGGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATGCCCCGCCGGGGGGTTTGAACGGAGTGGTCGCGATATCAGCCTTGACGCTCTGGAACGGGAAGTGATTAAAGAGGACAtt                                         <  2:36735/247‑1 (MQ=255)                                                                                                                                                                                                                   gAATGACGGCGAGGGCATTCGTACGGTGGTCT‑TTTTTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATGCCCCTCCGGGGCGTTTGAACGGATTGGTCGCGATATCAGCCTTGACGCTCTGGAACGGGAAGTGATGAAAGATGACATTTTCTTTCGc                                <  2:300296/251‑1 (MQ=255)                                                                                                                                                                                                                                      cGTACGGTGGTCT‑TTTTTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATGCCCCTCCGGGGCGTTTGAACGGATTGGTCGCGATATCAGCCTTGACGCTCTGGAACGGGAAGTGATGAAAGATGACATTTTCTTTCGCACGTCCGGCGGCGGCGTc              >  2:197835/1‑250 (MQ=255)                                                                                                                                                                                                                                        tACGGTGGTCT‑TTTTTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATGCCCCTCCGGGGCGTTTGAACGGATTGGTCGCGATATCAGCCTTGACGCTCTGGAACGGGAAGTGATGAAAGATGACATTTTCTTTCGCACGTCCGGCGGCGGCGTCACg           >  2:148647/1‑251 (MQ=255)                                                                                                                                                                                                                                          gggtggtCT‑TTTTTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCGAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATGCACCTCCGGGGCGTTTGAACGGATTGGTCGCGATATCAGCCTTGACGCTCTGGAACGGGAAGTGATGAAATATGACATTTTCTTTCGCACGTCCGGAGGCGGCGTCACGcc         <  2:67612/250‑2 (MQ=255)                                                                                                                                                                                                                                               tctt‑ttttt‑AAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATGCCCCTCCGGGGCGTTTGAACGGATTGGTCGCGATATCAGCCTTGACGCTCTGGAACGGGAAGTGATGAAAGATGACATTTTCTTTCGCACGTCCGGCGGCGGCGTCACGCTTTCTgg   >  1:278395/4‑251 (MQ=255)                                                                                                                                                                                                                                                  t‑ttttt‑AAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATGCCCCTCCGGGGCGTTTGAACGGATTGGTCGCGATATCAGCCTTGACGCTCTGGAACGGGAAGTGATGAAAGATGACATTTTCTTTCGCACGTCCGGCGGCGGCGTCACGCTTTCTggc  <  1:214202/249‑1 (MQ=255)                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                    GTGAACGCCGACACGTTGCGGGAAGCGCAACAGCGCCCACAAGATTATGCCGGGCTGGTGGTGCGCGTTGCCGGATACAGCGCCTTCTTTGTCGAACTGTCGAAGGAGATCCAGGATGACATCATCCGCCGGACAGCGCATCAGCTGTAACGTTGTGGAAACGCGCCGCGATGATGTGGCGCGCATTTTCAACATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGTCT‑TTTTT‑AAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATGCCCCTCCGGGGCGTTTGAACGGATTGGTCGCGATATCAGCCTTGACGCTCTGGAACGGGAAGTGATGAAAGATGACATTTTCTTTCGCACGTCCGGCGGCGGCGTCACGCTTTCTGGC  >  NC_000913/4146143‑4146632

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