Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,575,673 G→A intergenic (+9/+48) yhgN → / ← gntU MarC family putative inner membrane protein YhgN/low‑affinity gluconate transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,575,6730GA86.7% 25.5 / ‑1.2 15intergenic (+9/+48)yhgN/gntUMarC family putative inner membrane protein YhgN/low‑affinity gluconate transporter
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (5/8);  total (6/9)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.31e-01

CGCTGATGTTGTTGTCTCATCAGTACCCGAATCAGATGGGGCATCTGGTGATTGCTCTGCTGCTGGCCTGGGGCGGCACCTTTGTCATCCTGCTACAGTCTTCGCTATTTTTACGTCTGCTGGGCGAGAAAGGGGTGAACGCACTTGAACGCCTGATGGGATTGATTCTGGTGATGATGGCAACCCAGATGTTCCTCGACGGCATTCGAATGTGGATGAAGGGGTAACGAGGAGAGTTTTTCGAAAAGAACCCGGTAGCGTCACAACTACCGGGCGAACAAACTTAACTCAACAGCTGGAACGCAATCATCCCAACGATGGCACCGACAGTGCCGAGGATGGTTTCCATCATGGTCCAGGTTTTCAGCGTTTCGGCTTCGGTCGCGCCGGTAAATTTACCGAACAACCAGAAACCGGCGTCGTTAACGTGGCTGACAACAATC  >  NC_000913/3575438‑3575880
                                                                                                                                                                                                                                           |                                                                                                                                                                                                               
cGCTGATGTTGTTGTCTCATCAGTACCCGAATCAGATGGGGCATCTGGTGATTGCTCTGCTGCTGGCCTGGGGCGGCACCTTTGTCATCCTGCTACAGTCTTCGCTATTTTTACGTCTGCTGGGCGAGAAAGGGGTGAACGCACTTGAACGCCTGATGGGATTGATTCTGGTGATGATGGCAACCCAGATGTTCCTCGACGGCATTCGAATGTGGATGAAGGGGTAACGAGGAGAATTTTTCGAAAAGAAc                                                                                                                                                                                                  >  1:114259/1‑251 (MQ=255)
                                                                      gggCGGCACCTTTGTCATCCTGCTACAGTCTTCGCTATTTTTACGTCTGCTGGGCGAGAAAGGGGTGAACGCACTTGAACGCCTGATGGGATTGATTCTGGTGATGATGGCAACCCAGATGTTCCTCGACGGCATTCGAATGTGGATGAAGGGGTAACGAGGAGAATTTTTCGAAAAGAACCCGGTAGCGTCACAACTACCGGGCGAACAAACTTAACTCAACAGCTGGAACGCAATCATCCCAACGATgg                                                                                                                            <  1:123168/251‑1 (MQ=255)
                                                                      gggCGGCACCTTTGTCATCCTGCTACAGTATTCGCTATTTTTACGTCTGCTGGGCGAGAAAGGGGTGAACGCACTTGAACGCCTGATGGGATTGATTCTGGTGATGATGGCAACCCAGATGTTCCTCGACGGCATTCGAATGTGGATGAAGGGGTAACGAGGAGAATTTTTCGAAAAGAACCCGGTAGCGTCACAAc                                                                                                                                                                                  >  1:169413/1‑197 (MQ=255)
                                                                        gcggcACCGTTGTCATCCTTCTACAGTCTGCGTTATTTGTACGTCTGCTGGGCGAGAAAGGGGTGAACACTCTTTACCGCCTGATGTGATTGATTCTGGTGATTATGGCAATTCTGATTGTCGTCGACGGCTTTCGTATGTGGATGAAGGGGTAACGAGGAGAATTTTTCGaaaan                                                                                                                                                                                                     <  1:34995/176‑2 (MQ=255)
                                                                                               cAGTCTTCGCTATTTTTACGTCTGCTGGGCGAGAAAGGTGTGATCGCACTTGAACGCCTGATGGGATTGATTCTGGAGATGATGTCAACCCAGATGTTCCTCGACGGCATTCGAATGTGGATGAAGGGGTAACGAGGAGAGTTTTTCGAAAAGAACCCGGTAGCGTCACAACTACCGGGCGAACAGACTTAACTCATCATCTGGAACGCAATCTTCCCAACGATGGCACCGACAGTGCCGAGGATGGTTTc                                                                                                   <  1:144537/251‑1 (MQ=255)
                                                                                                 gTCTTCGGTGTTTGTACGTCTGTTGGGGGAGAAAGGGGAGAACGCACTTGAACTTTTGATGGGGTTGGTTTTGGTGATGATGCTAAAGCATATGTGCCGCGAGTGTATTCGAATGTGGATGAAGGGGTAACGAGGTGAATTTTTCGAAAAGAAc                                                                                                                                                                                                  <  1:49596/154‑1 (MQ=255)
                                                                                                      cGCTATTTTTACGTCTGCTGGGCGCGAACGGGGTGAACGCACTTGAACGCCTGATGGGATTGATTCTGGTGATGATGGCAACCCAGATGTTCCTCGACGGCATTCGAATGTGGATGAAGGGGTAACGAGGAGAGTTTTTCGAAAAGAACCCGGTAGCGTCACAACTACCGGGCGAACAAACTTAACTCAACAGCTg                                                                                                                                                   >  1:115695/1‑196 (MQ=255)
                                                                                                                                        gAACGCACTTGAACGCCTGATGGGCTTGATTCTGGTGATATTGGCAACCCAGATGTTCATCGGCGGCATTCGAATGTGGCAGAAGGGCAACCGAGGAGAATTTTTCGAAAATCACCCGGTAGCGTCCTAACTCCCGGGCGAACAAACATAACTCAACAGCTGGAACGCAATCATCCCAACGATGGCACCGACAGTGCCGCGGATGGTTTCACTCATGGTCAAGGTTTTCAGCGTTTCGGCTTCGGTCgcgc                                                          >  1:54670/1‑251 (MQ=255)
                                                                                                                                                             gggATTGATTCTGGTGATGATGGCAACCCAGATGTTCCTCGACGGCATTCGAATGTGGATGAAGGGGTAACGAGGAGAATTTTTCGAAAAGAACCCGGTAGCGTc                                                                                                                                                                                       >  1:21512/1‑105 (MQ=255)
                                                                                                                                                                attgattCTGGTGATGATGGCAACCCAGATGTTCCTCGACGGCAGTCGAATGTGGATGAAGTGGTAACGAGGAGAATTTTTCGAAAAGAACCCGGTAGCGTCACAACTACCGGGCGAAc                                                                                                                                                                      <  1:294717/119‑1 (MQ=255)
                                                                                                                                                                   gnntCTGGTNNTGATGGCAACCCAGATGTTCCTCGACGGCATTCGAATGTGGATGAAGGGGTAACGAGGAGAATTTTTCGAAAAGACCCCGGTAGCGTc                                                                                                                                                                                       <  2:21512/99‑1 (MQ=255)
                                                                                                                                                                                       cccAGATGTTCCTCGACGGCATTCGAATGTGGATGAAGGGGTAACGAGGAGAATTTTTCGAAAAGac                                                                                                                                                                                                   >  1:112805/1‑66 (MQ=255)
                                                                                                                                                                                           gATGTTCCTCGACGGCATTCGAATGTGGATGAAGGGTTAACGAGGAGAAGTTTTCGAAAAGAACCCGGTAGCGTCACAACTACCGGGc                                                                                                                                                                          <  1:128355/88‑1 (MQ=255)
                                                                                                                                                                                                         gCATTCGAATGGGGATGAAGGGGTNACTAGGAGAANTTTNNGAAANGNANCCGGTAGCGTCACAACTACCGGGCGAACANACTTAACTCAACAGCTGGAACGCNNTCNNNCCNNCGATg                                                                                                                             <  2:240285/119‑1 (MQ=255)
                                                                                                                                                                                                                                      ggagAATTTTTCGAAAAGAACCCGGTAGCGTCACAACTACCGGGCGAACAAACTTAACTCAACAGCTGGAACGCAATCATCCCAACGAGGGCACCGACAGTGCCGAGGATGGTTTCCATCATGGTCCATGTTTTCAGCGTTTCGGCTTCGGTCGCGCCgg                                                       <  1:249245/160‑1 (MQ=255)
                                                                                                                                                                                                                                           gTTTTTCGAAAAGAACCCGGTAGCGTCACAACTACCGGGCGAACAATCTTAACTCAACAGCTGGAACGCAATCATCCCAGCGATGGCACCGACAGTGCCGATGATTGTTTCCATCATGGTCCAGGTTTTCATCTTTTCGTCGTCGGTCGCGCCGGTAAATTTACCGAACAACCAGAAACCTGCTTCGTTAACGAGGCTGACAACAATc  <  1:273005/208‑1 (MQ=255)
                                                                                                                                                                                                                                           |                                                                                                                                                                                                               
CGCTGATGTTGTTGTCTCATCAGTACCCGAATCAGATGGGGCATCTGGTGATTGCTCTGCTGCTGGCCTGGGGCGGCACCTTTGTCATCCTGCTACAGTCTTCGCTATTTTTACGTCTGCTGGGCGAGAAAGGGGTGAACGCACTTGAACGCCTGATGGGATTGATTCTGGTGATGATGGCAACCCAGATGTTCCTCGACGGCATTCGAATGTGGATGAAGGGGTAACGAGGAGAGTTTTTCGAAAAGAACCCGGTAGCGTCACAACTACCGGGCGAACAAACTTAACTCAACAGCTGGAACGCAATCATCCCAACGATGGCACCGACAGTGCCGAGGATGGTTTCCATCATGGTCCAGGTTTTCAGCGTTTCGGCTTCGGTCGCGCCGGTAAATTTACCGAACAACCAGAAACCGGCGTCGTTAACGTGGCTGACAACAATC  >  NC_000913/3575438‑3575880

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