Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,146,384	A→G	V71A (GTA→GCA)	yceF ←	m(7)GTP pyrophosphatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,146,384	0	A	G	98.6%	223.5 / ‑4.6	70	V71A (GTA→GCA)	yceF	m(7)GTP pyrophosphatase
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (37/32);  total (37/33)
Fisher's exact test for biased strand distribution p-value = 4.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.26e-01

CTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAGATGCCCATTCGCCGAATTAAACAACGCCAGTCCAGTATAGAAGGTGACGATATTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATACCTGGTCTGAACCAATAATTAAATGATCCGGATAAC‑GTGACGCCAGAGATTGCGCTTTTTCTTGTGCCAGTCGAAGCACCAACTGT‑CGCGGTGATTCGTCGCTGCGTG‑GGGTTTCGTCGACCTCTG‑GTGCTGCACATTCAAAAGAGATTTGCAGTTTTTCCAGAAGAGCGCGACGCCATGGCGAGGTGGAGGCTAAAATAAGTTTAGGCATATTTTTTTCCATCAGATATAGCGTATTGATGA  >  NC_000913/1146138‑1146627                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                        ctctTAAAGCTACCCGCGCAGTGCAGGGTATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAGATGCCCATTCGCCGAATTAAACAACGCCAGTCCAGTATAGAAGGTGACGATATTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCCTCAAGAACACATGCCTg                                                                                                                                                                                                                                                     >  2:415865/1‑251 (MQ=255)     tAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAGATGCCCATTCGCCGAATTAAACAACGCCAGTCCAGTATAGAAGGTGACGATATTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATGCCTGGTc                                                                                                                                                                                                                                                  <  2:238072/250‑1 (MQ=255)       aaGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAGATGCCCATTCGCCGAATTAAACAACGCCAGTCCAGTATAGAAGGTGAAGATATTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATGCATGGTCTg                                                                                                                                                                                                                                                >  1:559388/1‑250 (MQ=255)            tcccGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAGATGCCCATTCGCCGAATTAAACAACGCCAGTCCAGTATAGAAGGTGACGATATTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATGCCTGGTCTGAAcc                                                                                                                                                                                                                                            <  1:540028/248‑1 (MQ=255)               cgcgCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAGATGCCCATTCGCCGAATTAAACAACGCCAGTCCAGTATAGAAGGTGACGATATTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATGCCt                                                                                                                                                                                                                                                      <  1:275423/236‑1 (MQ=255)               cgcgCAGTGAAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAGATGCCCATTCGCCGAATTAAACAACGCCAGTCCAGTATAGAAGGTGACGATATTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGCTTTCGCCATCAAGAACACATGCCt                                                                                                                                                                                                                                                      >  2:275423/1‑236 (MQ=255)                  gcAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAGATGCCCATTCGCCGAATTAAACAACGCCAGTCCAGTATAGAAGGTGACGATATTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATGCCTGGTCTGAACCAATAATTa                                                                                                                                                                                                                                    <  1:41258/251‑1 (MQ=255)                     gtgcagtgtaatctctTTATGCACGTAATTATCACTCGCCGCCTCGTTCAGGTGTCGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAGATGCCCATTCGCCGAATTAAACAACGCCAGTCCAGTCTAGAAGGTGACGATATTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTGCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATGCCTGGTCTGAACCAATAATTaa                                                                                                                                                                                                                                   <  2:513579/237‑1 (MQ=255)                     gTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGAGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAGGTGCCCATTCGCCGAATTAAACAACGCCAGTCCAGTATAGAAGGTGACGATATTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATGCCTGGTCTGAACCAATAATTAAAt                                                                                                                                                                                                                                 >  1:586744/1‑251 (MQ=255)                                gCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAGATGCCCATTCGCCGAATTAAACAACGCCAGTCCAGTATAGAAGGTGACGATATTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATGCCTGGTCTGAACCAATAATTAAATGATCCGGATa                                                                                                                                                                                                                       <  1:440109/250‑1 (MQ=255)                                   ctTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAGATGCCCATTCGCCGAATTAAACAACGCCAGTCCAGTATAGAAGGTGACGATATTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATGCCTGGTCTGAACCAATAATTAAATGATCCGGATAAC‑gg                                                                                                                                                                                                                  >  1:651471/1‑250 (MQ=255)                                       aCGCACGTAATTATCAAGCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAGATGCCCATTCGCCGAATTAAACAACGCCAGTCCAGTATAGAAGGTGACGATATTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATGCCTGGTCTGAACCAATAATTAAATGATCCGGATAAC‑GTGACg                                                                                                                                                                                                              <  1:103339/251‑1 (MQ=255)                                                 ttATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAGATGCCCATTCGCCGAATTAAACAACGCCAGTCCAGTATAGAAGGTGACGATATTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATGCCTGGTCTGAACCAATAATTAAATGATCCGGATAAC‑GTGACGCCAGAGATTg                                                                                                                                                                                                    >  1:278388/1‑251 (MQ=255)                                                  tATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAGATGCCCATTCGCCGAATTAAACAACGCCAGTCCAGTATAGAAGGTGACGATATTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATGCCTGGTCTGAACCAATAATTAAATGATCCGGATAAC‑GTGACGCCAGAGATTgc                                                                                                                                                                                                   >  1:306830/1‑251 (MQ=255)                                                     ccatcgccGCCTCGCACAGGTGACGGAAATGGGCGTCAAAAGGCTCCACCTCGGTTTGCAGATGCCCATTCGCCGAATTTAACAACGCCAGTCCAGTATAGAAGGTGACGATATTGCCGCTGGCTGTGCGTAATTGCAGACGGTCATTTTCTTCCGTTAACGTTTTACCGGTGATTTCGCCATCAATAACACATGCCTGGTCTGAAACAATAATTAAATGAGCCGGATAAC‑GTGACGCCAGAGATTGCGCt                                                                                                                                                                                                <  2:642354/245‑1 (MQ=255)                                                     cAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAGATGCCCATTCGCCGAATTAAACAACGCCAGTCCAGTATAGAAGGTGACGATATTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATGCCTGGTCTGAACCAATAATTAAATGATCCGGATAAC‑GTGACGCCAGAGATTGCGCt                                                                                                                                                                                                >  1:423914/1‑251 (MQ=255)                                                        tctcCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAGATGCCCATTCGCCGAATTAAACAACGCCAGTCCAGTATAGAAGGTGACGATATTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATGCCTGGTCTGAACCAATAATTAAATGAt                                                                                                                                                                                                                              <  2:479159/219‑1 (MQ=255)                                                        tctcCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAGATGCCCATTCGCCGAATTAAACAACGCCAGTCCAGTATAGAAGGTGACGATATTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATGCCTGGTCTGAACCAATAATTAAATGAt                                                                                                                                                                                                                              >  1:479159/1‑219 (MQ=255)                                                           tcGCCTCGATCAGATGACAGCAATGCACGGCAACATGCTCAACTTCTGTTTGCAGATGCCCATTCGCCGTATTAAACAACGCCAGTCCAGGATAGTAGGTGACGCTATTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATATCGCCATCACGAACTCATGCCTGGTCTGAACCAATAATTAATTGATCCGGATACA‑GTGACGCCAGAGATTACGCTTTTTCt                                                                                                                                                                                          >  2:465039/2‑251 (MQ=255)                                                           ccGCCTCTCTCAGGTGACGGAAATGGACGTCAAAATGCCCCACGACTGTTTGCAGATGCCCATTCGCCGAATTAAACAACGCCAGTCCAGTATAGAAGTTGACGAAATGGCCGCTTGCTTTGCGTAATTGCAGACGGGCATTTTCTTCCGTTATCGGTTTACCGGTGATTTCGCCATCAAGAACGCATGCCTGGGATGAACCAATAATTAAATGATCCGGATAAC‑TTGGCGCCAGAGATTGCGCtttt                                                                                                                                                                                             <  2:54476/248‑1 (MQ=255)                                                           ccGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAGATGCCCATTCGCCGAATTAAACAACGCCAGTCCAGTATAGAAGGTGACGCTATTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCCGTGAACGGTTTACCGGTGATTTCGCCATCAAGAACACATGCCTGGTCTGAACCATTAATTAAATGATCCGGATAAC‑GTGACGCCAGAGATTGCGCTTTTTCt                                                                                                                                                                                          <  2:181823/251‑1 (MQ=255)                                                           ccGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAGATGCCCATTCGCCGAATTAAACAACGCCAGTCCAGTATAGAAGGTGACGATATTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATGCCTGGTCTGAACCAATAATTAAATGATCCGGATAAC‑GTGACGCCAGAGATTGCGCttt                                                                                                                                                                                              <  1:103958/247‑1 (MQ=255)                                                                      aGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAGATGCCCATTCGCCGAATTAAACAACGCCAGTCCAGTATAGAAGGTGACGATATTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATGCCTGGTCTGAACCAATAATTAAATGATCCGGATAAC‑GTGACGCCAGAGATTGCGCTTTTTCTTGTGCCAGTCg                                                                                                                                                                               >  1:108196/1‑251 (MQ=255)                                                                              gAAATGGACGTCAAAAGGCGTCACTTCTGTTTGCAGATGTCCAGTCGCCGAATTAAACAACGCCAGTCCAGTATAGAAGGTGACGATATTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATGCCTGGTCTGAACCAATAATTAAATGATCCGGATAAC‑GTGACGCCAGAGATTGCGCTGTTTCTTGTGCAAGGCGAAGCACCa                                                                                                                                                                       <  2:70453/251‑1 (MQ=255)                                                                                               gCTCCACTTCTGTTTGCAGATGCCCATTCGCCGAATTAAACAACGCCAGTCCAGTATAGAAGGTGACGATATTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATGCCTGGTCTGAACCAATAATTAAATGATCCGGATAAC‑GTGACGCGAGAGATTGCGCTTTTTCTTGTGCCAGTCGAAGCACCAACTGT‑CGCGGTGATtcg                                                                                                                                                     >  1:578202/1‑251 (MQ=255)                                                                                                      ttCTGTTTGCAGATGCCCATACGCCGAATTAAACAACGCCAGTCCAGTATAGTAGGTGACGATATTGCCGCTGGCTTTGTGTAAGTGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACGCATGCCTTGTCTGAACCAATAATTAAATGATCCGGATAAC‑GTGACGCCAGAGATTGCGCTTTTTCTTGTGCCAGTCGAAGCACCAACTGT‑CGCGGTGATTCGTCGCt                                                                                                                                                <  2:538189/249‑1 (MQ=255)                                                                                                        cTGTTTGCAGATGCCCATTCGCCGAATTAAACAACGCCAGTCCAGTATAGAAGGTGACGATATTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATGCCTGGTCTGAACCAATAATTAAATGATCCGGATAAC‑GTGACGCCAGAGATTGCGCTTTTTCTTGTGCCAGTCGAAGCACCAACTGT‑CGCGGTGATTCGTCGCTGCGt                                                                                                                                            >  1:434416/1‑251 (MQ=255)                                                                                                         tgtttaaagctGCCCATTCGCAGAATTAAACACCGCCAGTCCAGTATAGAAGGGGACGATATTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATGCCTGGTATGAACCAATAATTAAATGATCCGGATAAC‑GTGACGCCAGAGCTTGCGCTTTTTCTTGTGCCAGTCGAAGCACCAACTGT‑CGCGGTGATTCGTCGCTGCGTg                                                                                                                                           >  2:425227/11‑251 (MQ=255)                                                                                                          gtgtGCAGATGACCATTCGCCGAATTAATCAACGCCAGTCCAGTATAGAAGGTGACGATAGTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATGCCTGGTCTGAACCAATAATTAAATGATCCGGATAAC‑GTGACGCCAGAGATTGCGCTTTTTCTTGTGCCAGTCGAAGCACCAACTGT‑CGCGGTGATTCGTCGCTGCg                                                                                                                                             <  1:47371/245‑1 (MQ=255)                                                                                                          gTTTGCAGATGCCCATTCGCCGAATTAAACAACGCCAGTCCAGTATAGAAGGTGACGATATTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATGCCTGGTCTGAACCAATAATTAAATGATCCGGATAAC‑GTGACGCCAGAGATTGCGCTTTTTCTTGTGCCAGTCGAAGCACCAACTGT‑CGCGGTGATTCGTCGCTGCGTg                                                                                                                                           <  1:483538/250‑1 (MQ=255)                                                                                                                    gCCCATTCGCCGAATTAAACAACGCCAGTCCAGTATAGAAGGTGACGATATTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATGCCTGGTCTGAACCAATAATTAAATGATCCGGATAAC‑GTGACGCCAGAGATTGCGCTTTTTCTTGTGCCAGTCGAAGCACCAACTGT‑CGCGGTGATTCGTCGCTGCGTG‑GGGTTtcgtc                                                                                                                                <  2:278388/250‑1 (MQ=255)                                                                                                                          tCGCCGAATTAAACAACGCCAGTCCAGTATAGAAGGTGACGATATTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATACCTGGTCTGAACCAATAATTAAATGATCCGGATAAC‑GTGACGCCAGAGATTGCGCTTTTTCTTGTGCCAGTCGAAGCACCAACTGT‑CGCGGTGATTCGTCGCTGCGTG‑GGGTTTCGTCGACctc                                                                                                                          <  1:327061/250‑1 (MQ=255)                                                                                                                           cGCCGAATTAAACAACGCCAGTCCAGTATAGAAGGTGACGATATTGCCGCTGGCTTTGCGTAAATGCAGACGGGCATGTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACTGATGCCTGGTCTGAACCAATAAGTAAAGGATCCGGAGAAC‑GTGACGCCAGAGATTGCGCTGTTTCTTGTGCCAGTCGAAACGCCAAC‑GTGCGCGGGGATGCGTCGCTGCGTG‑TGGTTTCGTCGACctct                                                                                                                         <  2:332606/250‑1 (MQ=255)                                                                                                                                aaTTAAACAACGCCAGTCCAGTATAGAAGGGGACGATATTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCACGAACACATGCCTGGTCTGAACCAATAATTAAATGATCCGGATAAC‑GTGACGCCAGAGATTGCGCTTTTTCTTGTGCCAGTCGAAGCACCAACTGT‑CGCGGTGATTCGTCGCTGCGTG‑GGGTTTCGTCGACCTCTG‑Gtgat                                                                                                                  >  2:132554/1‑249 (MQ=255)                                                                                                                                 aTTAAACAACGCCAGTCCAGTATAGAAGGAGACGATATTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATGCCTGGTCTGAACCAATAATTAAATGATCCGGATAAC‑GTGACGCCAGAGATTGCGCTTTTTCTTGTGCCAGTCGAAGCACCAACTGT‑CGCGGTGATTCGTCGCTGCGTG‑GGGTTTCGTCGACCTCTG‑Gagctg                                                                                                                 >  1:459611/1‑251 (MQ=255)                                                                                                                                   tAAACAACGCCAGTCCAGTATAGAAGGTGACGATATTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATGCCTGGTCTGAACCAATAATTAAATGATCCGGATAAC‑GTGACGCCAGAGATTGCGCTTTTTCTTGTGCCAGTCGAAGCACCAACTGT‑CGCGGTGATTCGTCGCTGCGTG‑GGGTTTCGTCGACCTCTG‑Gtgctgc                                                                                                                >  1:403118/1‑250 (MQ=255)                                                                                                                                     aacaaaGCCAGTACAGTCTAGAAGGCGACGATATTGCAGCTGGCTTTGCGTAATTGCAGACGGGCATTTTGATCCGATAACGGATTACCGGTGATTTCGAACTCAAGACCACATGCCTGTTCTGAACAAATAATTATATGATCCGGCTAAC‑GTAACGCCAGAGATTGCGATTTGTCTAGTGCCAGCAGACGAACCAACTGT‑CGAGGTGCTTCGTCGCTGC‑TGTGGGTTTAGTCGACCTCTA‑GTGCTGcac                                                                                                              >  2:328048/1‑250 (MQ=255)                                                                                                                                         acGCCAGTCCAGTATAGAAGGTGACGATATTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATGCCTGGTCTGAACCAATAATTAAATGATCCGGATAAC‑GTGACGCCAGAGATTGCGCTTTTTCTTGTGCCAGTCGAAGCACCAACTGT‑CGCGGTGATTCGTCGCTGCGTG‑GGGTTTCGTCGACCTCTG‑GTGCTGCACATTCa                                                                                                         <  1:530179/251‑1 (MQ=255)                                                                                                                                          cGCCAGTCCAGAATAGAAGGTGACGATATTGACGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATGCCTGGTCTGAACCAATAATTAAATGATCCGGATAAC‑GTGACGCCAGAGATTGCGCTTTTTCTTGTGCAAGTCGAAGCACCAACTGT‑CGCGGTGATTCGTCGCTGCGTG‑GGGTTTCGTCGACCTCTG‑GAGCTGCACATTCa                                                                                                         >  2:446824/1‑250 (MQ=255)                                                                                                                                            acacaccaaaataGAAGATGACGAAATTGACTCTGTCTGTGCGTAATTGCAGACGGGCATTTTCTTCCGTTCACGGCTTACCGGTGATTTCGCCATCAAGAACCCATGCCTGGTCTGACCCAATAATTAAATTATC‑AGTTAACAGATAAGCACGAGATTTCGCTTTTTCTTGTGTCAGTAGAACCACCAACTGA‑CGCGGTAGTTCGTCGCTGCGTG‑GGGTTTCGGCCACCT‑TGTGTGCTGCACATTCaaa                                                                                                       >  2:651702/11‑250 (MQ=255)                                                                                                                                                    gtATAGAAGGTGACGATATTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATGCCTGGTCTGAACCAATAATTAAATGATCCGGATAAC‑GTGACGCCAGAGATTGCGCTTTTTCTTGTGCCAGTCGAAGCACCAACTGT‑CGCGGTGATTCGTCGCTGCGTG‑GGGTTTCGTCGACCTCTG‑GTGCTGCACATTCAAAAGAGATTTg                                                                                              >  1:574829/1‑251 (MQ=255)                                                                                                                                                      ataGAAGGTGACGATATTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATGCCTGGTCTGAACCAATAATTAAATGATCCGGATAAC‑GTGACGCCAGAGATTGCGCTTTTTCTTGTGCCAGTCGAAGCACCAACTGT‑CGCGGTGATTCGTCGCTGCGTG‑GGGTTTCGTCGACCTCTG‑GTGCTGCACATTCAAAAGAGATTTGc                                                                                             >  1:298088/1‑250 (MQ=255)                                                                                                                                                      ataGAAGGTGACGATATTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATGCCTGGTCTGAACCAATAATTAAATGATCCGGATAAC‑GTGACGCCAGAGATTGCGCTTTTTCTTGTGCCAGTCGAAGCACCAACTGT‑CGCGGTGATTCGTCGCTGCGTG‑GGGTTTCGTCGACCTCTG‑GTGCTGCACATTCAAAAGAGATTTGc                                                                                             >  1:421953/1‑250 (MQ=255)                                                                                                                                                       tagcaGGTGACGATATTGACGCTGGATTTGCGTAATTGCAGACGGGCATTTTATTCAGTAAGCGGTTTACCGGTGATTTCGCCATCAAGAACACATGCCTGGTCTGAACCAATAATTAAATGATCCGTATAAC‑GTGACGCCAGAGATAGCGCTTTTTCTTGTGCCAGTCGAAGCACAAACTGT‑CGCGGCGATTCGTCGCTGCGTG‑GGGTTTCGTCGACATCTG‑GTGCTGCACATTCAAACGAGCTTTGc                                                                                             >  1:401815/5‑249 (MQ=255)                                                                                                                                                        aGAAGGTGACGATATTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATGCCTGGTCTGAACCAATAATTAAATGATCCGGATAAC‑GTGAc                                                                                                                                                                                                               <  1:136523/137‑1 (MQ=255)                                                                                                                                                        aGAAGGTGACGATATTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATGCCTGGTCTGAACCAATAATTAAATGATCCGGATAAC‑GTGAc                                                                                                                                                                                                               >  2:136523/1‑137 (MQ=255)                                                                                                                                                         gAAGGTGACGATATTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATGCCTGGTCTGAACCAATAATTAAATGATCCGGATAAC‑GTGACGCCAGAGATTGCGCTTTTTCTTGTGCCAGTCGAAGCACCAACTGT‑CGCGGTGATTCGTCGCTGCGTG‑GGGTTTCGTCGACCTCTG‑GTGCTGCACATTCAAAAGAGATTTGCAGtt                                                                                         <  1:425227/251‑1 (MQ=255)                                                                                                                                                              tGACGATATTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATGCCTGGTCTGAACCAATAATTAAATGATCCGGATAAC‑GTGACGCCAGAGATTGCGCTTTTTCTTGTGCCAGTCGAAGCACCAACTGT‑CGCGGTGATTCGTCGCTGCGTG‑GGGTTTCGTCGACCTCTG‑GTGCTGCACATTCAAAAGAGATTTGCAGTTTTTcc                                                                                    >  1:362368/1‑251 (MQ=255)                                                                                                                                                                 cGATATTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATGCCTGGTCTGAACCAATAATTAAATGATCCGGATAAC‑GTGACGCCAGAGATTGCGCTTTTTCTTGTGCCAGTCGAAGCACCAACTGT‑CGCGGTGATTCGTCGCTGCGTG‑GGGTTTCGTCGACCTCTG‑GTGCTGCACATTCAAAAGAGATTTGCAGTTTTTCCaga                                                                                 <  2:298088/251‑1 (MQ=255)                                                                                                                                                                                       tAATTGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATGCCTGGTCTGAACCAATAATTAAATGATCCGGATAAC‑GTGACGCCAGAGATTGCGCTTTTTCTTGTGCCAGTCGAAGCACCAACTGT‑CGCGGTGATTCGTCGCTGCGTG‑GGGtt                                                                                                                                     <  1:351027/178‑1 (MQ=255)                                                                                                                                                                                       tAATTGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATGCCTGGTCTGAACCAATAATTAAATGATCCGGATAAC‑GTGACGCCAGAGATTGCGCTTTTTCTTGTGCCAGTCGAAGCACCAACTGT‑CGCGGTGATTCGTCGCTGCGTG‑GGGtt                                                                                                                                     >  2:351027/1‑178 (MQ=255)                                                                                                                                                                                           tGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATGCCTGGTCTGAACCAATAATTAAATGATCCGGATAAC‑GTGACGCCAGAGATTGCGCTTTTTCTTGTGCCAGTCGAAGCACCAACTGT‑CGCGGTGATTCGTCGCTGCGTG‑GGGGTTCGCCGACCTCTG‑GAGCTGCACATTCAAAAGAGATTTGCAGTTTTTCCAGAAGAGCGCGACGCCATGGCGAgggggg                                                       >  2:28982/1‑247 (MQ=255)                                                                                                                                                                                            gCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATGCCTGGTCTGAACCAATAATTAAATGATCCGGATAAC‑GTGACGCCAGAGATTGCGCTTTTTCTTGTGCCAGTCGAAGCACCAACTGT‑CGCGGTGATTCGTCGCTGCGTG‑GGGTTTCGTCGACCTCTG‑GTGCTGCACATTCAAAAGAGATTTGCAGTTTTTCCAGAAGAGCGCGACGCCATGGCGAGGTgg                                                        >  1:235544/1‑249 (MQ=255)                                                                                                                                                                                             cAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATGCCTGGTCTGAACCAATAATTAAATGATCCGGATAAC‑GTGACGCCAGAGATTGCGCTTTTTCTTGTGCCAGTCGAAGCACCAACTGT‑CGCGGTGATTCGTCGCTGCGTG‑GGGTTTCGTCGACCTCTG‑GTGCTGCCCATTCAAAAGAGATTTGCAGTTTTTCCAGAAGCGCGCGACGCCATGGCGAGGTGGAgg                                                     >  2:139590/1‑251 (MQ=255)                                                                                                                                                                                             cAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATGCCTGGTCTGAACCAATAATTAAATGATCCGGATAAC‑GTGACGCCAGAGATTGCGCTTTTTCTTGTGCCAGTCGAAGCACCAACTGT‑CGCGGTGATTCGTCGCTGCGTG‑GGGTTTCGTCGACCTCTG‑GTGCTGCACATTCAAAAGAGATTTGCAGTTTTTCCAGAAGAGCGCGACGCCATGGCGAGGTGGAgg                                                     <  1:47811/251‑1 (MQ=255)                                                                                                                                                                                              aGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATGCCTGGTCTGAACCAATAATTAAATGATCCGGATAAC‑GTGACGCCAGAGATTGCGCTTTTTCTTGTGCCAGTCGAAGCACCAACTGT‑CGCGGTGATTCGTCGCTGCGTG‑GGGTTTCGTCGACCTCTG‑GTGCTGCACATTCAAAAGAGATTTGCAGTTTTTCCAGAAGAGCGCGACGCCATGGCGAGGTGGAGGc                                                    >  1:241694/1‑251 (MQ=255)                                                                                                                                                                                                  gggCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATGCCTGGTCTGAACCAATAATTAAATGATCCGGATAAC‑GTGACGCCAGAGATTGCGCTTTTTCTTGTGCCAGTCGAAGCACCAACTGT‑CGCGGTGATTCGTCGCTGCGTG‑GGGTTTCGTCGACCTCTG‑GTGCTGCACATTCAAAAGAGATTTGCAGTTTTTCCAGAAGAGCGCGACGCCATGGCGAGGTGGAGGCTaa                                                 >  1:275667/1‑250 (MQ=255)                                                                                                                                                                                                   ggCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACCCATGCCTGGTCTTAAACAATAATTAAATGATCCGGATAAC‑GTGACGCCAGAGATTGCGCTTTTTCTTGTGCCAGTCGAAGCACCAACTGT‑CGCGGTGATTCGTCGCTGCGTG‑GGGTTTCGTCGACCTCTG‑GTGCTGCACATTCAAAAGAGATTTGCAGTTTGTCCAGAAGAGCGCGACGCCATGGCGAGGTGGAGGCTaaaa                                               <  2:108196/251‑1 (MQ=255)                                                                                                                                                                                                   ggCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATGCCTGGTCTGAACCAATAATTAAATGATCCGGATAAC‑GTGACGCCAGAGATTGCGCTTTTTCTTGTGCCAGTCGAAGCACCAACTGT‑CGCGGTGATTCGTCGCTGCGTG‑GGGTTTCGTCGACCTCTG‑GTGCTGCACATTCAAAAGAGATTTGCAGTTTTTCCAGAAGAGCGCGACGCCATGGCGAGGTGGAGGCTaaaa                                               <  1:132554/251‑1 (MQ=255)                                                                                                                                                                                                            ttcCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATGCCTGGTCTGAACCAATAATTAAATGATCCGGATAAC‑GTGACGCCAGAGATTGCGCTTTTTCTTGTGCCAGTCGAAGCACCAACTGT‑CGCGGTGATTCGTCGCTGCGTG‑GGGTTTCGTCGACCTCTG‑GTGCTGCACATTCAAAAGAGATTTGCAGTTTTTCCAGAAGAGCGCGACGCCATGGCGAGGTGGAGGCTAAAATAAGttt                                        >  1:240442/1‑249 (MQ=255)                                                                                                                                                                                                             tcCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATGCCTGGTCTGAACCAATAATTAAATGATCCGGATAAC‑GTGACGCCAGAGATTGCGCTTTTTCTTGTGCCAGTCGAAGCACCAACTGT‑CGCGGTGATTCGTCGCTGCGTG‑GGGTTTCGTCGACCTCTG‑GTGCTGCACATTCAAAAGAGATTTGCAGTTTTTCCAGAAGAGCGCGACGCCATGGCGAGGTGGAGGCTAAAATAAGTTTa                                       <  1:138584/249‑1 (MQ=255)                                                                                                                                                                                                                gctAACGGGTTACCGGTGAATTCGCCATCAAGAACACATGCCTGGTCTGAACCAATAATTAAATGAGCCGGATAAC‑GTTACGCCAGAGATTTCGCTTTTTCTTGTGCCAGTCGAAGCACCTATTGT‑CGCGGTGATTCGACGCTGCGTG‑GGGTTTCGTCGACCTCTG‑GTGCTGCACGTTCAAAAGAGATTGGCAGTTTTTCCAGAAGAGCGCGACTCCATGTGTAGGTGGATGCTAAAATAAGTTTAGGCat                                  <  2:44304/249‑1 (MQ=255)                                                                                                                                                                                                                   aaCGGTTTACCGGTGATTTCGCCAGCAAGAACACATGCCTGGTCTTAACCAATAATTAAATGATCCGGATAAC‑GTGACGCCAGAGATTGCGCTTTTTCTTGTGCCAGTCGAAGCACCAACTGT‑CGCGGTGATTCGTCGCTGCGTG‑GGGTTTCGTCGACCTCTG‑GTGCTGCACATTCAAAAGAGATTTGCAGTTTTTCCAGAAGAGCGCGACGCCATGGCGAGGTGGAGGCTAAAATAAGTTTAGGCATAtt                               <  2:650803/251‑1 (MQ=255)                                                                                                                                                                                                                       gTTTACCGGTGCGTTCGCCACCACGAACACATGCCTGGTCTGAACCAATAATTAAATGAACCCGACAAC‑GTGGCGCCCGAGCTTGCGTATTTTCTTTAGCCAGTCGAAGCAACAGCTGT‑CGCGGTGATCCGGCGCTGCGTG‑GGGTTTCGTCGACCTGTG‑GTGCTGCACATTCAAAAGAGATTTGCAGTTTTTACAGAAGCGCGCGACGCCATCGCGAGGTGGAGGCTAAAATAAGTTTAGGCATAtttttt                           <  1:249556/251‑1 (MQ=255)                                                                                                                                                                                                                       gTTTACCGGTGATTTCGCCATCAAGAACACATGCCTGGTCTGAACCAATAATTAAATGATCCGGATAAC‑GTGACGCCAGAGATTGCGCTTTTTCTTGTGCCAGTCGAAGCACCAACTGT‑CGCGGTGATTCGTCGCTGCGTG‑CGGTTTCGTCGTCCTCTT‑GTGCTGCACATTCAAAAGAGATGTGCAGTTTTTCCAGAAGAGCGCGACGCCATgg                                                                <  1:323219/214‑1 (MQ=255)                                                                                                                                                                                                                         ttACCGGTGATTTCGCCATCAAGAACACATGCCTGGTCTGAACCAATAATTAAATGATCCGGATAAC‑GTGACGCCAGAGATTGCGCTTTTTCTTGTGCCAGTCGAAGCACCAACTGT‑CGCGGTGATTCGTCGCTGCGTG‑GGGTTTCGTCGACCTCTG‑GTGCTGCACATTCAAAAGAGATTTGCAGTTTTTCCAGAAGAGCGCGACGCCATGGCGAGGTGGAGGCTAAAATAAGTTTAGGAATATTTTTTTc                         >  2:444140/1‑251 (MQ=255)                                                                                                                                                                                                                                  aTTTCGCCATCAAGAACACATGCCTGGTCTGAACCAATAATTAAATGATCCGGATAAC‑GTGACGCCAGAGATTGCGCTTTTTCTTGTGCCAGTCGAAGCACCAACTGT‑CGCGGTGATTCGTCGCTGCGTG‑GGGTTTCGTCGACCTCTG‑GTGCTGCACATTCAAAAGAGATTTGCAGTTTTTCCAGAAGAGCGCGACGCCATGGCGAGGTGGAGGCTAAAATAAGATTAGGCATATTTTTTTCCATCAGat                 >  1:102935/1‑250 (MQ=255)                                                                                                                                                                                                                                       gCCATCAAGAACACATGCCTGGTCTGAACCAATAATTAAATGATCCGGATAAC‑GTGACGCCAGAGATTGCGCTTTTTCTTGTGCCAGTCGAAGCACCAACTGT‑CGCGGTGATTCGTCGCTGCGTG‑GGGTTTCGTCGACCTCTG‑GTGCTGCACATTCAAAAGAGATTTGCAGTTTTTCCAGAAGAGCGCGACGCCATGGCGAGGTGGAGGCTAAAATAAGTTTAGGCATATTTTTTTCCATCAGATATAGCg           >  2:154452/1‑251 (MQ=255)                                                                                                                                                                                                                                                aaCACATGCCTGGTCTGAACCAATAATTAAATGATCCGGATAAC‑GTGACGCCAGAGATTGCTCTTTTTCTTGTGCCAGTCGAAGCACCAACTGT‑CGCGGTGTTTCGTCGCTGCGTG‑GGGTTTCGTCGACCTCTG‑GTGCTGCACATTCAAAAGAGATTTGCAGTTTTTCCAGAAGAGCGCGACGCCATGGCGAGGTGGAGGCTAAAATAAGTTTAGGCATATTTTTTGCCATCATATATAGCGTATtgatga  <  2:341682/251‑1 (MQ=255)                                                                                                                                                                                                                                                aaCACATGCCTGGTCTGAACCAATAATTAAATGATCCGGATAAC‑GTGACGCCAGAGATTGCGCTTTTTCTTGTGCCAGTCGAAGCACCAACTGT‑CGCGGTGATTCGTCGCTGCGTG‑GGGTTTCGTCGACCTCTG‑GTGCTGCACATTCAAAAGAGATTTGCAGTTTTTCCAGAAGAGCGCGACGCCATGGCGAGGTGGAGGCTAAAATAAGTTTAGGCATATTTTTTTCCATCAGATATAGCGTATtgatg   >  1:382223/1‑250 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                        CTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAGATGCCCATTCGCCGAATTAAACAACGCCAGTCCAGTATAGAAGGTGACGATATTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATACCTGGTCTGAACCAATAATTAAATGATCCGGATAAC‑GTGACGCCAGAGATTGCGCTTTTTCTTGTGCCAGTCGAAGCACCAACTGT‑CGCGGTGATTCGTCGCTGCGTG‑GGGTTTCGTCGACCTCTG‑GTGCTGCACATTCAAAAGAGATTTGCAGTTTTTCCAGAAGAGCGCGACGCCATGGCGAGGTGGAGGCTAAAATAAGTTTAGGCATATTTTTTTCCATCAGATATAGCGTATTGATGA  >  NC_000913/1146138‑1146627

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