Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,148,214 T→A 9.1% intergenic (+17/+56) mscMB → / ← PP_5695 miniconductance mechanosensitive channel protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,148,2140TA 9.1% 156.8 / 5.9 55intergenic (+17/+56)mscMB/PP_5695miniconductance mechanosensitive channel protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (23/26);  new base A (3/2);  total (26/29)
Fisher's exact test for biased strand distribution p-value = 6.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01

GCGGTGTTGCCGGAGTTCGGCTTGAGCCTGTATCAGCAGCCGAGTGGCAATGACATGCGCGTGGGGTTGGCCGGGCGTACGGCGGCGGAGCCGATGCCGCGACACCTTGAAGAAATGAGCGAGGCCTGAGATTGCCGGGGCTGCCTTGCAGCCCCTTACGAGCCAAGGCCGCTCTGCAAAGGAGCGGCGGGATGATGGTGTTTAAGGGGTTTTCATCATGCCCAATTCAGCATCATCCATCAGCGCCTTGGCCAAGGCACAGAGGTAGTGCGAAGCCCAGATTAGCT  >  NC_002947/5148069‑5148355
                                                                                                                                                 |                                                                                                                                             
gCGGTGTTGCCGGAGTTCGGCTTGAGCCTGTATCAGCAGCCGAGTGGCAATGACATGCGCGTGGGGTTGGCCGGGCGTACGGCGGCGGAGCCGATGCCGCGACACCTTGAAGAAATGAGCGAGGCCTGAGATTGCCGGGGCTGCCTTg                                                                                                                                             <  2:190248/148‑1 (MQ=255)
 cGGTGTTGCCGGAGTTCGGCTTGAGCCTGTATCAGCAGCCGAGTGGCAATGACATGCGCGTGGGGTTGGCCGGGCGTACGGCGGCGGAGCCGATGCCGCGACACCTTGAAGAAATGAGCGAGGCCTGAGATTGCCGGGGCTGCCTTGCa                                                                                                                                           <  6:107607/149‑1 (MQ=255)
 cGGTGTTGCCGGAGTTCGGCTTGAGCCTGTATCAGCAGCCGAGTGGCAATGACATGCGCGTGGGGTTGGCCGGGCGTACGGCGGCGGAGCCGATGCCGCGACACCTTGAAGAAATGAGCGAGGCCTGAGATTGCCGGGGCTGCCTTGCa                                                                                                                                           <  5:248081/149‑1 (MQ=255)
      ttGCCGGAGTTGGGCTTGAGCCTGTATCAGCAGCCGAGTGGCAATGACATGCGCGTGGGGTTGGCCGGGCGTACGGCGGCGGAGCCGATGCCGCGACACCTTGAAGAAATGAGCGAGGCCTGAGATTGCCGGGGCTGCCTTGCAGccc                                                                                                                                       <  4:43046/148‑1 (MQ=255)
        gCCGGAGTTCGGCTTGAGCCTGTATCAGCAGCCGAGTGGCAATGACATGCGCGTGGGGTTGGCCGGGCGTACGGCGGCGGAGCCGATGCCGCGACACCTTGAAGAAATGAGCGAGGCCTGAGATTGCCGGGGCTGCAAGGCAGccccg                                                                                                                                     <  6:251392/148‑2 (MQ=255)
        gCCGGAGTTCGGCTTGAGCCTGTATCAGCAGCCGAGTGGCAATGACATGCGAGTGGGGTTGGCCGGGCGTACGGCGGCGGAGCCGATGCCGCGACACCTTGAAGAAATGAGCGAGGCCTGAGATTGCCGGGGCTGCCTAGCAGCCCCt                                                                                                                                     >  8:71808/1‑148 (MQ=255)
           ggAGTTCGGCTTGAGCCTGTATCAGCAGCCGAGTGGCAATGACATGCGCGTGGGGTTGGCCGGGCGTACGGCGGCGGAGCCGATGCCGCGACACCTTGAAGAAATGAGCGAGGCCTGAGATTGCAGGGGCTGCAAGGCAGccccggcaa                                                                                                                                 >  3:131778/1‑144 (MQ=17)
                 cGGCTTGAGCCTGTATCAGCAGCCGAGTGGCAATGACATGCGCGTGGGGTTGGCCGGGCGTACGGCGGCGGAGCCGATGCCGCGACACCTTGAAGAAATGAGCGAGGCCTGAGATTGCCGGGGCTGCCTTGCAGCCCCTTa                                                                                                                                   >  4:176934/1‑141 (MQ=255)
                 cGGCTTGAGCCTGTATCAGCAGCCGAGTGGCAATGACATGCGCGTGGGGTTGGCCGGGCGTACGGCGGCGGAGCCGATGCCGCGACACCTTGAAGAAATGAGCGAGGCCTGAGATTGCCGGGGCTGCCTTGCAGCCCCTTa                                                                                                                                   <  3:176934/141‑1 (MQ=255)
                    cTTGAGCCTGTATCAGCAGCCGAGTGGCAATGACATGCGCGTGGGGTTGGCCGGGCGTACGGCGGCGGAGCCGATGCCGCGACACCTTGAAGAAATGAGCGAGGCCTGAGATTGCCGGGGCTGCCTTGCAGCCCCTTACGAGCCAAGGc                                                                                                                        >  5:6984/1‑149 (MQ=255)
                           cTGTATCAGCAGCCGAGTGGCAATGACATGCGCGTGGGGTTGGCCGGGCGTACGGCGGCGGAGCCGATGCCGCGACACCTTGAAGAAATGAGCGAGGCCTGAGATTGCCGGGGCTGCCTTGCAGccc                                                                                                                                       <  1:328213/127‑1 (MQ=255)
                           cTGTATCAGCAGCCGAGTGGCAATGACATGCGCGTGGGGTTGGCCGGGCGTACGGCGGCGGAGCCGATGCCGCGACACCTTGAAGAAATGAGCGAGGCCTGAGATTGCCGGGGCTGCCTTGCAGCCCCTTACGAGCCAAGGCCGCTCTg                                                                                                                 >  4:214832/1‑149 (MQ=255)
                              tATCAGCAGCCGAGTGGCAATGACATGCGCGTGGGGTTGGCCGGGCGTACGGCGGCGGAGCCGATGCCGCGACACCTTGAAGAAATGAGCGAGGCCTGAGATTGCCGGGGCTGCCTTGCAGCCCCTTACGAGCCAAGGCCGCTCTGCaa                                                                                                              <  3:316259/149‑1 (MQ=255)
                                            tGGCAATGACATGCGCGTGGGGTTGGCCGGGCGTACGGCGGCGGAGCCGATGCCGCGACACCTTGAAGAAATGAGCGAGGCCTGAGATTGCCGGGGCTGCCTTGCAGCCCCTTACGAGCCAAGGCCGCTCTGCAAAGGAGCGGCGGgat                                                                                                <  1:9524/149‑1 (MQ=255)
                                                aaTGACATGCGCGTGGGGTTGGCCGGGCGTACGGCGGCGGAGCCGATGCCGCGACACCTTGAAGAAATGAGCGAGGCCTGAGATTGCCGGGGCTGCCTTGCAGCCCCTTACGAGCCAAGGCCGCTCTGCAAAGGAGCGGCGGGatgatg                                                                                            <  7:71808/149‑1 (MQ=255)
                                                aaTGACATGCGCGTGGGGTTGGCCGGGCGTACGGCGGCGGAGCCGATGCCGCGACACCTTGAAGAAATGAGCGAGGCCTGAGATTGCCGGGGCTGCCTTGCAGCCCCTTACGAGCCAAGGCCGCTCTGCAAAGGAGCGGCGGGatgatg                                                                                            <  7:58442/149‑1 (MQ=255)
                                                      aTGCGCGTGGGGTTGGCCGGGCGTACGGCGGCGGAGCCGATGCCGCGACACCTTGAAGAAATGAGCGAGGCCTGAGATTGCCGGGGCTGCCTTGCAGCCCCTTACGAGCCAAGGCCGCTCTGCAAAGGAGCGGCGGGATGATGGTGtt                                                                                       >  3:25875/1‑148 (MQ=255)
                                                        gcgcGTGGGGTTGGCCGGGCGTACGGCGGCGGAGCCGATGCCGCGACACCTTGAAGAAATGAGCGAGGCCTGAGATTGCCGGGGCTGCCTTGCAGCCCCTTACGAGCCAAGGCCGCTCTGCAAAGGAGcggcgg                                                                                                   >  8:115236/1‑134 (MQ=255)
                                                        gcgcGTGGGGTTGGCCGGGCGTACGGCGGCGGAGCCGATGCCGCGACACCTTGAAGAAATGAGCGAGGCCTGAGATTGCCGGGGCTGCCTTGCAGCCCCTTACGAGCCAAGGCCGCTCTGCAAAGGAGcggcgg                                                                                                   <  7:115236/134‑1 (MQ=255)
                                                         cgcgTGGGGTTGGCCGGGCGTACGGCGGCGGAGCCGATGCCGCGACACCTTGAAGAAATGAGCGAGGCCTGAGATTGCCGGGGCTGCCTTGCAGCCCCTTACGAGCCAAGGCCGCTCTGCAAAGGAGCGGCGGGATGATGGTGTTTa                                                                                     >  6:179079/1‑147 (MQ=255)
                                                         cgcgTGGGGTTGGCCGGGCGTACGGCGGCGGAGCCGATGCCGCGACACCTTGAAGAAATGAGCGAGGCCTGAGATTGCCGGGGCTGCCTTGCAGCCCCTTACGAGCCAAGGCCGCTCTGCAAAGGAGCGGCGGGATGATGGTGTTTa                                                                                     <  5:179079/147‑1 (MQ=255)
                                                          gcgTGGGGTTGGCCGGGCGTACGGCGGCGGAGCCGATGCCGCGACACCTTGAAGAAATGAGCGAGGCCTGAGATTGCCGGGGCTGCCTTGCAGCCCCTTACGAGCCAAGGCCGCTCTGCAAAGGAGCGGCGGGATGATGGTGTTTAAgg                                                                                  >  7:20525/1‑149 (MQ=255)
                                                               gggTTGGCCGGGCGTACGGCGGCGGAGCCGATGCCGCGACACCTTGAAGAAATGAGCGAGGCCTGAGATTGCCGGGGCTGCCTTGCAGCCCCTTACGAGCCAAGGCCGCTCTGCAAAGGAGCGGCGGGATGATGGTGTTTAAGGGGttt                                                                             >  8:149727/1‑149 (MQ=255)
                                                                ggTTGGCCGGGCGTACGGCGGCGGAGCCGATGCCGCGACACCTTGAAGAAATGAGCGAGGCCTGAGATTGCCGGGGCTGCCTTGCAGCCCCTTACGAGCCAAGGCCGCTCTGCAAAGGAGCGGCGGGATGATgg                                                                                           >  6:250102/1‑134 (MQ=255)
                                                                 gTTGGCCGGGCGTACGGCGGCGGAGCCGATGCCGCGACACCTTGAAGAAATGAGCGAGGCCTGAGATTGCCGGGGCTGCCTTGCAGCCCCTTACGAGCCAAGGCCGCTCTGCAAAGGAGCGGCGGGATGATGGTGTTTAAGGGGTTTtc                                                                           <  6:6984/149‑1 (MQ=255)
                                                                     gCCGGGCGTACGGCGGCGGAGCCGATGCCGCGACACCTTGAAGAAATGAGCGAGGCCTGAGATTGCCGGGGCTGCCTTGCAGCCCCTTACGAGCCAAGGCCGCTCTGCAAAGGAGCGGCGGGATGATGGTGTTTAAGGGGTTTtcatc                                                                        <  7:149727/148‑1 (MQ=255)
                                                                         ggCGTACGGCGGCGGAGCCGATGCCGCGACACCTTGAAGAAATGAGCGAGGCCTGAGATTGCCGGGGCTGCCTTGCAGCCCCTTACGAGCCAAGGCCGCTCTGCAAAGGAGCGGCGGGATGATGGTGTTTAAGGGGTTTTCATCATGcc                                                                   <  7:157085/149‑1 (MQ=255)
                                                                                ggcggcggAGCCGATGCCGCGACACCTTGAAGAAATGAGCGAGGCCTGAGATTGCCGGGGCTGCCTTGCAGCCCCTTACGAGCCAAGGCCGCTCTGCAAAGGAGCGGCGGGATGATGGTGTTTAAgggg                                                                                >  8:55452/1‑129 (MQ=255)
                                                                                ggcggcggAGCCGATGCCGCGACACCTTGAAGAAATGAGCGAGGCCTGAGATTGCCGGGGCTGCCTTGCAGCCCCTTACGAGCCAAGGCCGCTCTGCAAAGGAGCGGCGGGATGATGGTGTTTAAgggg                                                                                <  7:55452/129‑1 (MQ=255)
                                                                                 gcggcggAGCCGATGCCGCGACACCTTGAAGAAATGAGCGAGGCCTGAGATTGCCGGGGCTGCCTTGCAGCCCCTTACGAGCCAAGGCCGCTCTGCAAAGGAGCGGCGGGATGATGGTGTTTAAGGGGTTTTCATCATGCCCAATTCAg                                                           >  3:203795/1‑149 (MQ=255)
                                                                                    gcggAGCCGATGCCGCGACACCTTGAAGAAATGAGCGAGGCCTGAGATTGCCGGGGCTGCCTTGCAGCCCCTTACGAGCCAAGGCCGCTCTGCAAAGGAGCGGCGGGAAGATGGTGATTAAGGGGTGTTCATCATGACAAATTCAGca                                                         >  8:225957/1‑148 (MQ=255)
                                                                                      ggAGCCGATGCCGCGACACCTTGAAGAAATGAGCGAGGCCTGAGATTGCCGGGGCTGCCTTGCAGCCCCTTACGAGCCAAGGCCGCTCTGCAAAGGAGCGGCGGGATGATGGTGTTTAAGGGGTTTTCATCATGCCCAATTCAGcatca                                                      >  7:220465/1‑149 (MQ=255)
                                                                                                 cgcgACACCTTGAAGAAATGAGCGAGGCCTGAGATTGCCGGGGCTGCCTTGCAGCCCCTTACGAGCCAAGGCCGCTCTGCAAAGGAGCGGCGGGATGATGGTGTTTAAGGGGTTTTCATCATGCCCAATTCAGCATCATCCATCAgcgc                                           >  4:319962/1‑149 (MQ=255)
                                                                                                      cacCTTGAAGAAATGAGCGAGGCCTGAGATTGCCGGGGCTGCCTTGCAGCCCCTTACGAGCCAAGGCCGCTCTGCAAAGGAGCGGCGGGATGATGGTGTTTAAGGGGTTTTCATCATGCCCAATTCAGCATCATCCATCAGCGCCTTgg                                      <  4:25875/149‑1 (MQ=255)
                                                                                                          ttGAAGAAATGAGCGAGGCCTGAGATTGCCGGGGCTGCCTTGCAGCCCCTTACGAGCCAAGGCCGCTCTGCAAAGGAGCGGCGGGATGATGGTGTTTAAGGGGTTTTCATCATGCCCAATTCAGCATCATCCATCAGCGCCTTGGcc                                    >  5:251269/1‑147 (MQ=255)
                                                                                                          ttGAAGAAATGAGCGAGGCCTGAGATTGCCGGGGCTGCCTTGCAGCCCCTTACGAGCCAAGGCCGCTCTGCAAAGGAGCGGCGGGATGATGGTGTTTAAGGGGTTTTCATCATGCCCAATTCAGCATCATCCATCAGCGCCTTGGcc                                    >  3:57887/1‑147 (MQ=255)
                                                                                                             aagaaATGATAGAGGCATGAGATTGCCGGGGCTGCCTTGCAGCCCCTTACGAGACAAGGCCGCTCTGCAAAGGATCGGCGGGATGATGGTGTTTAAGGGGTTTTCATCATGCCCAATTCAGCATCATCCATCAGCGCCtttg                                      <  7:308800/142‑3 (MQ=255)
                                                                                                             aagaaATGAGCGAGGCCTGAGATTGCCGGGGCTGCCTTTCAGCCCCTTACGAGCCAAGGCCGCTCTGCAAAGGAGCGGCGTGATGATGGTGTTTAAGGTGTTTTCATCATGCCCAATACAGCATCATCCATAAGAGCCTTgt                                      >  8:308800/1‑141 (MQ=255)
                                                                                                              agaaATGAGCGAGGCCTGAGATTGCCGGGGCTGCCTTGCAGCCCCTTACGAGCCAAGGCCGCTCTGCAAAGGAGCGGCGGGATGATGGTGTTTAAGGGGTTTTCATCATGCCCAATTCAGCATCATCCATCAGCGCCTTGGCCAAGGc                               >  5:126848/1‑148 (MQ=255)
                                                                                                               gaaATGAGCGAGGCCTGAGATTGCCGGGGCTGCCTTGCAGCCCCTTACGAGCCAAGGCCGCTCTGCAAAGGAGCGGCGGGATGATGGTGTTTAAGGGGTTTTCATCATGCCCAATTCAGCATCATCCATCAGCGCCTTGGCCAAGGca                              >  4:126205/1‑148 (MQ=255)
                                                                                                                 aaTGAGCGAGGCCTGAGATTGCCGGGGCTGCCTTGCAGCCCCTTACGAGCCAAGGCCGCTCTGCAAAGGAGCGGCGGGATGATGGTGTTTAAGGGGTTTTCATCATGCCCAATTCAGCATCATCCATCAGCGCCTTGGCCAAGGCACag                           <  7:225957/149‑1 (MQ=255)
                                                                                                                  aTGCGCGCGGCCTGTGTATGCCGGGGCTGGCGTGCAGCCCCTTACGACCCAACGCCGCTCTGCAAAGGCGCGGCGGGCTGATGGTGTTTAAGGGGTTTTCATCATGCCCAATTCAGCATCATCCATCAGCGCCTTGGCCAAGGCACaga                          <  4:203795/149‑1 (MQ=255)
                                                                                                                        gAGGCCTGAGATTGCCGGGGCTGCCTTGCAGCCCCTTACGAGCCAAGGCCGCTCTGCAAAGGAGCGGCGGGATGATGGTGTTTAAGGGGTTTTCATCATGCCCAATTCAGCATCATCCATCAGCGCCTTGGCCAAGGCACAGAGGTAGt                    <  4:57887/149‑1 (MQ=255)
                                                                                                                        gAGGCCTGAGATTGCCGGGGCTGCCTTGCAGCCCCTTACGAGCCAAGGCCGCTCTGCAAAGGAGCGGCGGGATGATGGTGTTTAAGGGGTTTTCATCATGCCCAATTCAGCATCATCCATCAGCGCCTTGGCCAAGGCACAGAGGTAGt                    <  3:126205/149‑1 (MQ=255)
                                                                                                                              tGAGATTGCCGGGGCTGCCTTGCAGCCCCTTACGAGCCAAGGCCGCTCTGCAAAGGAGCGGCGGGATGATGGTGTTTAAGGGGTTTTCATCATGCCCAATTCAGCATCATCCATCAGCGCCTTGGCCAAGGcac                             >  7:131641/1‑134 (MQ=255)
                                                                                                                              tGAGATTGCCGGGGCTGCCTTGCAGCCCCTTACGAGCCAAGGCCGCTCTGCAAAGGAGCGGCGGGATGATGGTGTTTAAGGGGTTTTCATCATGCCCAATTCAGCATCATCCATCAGCGCCTTGGCCAAGGcac                             <  8:131641/134‑1 (MQ=255)
                                                                                                                              tGAGATTGCCGGGGCTGCCTTGCAGCCCCTTACGAGCCAAGGCCGCTCTGCAAAGGAGCGGCGGGATGATGGTGTTTAAGGGGTTTTCATCATGCCCAATTCAGCATCATCCATCAGCGCCTTGGCCAAGGCACAGAGGTAGTGCGAAg              >  1:80675/1‑149 (MQ=255)
                                                                                                                               gagaTTGCCGGGGCTGCCTTGCAGCCCCTTACGAGCCAAGGCCGCTCTGCAAAGGAGCGGCGGGATGATGGTGTTTAAGGGGTTTtcat                                                                         <  3:167364/89‑1 (MQ=255)
                                                                                                                                  tcgtaaggggCTGCAAGGCAGCCCCTTACGAGCCAAGGCCGCTCTGCAAAGGAGCGGCGGGATGATGGTGTTTAAGGGGTTTTCATCATGCCCAATTCAGCATCATCCATCAGCGCCTTGGCCAAGGCACAGAGATAGTGCGaa               <  4:14786/138‑1 (MQ=255)
                                                                                                                                  tcgtaaggggCTGCAAGGCAGCCCCTTACGAGCCAAGGCCGCTCTGCAAAGGAGCGGCGGGATGATGGTGTTTAAGGGGTTTTCATCATGCCCAATTCAGCATCATCCATCAGCGCCTTGGCCAAGGCACAGAGATAGTGCGaa               >  3:14786/7‑144 (MQ=255)
                                                                                                                                  aTTGCCGGGGCTGCCTTGCAGCCCCTTACGAGCCAAGGCCGCTCTGCAAAGGAGCGGCGGGATGATGGTGTTTAAGGGGTTTTCATCATGCCCAATTCAGCATCATCCATCAGCGCCTTGGCCAAGGCACAGAGGTAGTGCGAAGCCCa          <  6:126848/149‑1 (MQ=255)
                                                                                                                                   ttGCCGGGGCTGCCTTGCAGCCCCTTACGAGCCAAGGCCGCTCTGCAAAGGAGCGGCGGGATGATGGTGTTTAAGGGGTTTTCATCATGCCCAATTCAGCATCATCCATCAGCGCCTTGGCCAAGGCACAGAGGTAGTGCGAAGCCCAg         <  6:251269/149‑1 (MQ=255)
                                                                                                                                       aggggCTGCAAGGCAGCCCCTTACGAGCCAAGGCCGCTCTGCAAAGGAGCGGCGGGATGATGGTGTTTAAGGGGTTTTCATCATGCCCAATTCAGCATCATCCATCAGCGCCTTGGCCAAGGCACAGAGGTAGTGCGAAGCCCAGAt       >  8:91115/2‑147 (MQ=255)
                                                                                                                                             tGCCTTGCAGCCCCTTACGAGCCAAGGCCGCTCTGCAAAGGAGCGGCGGGATGATGGTGTTTAAGGGGTTTTCATCATGCCCAATTCAGCATCATCCATCAGCGCCTTGGCCAAGGCACAGAGGTAGTGCGAAGCCCAGATTAGCt  >  1:315185/1‑146 (MQ=255)
                                                                                                                                             tGCCTTGCAGCCCCTTACGAGCCAAGGCCGCTCTGCAAAGGAGCGGCGGGATGATGGTGTTTAAGGGGTTTTCATCATGCCCAATTCAGCATCATCCATCAGCGCCTTGGCCAAGGCACAGAGGTAGTGCGAAGCCCAGATTAGCt  <  2:315185/146‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                             
GCGGTGTTGCCGGAGTTCGGCTTGAGCCTGTATCAGCAGCCGAGTGGCAATGACATGCGCGTGGGGTTGGCCGGGCGTACGGCGGCGGAGCCGATGCCGCGACACCTTGAAGAAATGAGCGAGGCCTGAGATTGCCGGGGCTGCCTTGCAGCCCCTTACGAGCCAAGGCCGCTCTGCAAAGGAGCGGCGGGATGATGGTGTTTAAGGGGTTTTCATCATGCCCAATTCAGCATCATCCATCAGCGCCTTGGCCAAGGCACAGAGGTAGTGCGAAGCCCAGATTAGCT  >  NC_002947/5148069‑5148355

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