Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,906,163	G→A	G159D (GGC→GAC)	mntP →	Mn(2(+)) exporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,906,163	0	G	A	87.5%	29.0 / ‑2.8	16	G159D (GGC→GAC)	mntP	Mn(2(+)) exporter
Reads supporting (aligned to +/- strand):  ref base G (2/0);  new base A (6/8);  total (8/8)
Fisher's exact test for biased strand distribution p-value = 4.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.52e-01

GCGGGCGAATGATTATTGAGGGTTTTCGTGGCGCAGATGATGAAGATGAAGAGCCGCGCCGTCGACACGGTTTCTGGCTACTGGTAACCACCGCGATTGCCACCAGCCTGGATGCCATGGCTGTGGGTGTTGGTCTTGCTTTCCTGCAGGTCAACATTATCGCGACCGCATTGGCCATTGGTTGTGCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTGCCAGACTTCTGGCTTCGCACGCCAGGCAAACGGCGTCATCTGCAATAATGCGAC  >  NC_000913/1905929‑1906379                                                                                                                                                                                                                                           |                                                                                                                                                                                                                         gcgggcgAATGATTATTGAGGGGTTTCGTGGCGCAGATGATGAAGATGAAGAGCCGCGACGTCGACCCGGATTATGGCCACCGGTAAACCCCGCGATTGCCACAAACCTGGATGCCATGGCTGTGGGTGTTGGTCTTGCTTTCCTGCAGGGCAACATTATCCCGACCGCATTGGCCATTGGTTGTGCAACCTTGCTTATGTCAACATTAGGGGTGATGGTTGGTCGCTTTTTCGACTCCATTATTGGGaa                                                                                                                                                                                                           >  1:123698/1‑250 (MQ=255)                                                                                                                        gCTGTGGGTGTTGGTCTTGCTTTCCTGCAGGTCAACATTATCGCGACCGCATTGGCCATTGGTTGTGCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGACTCAATTATTGGGAAAAAAGCGGAAATTCTCGTTGGGCTGGTGCTGATCGGCATCGGCGGCCAGATcc                                                                                                                                                      <  1:275933/184‑1 (MQ=255)                                                                                                                                         tGCTTTCCTGCAGGTCAACATTATCGCGACCGCATTGGCCATTGGTTGTGCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGACTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTtgn                                                                    <  1:86058/249‑2 (MQ=255)                                                                                                                                                        cAACATTATCGCGACCGCATTGGCCATTGGTTGTGCAACCTTGATTATGTAAACATTAGGGCTGATAATTGGTCGATTTATCGACTCAACTATTAGGAAAAAAGCGGAACTTCTCGGAGGGCTGGCGCTGCTCGACATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGACTATTAAAGTCCGTCTGGCAGTAGAACACTTCTTTTGCTGCAAGTGTTTGcc                                                      >  1:240920/1‑248 (MQ=255)                                                                                                                                                           cATTATCGCGACCGCATTGGCCATTGGTTGTGCAACCTTGAGTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGACTCAATTATTGTGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGatcg                                                                                                                                                                        <  1:312770/131‑1 (MQ=255)                                                                                                                                                                               ncATTGGTTGTGCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCCGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTGCCAGACTTCTGGCTTCGCACGCCAGGc                             >  1:203596/2‑250 (MQ=255)                                                                                                                                                                                     gTTGTGCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCggn                                                                                                                                                                                                                         <  1:191457/56‑2 (MQ=255)                                                                                                                                                                                     gTTGTGCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGc                                                                                                                                                                                                                         >  2:191457/1‑56 (MQ=255)                                                                                                                                                                                       tgtgCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGACTCAAGTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTGCCAGACTTCTGGCTTCGCACGCCAGGCAAACGGCg                     <  1:19319/250‑1 (MQ=255)                                                                                                                                                                                           cAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGACTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGTTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTGc                                                       >  1:317002/1‑212 (MQ=255)                                                                                                                                                                                               cTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCGTTATCGACTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTgn                                                                                                                                                 <  1:81284/118‑2 (MQ=255)                                                                                                                                                                                                     tATGTNNACNNNAGNNATGATGGTTGGTCGCTTTATCGANTCAATTATTGGGAAAAAGGCGGAAATTCTngncn                                                                                                                                                                                       >  2:235754/1‑73 (MQ=255)                                                                                                                                                                                                     tATGTCCACATTAGGGATGATGGTTGGTCGCTTTATCGACTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGaa                                                                                                         >  1:276539/1‑152 (MQ=255)                                                                                                                                                                                                     tATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGACTCAATTATTGGGAAAAAAGCGGAAATTCTcggcg                                                                                                                                                                                       <  1:235754/74‑1 (MQ=255)                                                                                                                                                                                                          cAACATTAGGGATGATGGTTGGTCGCTTTATCGACTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCGCGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCGGTCGAACGCTTCTTTTGCTTCCAGTGTTTGCCAGACTTCTGGCGTCGCACGCCAGGCAAACGGCGTCATCTGCAATAATGCGAc  <  1:184888/250‑1 (MQ=255)                                                                                                                                                                                                           aaCATTAGGGATGATGGTTGGTCGCTTGATCGACTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATcc                                                                                                                                                      <  1:151480/101‑1 (MQ=255)                                                                                                                                                                                                                    gatgatGGTTGGTCGCTTTATCGACTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGAGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTGCCAGACTTCTGGCTTCGCACGCCAGGCAAACGGCGTCa                  >  1:163995/1‑224 (MQ=255)                                                                                                                                                                                                                                           |                                                                                                                                                                                                                         GCGGGCGAATGATTATTGAGGGTTTTCGTGGCGCAGATGATGAAGATGAAGAGCCGCGCCGTCGACACGGTTTCTGGCTACTGGTAACCACCGCGATTGCCACCAGCCTGGATGCCATGGCTGTGGGTGTTGGTCTTGCTTTCCTGCAGGTCAACATTATCGCGACCGCATTGGCCATTGGTTGTGCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTGCCAGACTTCTGGCTTCGCACGCCAGGCAAACGGCGTCATCTGCAATAATGCGAC  >  NC_000913/1905929‑1906379

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