Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,470,108 A→C 8.7% intergenic (‑255/+37) chiA ← / ← tufA periplasmic endochitinase/translation elongation factor EF‑Tu 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,470,1080AC 8.7% 147.0 / 5.9 46intergenic (‑255/+37)chiA/tufAperiplasmic endochitinase/translation elongation factor EF‑Tu 1
Reads supporting (aligned to +/- strand):  ref base A (26/16);  new base C (2/2);  total (28/18)
Fisher's exact test for biased strand distribution p-value = 6.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.63e-01

AAAAAGCCTGACAAATGAAGCATTTTAAAAACAGAAACATTCATATTTAAAATGTTAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCAGACCGTCGTCCATCGCGATCGGGTGGATCAGGGTAAC  >  NC_000913/3469963‑3470249
                                                                                                                                                 |                                                                                                                                             
aaaaaGCCTGACAAATGAAGCATTTTAAAAACAGAAACATTCATATTTAAAATGTTAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAgc                                                                                                                                            <  1:222645/149‑1 (MQ=255)
        tGACAAATGAAGCATTTTAAAAACAGAAACATTCATATTTAAAATGTTAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCttt                                                                                                                                     >  4:401987/1‑148 (MQ=255)
         gACAAATGAAGCATTTTAAAAACAGAAACATTCATATTTAAAATGTTAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCTTCGGCGCCCttt                                                                                                                                     <  2:336952/147‑1 (MQ=255)
         gACAAATGAAGCATTTTAAAAACAGAAACATTCATATTTAAAATGTTAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCTTCGGCGCCCttt                                                                                                                                     >  1:336952/1‑147 (MQ=255)
         gACAAATGAAGCATTTTAAAAACAGAAACATTCATATTTAAAATGTTAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCtttt                                                                                                                                    >  3:221600/1‑148 (MQ=255)
           cAAATGAAGCATTTTAAAAACAGAAACATTCATATTTAAAATGTTAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTttcaa                                                                                                                                 <  1:369683/149‑1 (MQ=255)
                      ttttAAAAACAGAAACATTCATATTTAAAATGTTAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTaa                                                                                                                      >  1:327547/1‑149 (MQ=255)
                              aCAGAAACATTCATATTTAAAATGTTAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACgtg                                                                                                              >  2:382174/1‑149 (MQ=255)
                                             tttAAAATGTTAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAAc                                                                                               >  2:363496/1‑149 (MQ=255)
                                              ttAAAATGTTAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACt                                                                                              >  1:378884/1‑149 (MQ=255)
                                                aaaaTGTTAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAAGTCAAAACTAATTAACGTGTAATTAGCCCAGAACtt                                                                                             <  4:383301/148‑1 (MQ=255)
                                                 aaaTGTTAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACttt                                                                                            >  2:77713/1‑148 (MQ=255)
                                                 aaaTGTTAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACttt                                                                                            <  1:77713/148‑1 (MQ=255)
                                                 aaaTGTTAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACttt                                                                                            >  4:358785/1‑148 (MQ=255)
                                                  aaTGTTAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTa                                                                                           >  4:269420/1‑148 (MQ=255)
                                                  aaTGTTAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAg                                                                                          <  1:325104/149‑1 (MQ=255)
                                                     gTTAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGcaa                                                                                       <  2:378884/149‑1 (MQ=255)
                                                      ttAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATATGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGcaac                                                                                      >  2:333022/1‑149 (MQ=255)
                                                      ttAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGcaac                                                                                      >  1:14174/1‑149 (MQ=255)
                                                                 tgaTATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCgcgc                                                                           <  2:327547/149‑1 (MQ=255)
                                                                  gaTATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCgcgc                                                                           >  3:75339/1‑148 (MQ=255)
                                                                      ttttAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACg                                                                      <  1:382174/149‑1 (MQ=255)
                                                                       tttAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACg                                                                      >  4:270456/1‑148 (MQ=255)
                                                                                 aataatTTATTCGTTCTGACAGTACGAATAAGATATTCAGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGgccgcc                                                           >  1:376291/1‑149 (MQ=255)
                                                                                      tttAGTCGTTCGGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCAc                                                      >  2:187094/1‑149 (MQ=255)
                                                                                             gTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCTAACGGTACGGCCGCCTTCACGGATTg                                                <  3:270456/148‑1 (MQ=255)
                                                                                                 tGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGaaa                                           <  1:345841/149‑1 (MQ=255)
                                                                                                  gACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGaaa                                           <  2:25770/148‑1 (MQ=255)
                                                                                                  gACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGaaa                                           >  1:25770/1‑148 (MQ=255)
                                                                                                      gTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCTTCGGCGCCCttt                                                                                                                                     >  3:397831/1‑54 (MQ=25)
                                                                                                      gTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCTTCGGCGCCCttt                                                                                                                                     <  4:397831/54‑1 (MQ=25)
                                                                                                      gTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCAg                                      >  1:240731/1‑149 (MQ=255)
                                                                                                       tACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCAGa                                     >  2:16029/1‑149 (MQ=255)
                                                                                                          gAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCAGACcg                                  >  1:325381/1‑149 (MQ=255)
                                                                                                            aTAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCAGACcgtc                                <  2:240731/149‑1 (MQ=255)
                                                                                                            aTAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGCTACGGCCGCCTTCACGGATTGCGAAACGCAGACcgt                                 >  3:270212/1‑148 (MQ=255)
                                                                                                                gATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCAGACCgtcgtc                             >  3:113974/1‑148 (MQ=255)
                                                                                                                gATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCAGACCgtcgtc                             >  4:218217/1‑148 (MQ=255)
                                                                                                                   atGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCAGACCGTCGTCCATc                         <  2:376291/149‑1 (MQ=255)
                                                                                                                          caacaaATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCAGACCGTCGTCCATCGCGATc                   >  4:144900/1‑148 (MQ=255)
                                                                                                                             caaATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCAGACCGTCGTCCATCGCGATCggg                >  3:111256/1‑148 (MQ=255)
                                                                                                                               aaTGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTa                                                                   <  3:181074/95‑1 (MQ=255)
                                                                                                                               aaTGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTa                                                                   >  4:181074/1‑95 (MQ=255)
                                                                                                                                   cAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCAGACCGTCGTCCATCGCGATCGGGTGGATc          <  3:358785/148‑1 (MQ=255)
                                                                                                                                           gcgcCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCAGACCGTCGTCCATCGCGATCGGGTGGATCAGGGTaac  >  4:164245/1‑148 (MQ=255)
                                                                                                                                           gcgcCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCAGACCGTCGTCCATCGCGATCGGGTGGATCAGGGTaac  <  3:144900/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                             
AAAAAGCCTGACAAATGAAGCATTTTAAAAACAGAAACATTCATATTTAAAATGTTAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCAGACCGTCGTCCATCGCGATCGGGTGGATCAGGGTAAC  >  NC_000913/3469963‑3470249

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