Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,867,576 G→T intergenic (‑154/‑152) ibpA ← / → yidQ small heat shock protein IbpA/DUF1375 domain‑containing putative lipoprotein YidQ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,867,5760GT97.1% 129.6 / ‑2.5 35intergenic (‑154/‑152)ibpA/yidQsmall heat shock protein IbpA/DUF1375 domain‑containing putative lipoprotein YidQ
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base T (21/13);  minor base . (1/0);  total (22/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 = 7.53e-01

GGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTGCAGTTTATGGTCTATTGCATAGACTGAGGGGGCAGCACGAAACGTTAAAGTGCGATAGCCGCCACAGTGCGACGTAATGATGATGTTATTTTCAGCAAGGGATCGCTATAACTCATCAGGCCAACACCCACGATCAACAGATGAATTTATTATGATAAGAAATGTGTTGTTAGCGTTCATGATATGCAGCGGAATGACATTACTCGGAGGATGCTCCAGCGTGATGTCCCATACCGGA  >  NC_000913/3867338‑3867814
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                              
ggctttggTTTTTTTCTAAGTTGTTAAACAAACGGTCTAATCCAAGAGCAGAGCGGTGAAGCGGGGATAAATCAAATTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGTCAGTTCAGGTAGCCAGAACTTCCAGAAGGCGCGTCCTCATCGGCTACTATGTAAAAATGGGTCTGAAAATGAGCTTTGCAAGCCTTATCAGCGTACTTTTTGTTTCAGTTTATGGt                                                                                                                                                                                                                                     <  2:394550/250‑1 (MQ=255)
ggctctggTTGTTTTCTAAGTGGTTAAACAAACGGTAAAATCAAATAGCAGAACGGTAAAGAGGGGATAAATCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCCAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTc                                                                                                                                                                                                                                    >  2:327860/1‑251 (MQ=255)
      ggTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTAt                                                                                                                                                                                                                                 <  1:121698/248‑1 (MQ=255)
      ggTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTGCCACAATGGACAGGTCATGTAGCCAGAACGCCCATCAGGCGCGTCCTCGTCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTg                                                                                                                                                                                                                               <  1:393232/250‑1 (MQ=255)
       gttgttTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTCAGTTTATGGTCTATTGCAt                                                                                                                                                                                                                            >  1:359624/1‑251 (MQ=255)
                 aaGTGGTTAAACAAACGGTCAAATCCAAGAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTCCTGAAACCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAATTGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGTTCTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagat                                                                                                                                                                                                                     <  2:407822‑M1/249‑6 (MQ=255)
                             aaacGGTAAAATCCAATAGCAGAAGGGTACAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTACTGAAATCAGCGAGAATGTCAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGACGCGTCCTCCTCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAACCTTATCAGCGGACTTTTTTTTTCAGTTTATGTTTTATTGCAGAGtagtctttcgggttgtta                                                                                                                                                                                                        >  2:90462‑M1/1‑232 (MQ=255)
                                            aTAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcaga                                                                                                                                                                                         >  2:196790‑M1/1‑217 (MQ=255)
                                                            gCGGGGATAAATCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcgg                                                                                                                                                                           >  1:141897‑M1/1‑201 (MQ=255)
                                                                    aaaTTAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTTGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctcc                                                                                                                                                                <  1:379746‑M1/251‑59 (MQ=255)
                                                                    aaaTCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctcc                                                                                                                                                                <  2:196882‑M1/251‑59 (MQ=255)
                                                                            gTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgccctttt                                                                                                                                                        >  1:50362‑M1/1‑185 (MQ=255)
                                                                                  cATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGgagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctacgcccttttcttat                                                                                                                                                   >  1:205301‑M1/1‑179 (MQ=255)
                                                                                        cAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttct                                                                                                                                            <  2:50362‑M1/251‑79 (MQ=255)
                                                                                            aGCTCCTGAAATCAGCGAGAATGGAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGGCCTCATCGGCTACGATGAAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGGCTATTGCATAGtagatttgcgggatgataattgcggcttgcagaaataaaaagggcggagatgctctccgccctgttcttatagcttcttgcc                                                                                                                                        >  1:236891‑M1/1‑169 (MQ=255)
                                                                                                          gCGAGAATGTAAGACCTTCCACAATGGCCAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGAATTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtga                                                                                                                           >  1:341294‑M1/1‑155 (MQ=255)
                                                                                                            gagaATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacg                                                                                                                        >  1:173836‑M1/1‑153 (MQ=255)
                                                                                                             agaATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgc                                                                                                                       >  1:146805‑M1/1‑152 (MQ=255)
                                                                                                                aTGTAAGACCGTCCACAATGGACAGGTCAGGTAGCCAGAACACCCAGCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGATCTTTTCAAGCCTTATCAGCGGACTGTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgcctt                                                                                                                    <  2:36213‑M1/251‑103 (MQ=255)
                                                                                                                     aGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGAGGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatccg                                                                                                               <  2:189352‑M1/251‑108 (MQ=255)
                                                                                                                                  tGGACAGGTCAGGGAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagacgatctccgcccttttcttatagattcttgcaggatgcggcgtgaaccccttttccggcctacaaaatcc                                                                                                  >  2:447875‑M1/1‑131 (MQ=255)
                                                                                                                                          tcaggtAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatccggcctacaaaatcatgaaaatt                                                                                          >  1:172143‑M1/1‑123 (MQ=255)
                                                                                                                                                                  gcgcGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtaaatttgcgggatgctaattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaaagccttatccgtctttcacaataatcaaaattcaataaagttaaaaattttcct                                                                    >  2:325466‑M1/1‑99 (MQ=255)
                                                                                                                                                                     cgTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatccggcctacaaaatcatgaaaattcaatacattgcaagattttcgtaggcc                                                               >  1:156507‑M1/1‑96 (MQ=255)
                                                                                                                                                                         cTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaatgccttatccggcctacaaaatcatgaaaattcaat                                                                                      <  1:295452‑M1/224‑133 (MQ=255)
                                                                                                                                                                          tCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggttgataattgcggattgcagaaataaaaagagcggagatgatctccgcccttttcttatagcatcttgccggatgcggcgtgaacgccttatcctgccttcaaaatcatggaaattaaataaattgcaaacatttctgaaccctgct                                                           >  2:263294‑M1/1‑91 (MQ=255)
                                                                                                                                                                            aTCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatccggcctacaaaataatgaaaattcaatacattgcaagattttcttaagcctgata                                                          >  1:159575‑M1/1‑89 (MQ=255)
                                                                                                                                                                                     gATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatccggcctacaaaatcatgaaaattcaatacattgcaagattttcgtaggcctgataagcgtagcgca                                               <  2:115820‑M1/251‑172 (MQ=255)
                                                                                                                                                                                        gTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatccggcctacaaaatcatgaaaattcaatacattgcaagattttcgtaggcctgataagcgtagcgcatc                                             >  1:353101‑M1/1‑77 (MQ=255)
                                                                                                                                                                                                       gAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatccggcctacaaaatcatgaaaattcaatacattgcaagattttcgtaggcctgataagcgtagcgcatcaggcacgctcgcatgg                             >  1:76939‑M1/1‑62 (MQ=255)
                                                                                                                                                                                                                        ccTTATCAGCGGACTTTTGTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagaggatctccgcccttttcttatagcttctttccggatgcggtgtgaacgccttatccggcctacaaaatcatgaaaattcaatacatggcaagattttcgtagtcctgataagcgtagcgcatcaggcacgctcgcatggttagcgccattaaata             <  2:383096‑M1/250‑206 (MQ=255)
                                                                                                                                                                                                                           tATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatccggcctacaaaatcatgaaaattcaatacattgcaagattttcgtaggcctgataagcgtagcgcatcaggcacgctcgcatggttagcgccattaaatatcg          >  1:207808‑M1/1‑42 (MQ=255)
                                                                                                                                                                                                                                 cGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggttgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatccggcctacaaaatcatgtaaattcaatacattgcaagattttcgtaggcctgataagcgtagcgcatcatgcacgctcgcatggttagcgccattaaatatcgttattcg   >  2:52646‑M1/1‑36 (MQ=255)
                                                                                                                                                                                                                                 cGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatccggcctacaaaatcatgaaaattcaatacattgcaagattttcttaggcctgataagcgtagcgcatcaggcacgctcgcatggttagcgccattaaatatcgatattcg   <  1:388918‑M1/251‑216 (MQ=255)
                                                                                                                                                                                                                                  ggACTTTTTTTTTCAGTTTATGGTCTATTGAATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggcgatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatcccgcctacaaaatcatgaaaattcaatacattgcaagattttcgtcggcctgataagcgtagcgcatcaggcacgctcgcatggttagcgccattaaatatcgatattcgc  >  1:126472‑M1/1‑35 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                              
GGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTGCAGTTTATGGTCTATTGCATAGACTGAGGGGGCAGCACGAAACGTTAAAGTGCGATAGCCGCCACAGTGCGACGTAATGATGATGTTATTTTCAGCAAGGGATCGCTATAACTCATCAGGCCAACACCCACGATCAACAGATGAATTTATTATGATAAGAAATGTGTTGTTAGCGTTCATGATATGCAGCGGAATGACATTACTCGGAGGATGCTCCAGCGTGATGTCCCATACCGGA  >  NC_000913/3867338‑3867814

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