Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,513,224 C→T 11.7% intergenic (+157/‑138) PP_mr18 → / → mraZ RNaseP_bact_a ncRNA/cell division/cell wall cluster transcriptional repressor MraZ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,513,2240CT11.7% 79.8 / 3.2 34intergenic (+157/‑138)PP_mr18/mraZRNaseP_bact_a ncRNA/cell division/cell wall cluster transcriptional repressor MraZ
Reads supporting (aligned to +/- strand):  ref base C (13/17);  new base T (2/2);  total (15/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.15e-01

TTTCCCTGCTTGAATCGATAGCAGATCCTCTCTGGTAATACCGAAAAAATCTTACTCTTAATAAATTACTTTAACTTCGAACTCCATTCGCTTGCGTGGGTTTGATTTTTCTCGTGTTTTTTTCTCTCTGGCAGTCCTTCCGCATTCCCTTGTTGCGCTAAATCTCGGTCCTGTAAGGCTTTTCCTTATGAACGTGCCTTGACGGTGCTGCGGGCGGATTCCTATAGTGTGCGCAAGTGGCAGAAAGTGGGATGAAGTGGGTTTTCTGACA  >  NC_002947/1513077‑1513347
                                                                                                                                                   |                                                                                                                           
tttcCCTGCTTGAATCGATAGCAGATCCTCTCTGGTAATACCGAAAAAATCTTACTCTTAATAAATTACTTTAACTTCGAACTCCATTCGCTTGCGTGGGTTTGATTTTTCTCGTGTTTTTTTCTCTCTGGCAGTCCTTCCGCATTccc                                                                                                                            >  3:14199/1‑149 (MQ=255)
      tGCTTGAATCGATAGCAGATCCTCTCTGGTAATACCGAAAAAATCTTACTCTTAATAAATTACTTTAACTTCGAACTCCATTCGCTTGCGTGGGTTTGATTTTTCTCGTGTTTTTTTCTCTCTGGCAGTCCTTCCGCATTCTCttgttg                                                                                                                      >  6:30081/1‑149 (MQ=255)
        cTTGAATCGATAGCAGATCCTCTCTGGTAATACCGAAAAAATCTTACTCTTAATAAATTACTTTAACTTCGAACTCCATTCGCTTGCGTGGGTTTGATTTTTCTCGTGTTTTTTTCTCTCTGGCAGTCCTTCCGCATTCCCTTGTTgcg                                                                                                                    >  6:204549/1‑149 (MQ=255)
                    gCAGATCCTCTCTGGTAATACCGAAAAAATCTTACTCTTAATAAATTACTTTAACTTCGAACTCCATTCGCTTGCGTGGGTTTGATTTTTCTCGTGTTTTTTTCTCTCTGGCAGTCCTTCCGCATTCTCTTGTTGCGCTAAATCTCgg                                                                                                         <  5:30081/148‑1 (MQ=255)
                    gCAGATCCTCTCTGGTAATACCGAAAAAATCTTACTCTTAATAAATTACTTTAACTTCGAACTCCATTCGCTTGCGTGGGTTTGATTTTTCTCGTGTTTTTTTCTCTCTGGCAGTCCTTCCGCATTCCCTTGTTGCGCTAAATCTCgg                                                                                                         <  5:204549/148‑1 (MQ=255)
                        aTCCTCTCTGGTAATACCGAAAAAATCTTACTCTTAATAAATTACTTTAACTTCGAACTCCATTCGCTTGCGTGGGTTTGATTTTTCTCGTGTTTTTTTCTCTCTGGCAGTCCTTCCGCATTCCCttgtt                                                                                                                       <  2:169978/130‑1 (MQ=255)
                        aTCCTCTCTGGTAATACCGAAAAAATCTTACTCTTAATAAATTACTTTAACTTCGAACTCCATTCGCTTGCGTGGGTTTGATTTTTCTCGTGTTTTTTTCTCTCTGGCAGTCCTTCCGCATTCCCttgtt                                                                                                                       >  1:169978/1‑130 (MQ=255)
                         tCCTCTCTGGTAATACCGAAAAAATCTTACTCTTAATAAATTACTTTAACTTCGAACTCCATTCGCTTGCGTGGGTTTGATTTTTCTCGTGTTTTTTTCTCTCTGGCAGTCCTTCCGCATTCCCTTGTTGCGCTAAATCTCGGTCCTGt                                                                                                   <  4:220235/149‑1 (MQ=255)
                            tctctGGTAATACCGAAAAAATCTTACTCTTAATAAATTACTTTAACTTCGAACTCCATTCGCTTGCGTGGGTTTGATTTTTCTCGTGTTTTTTTCTCTCTGGCAGTCCTTCCGCATTCCCTTGTTGCGCTAAATCTCGGTCCTGTAAg                                                                                                >  1:101781/1‑149 (MQ=255)
                                  gTAATACCGAAAAAATCTTACTCTTAATAAATTACTTTAACTTCGAACTCCATTCGCTTGCGTGGGTTTGATTTTTCTCGTGTTTTTTTCTCTCTGGCAGTCCTTCCGCATTCCCTTGTTGCGCTAAATCTCGGTCCTGTaa                                                                                                 <  3:12006/142‑1 (MQ=255)
                                  gTAATACCGAAAAAATCTTACTCTTAATAAATTACTTTAACTTCGAACTCCATTCGCTTGCGTGGGTTTGATTTTTCTCGTGTTTTTTTCTCTCTGGCAGTCCTTCCGCATTCCCTTGTTGCGCTAAATCTCGGTCCTGTaa                                                                                                 >  4:12006/1‑142 (MQ=255)
                                      tACCGAAAAAATCTTACTCTTAATAAATTACTTTAACTTCGAACTCCATTCGCTTGCGTGGGTTTGATTTTTCTCGTGTTTTTTTCTCTCTGGCAGTCCTTCCGCATTCCCTTGTTGCGCTAAATCTCGGTCCTGTAAGGCTTTTCCtt                                                                                      <  4:14199/149‑1 (MQ=255)
                                       aCCGAAAAAATCTTACTCTTAATAAATTACTTTAACTTCGAACTCCATTCGCTTGCGTGGGTTTGATTTTTCTCGTGTTTTTTTCTCTCTGGCAGTCCTTCCGCATTCCCTTGTTGCGCTAAATCTCGGTCCTGTAAGGCTTTTCCTTa                                                                                     <  6:236951/149‑1 (MQ=255)
                                        ccGAAAAAATCTTACTCTTAATAAATTACTTTAACTTCGAACTCCATTCGCTTGCGTGGGTTTGATTTTTCTCGTGTTTTTTTCTCTCTGGCAGTCCTTCCGCATTCCCTTGTTGCGCTAAATCTCGGTCCTGTAAGGCTTTTCCTTAt                                                                                    <  2:101781/149‑1 (MQ=255)
                                         cGAAAAAATCTTACTCTTAATAAATTACTTTAACTTCGAACTCCATTCGCTTGCGTGGGTTTGATTTTTCTCGTGTTTTTTTCTCTCTGGCAGTCCTTCAGCATTCCCTTGTTGCGCTAAATCTCGGTCCTGTAAGGCTTTTCCTTAt                                                                                    >  8:216406/1‑148 (MQ=255)
                                          gAAAAAATCTTACTCTTAATAAATTACTTTAACTTCGAACTCCATTCGCTTGCGTGGGTTTGATTTTTCTCGTGTTTTTTTCTCTCTGGCAGTCCTTCCGCATTCCCTTGTTGCGCTAAATCTCGGTCCTGTAAGGCTTTTCCTTATg                                                                                   >  3:211402/1‑148 (MQ=255)
                                                         ttAATAAATTACTTTAACTTCGAACTCCATTCGCTTGCGTGGGTTTGATTTTTCTCGTGTTTTTTTCTCTCTGGCAGTCCTTCCGCATTCCCTTGTTGCGCTAAATCTCGGTCCTGTAAGGCTTTTCCTTATGAACGTGCCTTGACg                                                                     <  4:211402/147‑1 (MQ=255)
                                                         ttAATAAATTACTTTAACTTCGAACTCCATTCGCTTGCGTGGGTTTGATTTTTCTCGTGTTTTTTTCTCTCTGGCAGTCCTTCCGCATTCCCTTGTTGCGCTAAATCTCGGTCCTGTAAGGCTTTTCCTTATGAACGTGCCTTGACGGt                                                                   >  1:180597/1‑149 (MQ=255)
                                                           aataaATTACTTTAACTTCGAACTCCATTCGCTTGCGTGGGTTTGATTTTTCTCGTGTTTTTTTCTCTCTGGCAGTCCTTCCGCATTCCCTTGTTGCGCTAAATCTCGGTCCTGTAAGGCTTTTCCTTATGAACGTGCCTTGACgg                                                                    <  7:216406/146‑1 (MQ=255)
                                                            ataaATTACTTTAACTTCGAACTCCATTCGCTTGCGTGGGTTTGATTTTTCTCGTGTTTTTTTCTCTCTGGCAGTCCTTCCGCATTCCCTTGTTGCGCTAAATCTCGGTCCTGTAAGGCTTTTCCTTATGAACGTGCCTTGACGGtgct                                                                <  2:180597/149‑1 (MQ=255)
                                                                aTTACTTTAACTTCGAACTCCATTCGCTTGCGTGGGTTTGATTTTTCTCGTGTTTTTTTCTCTCTGGCAGTCCTTCCGCATTCCCTTGTTGCGCTAAATCTCGGTCCTGTAAGGCTTTTCCTTATGAACGTGCCTTGACGGTGCTgcgg                                                            >  3:167986/1‑149 (MQ=255)
                                                                         aCTTCGAACTCCATTCGCTTGCGTGGGTTTGATTTTTCTCGTGTTTTTTTCTCTCTGGCAGTCCTTCCGCATTCCCTTGTTGCGCTAAATCTCGGTCCTGTAAGGCTTTTCCTTATGAACGTGCCTTGACGGTGCTGCGGGCGGATTc                                                    <  4:167986/148‑1 (MQ=255)
                                                                          catCGACCTCCATTCGCTTGCGTGGGGTTGATTTTTCTCGTGTTTTTTTCTCTCTGGCAGTCCTTCCGCATTCCCTTGTTGCGCTAAATCTCGGTCCTGTAAGGCTTTTCCTTATGAACGTGCCTTGACGGTGCTGCGGGCGGATTCCt                                                  <  2:203347/147‑1 (MQ=255)
                                                                          cTTCGAACTCCATTCGCTTGCGTGGGTTTGATTTTTCTCGTGTTTTTTTCTCTCTGGCAGTCCTTCCGCATTCCCTTGTTGCGCTAAATCTCGGTCCTGTAAGGCTTTTCCTTATGAACGTGCCTTGACGGTGCTGCGGGCGGATTcc                                                   >  1:246189/1‑148 (MQ=255)
                                                                           ttCGAACTCCATTCGCTTGCGTGGGTTTGATTTTTCTCGTGTTTTTTTCTCTCTGGCAGTCCTTCCGCATTCTCTTGTTGCGCTAAATCTCGGTCCTGTAAGGCTTTTCCTTATGAACGTGCCTTGACGGTGCTGCGGGCGGATTCCta                                                 >  4:20779/1‑149 (MQ=255)
                                                                                  tCCATTCGCTTGCGTGGGTTTGATTTTTCTCGTGTTTTTTTCTCTCTGGCAGTCCTTCCGCATTCCCTTGTTGCGCTAAATCTCGGTCCTGTAAGGCTTTTCCTTATGAACGTGCCTTGACGGTGCTGCGGGCGGATTCCTATAGtgtg                                          >  8:4332/1‑149 (MQ=255)
                                                                                         gCTTGCGTGGGTTTGATTTTTCTCGTGTTTTTTTCTCTCTGGCAGTCCTTCCGCATTCCCTTGTTGCGCTAAATCTCGGTCCTGTAAGGCTTTTCCTTATGAACGTGCCTTGACGGTGCTGCGGGCGGATTCCTATAGTGTGCGCAAg                                    <  7:4332/148‑1 (MQ=255)
                                                                                                           tttCTCGTGTTTTTTTCTCTCTGGCAGTCCTTCCGCATTCCCTTGTTGCGCTAAATCTCGGTCCTGTAAGGCTTTTCCTTATGAACGTGCCTTGACGGTGCTGCGGGCGGATTCCTATAGTGTGCGCAAGTGGCAGAAAGTGGGATg                   >  3:129811/1‑147 (MQ=255)
                                                                                                             tctcGTGTTTTTTTCTCTCTGGCAGTCCTTCCGCATTCTCTTGTTGCGCTAAATCTCGGTCCTGTAAGGCTTTTCCTTATGAACGTGCCTTGACGGTGCTGCGGGCGGATTCCTATAGTGTGCGCAAGTGGCAGAAAGTGGGATGAAg                <  3:20779/148‑1 (MQ=255)
                                                                                                                 gtgtTTTTTTCTCTCTGGCAGTCCTTCCGCATTCCCTTGTTGCGCTAAATCTCGGTCCTGTAAGGCTTTTCCTTATGAACGTGCCTTGACGGTGCTGCGGGCGGATTCCTATAGTGTGCGCAAGTGGCAGAAAGTGGGATGAAGTggg            <  2:246189/148‑1 (MQ=255)
                                                                                                                 gtgtTTTTTTCTCTCTGGCAGTCCTTCCGCATTCCCTTGTTGCGCTAAATCTCGGTCCTGTAAGGCTTTTCCTTATGAACGTGCCTTGACGGTGCTGCGGGCGGATTCCTATAGTGTGCGCAAGTGGCAGAAAGTGGGATGAAGTg              <  4:129811/146‑1 (MQ=255)
                                                                                                                            tctctGGCAGTCCTTCCGCATTCCCTTGTTGCGCTAAATCTCGGTCCTGTAAGGCTTTTCCTTATGAACGTGCCTTGACGGTGCTGCGGGCGGATTCCTATAGTGTGCGCAAGTGGCAGAAAGTGGGATGAAGTGGGTTTTCTGaca  >  2:102231/1‑147 (MQ=255)
                                                                                                                            tctctGGCAGTCCTTCCGCATTCCCTTGTTGCGCTAAATCTCGGTCCTGTAAGGCTTTTCCTTATGAACGTGCCTTGACGGTGCTGCGGGCGGATTCCTATAGTGTGCGCAAGTGGCAGAAAGTGGGATGAAGTGGGTTTTCTGaca  <  1:102231/147‑1 (MQ=255)
                                                                                                                                         ttCCGCATTCCCTTGTTGCGCTAAATCTCGGTCCTGTAAGGCTTTTCCTTATGAACGTGCCTTGACGGTGCTGCGGGCGGATTCCTATAGTGTGCGCAAGTGGCAGAAAGTgg                       >  5:175456/1‑113 (MQ=255)
                                                                                                                                         ttCCGCATTCCCTTGTTGCGCTAAATCTCGGTCCTGTAAGGCTTTTCCTTATGAACGTGCCTTGACGGTGCTGCGGGCGGATTCCTATAGTGTGCGCAAGTGGCAGAAAGTgg                       <  6:175456/113‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                           
TTTCCCTGCTTGAATCGATAGCAGATCCTCTCTGGTAATACCGAAAAAATCTTACTCTTAATAAATTACTTTAACTTCGAACTCCATTCGCTTGCGTGGGTTTGATTTTTCTCGTGTTTTTTTCTCTCTGGCAGTCCTTCCGCATTCCCTTGTTGCGCTAAATCTCGGTCCTGTAAGGCTTTTCCTTATGAACGTGCCTTGACGGTGCTGCGGGCGGATTCCTATAGTGTGCGCAAGTGGCAGAAAGTGGGATGAAGTGGGTTTTCTGACA  >  NC_002947/1513077‑1513347

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: