Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 702,766:1 +G 8.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.G 8.9% 165.5 / 15.3 56intergenic (+11/‑121)PP_23SD/PP_5SD23S ribosomal RNA/5S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base . (28/23);  new base G (3/2);  total (31/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.03e-01

ACTGGGATCTTGAATCCCCTAAAGGGCCGTCGAAGACTACGACGTTGATAGGTTGGGTGTGTAAGCGCTGTGAGGCGTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAGCTTCGCAAATATCGCATATCCGAATTCGCTGGGCTGTCCATCTGGACATTCTGGCAACAGAATTTCTTGACGACCATAGAGCATTGGAACC  >  NC_002947/702630‑702914
                                                                                                                                         |                                                                                                                                                    
aCTGGGATCTTGAATCCCCTAAAGGGCCGTCGAAGACTACGACGTTGATAGGTTGGGTGTGTAAGCGCTGTGAGGCGTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGCg                                                                                                                                         >  1:61331/1‑150 (MQ=31)
aCTGGGATCTTGAATCCCCTAAAGGGCCGTCGAAGACTACGACGTTGATAGGTTGGGTGTGTAAGCGCTGTGAGGCGTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGCg                                                                                                                                         >  1:242291/1‑150 (MQ=31)
aCTGGGATCTTGAATCCCCTAAAGGGCCGTCGAAGACTACGACGTTGATAGAAAGGGTGTGTAAGCGCTGTGAGGCGTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGCg                                                                                                                                         >  1:869423/1‑150 (MQ=9)
   gggATCTTGAATCCCCTAAAGGGCCGTCGAAGACTACGACGTTGATAGGTTGGGTGTGTAAGCGCTGTGAGGCGTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGCggtg                                                                                                                                      <  1:2958/150‑1 (MQ=32)
   gggATCTTGAATCCCCTAAAGGGCCGTCGAAGACTACGACGTTGATAGGTTGGGTGTGTAAGCGCTGTGAGGCGTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAACAATCTGCGTAGCAGATTGcttgg                                                                                                                                      <  2:283125/150‑5 (MQ=11)
     gATCTTGAATCCCCTAAAGGGCCGTCGAAGACTACGACGTTGATAGGTTGGGTGTGTAAGCGCTGTGAGGCGTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCGTAGCAGATTGcttgggt                                                                                                                                    <  2:711383/150‑7 (MQ=11)
     gATCTTGAATCCCCTAAAGGGCCGTCGAAGACTACGACGTTGATAGGTTGGGTGTGTAAGCGCTGTGAGGCGTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGCggtggt                                                                                                                                    >  1:498753/1‑150 (MQ=32)
                cccTAAAGGGCCGTCGAAGACTACGACGTTGATAGGTTGGGTGTGTAAGCGCTGTGAGGCGTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAg                                                                                                                         >  2:650993/1‑150 (MQ=32)
                     aaGGGCCGTCGAAGACTACGACGTTGATAGGTTGGGTGTGTAAGCGCTGTGAGGCGTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGACCCg                                                                                                                    >  1:568818/1‑150 (MQ=32)
                           cGTCGAAGACTACGACGTTGATAGGTTGGGTGTGTAAGCGCTGTGAGGCGTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAAtt                                                                                                              >  1:227097/1‑150 (MQ=32)
                              cGAAGACTACGACGTTGATAGGTTGGGTGTGTAAGCGCTGTGAGGCGTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGAACATATAACACCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGt                                                                                                           <  2:278108/150‑1 (MQ=17)
                               gAAGACTACGACGTTGATAGGTTGGGTGTGTAAGCGCTGTGAGGCGTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTg                                                                                                          >  2:860030/1‑150 (MQ=32)
                                  gACTACGACGTTGATAGGTTGGGTGTGTAAGCGCTGTGAGGCGTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCtt                                                                                                       <  2:242291/150‑1 (MQ=32)
                                    cTACGACGTTGATAGGTTGGGTGTGTAAGCGCTGTGAGGCGTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTa                                                                                                     <  2:568818/150‑1 (MQ=32)
                                      acgacgTTGATAGGTTGGGTGTGTAAGCGCTGTGAGGCGTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTaca                                                                                                   <  1:650993/150‑1 (MQ=32)
                                       cgacgTTGATAGGTTGGGTGTGTAAGCGCTGTGAGGCGTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTacac                                                                                                  >  1:25962/1‑150 (MQ=32)
                                              gATAGGTTGGGTGTGTAAGCGCTGTGAGGCGTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGAcc                                                                                                                      >  1:179871/1‑123 (MQ=33)
                                              gATAGGTTGGGTGTGTAAGCGCTGTGAGGCGTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGAcc                                                                                                                      <  2:179871/123‑1 (MQ=33)
                                              gATAGGTTGGGTGTGTAAGCGCTGTGAGGCGTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAGCt                                                                                           >  1:828981/1‑150 (MQ=32)
                                                   gTTGGGTGTGTAAGCGCTGTGAGGCGTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACATCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAGCTTCGCa                                                                                      >  2:968913/1‑150 (MQ=17)
                                                    ttGGGTGTGTAAGCGCTGTGCGGCGTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAGCTTCGCaa                                                                                     >  2:943042/1‑150 (MQ=17)
                                                         tgtgtAAGCGCTGTGAGGCGTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAGCTTCGCAAATATc                                                                                >  2:463261/1‑150 (MQ=32)
                                                           tgtAAGCGCTGTGAGGCGTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAGCTTCGCAAATATCGc                                                                              <  2:472206/150‑1 (MQ=32)
                                                             tAAGCGCTGTGAGGCGTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGCGGTGGTGaa                                                                                                                                 <  2:135667/97‑1 (MQ=34)
                                                             tAAGCGCTGTGAGGCGTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGCGGTGGTGaa                                                                                                                                 >  1:135667/1‑97 (MQ=34)
                                                                gcgcTGTGAGGCGTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAGCTTCGCAAATATCGCATATc                                                                         <  2:498753/150‑1 (MQ=32)
                                                                    tgtgAGGCGTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCGTAGCAGATTGcttg                                                                                                                                       >  1:420867/1‑81 (MQ=11)
                                                                    tgtgAGGCGTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCGTAGCAGATTGcttg                                                                                                                                       <  2:420867/84‑4 (MQ=11)
                                                                    tgtgAGGCGTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAGCTTCGCAAATCTCGCATATCCGaa                                                                     <  2:565293/150‑1 (MQ=17)
                                                                    tgtgAGGCGTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAGCTTCGCAAATATCGCATATCCGaa                                                                     >  1:714122/1‑150 (MQ=32)
                                                                    tgtgAGGCGTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAGCTTCGCAAATATCGCATATCCGaa                                                                     >  1:116978/1‑150 (MQ=32)
                                                                      tgAGGCGTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAGCCTCGCAAATATCGCATATCCGAAtt                                                                   <  2:752779/150‑1 (MQ=17)
                                                                           cGTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATACGTGCTTTACACAAGAGCTTAGCAAATATCGCATATCCGAATTCGCTg                                                              >  1:567037/1‑150 (MQ=17)
                                                                                     aCCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAGCTTCGCAAATATCGCATATCCGAATTCGCTGGGCTGTCCAt                                                    >  1:255921/1‑150 (MQ=32)
                                                                                      ccAATACTAATTGCCCGTGAGGCTTGACCGTATAACACCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAGCTTCGCAAATATCGCATATCCGAATTCGCTGGGCTGTCCATc                                                   >  1:1009219/1‑150 (MQ=17)
                                                                                       cAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGTTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAGCTTCGCAAATATCGCATATCCGAATTCGCTGGGCTGTCCATCt                                                  <  2:250819/150‑1 (MQ=17)
                                                                                       cAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAGCTTCGCAAATATCGCATATCCGAATTCGCTGGGCTGTCCATCt                                                  <  2:1129404/150‑1 (MQ=32)
                                                                                          tACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAGCTTCGCAAATATCGCATATCCGAATTCGCTGGGCTGTCCATCTGGa                                               >  1:442545/1‑150 (MQ=32)
                                                                                               aTTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAGCTTCGCAAATATCGCATATCCGAATTCGCTGGGCTGTCCATCTGGACATTc                                          <  2:442545/150‑1 (MQ=32)
                                                                                                    ccGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAGCTTCGCAAATATCGCATATCCGAATTCGCTGGGCTGTCCATCTGGACATTCTGGCa                                     >  2:1090358/1‑150 (MQ=32)
                                                                                                     cGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAGCTTCGCAAATATCGCATATCCGAATTCGCTGGGCTGTCCATCTGGACATTCTGGCaa                                    >  2:838724/1‑150 (MQ=32)
                                                                                                              tGACCGTATAACACCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAGCTTCGCAAATATCGCATATCCGAATTCGCTGGGCTGTCCATCTGGACATTCTGGCAACAGAATTTc                           <  2:1009219/150‑1 (MQ=17)
                                                                                                                aCCATATAACACCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAGCTTCGCAAATATCGCATATCCGAATTCGCTGGGCTGTCCATCTGGACATTCTGGCAACAGAATTTCtt                         <  2:20377/150‑1 (MQ=32)
                                                                                                                     ataACACCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGAACCGAAAATTc                                                                                                             >  2:549866/1‑61 (MQ=255)
                                                                                                                     ataACACCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGAACCGAAAATTc                                                                                                             <  1:549866/61‑1 (MQ=255)
                                                                                                                      taACACCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAGCTTCGCAAATATCGCATATCCGAATTCGCTGGGCTGTCCATCTGGACATTCTGGCAACAGAATTTCTTgacgac                   >  2:64669/1‑150 (MQ=32)
                                                                                                                       aaCACCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAGCTTCGCAAATATCGCATATCCGAATTCGCTGGGCTGTCCATCTGGACATTCTgg                                       <  1:198056/129‑1 (MQ=33)
                                                                                                                       aaCACCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAGCTTCGCAAATATCGCATATCCGAATTCGCTGGGCTGTCCATCTGGACATTCTgg                                       >  2:198056/1‑129 (MQ=33)
                                                                                                                            ccAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAGCTTCGCAAATATCGCATATCCGAATTCGCTGGGCTGTCCATCTGGACATTCTGGCAACAGAATTTCTTGACGACCATaga             <  1:838724/150‑1 (MQ=32)
                                                                                                                            caccgCAATCTGCGTAGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAGCTTCGCAAATATCGCATATCCGAATTCGCTGGGCTGTCCATCTGGACATTCTGGCAACAGAATTTCTTGACGACCATaga             >  1:532048/5‑150 (MQ=11)
                                                                                                                              ccgCAATCTGCGTAGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAGCTTCGCAAATATCGCATATCCGAATTCGCTGGGCTGTCCATCTGGACATTCTGGCAACAGAATTTCTTGACGACCATAGAGc           >  1:22417/3‑150 (MQ=2)
                                                                                                                              aaGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAGCTTCGCAAATATCGCATATCCGaa                                                                     <  1:661707/92‑1 (MQ=34)
                                                                                                                              aaGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAGCTTCGCAAATATCGCATATCCGaa                                                                     >  2:661707/1‑92 (MQ=34)
                                                                                                                                   aTCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAGCTTCGCAAATATCGCATATCCGAATTCGCTGGGCTGTCCATCTGGACATTCTGGCAACAGAATTTCTTGACGACCATAGAGCATTgg      >  2:953114/1‑150 (MQ=32)
                                                                                                                                    tCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAGCTTCGCAAATATCGCATATCCGAATTCGCTGGGCTGTCCATCTGGACATTCTGGCAACAGAATTTCTTGACGACCATAGAGCATTGGa     <  2:881489/150‑1 (MQ=32)
                                                                                                                                       gCTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAGCTTCGCAAATATCGCATATCCGAATTCGCTGGGCTGTCCATCTGGACATTCTGGCAACAGAATTTCTTGACGACCATAGAGCATTGGAacc  <  2:61331/150‑1 (MQ=32)
                                                                                                                                        cTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAGCTTCGCAAATATCGCATATCCGAAt                                                                    >  2:4942/1‑83 (MQ=32)
                                                                                                                                        cTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAGCTTCGCAAATATCGCATATCCGAAt                                                                    <  1:4942/83‑1 (MQ=32)
                                                                                                                                         |                                                                                                                                                    
ACTGGGATCTTGAATCCCCTAAAGGGCCGTCGAAGACTACGACGTTGATAGGTTGGGTGTGTAAGCGCTGTGAGGCGTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAGCTTCGCAAATATCGCATATCCGAATTCGCTGGGCTGTCCATCTGGACATTCTGGCAACAGAATTTCTTGACGACCATAGAGCATTGGAACC  >  NC_002947/702630‑702914

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