Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 702,766:1 +G 17.9% intergenic (+11/‑121) 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,7661.G17.9% 71.9 / 17.2 28intergenic (+11/‑121)PP_23SD/PP_5SD23S ribosomal RNA/5S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base . (12/11);  new base G (2/3);  total (14/14)
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

CTGGGATCTTGAATCCCCTAAAGGGCCGTCGAAGACTACGACGTTGATAGGTTGGGTGTGTAAGCGCTGTGAGGCGTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAGCTTCGCAAATATCGCATATCCGAATTCGCTGGGCTGTCCATCTGGACATTCTGGCAACAGAATTTCTTGACGACCATAGAGC  >  NC_002947/702631‑702905
                                                                                                                                        |                                                                                                                                           
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)
                                                                                                                                        |                                                                                                                                           
CTGGGATCTTGAATCCCCTAAAGGGCCGTCGAAGACTACGACGTTGATAGGTTGGGTGTGTAAGCGCTGTGAGGCGTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAGCTTCGCAAATATCGCATATCCGAATTCGCTGGGCTGTCCATCTGGACATTCTGGCAACAGAATTTCTTGACGACCATAGAGC  >  NC_002947/702631‑702905

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: