Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,182,066 C→A 12.9% intergenic (‑124/+30) 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,0660CA12.9% 105.2 / 7.7 45intergenic (‑124/+30)PP_4560/csbDribonuclease BN/stress response protein
Reads supporting (aligned to +/- strand):  ref base C (22/17);  new base A (3/3);  total (25/20)
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.49e-01

AGGTCGGGGAAAATCATAGCGTCTCCTTTCGCCGCTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTG‑AGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCGGTCACCTTGCCAACAC  >  NC_002947/5181926‑5182211
                                                                                                                                            |                                                                                                                                                  
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                              >  5:220790/1‑149 (MQ=255)
                                                                                                             ttACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTG‑AGCTCCTGCGCCTTGCCTTCGGCCCGCAGc                              >  1:43216/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                         <  6:171439/134‑1 (MQ=255)
                                                                                                                                 ctTCTCGGGTGAACCCGCGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTG‑AGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTc                         >  5:171439/1‑134 (MQ=255)
                                                                                                                                         gTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTG‑AGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCGGTCACCTTGCCAacac  >  5:173909/1‑149 (MQ=255)
                                                                                                                                            |                                                                                                                                                  
AGGTCGGGGAAAATCATAGCGTCTCCTTTCGCCGCTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTG‑AGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCGGTCACCTTGCCAACAC  >  NC_002947/5181926‑5182211

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: ‑