Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,146,110	C→T	E162E (GAG→GAA)	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,110	0	C	T	91.3%	55.5 / NA	24	E162E (GAG→GAA)	yceF	m(7)GTP pyrophosphatase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  major base T (8/13);  minor base A (2/0);  total (11/13)
Fisher's exact test for biased strand distribution p-value = 1.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGTAACGCGCGCTAATACAAGCATATAAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGCTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCC‑GCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAGATGCCCATTCGCCGAATTAAACAACGCCAGTCCAGTATAGAAGGTGACGATATTGCCGCTGGCTTTGCGTAATTGCAGACGGG  >  NC_000913/1145884‑1146333                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                 gCGTAACGCGCGCTAATACAAGCATATAAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAATcc                                                                                                                                                                                                           <  1:328422/250‑1 (MQ=255)                aTACAAGCATATAAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCt                                                                                                                                                                                            >  1:384669/1‑250 (MQ=255)                             aaacttaaaatggctgtttGGCCGGAGAAGGGGTTTACGCCGCACCAGGCCGACAGCACCGATGGCCGGATGCGACGTAACCTCACCCGGCGCAAAAAACTAGCCCATCTGCGGGTTTTTCCCCTCCCGACGTAACATCTGTCACAGCACGATTAGCGGTAAACCAACCAGCGTGTTAGTATCACGCCCCTATAAACGTTCAAACAGCGTAATGCCAAATCCTTCAATCTTAAAGTTACCC‑GCGCAGTGCa                                                                                                                                                                             <  2:128931/234‑1 (MQ=255)                                   atagtggacggatagCCTGATAAGACGTTTACGCAGCATCAGGCCACCAGCACCGATTGACGGATGCGACGTAACCACATGCGGCGCAAAACACTAACCCATCAGCGGGTTTTTCACTTCACGACGTAACATCTGACACAGCGCGATTAGCGGTAAAACAACCAGAGTGTTATGATCACGCCCCTCTAACCGATCAACCAGCGCAATGCCAAATCCTTCACTCTTAAAGCAACAC‑GCGCAGTGCAGGGGAt                                                                                                                                                                       >  2:22517/15‑251 (MQ=255)                                    tAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCC‑GCGCAGTGCAGGGGATg                                                                                                                                                                      >  1:47404/1‑251 (MQ=255)                                                 gCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCC‑GCGCAGTGCAGGGGATGCTCTTTACGCACg                                                                                                                                                         <  1:362486/251‑1 (MQ=255)                                                      aaaaGGCGTTTACGCCTCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAACTATCCCAACAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCC‑GCGCAGTGCAGGGGATGCTCTTTACGCACGTAAtt                                                                                                                                                    <  1:166974/249‑1 (MQ=255)                                                       tAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCC‑GCTCAGTGCAGGGGATGCTCTTTACGCACGTAAtt                                                                                                                                                    >  2:264510/1‑250 (MQ=255)                                                                tACGCCGCATCAGGCCGCTAGCACCGTTTCGCGGATGCGACGTAACACCATCCGGCGCAAAAAACTATCCCATCAGCGGTTTTTTCCCATACCGACGTACCATCTGACACAGCGCGATTAGCGGAAAACAAACCAGCGTGCTAGGATCAAGACTCTCCAAACAATCAAACAGCGAACTGCAAAATCCTTAACTATTAAAGC‑ACCCAGAGCCGTGCAGGAGAAGCGCTTTTCGCACGAAATTATCAATCTcc                                                                                                                                          >  2:119353/1‑251 (MQ=255)                                                                                gccAGCACCGATTGCCGTAGGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCTATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCC‑GCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGAc                                                                                                                          <  1:292527/251‑1 (MQ=255)                                                                                    gCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCC‑GCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGaa                                                                                                                      <  1:264510/251‑1 (MQ=255)                                                                                        tgATTGCCGGATGCGACGTAACAACATCCGGCGCAAAAAACTATCCCATCAGCGGGGTTTTCCCTTCCCGACGTAACATCTGACACAGCGAGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGACCCTCTCAACGTTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCC‑GCGCTGTGCAGGGGATGCTCTTTACGAACGTAATTATCAATCTCCGCCTCGCTCAGGAGACGGAAATgg                                                                                                                  >  2:39326/2‑251 (MQ=255)                                                                                              ccGGATGCGGCGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCC‑GCGCAGTGCAGGGGATGCTCTATACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGCAGGAAATGGACGTc                                                                                                             >  2:164014/1‑250 (MQ=255)                                                                                                   ggCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCC‑GCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAgg                                                                                                       <  1:206013/250‑1 (MQ=255)                                                                                                                       cgcAAAAAACTATCCCCTCAGCGGTTTTTTCCCTCCCCTCCGTTACATCTGCCACAGCTCGATTACCGGTAAACCAACCATCTCGTTAGGATTCCGCCCCTCTACACGTTCAAACAGCGTAATGCCATATCCCTCCCTCTTAAAGCTCCCC‑GCGCAGTTCAGGGGCTTCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTTACGGAAATGGACGTCaaacttctcctctactttttgcagat                                                                                   <  2:202767/251‑24 (MQ=255)                                                                                                                                                   ttcccttcccGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCC‑GCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAGATGCCCATTCGCCGAATTAAACAACGcca                                                        <  1:120561/250‑1 (MQ=255)                                                                                                                                                     cccttcacGACGTAACATCTGACACAGCGCGCTTCACAGGAACCCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGCTCAAACAGAGTAATGCCAAATACTTCACTCTTCAAGCTACCC‑GCGCAGTGCAGGAGATGCCCTTTACGCACGTACTTATCAATCTCCGCCTCGCTCCGGTGACGGAAATGGACGTCAAAAGGCACCACTTCTGTTTGCCTCTGCCCATTCGCCGAATTACACGACGccagtc                                                     >  2:112670/1‑251 (MQ=255)                                                                                                                                                          cccGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCC‑GCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAGATGCCCATTCGCCGAATTAAACAACGCCAGTCCAGta                                                <  1:133142/251‑1 (MQ=255)                                                                                                                                                            cGACGTAACATCTGACACAGCGCGATTAGAGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCC‑GCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAGATGCCCATTCGCCGAATTAAACAACGCCAGTCCAGtat                                               >  1:317221/1‑250 (MQ=255)                                                                                                                                                               cGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCC‑GCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAAtc                                                                                                                                             >  1:162288/1‑153 (MQ=255)                                                                                                                                                                         gACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCC‑GCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAGATGCCCATTCGCCGAATTAAACAACGCCAGTCCAGTATAGAAGGTGAc                                     <  1:193253/247‑1 (MQ=255)                                                                                                                                                                                  cgATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCC‑GCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAGATGCCCATTCGCCGAATTACACAACGCCAGTCCAGTATAGAAGGTGACGATATTGCCGCTg                        >  1:157185/1‑251 (MQ=255)                                                                                                                                                                                                     caGCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCC‑GCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAGATGCCCATTCGCCGAATTAAACAACGCCAGTCCAGTATAGAAGGTGACGATATTGCCGCTGGCTTTGCGTAATTGCAGAc     <  1:164014/251‑1 (MQ=255)                                                                                                                                                                                                        cGTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAATCCTTCATTCTTAAAGCTACCC‑GCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAGATGCCCATTCGCCGAATTAAACAACGCCAGTCCAGTATAGAAGGTGACGATATTGCCGCTGGCTTTGCGTAATTGCAGACggg  <  1:190723/251‑1 (MQ=255)                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                 GCGTAACGCGCGCTAATACAAGCATATAAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGCTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCC‑GCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAGATGCCCATTCGCCGAATTAAACAACGCCAGTCCAGTATAGAAGGTGACGATATTGCCGCTGGCTTTGCGTAATTGCAGACGGG  >  NC_000913/1145884‑1146333

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