Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,312,070 Δ1 bp intergenic (‑28/+31) alsB ← / ← alsR D‑allose ABC transporter periplasmic binding protein/DNA‑binding transcriptional repressor AlsR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,312,0690G.88.9% 93.1 / 8.0 27intergenic (‑27/+32)alsB/alsRD‑allose ABC transporter periplasmic binding protein/DNA‑binding transcriptional repressor AlsR
Reads supporting (aligned to +/- strand):  ref base G (2/1);  new base . (12/12);  total (14/13)
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.20e-01

AGATTGAAAATCGCCTTCTGAAGGAGAGGCAAAAATATCAACGCTGACGCCCAGTGTTTTTGCTTCATCTTCAATGCCTTTTTTCATATCTACCCAAAATGGGTTGGAGAGGGTTTTCAATACGACAGCATATTCGGCGGCAGCAAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAAGTTAGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATAATCGGCCAGTTTCGCTATCGGT  >  NC_000913/4311833‑4312311
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                  
gttatcgacgaccaggttgattttgcaacctgatgcgcaataggggcaaaccgtgacgacttttttcatcggtctcgctccagttaatcaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttagtcaaaactctgacgaaaaacaggctgtcgtcaggtttgtcgtgaGTTAAGTGGACAGTAACTCCGCGATgg                                                                                                                                                                                                                                     <  2:157274‑M2/27‑1 (MQ=255)
   atcgacgaccaggttgattgtgcaacctgatgcgcaataggggcaaaccgtgacgacttttttcatcggtctcgctccagttaatcaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatgtattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCg                                                                                                                                                                                                                                   <  2:6368‑M2/29‑1 (MQ=255)
                  gattttgcaacctgatgcgcaataggggcaaaccgtgacgacttttttcatcggtctcgctccagttaatcaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaacaacaggctgtcgttagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCAGAATTACCTTTTTCtt                                                                                                                                                                                                                   >  1:167722‑M2/207‑251 (MQ=255)
                            cctgatgcgcaataggggcaaaccgtgacgacttttttcatcggtctcgatccagttaatcaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAgcgc                                                                                                                                                                                                         >  1:179637‑M2/197‑251 (MQ=255)
                                                 accgtgacgacttttttcatcggtctcgctccagttaatcaaatcacgcatacgcgctctagactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagtttttacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAATAAATcaac                                                                                                                                                                                    <  1:215575‑M2/76‑1 (MQ=255)
                                                 accgtgacgactgttttcatcggtctcgctccagttaatcaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacatgttatgtggggtcattccctgatattacgggcactatgtattcaaaactctgacgaaaaacaggctgtcgtcagtgttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTGAGTTCAGCTCGCCTGGTGAGAAGAAATcaac                                                                                                                                                                                    <  2:210678‑M2/76‑1 (MQ=255)
                                                      gacgacttttttcatcggtctcgctccagttactcaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCGTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAAtt                                                                                                                                                                                <  2:200430‑M2/80‑1 (MQ=255)
                                                                           cgctccagttaatcaaatcacgcattcgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCATGGTGAGAAGAAAGCAACAATTGCGCCGGTTTTTTCCatattaa                                                                                                                                                          <  2:227219‑M2/102‑1 (MQ=255)
                                                                                      atcaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGt                                                                                                                                               >  2:171979‑M2/139‑251 (MQ=255)
                                                                                               cgcatacgcgctctcgactaccgtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactctttattcaaacctctgacgaacaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTAAGCGCGCCTGGTGAGAAGAAATCAACACTTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTg                                                                                                                                      >  1:126456‑M2/130‑251 (MQ=255)
                                                                                                  atacgcgctctcgactacagtatgcatcttttatgccacattttatgtgggggcattccctgatattccgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGCAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGCGAGAAGAAATCCACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGc                                                                                                                                   >  1:186990‑M2/127‑251 (MQ=255)
                                                                                                  atacgcgctctcgactacagtatgcatcttttatgcaacattttatgtggggtcattccctgatattacgggcactatttattccaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGc                                                                                                                                   >  1:261513‑M2/127‑251 (MQ=255)
                                                                                                       cgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTGGTTCAATGTTGAGCTGGGc                                                                                                                              <  1:79961‑M2/130‑1 (MQ=255)
                                                                                                            ccgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCGTTTTCTTATTTCAGCGCGCCTTGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATAGTAGCTTGTTCAATGTTGAGCTGGGCGACag                                                                                                                         <  2:270651‑M2/135‑1 (MQ=255)
                                                                                                             cgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGAc                                                                                                                           >  2:232933‑M2/116‑248 (MQ=255)
                                                                                                                  acagtatgcatcttttatgccacattttatgtggggtcattccctgatattgcgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTTATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGAc                                                                                                                     <  1:268695‑M2/139‑1 (MQ=255)
                                                                                                                                                        attccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTCGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGAAAGGTTGGTTGTAATATTCTTGcc                                                                             >  1:201050‑M2/73‑251 (MQ=255)
                                                                                                                                                          tccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCg                                                                            >  2:70570‑M2/71‑250 (MQ=255)
                                                                                                                                                                                            ctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCt                                         >  1:189308‑M2/37‑251 (MQ=255)
                                                                                                                                                                                                     aacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCa                                 >  1:107264‑M2/28‑250 (MQ=255)
                                                                                                                                                                                                            tgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGCTGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAataatat                         >  1:19758‑M2/21‑251 (MQ=255)
                                                                                                                                                                                                            ggtcgtcagttttgacttgaGGTAAGGGGACAGTAACTCCGCGATGGCTGAATTTCCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCTGTTTTTTCCATATTAATATTAGCTTGTTCAAGGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAAGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAagaatat                         <  2:205064‑M2/231‑1 (MQ=255)
                                                                                                                                                                                                             gtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAaataatatg                        <  1:185260‑M2/232‑2 (MQ=255)
                                                                                                                                                                                                                           cgtgaGGTTAGGGGACAGTAACTCCGCGATGGCGGAATGACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGGTTTTTCCATATTAATATTAGCTTGTTCGATGTTGAGCTGGGCGACTGAGACAAAAAAAGCGTCCAGCAATGTTGGTTGTGATATTCTTGCCGAGGCAGTACGACCTAATAACGGCGTTTCGGGGGCTGGTGAGCAAATAATATTATCGGCCAGTTTCg          <  1:291477‑M2/246‑1 (MQ=255)
                                                                                                                                                                                                                                  tAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATCgg                   <  1:6762/235‑1 (MQ=255)
                                                                                                                                                                                                                                  tAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGAGCGCCTGGCGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTATTATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATCgg                   >  2:6762/1‑235 (MQ=255)
                                                                                                                                                                                                                                   aaGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCggg  >  1:173024/1‑250 (MQ=255)
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                  
AGATTGAAAATCGCCTTCTGAAGGAGAGGCAAAAATATCAACGCTGACGCCCAGTGTTTTTGCTTCATCTTCAATGCCTTTTTTCATATCTACCCAAAATGGGTTGGAGAGGGTTTTCAATACGACAGCATATTCGGCGGCAGCAAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAAGTTAGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATAATCGGCCAGTTTCGCTATCGGT  >  NC_000913/4311833‑4312311

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