Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 537,582 C→G 12.2% intergenic (+109/‑44) PP_mr07 → / → rpoB P26 ncRNA/DNA‑directed RNA polymerase subunit beta

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947537,5820CG12.2% 114.4 / 6.0 42intergenic (+109/‑44)PP_mr07/rpoBP26 ncRNA/DNA‑directed RNA polymerase subunit beta
Reads supporting (aligned to +/- strand):  ref base C (16/20);  new base G (2/3);  total (18/24)
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.93e-01

TTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAGCTGGGGAACGCTGATGGCTTACTCATACACTGAGAAAAAACGTATCCGCAAGGACTTTAGCAAGTTGCCGGACGTCATGGATGTGCCTTACCTCCTGGCCATCCAGCTGGATTCGTAT  >  NC_002947/537434‑537730
                                                                                                                                                    |                                                                                                                                                    
ttGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCccgcc                                                                                                                                                     >  2:1000439/1‑150 (MQ=255)
  gAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCgg                                                                                                                                                   <  2:62344/150‑1 (MQ=255)
      gAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCgg                                                                                                                                                   <  1:674834/146‑1 (MQ=255)
      gAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCgg                                                                                                                                                   >  2:674834/1‑146 (MQ=255)
         aGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGc                                                                                                                                            <  2:376905/150‑1 (MQ=255)
          ggCTGATGGCTGGTGGCTTATGCCCCCCGACGTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGcc                                                                                                                                           >  1:37378/1‑150 (MQ=255)
           gCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTTCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTTCTCCCGGCGGGAGCAGTCa                                                                                                                                          <  1:1420578/150‑1 (MQ=255)
                   cTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCa                                                                                                                                  <  1:757667/150‑1 (MQ=255)
                    tggtggCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAg                                                                                                                                 <  1:204153/150‑1 (MQ=255)
                    tggtggCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGACTGCTCCCGCCGGGAGAAGCCAAACAAGCAg                                                                                                                                 <  2:875473/150‑1 (MQ=255)
                      gtggCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGaa                                                                                                                                              >  2:753714/1‑135 (MQ=255)
                      gtggCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGaa                                                                                                                                              <  1:753714/135‑1 (MQ=255)
                                 cACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAGCTgg                                                                                                                    <  1:934555/150‑1 (MQ=255)
                                 cACCGGCCTGTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAGCTgg                                                                                                                    <  1:1000439/150‑1 (MQ=255)
                                   ccGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTTCTCCCGGCGGgagca                                                                                                                                              <  1:51744/122‑3 (MQ=255)
                                   ccGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTTCTCCCGGCGGgagca                                                                                                                                              >  2:51744/1‑120 (MQ=255)
                                        cTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAGCTGGGGAACGc                                                                                                             >  1:908363/1‑150 (MQ=255)
                                        cTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCCCAAACCCAGGGGCTTGCACGATTTTCTTGCTGCTCCTGTGGGGTGAGTCCAAACATGCAGGTGACCAAGCTGGGGAACGc                                                                                                             <  2:520738/150‑1 (MQ=255)
                                                gTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAg                                                                                                                        <  1:761839/131‑1 (MQ=255)
                                                gTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAg                                                                                                                        >  2:761839/1‑131 (MQ=255)
                                                           ccGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAGCTGGGGAACGCTGATGGCTTACTCATacac                                                                                          >  2:69639/1‑150 (MQ=255)
                                                                     cGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCCACGGTATCGCCGTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAGCTGGGGAACGCTGATGGCTTACTCATACACTGAGaaaaaa                                                                                <  2:37378/150‑1 (MQ=255)
                                                                                     cTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAGCTGGGGAACGCTGATGGCTTACTCATACACTGAGAAAAAACGTATCCGCAAGGACt                                                                <  1:69639/150‑1 (MQ=255)
                                                                                      tGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAGc                                                                                                                       <  1:908769/94‑1 (MQ=255)
                                                                                      tGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAGc                                                                                                                       >  2:908769/1‑94 (MQ=255)
                                                                                           gccacCCACATGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAGCTGGGGAACGCTGATGGCTTACTCATACACTGAGAAAAAACGTATCCGCAAGGACTTTAGCa                                                          >  2:1161134/3‑150 (MQ=255)
                                                                                              acCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAGCTGGGGAACGCTGATGGCTTACTCATACACTGAGAAAAAACGTATCCGCAAGGACTTTAGCaa                                                         <  1:272410/148‑1 (MQ=255)
                                                                                              acCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAGCTGGGGAACGCTGATGGCTTACTCATACACTGAGAAAAAACGTATCCGCAAGGACTTTAGCaa                                                         >  2:272410/1‑148 (MQ=255)
                                                                                              acCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAGCTGGGGAACGCTGATGGCTTACTCATACACTGAGAAAAAACGTATCCGCAAGGACTTTAGCAAGt                                                       >  1:1181842/1‑150 (MQ=255)
                                                                                                   tgtgGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAGCTGGGGAACGCTGATGGCTTACTCATACACTGAGAAAAAACGTATCCGCAATTATTTTAGCAAGTTGCCg                                                  <  2:1195280/150‑1 (MQ=255)
                                                                                                             gcAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAGCTGGGGAACGCTGATGGCTTACTCATACACTGAGAAAAAACGTATCCGCAAGGACTTTAGCAAGTTGcc                                                   <  2:994435/139‑1 (MQ=255)
                                                                                                             gcAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAGCTGGGGAACGCTGATGGCTTACTCATACACTGAGAAAAAACGTATCCGCAAGGACTTTAGCAAGTTGcc                                                   >  1:994435/1‑139 (MQ=255)
                                                                                                                ctccAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAGCTGGGGAACGCTGATGGCTTACTCATACACTGAGAAAAAACGTATCCGCAAGGACTTTAGCAAGTTGCCGGACGTCATGGAtg                                     <  1:866619/148‑1 (MQ=255)
                                                                                                                 aCCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAGCTGGGGAACGCTGATGGCTTACTCATACACTGAGAAAAAACGTATCCGCAAGGACTTTAGCAAGTTGCCGGACGTCATGGAtgt                                    <  1:894858/150‑1 (MQ=255)
                                                                                                                  ccAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAGCTGGGGAACGCTGATGGCTTACTCATACACTGAGAAAAAACGTATCCGCAAGGACTTTAGCAAGTTGCCGGACGGCATGGAtgtg                                   >  1:130012/1‑150 (MQ=255)
                                                                                                                   cAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAGCTGGGGAACGCTGATGGCTTACTCATACACTGAGAAAAAACGTATCCACAAGGACTTTAGCAAGTTGCCGGACGGCATGTATTTGc                                  >  2:662499/1‑150 (MQ=255)
                                                                                                                        ggTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAGCTGGGGAACGCTGATGGCTTACTCATACACTGAGAAAAAACGTATCCGCAAGGACTTTAGCAAGTTGCCGGACGTCATGGATGTGCCTTAc                             <  1:662499/150‑1 (MQ=255)
                                                                                                                         gTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAGCTGGGGAACGCTGATGGCTTACTCATACACTGAGAAAAAACGTATCCGCAAGGACTTTAGCAAGTTGCCGGACGTCATGGATGTGCCTTAcc                            >  1:70658/1‑150 (MQ=255)
                                                                                                                              catgATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAGCTGGGGAACGCTGATGGCTTACTCATACACTGAGAAAAAACGTATCCGCAAGGACTTTAGCAAGTTGCCGGACGTCATGGATGTGCCTTACCTCCTg                       <  2:130012/147‑1 (MQ=255)
                                                                                                                                          gcttctCCCGGCGGGAGCAGCCAAACAAGCAGGTGACCAAGCTGGGGAACGCTGATGGCTTACTCATACACTGAGAAAAAACGTATCCGCaa                                                                     >  1:911643/1‑92 (MQ=255)
                                                                                                                                          gcttctCCCGGCGGGAGCAGCCAAACAAGCAGGTGACCAAGCTGGGGAACGCTGATGGCTTACTCATACACTGAGAAAAAACGTATCCGCaa                                                                     <  2:911643/92‑1 (MQ=255)
                                                                                                                                            tgctCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAGCTGGGGAACGCTGATGGCTTACTCATACACTGAGAAAAAACGTATCCGCAAGGACTTTAGCAAGTTGCCGGACGTCATGGATGTGCCTTACCTCCTGGCCATCCAGCTGGa         >  1:594318/1‑150 (MQ=255)
                                                                                                                                            tgctCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAGCTGGGGAACGCTGATGGCTTACTCATACACTGAGAAAAAACGTATCCGCAAGGACTTTAGCAAGTTGCCGGACGTCATGGATGTGCCTTACCTCCTGGCCATCCAGCTGGa         >  1:1381285/1‑150 (MQ=255)
                                                                                                                                                 ccgccgGGAGAAGCCAAACAAGCAGGTGACCAAGCTGGGGAACGCTGATGGCTTACTCATACACTGAGAAAAAACGTATCCGCAAGGACTTTAGCAAGTTGCCGGACGTCATgg                                        >  2:1392198/1‑114 (MQ=255)
                                                                                                                                                 ccgccgGGAGAAGCCAAACAAGCAGGTGACCAAGCTGGGGAACGCTGATGGCTTACTCATACACTGAGAAAAAACGTATCCGCAAGGACTTTAGCAAGTTGCCGGACGTCATgg                                        <  1:1392198/114‑1 (MQ=255)
                                                                                                                                                   gccgGGAGAAGCCAAACAAGCAGGTGACCAAGCTGGGGAACGCTGATGGCTTACTCATACACTGAGAAAAAACGTATCCGCAAGGACTTTAGCAAGTTGCCGGACGTCATGGATGTGCCTTACCTCCTGGCCATCCAGCTGGATTCGTAt  >  2:301381/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAGCTGGGGAACGCTGATGGCTTACTCATACACTGAGAAAAAACGTATCCGCAAGGACTTTAGCAAGTTGCCGGACGTCATGGATGTGCCTTACCTCCTGGCCATCCAGCTGGATTCGTAT  >  NC_002947/537434‑537730

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