Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,346,888	G→T	14.2%	T31T (ACG→ACT)	ymiB →	uncharacterized protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,346,888	0	G	T	14.2%	91.3 / 4.4	33	T31T (ACG→ACT)	ymiB	uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (20/8);  new base T (3/2);  total (23/10)
Fisher's exact test for biased strand distribution p-value = 6.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGTTTTTTTTATATATTTATTTGTAATCCAGTTTTGGAAAAACGCCAGTTTTCAAACGAAAGTCAGTTAAAAAATCTGCCTGGATATAACGAAGGTAGAGCGGGGAAATAAACGGCCCATCCATGAGGAATGGGCCGTGAAAGGAGATTACGCGACCGGGCGCGCAATAATGCTGCGGGTTTCCATGCGGAC  >  NC_000913/1346756‑1347022                                                                                                                                     |                                                                                                                                       cGGTTATCTGCGGAGTCGGCGGTTATATTAAAGGAACATCATGACTTACATGTATTTAACGGCGGCTTAACGCTTTGTTTTTTTTATATATTTATTTGTAATCCTGTTTTGGATAAACGCCAGTTTTCAAACGAAAGTCa                                                                                                                                 <  2:940069/140‑1 (MQ=255) cGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGTTTTTTTTATATATTTATTTGTAATCCAGTTTTGGAAATACGCCAGTTTTCAAACGAAATTCa                                                                                                                                 >  1:1299889/1‑140 (MQ=255) cGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGTTTTTTTTATATATTTATTTGTAATCCAGTTTTGGAAAAACGCCAGTTTTCAAACGAAAGTCa                                                                                                                                 >  2:999661/1‑140 (MQ=255) cGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGTTTTTTTTATATATTTATTTGTAATCCAGTTTTGGAAAAACGCCAGTTTTCAAACGAAAGTCa                                                                                                                                 >  2:980716/1‑140 (MQ=255) cGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGTTTTTTTTATATATTTATTTGTAATCCAGTTTTGGAAAAACGCCAGTTTTCAAACGAAAGTCa                                                                                                                                 >  2:979716/1‑140 (MQ=255) cGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGTTTTTTTTATATATTTATTTGTAATCCAGTTTTGGAAAAACGCCAGTTTTCAAACGAAAGTCa                                                                                                                                 >  2:565540/1‑140 (MQ=255) cGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGTTTTTTTTATATATTTATTTGTAATCCAGTTTTGGAAAAACGCCAGTTTTCAAACGAAAGTCa                                                                                                                                 >  2:1457152/1‑140 (MQ=255)        cAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGTTTTTTTTATATATTTATTTGTAATCCAGTTTTGGAAAAACGCCAGTTTTCAAACGAAAGTCAGTTaaaa                                                                                                                          >  1:1678559/1‑140 (MQ=255)             gAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGTTTTTTTTATATATTTATTTGTAATCCAGTTTTGGAAAAACGCCAGTTTTCAAACGAAAGTCAGTTAAAAAATCt                                                                                                                     >  1:48482/1‑140 (MQ=255)                  cgcgGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAATCGCATTGTTTTGTTTATAGATTTATTTGTAATCCAGTTTTGGAAAAACGCCAGTTTTCAAACTAAAGTCAGTTAAAAAATCTGCCTg                                                                                                                >  2:333398/1‑140 (MQ=255)                      gTTAAATAAAAGGAACAAAATGACATACATGAAATTAACGGCGGCTAAACGCATTGTTTTTTTTATATATTTATTTGAAATCCAGTATTGGAAAAACGCCAGTTTTCAAACGAAAGTCCGTGAAAAAATCTGCCTGGata                                                                                                            >  2:287854/1‑140 (MQ=255)                          aataaaAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGTTTTTTTTATATATTTATTTGTAATCCAGTTTTGGAAAAACGCCAGTTTTCAAACGAAAGTCAGTTAAAAAATCTGCCTGGATATAAc                                                                                                        >  1:1503773/1‑140 (MQ=255)                                 ggAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGTTTTTTTTATATATTTATTTGTAATCCAGTTTTGGAAAAACGCCAGTTTTCAAACGAAAGTCAGTTAAAAAATCTGCCTGGATATAACGAAGGTa                                                                                                 >  1:846201/1‑140 (MQ=255)                                      aacaTGACTTACATGTAATTAACGGCGGCTAAACGCATTGTTTTTTTTATATATTTATTTGTAATCCAGTTTTGGAAAAACGCCAGTTTTCAAACGAAAGTCAGTTAAAAAATCTGCCTGGATATAACGAAGGTAGAGCg                                                                                            <  2:873967/140‑1 (MQ=255)                                             cTTACATGAAATTAACGGCGGCTAAACGCATTGTTTTTTTTATATTTTTATTAGTAATCCAGTTTTGGAAAAACGCCAGTTTTCAAACGAAAGTCAGATAAAAAATCTGCCTGGATATACCGATTGTAGAGCGGGGCAAt                                                                                     >  2:1104708/1‑140 (MQ=255)                                             cTTACATGAAATTAACGGCGGCTAAACGCATTGTTTTTTTTATATATTTATTTGTAATCCAGTTTTGGAAAAACGCCAGTTTTCAAACTAAAGTCAGTTAAAAAATCTGCCTGGATATAACGAAGGTAGAGCGGGGAAAt                                                                                     >  1:1108514/1‑140 (MQ=255)                                             cTTACATGAAATTAACGGCGGCTAAACGCATTGTTTTTTTTATATATTTATTTGTAATCCAGTTTTGGAAAAACGCCAGTTTTCAAACTAAAGTCAGTTAAAAAATCTGCCTGGATATAACGAAGGTAGAGCGGGGAAAt                                                                                     >  1:1108319/1‑140 (MQ=255)                                             cTTACATGAAATTAACGGCGGCTAAACGCATTGTTTTTTTTATATATTTATTTGTAATCCAGTTTTGGAAAAACGCCAGTTTTCAAACGAAAGTCAGTTAAAAAATCTGCCTGGATATAACGAAGGTAGAGCGGGGAAAt                                                                                     >  2:1324835/1‑140 (MQ=255)                                                    gAAATTAACGGCGGCTAAACGCATTGTTTTTTTTATATATTTATTTGTAATCCAGTTTTGGAAAAACGCCAGTTTTCAAACGAAAGTCAGTTAAAAAATCTGCCTGGATATAACGAAGGTAGAGCGGGGAAATAAACGGc                                                                              <  1:1457152/140‑1 (MQ=255)                                                       aTTAACGGCGGCTAAACGCATTGTTTTTTTTATATAATTATTTGTAATCCAGTTTTGGAAAAACGCCAGTTTTCAAAAGAAATTCAGTTAAAAAATATGCCTGGATATCACGAAGGTAGAGAGGGGACATAAACGGCcca                                                                           >  1:90144/1‑140 (MQ=255)                                                       aTTAACGGCGGCTAAACGCATTGTTTTTTGTATATATTTATTTGTAATCCAGTTTTGGAAAAACGCCAGTTTTCAAACGAAAGTCAGTGAAAAAATCTGCCTGGATATAACGAAGGTAGAGCGGGGAAATAAACGGCcca                                                                           >  1:1263332/1‑140 (MQ=255)                                                            cggcggCTAAACGCATTGTTTTTTTTATATATTTATTTGTAATCCAGTTTTGGAAAAACGCCAGTTTTCAAACGAAAGTCAGTTAAAAAATCTGCCTGGATATAACGAAGGTAGAGCGGGGAAATAAACGGCccatccat                                                                      >  1:947452/1‑140 (MQ=255)                                                            cgcccgcTAAACGCATTGTTTTTTTTATATATTTATTTGAAATCCAGTTTTGCAAAAACGCCAGTTTTCAAACGAAAGTCAGTTAAAAAATCTGCCTGGATATAACGAAGGTAGAGCGGGGAAATAAACGGCccatccat                                                                      <  1:1413679/135‑1 (MQ=255)                                                                            tGTTTTTTTTATATATTTATTTGTAATCCAGTTTTGGAAAAACGCCAGTTTTCAAACGAAAGTCAGTTAAAAAATCTGCCTGGATATAAc                                                                                                        >  1:876891/1‑90 (MQ=255)                                                                            tGTTTTTTTTATATATTTATTTGTAATCCAGTTTTGGAAAAACGCCAGTTTTCAAACGAAAGTCAGTTAAAAAATCTGCCTGGATATAAc                                                                                                        <  2:876891/90‑1 (MQ=255)                                                                            tGTTTTTTTTATATATTTATTTGTAATCCAGTTTTGGAAAAACGCCAGTTTTCAAACGAAAGTCAGTTAAAAAATCTGCCTGGATATAACGAAGGTAGAGCGGGGAAATAAACGGCCCATCCATGAGGAATGGGCCGTGa                                                      >  1:682213/1‑140 (MQ=255)                                                                                      atatatTTATTTGTAATCCAGTTTTGGAAAAACGCCAGTTTTCAAACGAAAGTCAGTTAAAAAATCTGCCTGGATATAACGAAGGTAGAGCGGGGAAATAAACGGCCCATCCATGAGGAATGGGCCGTGAAAGGAGATTa                                            >  2:1550840/1‑140 (MQ=255)                                                                                      atatatTTATTTGTAATCCAGTTTTGGAAAAACGCCAGTTTTCAAACGAAAGTCAGTTAAAAAATCTGCCTGGATATAACGAAGGTAGAGCGGGGAAATAAACGGCCCATCCATGAGGAATGGGCCGTGAAAGGAGATTa                                            >  2:1299210/1‑140 (MQ=255)                                                                                               tttGTAATCCAGTTTTGGAAAAACGCCAGTTTTCAAACTAAAGTCAGTTAAAAAATCTGCCTGGATATAACGAAGGTAGAGCGGGGAAATAAACGGCCCATCCATGAGGAATGGGCCGTGAAAGGAGATTACGCGACCgg                                   <  2:1108319/140‑1 (MQ=255)                                                                                               tttGTAATCCAGTTTTGGAAAAACGCCAGTTTTCAAACTAAAGTCAGTTAAAAAATCTGCCTGGATATAACGAAGGTAGAGCGGGGAAATAAACGGCCCATCCATGAGGAATGGGCCGTGAAAGGAGATTACGCGACCgg                                   <  2:1108514/140‑1 (MQ=255)                                                                                                                                cAAACGAAAGTCAGTTAAAAAATCTGCCTGGATATAACGAAGGTAGAGCGGGGAAATAAACGGCCCATCCATGAGGAATGGGCCGTGAAAGGAGATTACGCGACCGGGCGCGCAATACTGCTGCGGGTTTCCATGCGGAc  <  2:1025877/140‑1 (MQ=255)                                                                                                                                cAAACGAAAGTCAGTTAAAAAATCTGCCTGGATATAACGAAGGTAGAGCGGGGAAATAAACGGCCCATCCATGAGGAATGGGCCGTGAAAGGAGATTACGCGACCGGGCGCGCAATAATGCTGCGGGTTTCCATGCGGAc  <  2:1678559/140‑1 (MQ=255)                                                                                                                                cAAACGAAAGTCAGTTAAAAAATCTGCCTGGATATAACGAAGGTAGAGCGGGGAAATAAACGGCCCATCCATGAGGAATGGGCCGTGAAAGGAGATTACGCGACCGGGCGCGCAATAATGCTGCGGGTTTCCATGCGGAc  <  2:940120/140‑1 (MQ=255)                                                                                                                                     |                                                                                                                                       CGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGTTTTTTTTATATATTTATTTGTAATCCAGTTTTGGAAAAACGCCAGTTTTCAAACGAAAGTCAGTTAAAAAATCTGCCTGGATATAACGAAGGTAGAGCGGGGAAATAAACGGCCCATCCATGAGGAATGGGCCGTGAAAGGAGATTACGCGACCGGGCGCGCAATAATGCTGCGGGTTTCCATGCGGAC  >  NC_000913/1346756‑1347022

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑