Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,472,204 T→A 17.3% K95M (AAG→ATG) ‡ rlmH ← 23S rRNA (pseudouridine(1915)‑N(3))‑methyltransferase RlmH

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,472,2040TA17.3% 66.2 / 9.2 29K95M (AAG→ATG) ‡rlmH23S rRNA (pseudouridine(1915)‑N(3))‑methyltransferase RlmH
Reads supporting (aligned to +/- strand):  ref base T (10/14);  new base A (2/3);  total (12/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.95e-01

TCAGCGGCGACAGCGACCAACGTTGCTCGGCGCGCGCGCAAACCTCAGGCGCCAGGCCCTCCGGGCCGCCCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGATCAGGCGGGCGACGTCGGCATTCTTGCCACGGGTATTCAGCGGGATTTCCACCAGCTCCAGCGACAGCTCGGC  >  NC_002947/5472058‑5472352
                                                                                                                                                  |                                                                                                                                                    
tCAGCGGCGACAGCGACCAACGTTGCTCGGCGCGCGCGCAAACCTCAGGCGCCAGGCCCTCCGGGCCGCCCACAATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGc                                                                                                                                                   >  2:307683/1‑150 (MQ=255)
    cggcgACAGCGACCAACGTTGCTCGGCGCGCGCGCAAACCTCAGGCGCCAGGCCCTCCGGGCCGCCCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCAtt                                                                                                                                               <  2:40073/150‑2 (MQ=255)
           aGCGACCAACGTTGCTCGGCGCGCGCGCAAACCTCAGGCGCCAGGCCCTCCGGGCCGCCCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCAAGCCATg                                                                                                                                               <  1:4129/143‑1 (MQ=255)
           aGCGACCAACGTTGCTCGGCGCGCGCGCAAACCTCAGGCGCCAGGCCCTCCGGGCCGCCCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCAAGCCATg                                                                                                                                               >  2:4129/1‑143 (MQ=255)
             cGACCAACGTTGCTCGGCGCGCGCGCAAACCTCAGGCGCCAGGCCCTCCGGGCCGCCCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCAAGCCATGGagtaccga                                                                                                                                      <  1:188495/150‑8 (MQ=255)
                ccAACGTTGCTCGGCGCGCGCGCAAACCTCAGGCGCCAGGCCCTCCGGGCCGCCCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGt                                                                                                                                   >  1:508356/1‑150 (MQ=255)
                        gCTCGGCGCGCGCGCAAACCTCAGGCGCCAGGCCCTCCGGGCCGCCCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGc                                                                                                                           <  2:39998/150‑1 (MQ=255)
                           cGGCGCGCGCGCAAACCTCAGGCGCCAGGCCCTCCGGGCCGCCCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGtt                                                                                                                        <  1:658875/150‑1 (MQ=255)
                                  cgcgcAAACCTCAGGCGCCAGGCCCTCCGGGCCGCCCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAgg                                                                                                                 <  2:447777/150‑1 (MQ=255)
                                   gcgcAAACCTCAGGCGCCAGGCCCTCCGGGCCGCCCACCATCAAATTCACCGTGCGCGCATCCAAAGGACTGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGc                                                                                                                <  1:48998/150‑1 (MQ=255)
                                         aCCTCAGGCGCCAGGCCCTCCGGGCCGCCCACCATCACATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATTGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCGAGGCTGCACt                                                                                                          >  1:646646/1‑150 (MQ=255)
                                            tCAGGCGCCAGGCCCTCCGGGCCGCCCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTTCTCTGGCTGCACTTTg                                                                                                       <  2:669673/150‑1 (MQ=255)
                                              aGGCGCCAGGCCCTCCGGGCCGCCCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGtt                                                                                                                        <  1:132859/131‑1 (MQ=255)
                                              aGGCGCCAGGCCCTCCGGGCCGCCCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGtt                                                                                                                        >  2:132859/1‑131 (MQ=255)
                                                                      caccaTCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCa                                                                                                                   <  2:782844/112‑1 (MQ=255)
                                                                      caccaTCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCa                                                                                                                   >  1:782844/1‑112 (MQ=255)
                                                                         caTCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATg                                                                                              >  1:757481/1‑130 (MQ=255)
                                                                         caTCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATg                                                                                              <  2:757481/130‑1 (MQ=255)
                                                                           tCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGATCa                                                                        >  1:143904/1‑150 (MQ=255)
                                                                            cAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGATCAg                                                                       <  2:494900/150‑1 (MQ=255)
                                                                                                  aGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGccc                                                                                     >  2:270204/1‑114 (MQ=255)
                                                                                                  aGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGccc                                                                                     <  1:270204/114‑1 (MQ=255)
                                                                                                        cAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATg                                                                                                                            <  1:671464/69‑1 (MQ=255)
                                                                                                        cAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATg                                                                                                                            >  2:671464/1‑69 (MQ=255)
                                                                                                          gCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGATCAGGCGGGCGACGTCGGCATTCTTGCCACGGGt                                         <  1:790935/150‑1 (MQ=255)
                                                                                                                                   tCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGATCAGGCGGGCGACGTCGGCATTCTTGCCACGGGTATTCAGCGGGATTTCCACCAGCTcc                >  2:714646/1‑150 (MQ=255)
                                                                                                                                          tCCATGGCAAGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGATCAGGCGGGCGACGTCGGCATTCTTGCCACGGGTATTCAGCGGGATTTCCACCAGCTCCAGc             >  1:639381/1‑146 (MQ=255)
                                                                                                                                          tCCATGGCAAGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGATCAGGCGGGCGACGTCGGCATTCTTGCCACGGGTATTCAGCGGGATTTCCACCAGCTCCAGc             <  2:639381/146‑1 (MQ=255)
                                                                                                                                                 cTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGATCAGGCGGGCGACGTCGGCATTCTTGCCACGGGTATTCAGCGGGATTTCCACCAGCTCCAGCGACAGCTcggc  <  2:719040/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
TCAGCGGCGACAGCGACCAACGTTGCTCGGCGCGCGCGCAAACCTCAGGCGCCAGGCCCTCCGGGCCGCCCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGATCAGGCGGGCGACGTCGGCATTCTTGCCACGGGTATTCAGCGGGATTTCCACCAGCTCCAGCGACAGCTCGGC  >  NC_002947/5472058‑5472352

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