Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,182,059 G→T 11.5% intergenic (‑117/+37) PP_4560 ← / ← csbD ribonuclease BN/stress response protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,182,0590GT11.5% 104.4 / 7.3 44intergenic (‑117/+37)PP_4560/csbDribonuclease BN/stress response protein
Reads supporting (aligned to +/- strand):  ref base G (21/17);  new base T (2/3);  total (23/21)
Fisher's exact test for biased strand distribution p-value = 6.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.38e-01

AGCCGCGCAGGTCGGGGAAAATCATAGCGTCTCCTTTCGCCGCTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTG‑AGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCGGTCACCT  >  NC_002947/5181918‑5182202
                                                                                                                                             |                                                                                                                                                
aGCCGCGCAGGTCGGGGAAAATCATAGCGTCTCCTTTCGCCGCTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGtga                                                                                                                                           <  6:380112/149‑3 (MQ=255)
        aGGTCGGGGAAAATCATAGCGTCTCCTTTCGCCGCTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGACGGCCTCTTCGCGGGTTCACCCgag                                                                                                                                    >  1:48889/1‑148 (MQ=255)
             gtggAAAATCATAGCGTCTCCTTTCGCCGCTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCTCGGGTGAACCCGCGAAGAgg                                                                                                                              <  1:78888/147‑1 (MQ=255)
             ggggAAAATCATAGCGTCTCCTTTCGCCGCTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCTCGGGTGAACCCGCGAagag                                                                                                                               <  6:346294/148‑1 (MQ=255)
              gggAAAATCATAGCGTCTCCTTTCGCCGCTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGc                                                                                                                             <  3:313532/149‑1 (MQ=255)
              gggAAAATCATAGCGTCTCCTTTCGCCGCTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGc                                                                                                                             >  2:258557/1‑149 (MQ=255)
                      cATAGCGTCTCCTTTCGCCGCTTGATCGTTACACAGGCAAGTCTAATAGATCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCg                                                                                                                     >  1:28713/1‑149 (MQ=255)
                          gCGTCTCCTTTCGCCGCTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTGCACCCGAGAAGAGGCCGGCACCGCCtc                                                                                                                 <  6:373300/149‑1 (MQ=255)
                                ccTTTCGCCGCTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCAGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCtcaagtc                                                                                                            >  7:65554/1‑148 (MQ=255)
                                     cgccgcTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAATAATTTGCTAACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGATGAGGCCGGCACCGCCTCAAGTCAAGGc                                                                                                       <  1:15165/148‑1 (MQ=255)
                                     cgccgcTTGATCGTTACACAGCCAAGTCTAATAGCCCGTGCCCGAGCTGGGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGc                                                                                                       <  6:306845/148‑1 (MQ=255)
                                      gccgcTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGctt                                                                                                     >  7:237495/1‑149 (MQ=255)
                                         gcTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGcttctt                                                                                                  >  1:31580/1‑149 (MQ=255)
                                               tCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCg                                                                                             >  2:18913/1‑148 (MQ=255)
                                                          gCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGagaaga                                                                                                                                >  7:148582/1‑102 (MQ=255)
                                                          gCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGagaaga                                                                                                                                <  8:148582/102‑1 (MQ=255)
                                                              gTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTcc                                                                             >  1:236798/1‑149 (MQ=255)
                                                               tCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTccc                                                                            >  8:272822/1‑149 (MQ=255)
                                                               tCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTccc                                                                            >  1:194249/1‑149 (MQ=255)
                                                                  aaTAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACTAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTcc                                                                             <  8:291041/145‑1 (MQ=255)
                                                                  aaTAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTcc                                                                             >  7:291041/1‑145 (MQ=255)
                                                                       aCCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACctg                                                                    >  2:234667/1‑149 (MQ=255)
                                                                         cGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACctgct                                                                  >  3:196419/1‑149 (MQ=255)
                                                                                        gtgctCCCGGCCGTGGTTTAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTg                                                    <  2:48889/147‑1 (MQ=255)
                                                                                        ctgcgCCCGGCCGTGGTTGAAGAATTTGCTTCCGAATAAGTGCCGGCCTCTTCGCGGGTGCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTG‑a                                                  <  2:307377/149‑1 (MQ=255)
                                                                                              ccGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTG‑AGCTc                                              <  2:28713/147‑1 (MQ=255)
                                                                                                   cGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTG‑AGCTCCTGCgcc                                       >  3:276225/1‑149 (MQ=255)
                                                                                                     tGGTTTAAGAATTTGCTAACTAATACGTGCCGGCCTCTTCCCGGGTTCACCCTCGAAGTGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTG‑AGCTCCTGCgcctt                                     <  4:196419/149‑1 (MQ=255)
                                                                                                       gTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTG‑AGCTCCTGCgccttgc                                   <  1:258557/149‑1 (MQ=255)
                                                                                                        ttGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTG‑AGCTCCTGCgccttgcc                                  >  5:30424/1‑149 (MQ=255)
                                                                                                         tGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTG‑AGCTCCTGCgccttgcct                                 >  2:56090/1‑149 (MQ=255)
                                                                                                          gaagaaTTTGCTTAAGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACAGCATCCTTCACATTTCCCATGACCTGCTGCGCCACGCATTTG‑ACCTCCTCCgccttgcctt                                >  2:91936/1‑149 (MQ=255)
                                                                                                           aagaaTTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTgcgc                                                              >  5:19886/1‑119 (MQ=255)
                                                                                                            agaaTTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTG‑AGCTCCTGCGCCTTGCCTTCg                              <  2:31580/149‑1 (MQ=255)
                                                                                                                 ttGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTG‑AGCTCCTGCGCCTTGCCTTCGGccc                          <  4:276225/148‑1 (MQ=255)
                                                                                                                     ttACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTG‑AGCTCCTGCGCCTTGCCTTCGGCCCGCAGc                     >  1:43216/1‑149 (MQ=255)
                                                                                                                     ttACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTG‑AGCTCCTGCGCCTTGCCTTCGGCCCGCAGc                     >  5:220790/1‑149 (MQ=255)
                                                                                                                       aCGACTAAGTCCCGGCCGCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTG‑AGCTCCTGCGCCTTGCCTTCGGCCCGCAGCtt                   <  4:119400/149‑1 (MQ=255)
                                                                                                                         gAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTG‑AGCTCCTGCGCCTTGCCTTCGGCCCGCAGCtttt                 <  1:234667/149‑1 (MQ=255)
                                                                                                                            tAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGAAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCg               <  1:56090/149‑1 (MQ=255)
                                                                                                                             aaGTGCCGTCCTCTTCGCGGGGTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTG‑AGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGtt             <  6:264487/149‑1 (MQ=255)
                                                                                                                                    ggCCTCTTCTCGGGTGAACCCGCGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTG‑AGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCGGt       >  7:267719/1‑148 (MQ=255)
                                                                                                                                        tcttctCGGGTGAACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTG‑AGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCGGTCACCt  >  6:73346/1‑149 (MQ=255)
                                                                                                                                         ctTCTCGGGTGAACCCGCGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTG‑AGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTc                >  5:171439/1‑134 (MQ=255)
                                                                                                                                         ctTCTCGGGTGAACCCGCGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTG‑AGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTc                <  6:171439/134‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                
AGCCGCGCAGGTCGGGGAAAATCATAGCGTCTCCTTTCGCCGCTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTG‑AGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCGGTCACCT  >  NC_002947/5181918‑5182202

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