Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 702,766 C→A 14.6% 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,7660CA14.6% 120.7 / 11.7 48intergenic (+14/‑118)PP_23SD/PP_5SD23S ribosomal RNA/5S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base C (20/21);  new base A (3/4);  total (23/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 = 9.38e-01

ACTGGGATCTTGAATCCCCTAAAGGGCCGTCGAAGACTACGACGTTGATAGGTTGGGTGTGTAAGCGCTGTGAGGCGTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAGCTTCGCAAATATCGCATATCCGAATTCGCTGGGCTGTCCATCTGGACATTCTGGCAACAGAATTTCTTGACGACCATAGAGCATTGGAAC  >  NC_002947/702627‑702910
                                                                                                                                           |                                                                                                                                                
aCTGGGATCTTGAATCCCCTAAAGGGCCGTCGAAGACTACGACGTTGATAGGTTGGGTGTGTAAGCGCTGTGAGGCGTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGc                                                                                                                                         <  1:181987/149‑1 (MQ=11)
 cTGGGATCTTGAATCCCCTAAAGGGCCGTCGAAGACTACGACGTTGATAGGTTTGGTGTGTAAGCGCTGTGAGGCGTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCGTAGCAGATTGct                                                                                                                                        <  2:203969/149‑2 (MQ=2)
    ggATCTTGAATCCCCTAAAGGGCCGTCGAAGACTACGACGTTGATAGGTTGGGTGTGTAAGCGCTGTGAGGCGTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGCgg                                                                                                                                       <  1:226579/147‑1 (MQ=32)
    ggATCTTGAATCCCCTAAAGGGCCGTCGAAGACTACGACGTTGATAGGTTGGGTGTGTAAGCGCTGTGAGGCGTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGCgg                                                                                                                                       <  6:237816/147‑1 (MQ=32)
    ggATCTTGAATCCCCTAAAGGGCCGTCGAAGACTACGACGTTGATAGGTTGGGTGTGTAAGCGCTGTGAGGCGTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCGTAGCAGATTGctt                                                                                                                                       <  1:68832/147‑3 (MQ=2)
    ggATCTTGAATCCCCTAAAGGGCCGTCGAAGACTACGACGTTGATAGGTTGGGTGTGTAAGCGCTGTGAGGCGTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCGTAGCAGATTGctt                                                                                                                                       <  2:49990/147‑3 (MQ=2)
           gAATCCCCTAAAGGGCCGTCGAAGACTACGACGTTGATAGGTTGGGTGTGTAAGCGCTGTGAGGCGTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGCggtggt                                                                                                                                   <  7:46043/144‑1 (MQ=32)
                 ccTAAAGGGCCGTCGAAGACTACGACGTTGATAGGTTGGGTGTGTAAGCGCTGTGAGGCGTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGa                                                                                                                           >  4:195996/1‑146 (MQ=32)
                 ccTAAAGGGCCGTCGAAGACTACGACGTTGATAGGTTGGGTGTGTAAGCGCTGTGAGGCGTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGa                                                                                                                           <  3:195996/146‑1 (MQ=32)
                  cTAAAGGGCCGTCGAAGACTACGACGTTGATAGGTTGGGTGTGTAAGCGCTGTGAGGCGTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGa                                                                                                                       >  7:39787/1‑149 (MQ=32)
                    aaaGGGCCGTCGAAGACTACGACGTTGATAGGTTGGGTGTGTAAGCGCTGTGAGGCGTTGAGCTAAGCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGaaa                                                                                                                         <  5:133115/145‑1 (MQ=17)
                    aaaGGGCCGTCGAAGACTACGACGTTGATAGGTTGGGTGTGTAAGCGCTGTGAGGCGTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGaaa                                                                                                                         >  6:133115/1‑145 (MQ=32)
                       gggCCGTCGAAGACTACGACGTTGATAGGTTGGGTGTGTAAGCGCTGTGAGGCGTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGACCCg                                                                                                                   >  2:210987/1‑148 (MQ=32)
                          ccGTCGAAGACTACGACGTTGATAGGTTGGGTGTGTAAGCGCTGTGAGGCGTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCGTAGCAGATTGctt                                                                                                                                       <  3:81180/125‑3 (MQ=2)
                          ccGTCGAAGACTACGACGTTGATAGGTTGGGTGTGTAAGCGCTGTGAGGCGTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCGTAGCAGATTGctt                                                                                                                                       >  4:76947/1‑123 (MQ=2)
                          ccGTCGAAGACTACGACGTTGATAGGTTGGGTGTGTAAGCGCTGTGAGGCGTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCGTAGCAGATTGctt                                                                                                                                       >  4:81180/1‑123 (MQ=2)
                               gAAGACTACGACGTTGATAGGTTGGGTGTGTAAGCGCTGTGAGGCGTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGt                                                                                                          <  4:115415/149‑1 (MQ=32)
                                     tACGACGTTGATAGGTTGGGTGTGTAAGCGCTGTGAGGCGTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTAGTGCTTTa                                                                                                    >  2:21775/1‑149 (MQ=17)
                                          cgTTGATAGGTTGGGTGTGTAAGCGCTGTGAGGCGTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACaa                                                                                               <  6:186425/149‑1 (MQ=32)
                                           gTTGATAGGTTGGGTGTGTAAGCGCTGTGAGGCGTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAag                                                                                              <  4:8024/149‑1 (MQ=32)
                                              gATAGGTTGGGTGTGTAAGCGCTGTGAGGCGTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAGc                                                                                           >  2:130005/1‑149 (MQ=32)
                                                  ggatGGGTGTGTAAGCGCTGTGAGGCGTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAAATCGTGCTTTACACAAGAGCTTc                                                                                        <  4:227053/145‑1 (MQ=17)
                                                               aGCGCTGTGAGGCGTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGAcc                                                                                                                     >  3:275442/1‑106 (MQ=32)
                                                               aGCGCTGTGAGGCGTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGAcc                                                                                                                     <  4:275442/106‑1 (MQ=32)
                                                                       gAGGCGTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAGCTTCGCAAATATCGCATATCCGAAtt                                                                  >  7:49705/1‑149 (MQ=32)
                                                                        aGGCGTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAGCTTCGCAAATATCGCATATCCGAATTc                                                                 >  6:307533/1‑149 (MQ=32)
                                                                                     aCCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAGCTTCGCAAATATCGCATATCCGAATTCGCTGGGCTGTCCa                                                    >  3:270857/1‑149 (MQ=32)
                                                                                         aTACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAGCTTc                                                                                        <  5:324531/109‑1 (MQ=32)
                                                                                         aTACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAGCTTc                                                                                        >  6:324531/1‑109 (MQ=32)
                                                                                          tACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAGCTTCGCAAATATCGCATATCCGAATTCGCTGGGCTGTCCATCTg                                                >  1:119142/1‑148 (MQ=32)
                                                                                            cTAATTTCCCGTGAGGCTTGACCATATAACACCCAATCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAGCTTCGCAAATATCGCATATCCGAATTCGCTGGGCTGTCCATCTGGa                                              <  5:307533/148‑1 (MQ=17)
                                                                                             tAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAGCTTCGCAAATATCGCATATCCGAATTCGCTGGGCTGTCCATCTGGACa                                            >  7:212514/1‑149 (MQ=32)
                                                                                             tAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAGCTTCGCAAATATCGCATATCCGAATTCGCTGGGCTGTCCATCTGGACa                                            >  7:50571/1‑149 (MQ=32)
                                                                                                           gCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAGCTTCGCAAATATCGCATATCCGAATTCGCTGGGCTGTCCATCTGGACATTCTGGCAACAGaa                              >  4:102225/1‑149 (MQ=32)
                                                                                                           gCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAGCTTCGCAAATATCGCATATCCGAATTCGCTGGGCTGTCCATCTGGACATTCTGGCAACAGaa                              <  2:119142/149‑1 (MQ=32)
                                                                                                                aCCATATAACACCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAGCTTCGCAAATATTGCATATCCGAATTCGCTGGGCTGTCCATCTGGACATTCTGGCAACAGAATTTc                          <  1:130005/148‑1 (MQ=17)
                                                                                                                aCCATATAACACCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAGCTTCGCAAATATCGCATATCCGAATTCGCTGGGCTGTCCATCTGGACATTCTGGCAACAGAATTTCt                         >  6:210467/1‑149 (MQ=32)
                                                                                                                aCCATATAACACCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAGCTTCGCAAATATCGCATATCCGAATTCGCTGGGCTGTCCATCTGGACATTCTGGCAACAGAATTTCt                         >  2:244401/1‑149 (MQ=32)
                                                                                                                   atatAACACCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAGCTTCGCAAATATCGCATATCCGAATTCGCTGGGCTGTCCATCTGGACATTCTGGCAACAGAATTTCTTga                      <  1:210987/149‑1 (MQ=32)
                                                                                                                       aaCACCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAGCTTCGCAAATATCGCATATCCGAATTCGCTGGGCTGTCCATCTGGACATTCTGGCAACAGAATTTCTTgacgac                  >  7:94192/1‑149 (MQ=32)
                                                                                                                         cacCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAGCTTCGCAAATATCGCATATCCGAATTCGCTGGGCTGTCCATCTGGACATTCTGGCAACAGAATTTCTTGACGACCa                <  8:212514/149‑1 (MQ=32)
                                                                                                                            ccAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAGCTTCGCAAATATCGCATATCCGAATTCGCTGGGCTGTCCATCTGGACATTCTGGCAACAGAATTTCTTGACGACCa                >  8:34885/1‑146 (MQ=32)
                                                                                                                            caccgCAATCTGCGTAGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAGCTTCGCAAATATCGCATATCCGAATTCGCTGGGCTGTCCATCTGGACATTCTGGCAACAGAATTTCTTGACGACCa                >  7:127521/5‑146 (MQ=255)
                                                                                                                                 cAATCTGCTACGCAGATTGCGGTTGTGAAGAGGAAAGACCAGAAAATTCGACCTTAACACCATAGCGTCGCAAAACTCGCACATCCGACATCGCTGGGCTGGCCATGTGGCCATTCTGGCAACAGAATTTCTTGACGACCATAGAGCAt        <  6:170289/149‑1 (MQ=11)
                                                                                                                                 cAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAGCTTCGCAAATATCGCATATCCGAATTCGCTGGGCTGTCATCTGGACATTCTGGCAACAGAATTTCTTGACGACCATAGAGCAtt       <  8:39787/149‑1 (MQ=17)
                                                                                                                                  aaTCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAGCTTCGCAAATATCGCATATCCGAATTCGCTGGGCTGTCCATCTGGACATTCTGGCAACAGAATTTCTTGACGACCATAGAGCAtt       <  8:49705/149‑1 (MQ=32)
                                                                                                                                        cTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAGCTTCGCAAATATCGCATATCCGAATTCGCTGGGCTGTCCATCTGGACATTCTGGCAACAGAATTTCTTGACGACCATAGAGCATTGGa    <  8:94192/146‑1 (MQ=31)
                                                                                                                                        cTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAGCTTCGCAAATATCGCATATCCGAATTCGCTGGGCTGTCCATCTGGACATTCTGGCAACAGAATTTCTTGACGACCATAGAGCATTGGAac  >  3:189226/1‑148 (MQ=31)
                                                                                                                                           |                                                                                                                                                
ACTGGGATCTTGAATCCCCTAAAGGGCCGTCGAAGACTACGACGTTGATAGGTTGGGTGTGTAAGCGCTGTGAGGCGTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAGCTTCGCAAATATCGCATATCCGAATTCGCTGGGCTGTCCATCTGGACATTCTGGCAACAGAATTTCTTGACGACCATAGAGCATTGGAAC  >  NC_002947/702627‑702910

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