Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,332,076 C→T 63.6% intergenic (+97/‑196) PP_16SE → / → PP_23SE 16S ribosomal RNA/23S ribosomal RNA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,332,0760CT63.6% ‑1.1 / 19.8 22intergenic (+97/‑196)PP_16SE/PP_23SE16S ribosomal RNA/23S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base C (4/4);  new base T (5/9);  total (9/13)
Fisher's exact test for biased strand distribution p-value = 6.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.23e-01

TAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGG  >  NC_002947/1331949‑1332169
                                                                                                                               |                                                                                             
tAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGAATg                                                                                       >  2:236956/1‑136 (MQ=17)
     gAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGAt                                                                                  >  2:213034/1‑136 (MQ=17)
     gAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGAt                                                                                  >  2:1256/1‑136 (MQ=17)
          tGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTg                                                                             >  1:333309/1‑136 (MQ=17)
           gCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTgg                                                                            <  1:392656/136‑1 (MQ=17)
            cGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGc                                                                           <  2:333309/136‑1 (MQ=17)
                      aCCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGATTTCTGGCTTTTGTCAGa                                                                 >  2:247141/1‑136 (MQ=32)
                                tCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTa                                                       >  2:344859/1‑136 (MQ=17)
                                        aTCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCg                                               >  2:317515/1‑136 (MQ=32)
                                        aTCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCg                                               <  1:317515/136‑1 (MQ=32)
                                          gtgCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGa                                             <  1:186184/134‑1 (MQ=17)
                                              ttgctgATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATAtg                                         <  1:263206/135‑1 (MQ=17)
                                              ctgctgATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATAtg                                         <  1:87924/136‑1 (MQ=32)
                                              ctgctgATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATAtg                                         <  1:190655/136‑1 (MQ=32)
                                              atgctgATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATAtg                                         <  1:73733/135‑1 (MQ=17)
                                                             cacaCGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGATTTCTGACTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACt                          >  2:417077/1‑136 (MQ=17)
                                                              acacGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTg                         <  1:213034/136‑1 (MQ=17)
                                                              acacGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTg                         <  1:1256/136‑1 (MQ=17)
                                                                acGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGaa                       <  2:289784/136‑1 (MQ=32)
                                                                 cGAATTGATTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAac                      <  2:75545/136‑1 (MQ=17)
                                                                      tGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAg                 <  2:220478/136‑1 (MQ=17)
                                                                                     cGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCtggtgg  >  1:276421/1‑136 (MQ=31)
                                                                                                                               |                                                                                             
TAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGG  >  NC_002947/1331949‑1332169

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