Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	6,574	A→G	G462G (GGT→GGC)	yaaJ ←	putative transporter YaaJ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	6,574	0	A	G	95.8%	61.6 / NA	24	G462G (GGT→GGC)	yaaJ	putative transporter YaaJ
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (7/16);  minor base T (0/1);  total (7/17)
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.96e-01

GCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATATCCAGTCCTTGCAGGAAATTTATGCCGACTTTAGCAAAAAATGAGAATGAGTTGATCGATAGTTGTGATTACTCCTGCGAAACATCATCCCACGCGTCCGGAGAAAGCTGGCGACCGATATCCGGATAACGCAATGGATCAAACACCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGGCTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTTTAG  >  NC_000913/6363‑6812                                                                                                                                                                                                                    |                                                                                                                                                                                                                                               gctgGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATATCCAGTCCTTGCAGGAAATTTATGCCGACTTTAGCAAAAAAAGAGAATGAGTTGATCGATAGTTGTGATTACTCCTGCGAAACATCATCCCACGCGTCCGGAGAAAGCTGGCGGCCGATATCCGGATAACGCAACGGATCAAACACCGGgcg                                                                                                                                                                                                          <  1:232199/250‑1 (MQ=255)               aaCAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATATCCAGTCCTTGCAGGAAATTTATGCCGACTTTAGCAAAAAAAGAGAATGAGTTGATCGATAGTTGTGATTACTCCTGCGAAACATCATCCCACGCGTCCGGAGAAAGCTGGCGGCCGATAGCCGGATAACGCAACGGATCAAACACCGGGCGCACGCCGAGTTTAc                                                                                                                                                                                            <  2:89127/250‑1 (MQ=255)                                    cgcgTGGTGGTCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATATCCAGTCCTTGCAGGAAATTTATGCCGACTTTAGCAAAAAAAGAGAATGAGTTGATCGATAGTTGTGATTACTCCTGCGAAACATCATCCCACGCGTCCGGAGAAAGCTGGCGGCCGATATCCGGATAACGCAACGGATCAAACACCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGc                                                                                                                                                                      <  2:131583/251‑1 (MQ=255)                                                    ggCTTTGGTAATCAAGCGGTTTCGCAGGTGAAATAAGAATCAGCATATCCAGTCCTTGCAGGAAATTTATGCGGACTTTAGCAAAAAAAGAGAATGAGTTGATCGATAGTTGTGATTACTCCTTCGAAACATCATCCCACGCGTCCGGAGAAAGCTGGCGGCCGATATCCGGATAACGCAACGGATCAAACACCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACcaca                                                                                                                                                      <  1:71937/251‑1 (MQ=255)                                                             aagcAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATATCCAGTCCTTGCAGGAAATTTATGCCGACTTTAGCAAAAAAAGAGAATGAGTTGATCGATAGTTGTGATTACTCCTGCGAAACATCATCCCACGCGTCCGGAGAAAGCTGGCGTCCGATATCCGGATAACGCAACGGATCAAACACCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACATGCGAGAGc                                                                                                                                             <  2:103399/248‑1 (MQ=255)                                                                       tttCGCAGGTGAAATAAGAATCAGCATATCCAGTCCTTGCAGGAAATTTATGCCGACTTTAGCAAAAAAAGAGAATGAGTTGATCGATAGTTGTGATTACTCCTGCGAAACATCATCCCACGCGTCCGGAGAAAGCTGGCGGCCGATATCCGGATAACGCAACGGATCAAACACCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATg                                                                                                                                    <  1:138093/250‑1 (MQ=255)                                                                         tCGCAGGTGAAATAAGAATCAGCATATCCAGTCCTTGCAGGAAATTTATGCCGACTTTAGCAAAAAAAGAGAATGAGTTGATCGATAGTTGTGATTACTCCTGCGAAACATCATCCCACGCGTCCGGAGAAAGCTGGCGGCCGATATCCGGATAACGCAACGGATCAAACACCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGc                                                                                                                                  <  1:86065/250‑1 (MQ=255)                                                                                 gAAATAAGAATCAGCATATCCAGTCCTTGCAGGAAATTTATGCCGACTTTAGCAAAAAAAGAGAATGAGTTGATCGATAGTTGTGATTACTCCTGCGAAACATCATCCCACGCGTCCGGAGAAAGCTGGCGGCCGATATCCGGATAACGCAACGGATCAAACACCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCCAAtt                                                                                                                         >  1:139765/1‑251 (MQ=255)                                                                                       aGAATCAGCATATCCAGTCCTTGCAGGAAATTTATGCCGACTTTAGGAAAAAAAGAGAATGAGTGGATCGATAGTGGTGATTACTCCTGCGAAACATCAGCCCACGCGTCCGGAGAAAGCTGGCGGCCGATATCCGGATAACTCAACGGATCAAACACCGGTTTCATGCCGAGTGTAGTCTGTCGTAGATAATCACTGGCAATTGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTAAt                                                                                                                   <  1:155705/251‑1 (MQ=255)                                                                                              gCATATCCAGTCCTTGCAGGAAATTTATGCCGACTTTAGCAAAAAAAGAGAATGAGTTGATCGATAGTTGTGATTACTCCTGCGAAACATCATCCCACGCGTCCGGAGAAAGCTGGCGGCCGATATCCGGATAACGCCACGGATCAAACACCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTAAGAGcgga                                                                                                             >  1:8663/1‑247 (MQ=255)                                                                                              gCATATCCAGTCCGTGCAGGAAATTTATGCCGACTTTAGCAAAAAAAGAGAATGAGTTGATCGATAGTTGTGATTACTCCTGCGAAACATCATCCCACGCGTCCGGAGAAAGCTGGCGGCCGATATCCGGATAACGCAACGGATCAAACACCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTa                                                                                                                     <  2:8663/242‑1 (MQ=255)                                                                                                  atCCAGTCCTTGCAGGAAATTTATGCCGACTTTAGCAAAAAAAGAGAATGAGTTGATCGATAGTTGTGATTACTCCTGCGAAACATCATCCCACGCGTCCGGAGAAAGCTGGCGGCCGATATCCGGATAACGCAACGGATCAAACACCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAg                                                                                                         <  1:26587/250‑1 (MQ=255)                                                                                                     cAGTCCTTGCAGGAAATTTATGCCGACTTTAGCAAAACAAGAGAATGAGTTGATCGATAGTTGTGATTACTCCTGCGAAACATCATCCCACGCGTCCGGAGAAAGCTGGCGGCCGAAATCCGGATAACGCAACGGATCAACCACCGGGCGCACGCCGAGCTTACGCTGCCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGAAGTAAAATGGCGGTCAAATTGGTAAAAGCCATGCCGGcc                                                                                                      >  1:135989/1‑250 (MQ=255)                                                                                                                                     cAAAAAAAGAGAATCAGTTGATCGATAGTTCTGATTACTCCTGCGAAACATCATCCCACGCGTCCGGAGAAAGCTGGCGGCCGATATCCGGATAACGCAACGGATCAAACACCGGGCGCACGCCGAGTTTCCGCTGGCGTCGATAATCACTGGCAATGGTATGAACCACAGGCGCGAGCAGTAAAATGGCGGTCAAATTGGTAATCGCCAGGAAGGCCATTATGGTATCTGCCAGTTGCCACATCAGCgg                                                                      >  1:188198/1‑250 (MQ=255)                                                                                                                                      aaaaaaAGAGAATGAGTTGACAGATAGTTGTGATTACTCCTGCGCAACATCATCCCAAGCGTCCGGAGAAAGCTGACCGCCGATATCCGGATAACGCAACGGATCAAACCCCGGGAGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTCTGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTTGGAAAAGAAATGCAGGCCATTATCATATcatcccgttgcccaaacatcgga                                                                     >  2:139078/1‑228 (MQ=255)                                                                                                                                          aaagagaatgagTTGATCGATAGTTGTGATTACTCCTGCGAAACATCATCCCACGCGTCCGGAGAAAGCTGGCGGCCGATATCCGGATAACGCAACGGATCAAACACCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAgact                                                                <  1:149819/248‑4 (MQ=255)                                                                                                                                                        gATCGATAGTTGTGATTACTCCTGCGAAACATCATCCCACGCGTCCGGAGAAAGCTGGCGGCCGATATCCGGATAACGCAACGGATCAAACACCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGACTTAGCAAGGTgccg                                                   <  1:69069/250‑1 (MQ=255)                                                                                                                                                                  tgtgATTACTCCTGCGAAACATCATCCCACGCGTCCGGAGAAAGCTGGCGGCCGATATCCGGATAACGCAACGGATCAAACACCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCg                                         <  1:77783/250‑1 (MQ=255)                                                                                                                                                                          cTCCTGCGAAACATCATCCCACGCGTCCGGAGAAAGCTGGCGGCCGATATCCGGATAACGCAACGGATCAAACACCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGAGAAgg                                >  1:69806/1‑251 (MQ=255)                                                                                                                                                                                               cgcgTCCGGAGAAAGCTGGCGGCCGATATCCGGATAACGCAACGGATCAAACACCGGGCGCACGCCGAGTTTACGCTGGCGTAGATACTCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACCTCAGCGGAAGACTTAGCAAGGGGCCGCCGATGAACGTTGCGAAGGGGCAGATcaccaaccacc              >  2:33471/1‑238 (MQ=255)                                                                                                                                                                                                 cgTCCGGAGAAAGCTGGCGGCCGATATCCGGATAACGCAACGGATCAAACACCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATc         <  1:78559/251‑1 (MQ=255)                                                                                                                                                                                                        agaAAGCTGGCGGCCGATATCCGGATAACGCAACGGATCAAACACCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTTTAg  <  2:139765/251‑1 (MQ=255)                                                                                                                                                                                                        agaAAGCTGGCGGCCGATAACCGGAGAACGCAACGGTTCAAACACCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCCATTATGATATCGGCCAGTTTCCACATCAGCGGAAGACTGAGCAAGGTGCCGCCGATGTGCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTGTAg  <  2:188198/251‑1 (MQ=255)                                                                                                                                                                                                         gaAAGCTGGCGGCCGATATCCGGATAACGCAACGGATCAAACACCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCATGCCATTATGatat                                                                                            <  1:191265/160‑1 (MQ=255)                                                                                                                                                                                                                    |                                                                                                                                                                                                                                               GCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATATCCAGTCCTTGCAGGAAATTTATGCCGACTTTAGCAAAAAATGAGAATGAGTTGATCGATAGTTGTGATTACTCCTGCGAAACATCATCCCACGCGTCCGGAGAAAGCTGGCGACCGATATCCGGATAACGCAATGGATCAAACACCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGGCTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTTTAG  >  NC_000913/6363‑6812

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