Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 702,769 Δ1 bp 17.2% intergenic (+14/‑118) PP_23SD → / → PP_5SD 23S ribosomal RNA/5S ribosomal RNA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947702,7690C.17.2% 68.5 / 13.4 29intergenic (+14/‑118)PP_23SD/PP_5SD23S ribosomal RNA/5S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base C (12/11);  new base . (2/3);  total (14/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.78e-01

CTGGGATCTTGAATCCCCTAAAGGGCCGTCGAAGACTACGACGTTGATAGGTTGGGTGTGTAAGCGCTGTGAGGCGTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAGCTTCGCAAATATCGCATATCCGAATTCGCTGGGCTGTCCATCTGGACATTCTGGCAACAGAATTTCTTGACGACCATAGAGCATTGGAACCAC  >  NC_002947/702631‑702916
                                                                                                                                           |                                                                                                                                                   
cTGGGATCTTGAATCCCCTAAAGGGCCGTCGAAGACTACGACGTTGATAGGTTGGGTGTGTAAGCGCTGTGAGGCGTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCGTAGCAGATTGctt                                                                                                                                          <  1:159857/150‑3 (MQ=2)
  gggATCTGGAATCCCCTAAAGGGCCGTCGAAGACTACGACGTTGATAGGGTGGGTGTGTAAGCGCTGTGAGGCGTTGAGCTAACCAACACTAATTGGCCGTGAGGCTTGACCATATAACACCCCCACAATCTGCGTAGCAGATTGcttgg                                                                                                                                        <  2:65807/150‑5 (MQ=11)
        ttGAATCCCCTAAAGGGCCGTCGAAGACTACGACGTTGATAGGTTGGGTGTGTAAGCGCTGTGAGGCGTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAg                                                                                                                                  >  2:398903/1‑150 (MQ=32)
         tGAATCCCCTAAAGGGCCGTCGAAGACTACGACGTTGATAGGTTGGGTGTGTAAGCGCTGTGAGGCGTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGa                                                                                                                                 <  1:398903/150‑1 (MQ=32)
                                                   ttGGGTGTGTAAGCGCTGTGAGGCGTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCGTAGCAGATTGct                                                                                                                                           <  2:466191/98‑2 (MQ=11)
                                                   ttGGGTGTGTAAGCGCTGTGAGGCGTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCGTAGCAGATTGct                                                                                                                                           >  1:466191/1‑97 (MQ=11)
                                                      ggTGTGTAAGCGCTGTGAGGCGTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAGCTTCGCAAata                                                                                    >  2:39947/1‑150 (MQ=32)
                                                            tAAGCGCTGTGAGGCGTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGAc                                                                                                                                <  1:281808/100‑1 (MQ=34)
                                                            tAAGCGCTGTGAGGCGTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGAc                                                                                                                                >  2:584789/1‑100 (MQ=33)
                                                            tAAGCGCTGTGAGGCGTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGAc                                                                                                                                <  1:584789/100‑1 (MQ=34)
                                                            tAAGCGCTGTGAGGCGTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGAc                                                                                                                                >  2:281808/1‑100 (MQ=34)
                                                                       aGGCGTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAGCTTCGCAAATATCGCATATCCGAATTCg                                                                   >  1:40196/1‑150 (MQ=32)
                                                                         gCGTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAACAATCTGCGTAGCAGATTGctt                                                                                                                                          >  1:401231/1‑75 (MQ=11)
                                                                         gCGTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAACAATCTGCGTAGCAGATTGctt                                                                                                                                          <  2:401231/77‑3 (MQ=11)
                                                                               aGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAGCTTCGCAAATATCGCATATCCGAATTCGCTGGGCTg                                                           >  2:429158/1‑150 (MQ=32)
                                                                                    aCCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTg                                                                                                            >  1:103030/1‑96 (MQ=34)
                                                                                    aCCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTg                                                                                                            <  2:103030/96‑1 (MQ=34)
                                                                                                   ccGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAGCTTCGCAAATATCGCATATCCGAATTCGCTGGGCTGTCCATCt                                                    <  1:540823/137‑1 (MQ=33)
                                                                                                   ccGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAGCTTCGCAAATATCGCATATCCGAATTCGCTGGGCTGTCCATCt                                                    >  2:540823/1‑137 (MQ=33)
                                                                                                      tGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAGCTTCGCAAATATCGCATATCCGAATTCGCTGGGCTGTCCATCTGGACATTCTGGCAACa                                    <  2:40196/150‑1 (MQ=32)
                                                                                                          gCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAGCTTCGCAAATATCGCATATCCGAATTCGCTGGGCTGTCCATCTGGACATTCTGGCAACAGAAt                                >  2:235348/1‑150 (MQ=32)
                                                                                                               aCCATATAACACCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAGCTTCGCAAATATCGCATATCCGAAt                                                                      <  2:429069/107‑1 (MQ=33)
                                                                                                               aCCATATAACACCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAGCTTCGCAAATATCGCATATCCGAATTCGCTGGGCTGTCCATCTGGACATTCTGGCAACAGAATTTCtt                           <  1:429158/150‑1 (MQ=32)
                                                                                                               aCCATATAACACACAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAGCTTCGCAAATATCGCATATCCGAAt                                                                      >  1:429069/1‑107 (MQ=17)
                                                                                                                ccATATAACACCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAGCTTCGCAAATATCGCATATCCGAATTCGCTGGGCTGTCCATCTGGACATTCTGGCAACAGAATTTCTTg                          >  2:512825/1‑150 (MQ=32)
                                                                                                                         acCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAGCTTCGCAAATATCGCATATCCGAATTCGCTGGGCTGTCCATCTGGACATTCTGGCAACAGAATTTCTTGACGACCATa                 >  1:577929/1‑150 (MQ=32)
                                                                                                                         acCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAGCTTCGCAAATATCGCATATCCGAATTCGCTGGGCTGTCCATCTGGACATTCTGGCAACAGAATTTCTTGACGACCATa                 <  1:512825/150‑1 (MQ=32)
                                                                                                                             aaGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAGCTTCGCAAATATCGCATATCCGAATTCGCTGGGCTGTCCATCTGGACATTCTGGCAACAGAATTTCTTGACGACCATAGAGc             <  1:235348/150‑1 (MQ=32)
                                                                                                                                         tACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAGCTTCGCAAATATCGCATATCCGAATTCGCTGGGCTGTCCATCTGGACATTCTGGCAACAGAATTTCTTGACGACCATAGAGCATTGGAaccac  <  1:374208/150‑1 (MQ=11)
                                                                                                                                           |                                                                                                                                                   
CTGGGATCTTGAATCCCCTAAAGGGCCGTCGAAGACTACGACGTTGATAGGTTGGGTGTGTAAGCGCTGTGAGGCGTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAGCTTCGCAAATATCGCATATCCGAATTCGCTGGGCTGTCCATCTGGACATTCTGGCAACAGAATTTCTTGACGACCATAGAGCATTGGAACCAC  >  NC_002947/702631‑702916

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