Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,029,019 A→G G317G (GGT→GGC)  rlmI ← 23S rRNA m(5)C1962 methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,029,0190AG96.2% 174.2 / ‑1.9 52G317G (GGT→GGC) rlmI23S rRNA m(5)C1962 methyltransferase
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (23/27);  total (24/28)
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.46e-01

CACGACCGGCATCAATTGCGGCATCCGCGATGATTTTCTGAAATAAATCGCTGGTCATCAGACCGGAACAGGAGAAAGTCAGGAGAATACCGCCTTCATTCAGCAGCTGAATCGCCAGCATGTTGATGTCTTTATAACCACGACACGCGCCCATCAACTGGCTTTTATTCTCAACAAACTTCGGCGGGTCCATCACGATAACGTCAAATTTTTCACCGCGATCGCGATAAGTACGCAGCAATTTAAAGACATCATCACGGACAAACTC  >  NC_000913/1028883‑1029150
                                                                                                                                        |                                                                                                                                   
cACGACCGGCATCAATTGCGGCATCCGCGATGATTTTCTGAAATAAATCGCTGGTCATCAGACCGGAACAGGAGAAAGTCAGGAGAATACCGCCTTCATTCAGCAGCTGAATCGCCAGCATGTTGATGTCTTTATAGCCa                                                                                                                                  >  1:294359/1‑140 (MQ=255)
cACGACCGGCATCAATTGCGGCATCCGCGATGATTTTCTGAAATAAATCGCTGGTCATCAGACCGGAACAGGAGAAAGTCAGGAGAATACCGCCTTCATTCAGCAGCTGAATCGCCAGCATGTTGATGTCTTTATAGCCa                                                                                                                                  >  2:578618/1‑140 (MQ=255)
  cGACCGGCATCAATTGCGGCATCCGCGATGATTTTCTGAAATAAATCGCTGGTCATCAGACCGGAACAGGAGAAAGTCAGGAGAATACCGCCTTCATTCAGCAGCTGAATCGCCAGCATGTTGATGTCTTTATAACCACg                                                                                                                                >  2:903933/1‑140 (MQ=255)
   gACCGGCATCAATTGCGGCATCCGCGATGATTTTCTGAAATAAATCGCTGGTCATCAGACCGGAACAGGAGAAAGTCAGGAGAATACCGCCTTCATTCAGCAGCTGAATCGCCAGCATGTTGATGTCTTTATAGCCACGa                                                                                                                               <  1:325963/140‑1 (MQ=255)
     ccGGCATCAATTGCGGCATCCGCGATGATTTTCTGAAATAAATCGCTGGTCATCAGACCGGAACAGGAGAAAGTCAGGAGAATACCGCCTTCATTCAGCAGCTGAATCGCCAGCATGTTGATGTCTTTATAGCCACGaca                                                                                                                             >  1:1052658/1‑140 (MQ=255)
       ggCATCAATTGCGGCATCCGCGATGATTTTCTGAAATAAATCGCTGGTCATCAGACCGGAACAGGAGAAAGTCAGGAGAATACCGCCTTCATTCAGCAGCTGAATCGCCAGCATGTTGATGTCTTTATAGCCACGACAcg                                                                                                                           <  2:1128418/140‑1 (MQ=255)
        gCATCAATTGCGGCATCCGCGATGATTTTCTGAAATAAATCGCTGGTCATCAGACCGGAACAGGAGAAAGTCAGGAGAATACCGCCTTCATTCAGCAGCTGAATCGCCAGCATGTTGATGTCTTTATAGCCACGACAcgc                                                                                                                          >  1:487446/1‑140 (MQ=255)
           tCAATTGCGGCATCCGCGATGATTTTCTGAAATAAATCGCTGGTCATCAGACCGGAACAGGAGAAAGTCAGGAGAATACCGCCTTCATTCAGCAGCTGAATCGCCAGCATGTTGATGTCTTTATAGCCACGACACGCGcc                                                                                                                       <  2:839521/140‑1 (MQ=255)
             aaTTGCGGCATCCGCGATGATTTTCTGAAATAAATCGCTGGTCATCAGACCGGAACAGGAGAAAGTCAGGAGAATACCGCCTTCATTCAGCAGCTGAATCGCCAGCATGTTGATGTCTTTATAGCCACGACACGCGCCCa                                                                                                                     <  2:294359/140‑1 (MQ=255)
             aaTTGCGGCATCCGCGATGATTTTCTGAAATAAATCGCTGGTCATCAGACCGGAACAGGAGAAAGTCAGGAGAATACCGCCTTCATTCAGCAGCTGAATCGCCAGCATGTTGATGTCTTTATAGCCACGACACGCGCCCa                                                                                                                     >  1:863481/1‑140 (MQ=255)
                  cGGCATCCGCGATGATTTTCTGAAATAAATCGCTGGTCATCAGACCGGAACAGGAGAAAGTCAGGAGAATACCGCCTTCATTCAGCAGCTGAATCGCCAGCATGTTGATGTCTTTATAGCCACGACACGCGCCCATCAAc                                                                                                                <  2:1052658/140‑1 (MQ=255)
                   ggCATCCGCGATGATTTTCTGAAATAAATCGCTGGTCATCAGACCGGAACAGGAGAAAGTCAGGAGAATACCGCCTTCATTCAGCAGCTGAATCGCCAGCATGTTGATGTCTTTATAGCCACGACACGCGCCCATCAACt                                                                                                               >  2:1055398/1‑140 (MQ=255)
                   ggCATCCGCGATGATTTTCTGAAATAAATCGCTGGTCATCAGACCGGAACAGGAGAAAGTCAGGAGAATACCGCCTTCATTCAGCAGCTGAATCGCCAGCATGTTGATGTCTTTATAGCCACGACACGCGCCCATCAACt                                                                                                               <  2:204389/140‑1 (MQ=255)
                       tCCGCGATGATTTTCTGAAATAAATCGCTGGTCATCAGACCGGAACAGGAGAAAGTCAGGAGAATACCGCCTTCATTCAGCAGCTGAATCGCCAGCATGTTGATGTCTTTATAGCCACGACACGCGCCCATCAACTGGCt                                                                                                           >  1:284990/1‑140 (MQ=255)
                                  tttCTGAAATAAATCGCTGGTCATCAGACCGGAACAGGAGAAAGTCAGGAGAATACCGCCTTCATTCAGCAGCTGAATCGCCAGCATGTTGATGTCTTTATAGCCACGACACGCGCCCATCAACTGGCTTTTATTCTcaa                                                                                                <  2:581724/140‑1 (MQ=255)
                                            aaatCGCTGGTCATCAGACCGGAACAGGAGAAAGTCAGGAGAATACCGCCTTCATTCAGCAGCTGAATCGCCAGCATGTTGATGTCTTTATAGCCACGACACGCGCCCATCAACTGGCTTTTATTCTCAACAAACTTcgg                                                                                      <  2:863481/140‑1 (MQ=255)
                                              atCGCTGGTCATCAGACCGGAACAGGAGAAAGTCAGGAGAATACCGCCTTCATTCAGCAGCTGAATCGCCAGCATGTTGATGTCTTTATAGCCACGACACGCGCCCATCAACTGGCTTTTATTCTCAACAAACTTcggcg                                                                                    <  2:279633/140‑1 (MQ=255)
                                               tCGCTGGTCATCAGACCGGAACAGGAGAAAGTCAGGAGAATACCGCCTTCATTCAGCAGCTGAATCGCCAGCATGTTGATGTCTTTATAGCCACGACACGCGCCCATCAACTGGCTTTTATTCTCAACAAACTTcggcgg                                                                                   >  2:347006/1‑140 (MQ=255)
                                                cGCTGGTCATCAGACCTGAACAGGAGAAAGTCAGGAGAATACCGCCTTCATTCAGCAGCTGAATCGCCAGCATGTTGATGTCTTTATAGCCACGACACGCGCCCATCAACTGGCTTTTATTCTCAACAAACTTCGGCggg                                                                                  <  2:1045094/140‑1 (MQ=255)
                                                cGCTGGTCATCAGACCGGAACAGGAGAAAGTCAGGAGAATACCGCCTTCATTCAGCAGCTGAATCGCCAGCATGTTGATGTCTTTATAGCCACGACACGCGCCCATCAACTGGCTTTTATTCTCAACAAACTTCGGCggg                                                                                  <  1:856665/140‑1 (MQ=255)
                                                cGCTGGTCATCAGACCGGAACAGGAGAAAGTCAGGAGAATACCGCCTTCATTCAGCAGCTGAATCGCCAGCATGTTGATGTCTTTATAGCCACGACACGCGCCCATCAACTGGCTTTTATTCTCAACAAACTTCGGCggg                                                                                  <  2:284990/140‑1 (MQ=255)
                                                cGCTGGTCATCAGACCGGAACAGGAGAAAGTCAGGAGAATACCGCCTTCATTCAGCAGCTGAATCGCCAGCATGTTGATGTCTTTATAGCCACGACACGCGCCCATCAACTGGCTTTTATTCTCAACAAACTTCGGCggg                                                                                  <  1:578618/140‑1 (MQ=255)
                                                cGCTGGTCATCAGACCGGAACAGGAGAAAGTCAGGAGAATACCGCCTTCATTCAGCAGCTGAATCGCCAGCATGTTGATGTCTTTATAACCACGACACGCGCCCATCAACTGGCTTTTATTCTCAACAAACTTCGGCggg                                                                                  <  1:903933/140‑1 (MQ=255)
                                                  cTGGTCATCAGACCGGAACAGGAGAAAGTCAGGAGAATACCGCCTTCATTCAGCAGCTGAATCGCCAGCATGTTGATGTCTTTATAGCCACGACACGCGCCCATCAACTGGCTTTTATTCTCAACAAACTTCGGCGGGTc                                                                                <  2:993440/140‑1 (MQ=255)
                                                  cTGGTCATCAGACCGGAACAGGAGAAAGTCAGGAGAATACCGCCTTCATTCAGCAGCTGAATCGCCAGCATGTTGATGTCTTTATAGCCACGACACGCGCCCATCAACTGGCTTTTATTCTCAACAAACTTCGGCGGGTc                                                                                >  1:1041986/1‑140 (MQ=255)
                                                  cTGGTCATCAGACCGGAACAGGAGAAAGTCAGGAGAATACCGCCTTCATTCAGCAGCTGAATCGCCAGCATGTTGATGTCTTTATAGCCACGACACGCGCCCATCAACTGGCTTTTATTCTCAACAAACTTCGGCGGGTc                                                                                <  1:171853/140‑1 (MQ=255)
                                                     gTCATCAGACCGGAACAGGAGAAAGTCAGGAGAATACCGCCTTCATTCAGCAGCTGAATCGCCAGCATGTTGATGTCTTTATAGCCACGACACGCGCCCATCAACTGGCTTTTATTCTCAACAAACTTCGGCGGGTCCAt                                                                             >  1:944546/1‑140 (MQ=255)
                                                          caGACCGGAACAGGAGAAAGTCAGGAGAATACCGCCTTCATTCAGCAGCTGAATCGCCAGCATGTTGATGTCTTTATAGCCACGACACGCGCCCATCAACTGGCTTTTATTCTCAACAAACTTCGGCGGGTCCATCACGa                                                                        >  2:875141/1‑140 (MQ=255)
                                                           aGACCGGAACAGGAGAAAGTCAGGAGAATACCGCCTTCATTCAGCAGCTGAATCGCCAGCATGTTGATGTCTTTATAGCCACGACACGCGCCCATCAACTGGCTTTTATTCTCAACAAACTTCGGCGGGTCCATCACGAt                                                                       >  1:921449/1‑140 (MQ=255)
                                                            gACCGGAACAGGAGAAAGTCAGGAGAATACCGCCTTCATTCAGCAGCTGAATCGCCAGCATGTTGATGTCTTTATAGCCACGACACGCGCCCATCAACTGGCTTTTATTCTCAACAAACTTCGGCGGGTCCATCACGATa                                                                      <  2:921449/140‑1 (MQ=255)
                                                                ggAACAGGAGAAAGTCAGGAGAATACCGCCTTCATTCAGCAGCTGAATCGCCAGCATGTTGATGTCTTTATAGCCACGACACGCGCCCATCAACTGGCTTTTATTCTCAACAAACTTCGGCGGGGCCCTCACCATAAAGt                                                                  >  2:73292/1‑140 (MQ=255)
                                                                  aaCAGGAGAAAGTCAGGAGAATACCGCCTTCATTCAGCAGCTGAATCGCCAGCATGTTGATGTCTTTATAGCCACGACACGCGCCCATCAACTGGCTTTTATTCTCAACAAACTTCGGCGGGTCCATCACGATAACGTCa                                                                >  1:137139/1‑140 (MQ=255)
                                                                                  gagaATACCGCCTTCATTCAGCAGCTGAATCGCCAGCATGTTGATGTCTTTATAGCCACGACACGCGCCCATCAACTGGCTTTTATTCTCAACAAACTTCGGCGGGTCCATCACGATAACGTCAAATTTTTCACcgcgat                                                <  2:163125/140‑1 (MQ=255)
                                                                                   agaATACCGCCTTCATTCAGCAGCTGAATCGCCAGCATGTTGATGTCTTTATAGCCACGACACGCGCCCATCAACTGGCTTTTATTCTCAACAAACTTCGGCGGGTCCATCACGATAACGTCAAATTTTTCACcgcgatc                                               <  1:875141/140‑1 (MQ=255)
                                                                                     aaTACCGCCTTCATTCAGCAGCTGAATCGCCAGCATGTTGATGTCTTTATAGCCACGACACGCGCCCATCAACTGGCTTTTATTCTCAACAAACTTCGGCGGGTCCATCACGATAACGTCAAATTTTTCACcgcgatcgc                                             <  1:1055398/140‑1 (MQ=255)
                                                                                       tACCGCCTTCATTCAGCAGCTGAATCGCCAGCATGTTGATGTCTTTATAGCCACGACACGCGCCCATCAACTGGCTTTTATTCTCAACAAACTTCGGCGGGTCCATCACGATAACGTCAAATTTTTCACcgcgatcgcga                                           >  1:748998/1‑140 (MQ=255)
                                                                                       tACCGCCTTCATTCAGCAGCTGAATCGCCAGCATGTTGATGTCTTTATAGCCACGACACGCGCCCATCAACTGGCTTTTATTCTCAACAAACTTCGGCGGGTCCATCACGATAACGTCAAATTTTTCACcgccatcgcgg                                           >  2:201669/1‑139 (MQ=255)
                                                                                       tACCGCCTTCATTCAGCAGCTGAATCGCCAGCATGTTGATGTCTTTATAGCCACGACACGCGCCCATCAACTGGCTTTTATTCTCAACAAACTTCGGCGGGTCCATCACGATAACGGCAAATTTTTTACcgcgagcgcga                                           >  1:862414/1‑140 (MQ=255)
                                                                                               tcattcaGCAGCTGAATCGCCAGCATGTTGATGTCTTTATAGCCACGACACGCGCCCATCAACTGGCTTTTATTCTCAACAAACTTCGGCGGGTCCATCACGATAACGTCAAATTTTTCACcgcgatcgcgat                                          >  1:325796/1‑133 (MQ=255)
                                                                                               tcattcaGCAGCTGAATCGCCAGCATGTTGATGTCTTTATAGCCACGACACGCGCCCATCAACTGGCTTTTATTCTCAACAAACTTCGGCGGGTCCATCACGATAACGTCAAATTTTTCACcgcgatcgcgat                                          <  2:325796/133‑1 (MQ=255)
                                                                                               tcattcaGCAGCTGAATCGCCAGCATGTTGATGTCTTTATAGCCACGACACGCGCCCATCAACTGGCTTTTATTCTCAACAAACTTCGGCGGGTCCATCACGATAACGTCAAATTTTTCACCGCGATCGCGATAAGTACg                                   <  2:748998/140‑1 (MQ=255)
                                                                                               tcattcaGCAGCTGAATCGCCAGCATGTTGATGTCTTTATAGCCACGACACGCGCCCATCAACTGGCTTTTATTCTCAACAAACTTCGGCGGGTCCATCACGATAACGTCAAATTTTTCACCGCGATCGCGATAAGTACg                                   <  2:137139/140‑1 (MQ=255)
                                                                                                 attcaGCAGCTGAATCGCCAGCATGTTGGTGTCTTTATAGCCACGACACGCGCCCATCAACTGGCTTTTATTCTCAACAAACTTCGGCGGGGCCATCACGATAACGTCAAATTTTTTACCCCGATCGCGATAAATACgca                                 >  2:1109311/1‑140 (MQ=255)
                                                                                                  ttcaGCAGCTGAATCGCCAGCATGTTGATGTCTTTATAGCCACGACACGCGCCCATCAACTGGCTTTTATTCTCAACAAACTTCGGCGGGTCCATCACGATAACGTCAAATTTTTCACCGCGATCGCGATAAGTACgcag                                <  2:944546/140‑1 (MQ=255)
                                                                                                       cagcTGAATCGCCAGCATGTTGATGTCTTTATAGCCACGACACGCGCCCATCAACTGGCTTTTATTCTCAACAAACTTCGGCGGGTCCATCACGATAACGTCAAATTTTTCACCGCGATCGCGATAAGTACGCAGCAAtt                           <  1:347006/140‑1 (MQ=255)
                                                                                                          cTGAATCGCCAGCATGTTGATGTCTTTATAGCCACGACACGCGCCCATCAACTGGCTTTTATTCTCAACAAACTTCGGCGGGTCCATCACGATAACGTCAAATTTTTCACCGCGATCGCGATAAGTACGCAGCAATTTaa                        >  2:796889/1‑140 (MQ=255)
                                                                                                           tGAATCGCCAGCATGTTGATGTCTTTATAGCCACGACACGCGCCCATCAACTGGCTTTTATTCTCAACAAACTTCGGCGGGTCCATCACGATAACGTCAAATTTTTTACCCCGATCGCGGGAAGTACGCAGCAATTTaaa                       >  2:1126708/1‑140 (MQ=255)
                                                                                                            gAATCGCCAGCATGTTGATGTCTTTATAGCCACGACACGCGCCCATCAACTGGCTTTTATTCTCAACAAACTTCGGCGGGTCCATCACGATAACGTCAAATTTTTCACCGCGATCGCGGTAAGTACGCAGCAATTTAAAg                      <  1:1126708/140‑1 (MQ=255)
                                                                                                                      cATGTTGATGTCTTTATAGCCACGACACGCGCCCATCAACTGGCTTTTATTCTCAACAAACTTCGGCGGGTCCATCACGATAACGTCAAATTTTTCACCGCGATCGCGATAAGTACGCAGCAATTTAAAGACATCATCAc            >  1:1014951/1‑140 (MQ=255)
                                                                                                                          ttGATGTCTTTATAGCCACGACACGCGCCCATCAACTGGCTTTTATTCTCAACAAACTTCGGCGGGTCCATCACGATAACGTCAAATTTTTCACCGCGATCGCGATAAGTACGCAGCAATTTAAAGACATCATCACGGAc        <  2:487446/140‑1 (MQ=255)
                                                                                                                               gTCTTTATAGCCACGACACGCGCCCATCAACTGGCTTTTATTCTCAACAAACTTCGGCGGGTCCATCACGATAACGTCAAATTTTTCACCGCGATCGCGATAAGTACGCAGCAATTTAAAGACATCATCACGGACAAACt   <  1:73292/140‑1 (MQ=255)
                                                                                                                                tCTTTATAGCCACGACACGCGCCCATCAACTGGCTTTTATTCTCAACAAACTTCGGCGGGTCCATCACGATAACGGCAAAATTTTCACCGCGATCGCGCTAAATACGCAGCAATTTAAAGACATCATCACGGACAAACTc  >  1:650776/1‑140 (MQ=255)
                                                                                                                                        |                                                                                                                                   
CACGACCGGCATCAATTGCGGCATCCGCGATGATTTTCTGAAATAAATCGCTGGTCATCAGACCGGAACAGGAGAAAGTCAGGAGAATACCGCCTTCATTCAGCAGCTGAATCGCCAGCATGTTGATGTCTTTATAACCACGACACGCGCCCATCAACTGGCTTTTATTCTCAACAAACTTCGGCGGGTCCATCACGATAACGTCAAATTTTTCACCGCGATCGCGATAAGTACGCAGCAATTTAAAGACATCATCACGGACAAACTC  >  NC_000913/1028883‑1029150

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