Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,079,641 G→C 7.9% intergenic (‑376/+58) PP_2688 ← / ← PP_2689 hypothetical protein/endoribonuclease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,079,6410GC 7.9% 147.3 / 2.6 49intergenic (‑376/+58)PP_2688/PP_2689hypothetical protein/endoribonuclease
Reads supporting (aligned to +/- strand):  ref base G (21/24);  new base C (2/2);  total (23/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.04e-01

AAGGGCCGGCGCATCAAAAAAGGGCGCTTTTGACCCGGGCGCGGGTTTTGAACCGAAATGGGGCGCAGGCCGGTGCCTGGCGGGAGAAAAAGGGAGTATTCCGGGTGCATGACGTGGGGCGGAGGGGGCATGCATCGAGCGCCGCGCGGGCGGCACTCGATTGCGCGAACGCTGAATATC‑GTGAGACGGGCGCCTTCAGGGCGTCAGAGGCAGGCGTAGACCACCAGTTCCACCAGCATTTCTTCCCGGGCGAGGCCAGCCACCACTACGGCCGCGCGGTTGGGGTAGGGTGC  >  NC_002947/3079496‑3079788
                                                                                                                                                 |                                                                                                                                                    
aaGGGCCGGCGCATCAAAAAAGGGCGCTTTTGACCCGGGCGCGGGTTTTGAACCGAAATGGGGCGCAGGCCGGTGCCTGGCGGGAGAAAAAGGGAGTATTCCGGGTGCATGCCGTGGGGCGGAGGGGGCATGCATCGAGCGCCGCgcgg                                                                                                                                                   >  6:80114/1‑149 (MQ=255)
 aGGGCCGGCGCATCAAAAAAGGGCGCTTTTGACCCGGGCGCGGGTTTTGAACCGAAATGGGGCGCAGGCCGGTGCCTGGCGGGAGAAAAAGGGAGTATTCCGGGTGCATGACGTGGGGCGGAGGGGGCATGCATCGAGCGCCGCgcggg                                                                                                                                                  <  5:314329/149‑1 (MQ=255)
       ggCGCATCAAAAAAGGGCGCTTTTGACCCGGGCGCGGGTTTTGAACCGAAATGGGGCGCAGGCCGGTGCCTGGCGGGAGAAAAAGGGAGTATTCCGGGTGCATGACGTGGGGCGGAGGGGGCATGCATCGAGCGCCGCGCGGGCGGCAc                                                                                                                                            >  2:278888/1‑149 (MQ=255)
           cATCAAAAAAGGGCGCTTTTGACCCGGGCGCGGGTTTTGAACCGAAATGGGGCGCAGGCCGGTGCCTGGCGGGAGAAAAAGGGAGTATTCCGGGTGCATGACGTGGGGCGGAGGGGGCATGCATCGAGTGCCGCCCGCGCGGCGCTCGa                                                                                                                                        <  6:196041/149‑1 (MQ=255)
            aTCAAAAAAGGGCGCTTTTGACCCGGGCGCGGGTTTTGAACCGAAATGGGGCGCAGGCCGGTGCCTGGCGGGAGAAAAAGGGAGTATTCCGGGTGCATGACGTGGGGCGGAGGGGGCATGCATCGAGCGCCGCGCGGGCGGCACTCGAt                                                                                                                                       >  4:98754/1‑149 (MQ=255)
             tCAAAAAAGGGCGCTTTTGACCCGGGCGCGGGTTTTGAACCGAAATGGGGCGCAGGCCGGTGCCTGGCGGGAGAAAAAGGGAGTATTCCGGGTGCATGACGTGGGGCGGAGGGGGCATGCATCGAGCGCCGCGCGGGCGGCACTCGAtt                                                                                                                                      >  6:143324/1‑149 (MQ=255)
             tCAAAAAAGGGCGCTTTTGACCCGGGCGCGGGTTTTGAACCGAAATGGGGCGCAGGCCGGTGCCGGGCGGGAGAAAAAGGGAGTATTCCGGGTGCATGACGTGGGGCGGAGGGGGCATGCATCGAGCGCCGCGCGGGCGGCACTCGAtt                                                                                                                                      <  5:80114/149‑1 (MQ=255)
              cAAAAAAGGGCGCTTTTGACCCGGGCGCGGGTTTTGAACCGAAATGGGGCGCAGGCCGGTGCCTGGCGGGAGAAAAAGGGAGTATTCCGGGTGCATGACGTGGGGCGGAGGGGGCATGCATCGAGCGCCGCGCGGGCGGCACTCGATTg                                                                                                                                     >  6:310498/1‑149 (MQ=255)
               aaaaaaGGGCGCTTTTGACCCGGGCGCGGGTTTTGAACCGAAATGGGGCGCAGGCCGGTGCCTGGCGGGAGAAAAAGGGAGTATTCCGGGTGCATGACGTGGGGCGGAGGGGGCATGCATCGAGCGCCGCGCGGGCGGCACTCGATTgc                                                                                                                                    >  4:19170/1‑149 (MQ=255)
               aaaaaaGGGCGCTTTTGACCAGGGAGCGGGATTTGAACCGAAATGGGGCGCAGGCCGGAGCCTGTCGGGAGAAAAAGGGAGTATTGAGGTTGCATGACGTGGGGCGGAGCGGGCATGAATCGAGCGCAGCGCGGGCGGCACTCGATTgg                                                                                                                                    >  1:240194/1‑148 (MQ=255)
                 aaaaGGGCGCTTTTGACCCGGGCGCGGGTTTTGAACCGAAATGGGGCGCAGGCCGGTGCCTGGCGGGAGAAAAAGGGAGTATTCCGGGTGCATGACGTGGGGCGGAGGGGGCATGCATCGAGCGCCGCGCGGGCGGCACTCGATTgcgc                                                                                                                                  <  7:163802/149‑1 (MQ=255)
                       gcgcTTTTGACCCGGGCGCGGGTTTTGAACCGAAATGGGGCGCAGGCCGGTGCCTGGCGGGAGAAAAAGGGAGTATTCCGGGTGCATGACGTGGGGCGGAGGGGGCATGCATCGAGCGCCGCGCGGGCGGCACTCGATTGCGCGAACGc                                                                                                                            >  1:330803/1‑149 (MQ=255)
                          cTTTTGACCCGGGCGCGGGTTTTGAACCGAAATGGGGCGCAGGCCGGTGCCTGGCGGGAGAAAAAGGGAGTATTCCGGGTGCATGACGTGGGGCGGAGGGGGCATGTATCGAGCGCCGCGCGGGCGGCACTCGATTGCGCGAACGCTg                                                                                                                          >  5:214204/1‑148 (MQ=255)
                                aCCCGGGCGCGGGTTTTGAACCGAAATGGGGCGCAGGCCGGTGCCTGGCGGGAGAAAAAGGGAGTATTCCGGGTGCATGACGTGGGGCGGAGGGGGCATGCATCGAGCGCCGCGCGGGCGGCACTCGATTGCGCGAACGCTGAATATc                                                                                                                    <  6:281049/148‑1 (MQ=255)
                                 cccGGGCGCGGGTTTTGAACCGAAATGGGGCGCAGGCCGGTGCCTGGCGGGAGAAAAAGGGAGTATTCCGGGTGCATGACGTGGGGCGGAGGGGGCATGCATCGAGCGCCGCGCGGGCGGCACTCGATTGCGCGAACGCTGAATATC‑Gt                                                                                                                 <  3:19170/149‑1 (MQ=255)
                                 cccGGGCGCGGGTTTTGAACCGAAATGGGGCGCAGGCCGGTGCCTGGCGGGAGAAAAAGGGAGTATTCCGGGTGCATGACGTGGGGCGGAGGGGGCATGCATCGAGCGCCGCGCGGGCGGCACTCGATTGCGCGAACGCTGAATATC‑Gt                                                                                                                 <  1:252312/149‑1 (MQ=255)
                                          gggTTTTGAACCGAAATGGGGCGCAGGCCGGTGCCTGGCGGGAGAAAAAGGGAGTATTCCGGGTGCATGACGTGGGGCGGAGGGGGCATGCATCGAGCGCCGCGCGGGCGGCACTCGATTGCGCGAACGCTGAATATC‑GTGAGACGGgc                                                                                                        <  6:175777/149‑1 (MQ=255)
                                          gggTTTTGAACCGAAATGGGGCGCAGGCCGGTGCCTGGCGGGAGAAAAAGGGAGTATTCCGGGTGCATGACGTGGGGCGGAGGGGGCATGCATCGAGCGCCGCGCGGGCGGCACTCGATTGCGCGAACGCTGAATATC‑GTGAGACGGgc                                                                                                        >  2:221885/1‑149 (MQ=255)
                                          gggTTTTGAACCGAAATGGGGCGCAGGCCGGTGCCTGGCGGGAGAAAAAGGGAGTATTCCGGGTGCATGACGTGGGGCGGAGGGGGCATGCATCGAGCGCCGCGCGGGCGGCACTCGATTGCGCGAACGCTGAATATC‑GTGAGACGGgc                                                                                                        <  5:275952/149‑1 (MQ=255)
                                          gggTTTTGAACCGAAATGGGGCGCAGGCCGGTGCCTGGCGGGAGAAAAAGGGAGTATTCCGGGTGCATGACGTGGGGCGGAGGGGGCATGCAGCGAGCGCCGCGCGGGCGGCACTCGATTGCGCGAACGCTGAATATC‑GTGAGACGGgc                                                                                                        <  6:180401/149‑1 (MQ=255)
                                                 gAACCGAAATGGGGCGCAGGCCGGTGCCTGGCGGGAGAAAAAGGGAGTATTCCGGGTGCATGACGTGGGGCGGAGGGGGCATGCATCGAGCGCCGCGCGGGCGGCACTCGATTGCGCGAACGCTGAATATC‑GTGAGACGGGCGCCTTc                                                                                                  <  5:310498/148‑1 (MQ=255)
                                                  aaCCGAAATGGGGCGCAGGCCGGTGCCTGGCGGGAGAAAAAGGGAGTATTCCGGGTGCATGACGTGGGGCGGAGGGGGCATGCATCGAGCGCCGCGCGGGCGGCACTCGATTGCGCGAACGCTGAATATC‑GTGAGACGGGCGCCTTCa                                                                                                 <  6:39034/148‑1 (MQ=255)
                                                   aCCGAAATGGGGCGCAGGCCGGTGCCTGGCGGGAGAAAAAGGGAGTATTCCGGGTGCATGACGTGGGGCGGAGGGGGCATGCATCGAGCGCCGCGCGGGCGGCACTCGATTGCGCGAACGCTGAATATC‑GTGAGACGGGCGCCTTCAgg                                                                                               >  2:251475/1‑149 (MQ=255)
                                                    ccGAAATGGGGCGCAGGCCGGTGCCTGGCGGGAGAAAAAGGGAGTATTCCGGGTGCAAGACGTGGGGCGGAGGGGGCATGCATCGAGCGCCGCGCGGGCGGCACTCGATTGCGCGAACGCTGAATATC‑GTGAGACGGGCGCCTTCAggg                                                                                              >  1:18266/1‑149 (MQ=255)
                                                       aaaTGGGGCGCAGGCCGGTGCCTGGCGGGAGAAAAAGGGAGTATTCCGGGTGCATGACGTGGGGCGGAGGGGGCATGCATCGAGCGCCGCGCGGGCGGCACTCGATTGCGCGAACGCTGAATATC‑GTGAGACGGGCGCCTTCAGGGCGt                                                                                           <  5:330686/149‑1 (MQ=255)
                                                        aaTGGGGCGCAGGCCGGTGCCTGGCGGGAGAAAAAGGGAGTATTCCGGGTGCATGACGTGGGGCGGAGGGGGCATGCATCGAGCGCCGCGCGGGCGGCACTCGATTGCGCGAACGCTGAATATC‑GTGAGACGGGCGCCTTCAGGGCGTc                                                                                          <  4:295722/149‑1 (MQ=255)
                                                              gcgcAGGCCGGTGCCTGGCGGGAGAAAAAGGGAGTATTCCGGGTGCATGAGGTGGGGCGGAGGGGGCATGCATCGAGCGCCGCGCGGGCGGCTCTCGATTGCGCGAACGCTGAATATC‑GTGAGACGGGCGCCTTCAGGGCGTCAGaggc                                                                                    >  2:283254/1‑149 (MQ=255)
                                                                  aGGCCGGTGCCTGGCGGGAGAAAAAGGGAGTATTCCGGGTGCATGACGTGGGGCGGAGGGGGCATGCATCGAGCGCCGCGCGGGCGGCACTCGATTGCGCGAACGCTGAATATC‑GTGAGACGGGCGCCTTCAGGGCGTCAGaggcaggc                                                                                <  1:221885/149‑1 (MQ=255)
                                                                                      taaaaGGGTCTATTTCGGGTGCATGTCGTTGGGCGGAGGGGGCATGCATGGAGCGGCGCCCGGGCGGCACTGGATTGCGGGTACGCTGAATATCAGTGAGACGGGCGCCTTCAGGGCGTCAGAGGCAGGCGTAGACCACCAGTTccac                                                              <  1:251475/147‑1 (MQ=255)
                                                                                                         tGCATGACGTGGGGCGGAGGGGGCATGCATCGAGCGCCGCGCGGGCGGCACTCGATTGCGCGAACGCTGAATATC‑GTGAGACGGGCGCCTTCAGGGCGTCAGAGGCAGGCGTAGACCACCAGTTCCACCAGCATTTCTTCCCGGGCGAg                                         <  2:330803/149‑1 (MQ=255)
                                                                                                            aTGACGTGGGGCGGAGGGGGCATGCATCGAGCGCCGCGCGGGCGGCACTCGATTGCGCGAACGCTGAATATC‑GTGAGACGGGCGCCTTCAGGGCGTCAGAGGCAGGCGTAGACCACCAGTTCCACCAGCATTTCTTCCCGGGCGAGgcc                                      <  1:278888/149‑1 (MQ=255)
                                                                                                             tGACGTGGGGCGGAGGGGGCATGCATCGAGCGCCGCGCGGGCGGCACTCGATTGCGCGAACGCTGAATATC‑GTGAGACGGGCGCCTTCAGGGCGTCAGAGGCAGGCGTAGACCACCAGTTCCACCAGCATTTCTTCCCGGGCGAGgcca                                     <  6:23902/149‑1 (MQ=255)
                                                                                                             tGACGTGGGGCGGAGGGGGCATGCATCGAGCGCCGCGCGGGCGGCACTCGATTGCGCGAACGCTGAATATC‑GTGAGACGGGCGCCTTCAGGGCGTCAGAGGCAGGCGTAGACCACCAGTTCCACCAGCATTTCTTCCCGGGCGAGgc                                       >  5:248541/1‑147 (MQ=255)
                                                                                                             ggACGTGGGGCGGAGGGGGCATGCATCGAGCGCCGCGCGGGCGGCACTCGATTGCGCGAACGCTGAATATC‑GTGAGACGGGCGCCTTCAGGGCGTCAGAGGCAGGCGTAGACCACCAGTTCCACCAGCATTTCTTCCCGGGCGAGgcca                                     <  2:18266/148‑1 (MQ=255)
                                                                                                               aCGTGGGGCGGAGGGGGCATGCATCGAGCGCCGCGCGGGCGGCACTCGATTGCGCGAACGCTGAATATC‑GTGAGACGGgcgc                                                                                                      >  1:118192/1‑82 (MQ=255)
                                                                                                               aCGTGGGGCGGAGGGGGCATGCATCGAGCGCCGCGCGGGCGGCACTCGATTGCGCGAACGCTGAATATC‑GTGAGACGGgcgc                                                                                                      <  2:118192/82‑1 (MQ=255)
                                                                                                                  tGGGGCTGAGGGGGCATGTATCGAGCGCCGCGCGGGCGGCACTCGATTGCGCGAACGCTGAATATC‑GTGAGACGGGCGCCTTCAGGGCGTCAGAGGCAGGCGTAGACCACCAGTTCCACCAGCATTTCTTCCCGGGCGAGGCCAGccac                                <  6:214204/149‑1 (MQ=255)
                                                                                                                          aGGGGGCATGCATCGAGCGCCGCGCGGGCGGCACTCGATTGCGCGAACGCTGAATATC‑GTGAGACGGGCGCCTTCAGGGCGTCAGAGGCAGGCGTAGACCACCAGTTCCACCAGCATTTCTTCCCGGGCGAGGCCAGCCACCACTACgg                        >  3:16964/1‑149 (MQ=255)
                                                                                                                              ggCATGCATCGAGCGCCGCGCGGGCGGCACTCGATTGCGCGAACGCTGAATATC‑GTGAGACGGGCGCCTTCAGGGCGTCAGAGGCAGGCGTAGACCACCAGTTCCACCAGCATTTCTTCCCGGGCGAGGCCAGccacca                              <  5:218975/139‑1 (MQ=255)
                                                                                                                              ggCATGCATCGAGCGCCGCGCGGGCGGCACTCGATTGCGCGAACGCTGAATATC‑GTGAGACGGGCGCCTTCAGGGCGTCAGAGGCAGGCGTAGACCACCAGTTCCACCAGCATTTCTTCCCGGGCGAGGCCAGccacca                              >  6:218975/1‑139 (MQ=255)
                                                                                                                                   gcatCGAGCGCCGCGCGGGCGGCACTCGATTGCGTGAACGCTGAATATC‑GTGAGACGGGCGCCTTCAGGGCGTCAGAGGCAGGCGTAGACCACCAGTTCCACCAGCATTTCTTCCCGGGCGAGGCCAGCCACCACTACGGCCGCGCGGt               >  6:298305/1‑149 (MQ=255)
                                                                                                                                   gcatCGAGCGCCGCGCGGGCGGCACTCGATTGCGCGAACGCTGAATATC‑GTGAGACGGGCGCCTTCAGGGCGTCAGAGGCAGGCGTAGACCACCAGTTCCACCAGCATTTCTTCCCGGGCGAGGCCAGCCACCACTACGGCCGCGCGGt               >  8:129984/1‑149 (MQ=255)
                                                                                                                                      tCGAGTGCCGCCCGCGCGGCGCTCGATTGCGCGAACGCTGAATATC‑GTGAGACGGGCGCCTTCAGGGCGTCAGAGGCAGGCGTAGACCACCAGTTCCACCAGCATTTCTTCCCGGGCGAGGCCAGCCACCACTACGGCCGCGCGGTTgg            >  6:86826/1‑149 (MQ=255)
                                                                                                                                      tCGAGCGCCGCGCGGGCGGCACTCGATTGCGCGAACGCTGAATATC‑GTGAGACGGGCGCCTTCAGGGCGTCAGAGGCAGGCGTAGACCACCAGTTc                                                                 >  6:323852/1‑96 (MQ=255)
                                                                                                                                      tCGAGCGCCGCGCGGGCGGCACTCGATTGCGCGAACGCTGAATATC‑GTGAGACGGGCGCCTTCAGGGCGTCAGAGGCAGGCGTAGACCACCAGTTc                                                                 <  5:323852/96‑1 (MQ=255)
                                                                                                                                        gAGTGCCGCCCGCGCGGCGCTCGATTGCGCGAACGCTGAATATC‑GTGAGACGGGCGCCTTCAGGGCGTCAGAGGCAGGCGTAGACCACCAGTTCCACCAGCATTTCTTCCCGGGCGAGGCCAGCCACCACTACGGCCGCGCGGTTgggg          >  2:272372/1‑149 (MQ=255)
                                                                                                                                        gAGCGCCGCGCGGGAGGCACTCGATTTCGCGAACGCTGAATATC‑GTGAGACGGGCGCCTTCAGGGCGTCAGAGGCAGGCGTAGACCACCAGTTCCACCAGCATTTCTTCCCGGGCGAGGCCAGCCACCACTACGGCCGCGCGGTTggg           <  6:214742/148‑1 (MQ=255)
                                                                                                                                          gcgcCGCGCGGGCGGCACTCGATTGCGTGAACGCTGAATATC‑GTGAGACGGGCGCCTTCAGGGCGTCAGAGGCAGGCGTAGACCACCAGTTCCACCAGCATTTCTTCCCGGGCGAGGCCAGCCACCACTACGGCCGCGCGGTTGGGGTa        <  5:298305/149‑1 (MQ=255)
                                                                                                                                          gcgcCGCGCGGGCGGCACTCGATTGCGCGAACGCTGAATATC‑GTGAGACGTGCGCCTTCAGGGCGTCAGAGGCAGGCGTAGACCACCAGTTCCAGCAGCATTTCTTCCCGGGCGAGGCCAGCCACCACTACGGCCGCGCGGTTGGGGTa        <  7:129984/149‑1 (MQ=255)
                                                                                                                                          gcgcCGCGCGGGCGGCACTCGATTGCGCGAACGCTGAATATC‑GTGAGACGGGCGCCTTCAGGGCGTCAGAGGCAGGCGTAGACCACCAGTTCCACCAGCATTTCTTCCCGGGCGAGGCCAGCCACCACTACGGCCGCGCGGTTgggg          >  3:318895/1‑147 (MQ=255)
                                                                                                                                          gcgcCGCGCGGGCGGCAATCGATTGCGCGAACGCTGAATATC‑GTGAGACGGGCGCCTTCAGGGCGTCAGAGGCAGGCGTAGACCACCAGTTCCACCAGCATTTCTTCCCGGGCGAGGCCAGCCACCACTACGGCCGCGCGCTTgggg          >  7:148416/1‑147 (MQ=255)
                                                                                                                                                cgcgGGCGGCACTCGATTGCGCGAACGCTGAATATC‑GTGAGACGGGCGCCTTCAGGGCGTCAGAGGCAGGCGTAGACCACCAGTTCCACCAGCATTTCTTCCCGGGCGAGGCCAGCCACCACTACGGCCGCGCGGTTGGGGTAGGGTGc  <  3:98754/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
AAGGGCCGGCGCATCAAAAAAGGGCGCTTTTGACCCGGGCGCGGGTTTTGAACCGAAATGGGGCGCAGGCCGGTGCCTGGCGGGAGAAAAAGGGAGTATTCCGGGTGCATGACGTGGGGCGGAGGGGGCATGCATCGAGCGCCGCGCGGGCGGCACTCGATTGCGCGAACGCTGAATATC‑GTGAGACGGGCGCCTTCAGGGCGTCAGAGGCAGGCGTAGACCACCAGTTCCACCAGCATTTCTTCCCGGGCGAGGCCAGCCACCACTACGGCCGCGCGGTTGGGGTAGGGTGC  >  NC_002947/3079496‑3079788

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