Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,292,913 G→C 6.9% Y141* (TAC→TAG)  rsmI ← 16S rRNA C1402 2'‑O‑ribose methyltransferase, SAM‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,292,9130GC 6.9% 183.8 / 3.4 55Y141* (TAC→TAG) rsmI16S rRNA C1402 2'‑O‑ribose methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base G (24/27);  new base C (2/2);  total (26/29)
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.24e-01

GGGATTCGCCTAATACCGCAACGATATCTTCCAGGCTATCTAACAGACGGTGGGTAGATTCATAAAAAATCAGCGTGCGCGGCTCCGCTTCAATGGCTTTTAGCGCATCACGGCGGCCTTTTGATTTGGCAGGTAAAAAGCCTTCGTAACAGAAACGGTCAGAGGGTAAACCCGCTGCGCTTAACGCAGTGATAGCAGCACACGGCCCGGGTAGCGGCACCACGCGGATCCCCGCTTCACGGCAGGTACGCACCAGATGGTAGCCAGGATCGTTAATTAGCGGCGTTCCG  >  NC_000913/3292768‑3293057
                                                                                                                                                 |                                                                                                                                                
gggATTCGCCTAATACCGCAACGATATCTTCCAGGCTATCTAACAGACGGTGGGTAGATTCATAAAAAATCAGCGTGCGCGGCTCCGCTTCAATGGCTTTTAGCGCATCACGGCGGCCTTTTGATTTGGCAGGTAAAAAGCCTTCGTaa                                                                                                                                               <  4:337853/149‑1 (MQ=255)
gggAATCGCCTAATACCGCAACGATATCTTCCAGGCTATCTAACAGACGGTGGGTAGATTCATAAAAAATCAGCGTGCGCGGCTCCGCTTCAATGGCTTTTAGCGCATCACGGCGGCCTTTTGATTTGGCAGGTAAAAAGCCTTCGTaa                                                                                                                                               <  3:273749/149‑1 (MQ=255)
       gCCTAATACCGCAACGATATCTTCCAGGCTATCTAACAGACGGTGGGTAGATTCATAAAAAATCAGCGTGCGCGGCTCCGCTTCAATGGCTTTTAGCGCATCACGGCGGCCTTTTGATTTGGCAGGTAAAAAGCCTTCGTAACAGAAAc                                                                                                                                        >  7:226892/1‑149 (MQ=255)
        ccTAATACCGCAACGATATCTTCCAGGCTATCTAACAGACGGTGGGTAGATTCATAAAAAATCAGCGTGCGCGGCTCCGCTTCAATGGCTTTTAGCGCATCACGGCGGCCTTTTGATTTGGCAGGTAAAAAGCCTTCGTAACAGAAACg                                                                                                                                       >  3:116286/1‑149 (MQ=255)
               ccGCAACGATATCTTCCAGGCTATCTAACAGACGGTGGGTAGATTCATAAAAAATCAGCGTGCGCGGCTCCGCTTCAATGGCTTTTAGCGCATCACGGCGGCCTTTTGATTTGGCAGGTAAAAAGCCTTCGTAACAGAAACGGTCaga                                                                                                                                 <  5:83097/148‑1 (MQ=255)
                cGCAACGATATCTTCCAGGCTATCTAACAGACGGTGGGTAGATTCATAAAAAATCAGCGTGCGCGGCTCCGCTTCAATGGCTTTTAGCGCATCACGGCGGCCTTTTGATTTGGCAGGTAAAAAGCCTTCGTAACAGAAACGGTCAGAgg                                                                                                                               >  8:118045/1‑149 (MQ=255)
                   aaCGATATCTTCCAGGCTATCTAACAGACGGTGGGTAGATTCATAAAAAATCAGCGTGCGCGGCTCCGCTTCAATGGCTTTTAGCGCATCACGGCGGCCTTTTGATTTGGCAGGTAAAAAGCCTTCGTAACAGAAACGGTCAGAGGGt                                                                                                                             >  6:346532/1‑148 (MQ=255)
                            ttCCAGGCTATCTAACAGACGGTGGGTAGATTCATAAAAAATCAGCGTGCGCGGCTCCGCTTCAATGGCTTTTAGCGCATCACGGCGGCCTTTTGATTTGGCAGGTAAAAAGCCTTCCTAACAGAAACGGTCAGAGGGTAAACCCGCTg                                                                                                                   >  4:199860/1‑149 (MQ=255)
                            ttCCAGGCTATCTAACAGACGGTGGGTAGATTCATAAAAAATCAGCGTGCGCGGCTCCGCTTCAATGGCTTTTAGCGCATCACGGCGGCCTTTTGATTTGGCAGGTAAAAAGCCTTCCTAACAGAAACGGTCAGAGGGTAAACCCGCTg                                                                                                                   >  2:196386/1‑149 (MQ=255)
                            ttCCAGGCTATCTAACAGACGGTGGGTAGATTCATAAAAAATCAGCGTGATAGGCTCCGATTCAATGGATGTTAGCGCAAAACGGCGGCCTTATGATTTGGCATGTAAAAAGCCTTCGTAAAAGAAAAGGGCAGAGGGTAAACCCGCt                                                                                                                    >  6:96118/1‑148 (MQ=255)
                               cAGGCTATCTAACAGACGGTGGGTAGATTCATAAAAAATCAGCGTGCGCGGCTCCGCTTCAATGGCTTTTAGCGCATCACGGCGGCCTTTTGATTTGGCAGGTAAAAAGCCTTCGTAACAGAAACGGTCAGAGGGTAAACCCGCTgcgc                                                                                                                <  6:212080/149‑1 (MQ=255)
                                 ggCTATCTAACAGACGGTGGGTAGATTCATAAAAAATCAGCGTGCGCGGCTCCGCTTCAATGGCTTTTAGCGCATCACGGCGGCCTTTTGATTTGGCAGGTAAAAAGCCTTCGTAAc                                                                                                                                              >  2:189364/1‑117 (MQ=255)
                                   cTATCTAACAGACGGTGGGTAGATTCATAAAAAATCAGCGTGCGCGGCTCCGCTTCAATGGCTTTTAGCGCATCACGGCGGCCTTTTGATTTGGCAGGTAAAAAGCCTTCGTAACAGAAACGGTCAGAGGGTAAACCCGCTGCGCTTaa                                                                                                            <  3:160212/149‑1 (MQ=255)
                                    tATCTAACAGACGGTGGGTAGATTCATAAAAAATCAGCGTGCGCGGCTCCGCTTCAATGGCTTTTAGCGCATCACGGCGGCCTTTTGATTTGGCAGGTAAAAAGCCTTCGTAACAGAAACGGTCAGAGGGTAAACCCGCt                                                                                                                    <  8:223833/140‑1 (MQ=255)
                                    tATCTAACAGACGGTGGGTAGATTCATAAAAAATCAGCGTGCGCGGCTCCGCTTCAATGGCTTTTAGCGCATCACGGCGGCCTTTTGATTTGGCAGGTAAAAAGCCTTCGTAACAGAAACGGTCAGAGGGTAAACCCGCt                                                                                                                    >  7:223833/1‑140 (MQ=255)
                                     aTCTAACAGACGGTGGGTAGATTCATAAAAAATCAGCGTGCGCGGCTCCGCTTCAATGGCTTTTAGCGCATCACGGCGGCCTTTTGATTTGGCAGGTAAAAAGCCTTCGTAACAGAAACGGTCAGAGGGTAAACCCGCTGCGCTTAACg                                                                                                          >  2:311459/1‑149 (MQ=255)
                                      tCTAACAGACGGTGGGTAGATTCATAAAAAATCAGCGTGCGCGGCTCCGCTTCAATGGCTTTTAGCGCATCACGGCGGCCTTTTGATTTGGCAGGTAAAAAGCCTTCGTAACAGAAACGGTCAGAGGGTAAACCCGCTGCGCTTAACg                                                                                                          <  1:311459/148‑1 (MQ=255)
                                      tCTAACAGACGGTGGGTAGATTCATAAAAAATCAGCGTGCGCGGCTCCGCTTCAATGGCTTTTAGCGCATCACGGCGGCCTTTTGATTTGGCAGGTAAAAAGCCTTCGTAACAGAAACGGTCAGAGGGTAAACCCGCTGCGCTTAACg                                                                                                          >  6:178733/1‑148 (MQ=255)
                                       cTAACAGACGGTGGGTAGATTCATAAAAAATCAGCGTGCGCGGCTCCGCTTCAATGGCTTTTAGCGCATCACGGCGGCCTTTTGATTTGGCAGGTAAAAAGCCTTCGTAACAGAAACGGTCAGAGGGTAAACCCGCTGCGCTTAACg                                                                                                          >  8:389008/1‑147 (MQ=255)
                                          aCAGCCGGTGGGTCGATTCATAAAAAATCAGCGTGCGCGGCTCCGCTTCACTGGCTTTTCGCGCATCACGGCGGCCTTTTGATGTGGCAGGTAAAACGCCTTCGTAAGAGAAACGGTCAGAGGGTAAACCCGCTGCGCTTAACGCAg                                                                                                       <  3:188779/147‑1 (MQ=255)
                                            aGACGGTGGGTAGATTCATAAAAAATCAGCGTGCGCGGCTCCGCTTCAATGGCTTTTAGCGCATCACGGCGGCCTTTTGATTTGGCAGGTAAAAAGCCTTCGTAACAGAAACGGTCAGAGGGTAAACCCGCTGCGCTTAACGCAGTGAt                                                                                                   <  2:276847/149‑1 (MQ=255)
                                                 gTGGGTAGATTCATAAAAAATCAGCGTGCGCGGCTCCGCTTCAATGGCTTTTAGCGCATCACGGCGGCCTTTTGATTTGGCAGGTAAAAAGCCTTCGTAACAGAAACGGTCAGAGGGTAAACCCGCTGCGCTTAACGCAGTGATagcag                                                                                              >  4:197241/1‑149 (MQ=255)
                                                      tAGATTCATAAAAAATCAGCGTGCGCGGCTCCGCTTCAATGGCTTTTAGCGCATCACGGCGGCCTTTTGATTTGGCAGGTAAAAAGCCTTCGTAACAGAAACGGTCAGAGGGTAAACCCGCTGCGCTTAACGCAGTGATAGCAGCacac                                                                                         >  8:41993/1‑149 (MQ=255)
                                                          ttCATAAAAAATCAGCGTGCGCGGCTCCACTTCAATGGCTTTTAGCGCATCACGGCGGCCTTTTGATTTGGCAGGTAAATAGCCTTCGTAACAGAAACGGTCAGAGGGTAAACCCGCTGCGCTTAACGCAGTGATAGCAGCACACGGcc                                                                                     <  8:207767/149‑1 (MQ=255)
                                                            cATAAAAAATCAGCGTGCGCGGCTCCGCTTCAATGGCTTTTAGCGCATCACGGCGGCCTTTTGATTTGGCAGGTAAAAAGCCTTCGTAACAGAAACGGTCAGAGGGTAAACCCGCTGCGCTTAACGCAGTGATAGCAGCACACGGccc                                                                                    <  7:118045/148‑1 (MQ=255)
                                                             aTAAAAAATCAGCGTGCGCGGCTCCGCTTCAATGGCTTTTAGCGCATCACGGCGGCCTTTTGATTTGGCAGGTAAAAAGCCTTCGTAACAGAAACGGTCAGAGGGTAAACCCGCTGCGCTTAACGCAGTGa                                                                                                    <  3:56339/131‑1 (MQ=255)
                                                             aTAAAAAATCAGCGTGCGCGGCTCCGCTTCAATGGCTTTTAGCGCATCACGGCGGCCTTTTGATTTGGCAGGTAAAAAGCCTTCGTAACAGAAACGGTCAGAGGGTAAACCCGCTGCGCTTAACGCAGTGa                                                                                                    >  4:56339/1‑131 (MQ=255)
                                                             aTAAAAAATCAGCGTGCGCGGCTCCGCTTCAATGGCTTTTAGCGCATCACGGCGGCCTTTTGATTTGGCAGGTAAAAAGCCTTCGTAACAGAAACGGTCAGAGGGTAAACCCGCTGCGCTTAACGCAGTGa                                                                                                    <  2:174778/131‑1 (MQ=255)
                                                             aTAAAAAATCAGCGTGCGCGGCTCCGCTTCAATGGCTTTTAGCGCATCACGGCGGCCTTTTGATTTGGCAGGTAAAAAGCCTTCGTAACAGAAACGGTCAGAGGGTAAACCCGCTGCGCTTAACGCAGTGa                                                                                                    >  1:174778/1‑131 (MQ=255)
                                                              ttaaaaaTCACCGTGCGCGGCTCCGCTTCAATGGCTATTAGCGCATCACGGCGCCCTTTTGAATTGGCAGGAAAAAAGCCTACGTAACAGAAACCCTCAGAGGGTAAACCCGATCCGCTTAACGCAGTGATAGCAGCACAAGGCCCtgg                                                                                 >  2:17748/3‑149 (MQ=255)
                                                              tAAAAAATCAGCGTGCGCGGCTCCGCTTCAATGGCTTTTAGCGCATCACGGCGGCCTTTTGATTTGGCAGGTAAAAAGCCTTCGTAACAGAAACGGTCAGAGGGTAAACCCGCTGCGCTTAACGCAGTGATAGCAGCACACGGCCCggg                                                                                 >  8:323551/1‑149 (MQ=255)
                                                                        gCGTGCGCGGCTCCGCTTCAATGGCTTTTAGCGCATCACGGCGGCCTTTTGATTTGGCAGGTAAAAAGCCTTCGTAACAGAAACGGTCAGAGGGTAAACCCGCTGCGCTTAACGCAGTGATAGCAGCACACGGCCCGGGTAGCGGcacc                                                                       <  7:389008/149‑1 (MQ=255)
                                                                                 gCTCCGCTTCAATGGCTATTAGCGCATCACGGCGGCCTTTTGATTTGGCAGGTAAAAAGCCTTCGTAACAGAAACGGTCAGAGGGTAAAGGCGCTGCGCTTAACGCAGTGATAGCAGCACACGGCCCGGGTAGCGGCACCACGCGGATc                                                              <  1:17748/149‑1 (MQ=255)
                                                                                     cGCTTCAATGGCTTTTTGCGCATCACGGCGGCCTTTTGATTTGGCAGGTAAAAAGCCTTCGTAACAGAAACGGTCAGAGGGTAAACCCGCTGCGCTTAACGCAGTGATAGCAGCACACGGCCCGGGTAGCGGCACCACGCGGATCCCCg                                                          <  5:96118/149‑1 (MQ=255)
                                                                                     cGCTTCAATGGCTTTTAGCGCATCACGGCGGCCTTTTGATTTGGCAGGTAAAAAGCCTTCGTAACAGAAACGGTCAGAGGGTAAACCCGCTGCGCTTAACGCAGTGATAGCAGCACACGGCCCGGGTAGCGGCACCACGCGGATCCCCg                                                          >  6:53387/1‑149 (MQ=255)
                                                                                     cGCTTCAATGGCTTTTAGCGCATCACGGCGGCCTTTTGATTTGGCAGGTAAAAAGCCTTCCTAACAGAAACGGTCAGAGGGTAAACCCGCTGCGCTTAACGCAGTGATAGCAGCACACGGCCCGGGTAGCGGCACCACGCGGATccc                                                            <  1:196386/147‑1 (MQ=255)
                                                                                     cGCTTCAATGGCTTTTAGCGCATCACGGCGGCCTTTTGATTTGGCAGGTAAAAAGCCTTCCTAACAGAAACGGTCAGAGGGTAAACCCGCTGCGCCTAACGCAGTGATAGCAGCACACGGCCCGGGTAGCGGCACCACGCGGATccc                                                            <  3:199860/147‑1 (MQ=255)
                                                                                      gCTTCAATGGCTTTTAGCGCATCACGGCGGCCTTTTGATTTGGCAGGTAAAAAGCCTTCGTAACAGAAACGGTCAGAGGGTAAACCCGCTGCGCTTAACGCAGTGATAGCAGCACACGGCCCGGGTAGCGGCACCACGCGGATCCCCGc                                                         <  8:226892/149‑1 (MQ=255)
                                                                                      gCTTCAATGGCTTTTAGCGCATCACGGCGGCCTTTTGATTTGGCAGGTAAAAAGCCTTCGTAACAGAAACGGTCAGAGGGTAAACCCGCTGCGCTTAACGCAGTGATAGCAGCACACGGCCCGGGTAGCGGCACCACGCGGATCCCCGc                                                         <  4:116286/149‑1 (MQ=255)
                                                                                          cAATGGCTTTTAGCGCATCACGGCGGCCTTTTGATTTGGCAGGTAAAAAGCCTTCGTAACAGAAACGGTCAGAGGGTAAACCCGCTGCGCTTAACGCAGTGATAGCAGCACACGGCCCGGGTAGCGGCACCACGCGGATCCCCGCtt                                                       <  7:41993/147‑1 (MQ=255)
                                                                                           aaTGGCTTTTAGCGCATCACGGCGGCCTTTTGATTTGGCAGGTAAAAAGCCTTCGTAACAGAAACGGTCAGAGGGTAAACCCGCTGCGCTTAACGCAGTGATAGCAGCACACGGCCCGGGTAGCGGCACCACGCGGATCCCCGCTTCAc                                                    >  5:18995/1‑149 (MQ=255)
                                                                                           aaTGGCTTTTAGCGCATCACGGCGGCCTTTTGATTTGGCAGGTAAAAAGCCTTCGTAACAGAAACGGTCAGAGGGTAAACCCGCTGCGCTTAACGCAGTGATAGCAGCACACGGCCCGGGTAGCGGCACCACGCGGATCCCCGCTTCAc                                                    <  7:323551/149‑1 (MQ=255)
                                                                                             tGGCTTTTAGCGCATCACGGCGGCCTTTTGATTTGGCAGGTAAAAAGCCTTCGTAACAGAAACGGTCAGAGGGTAAACCCGCTGCGCTTAACGCAGTGATAGCAGCACACGGCCCGGGTAGCGGCACCACGCGGATCCCCGCTTCACgg                                                  >  6:112652/1‑149 (MQ=255)
                                                                                               gCTTTTAGCGCATCACGGCGGCCTTTTGATTTGGCAGGTAAAAAGCCTTCGTAACAGAAACGGTCAGAGGGTAAACCCGCTGCGCTTAACGCAGTGATAGCAGCACACGGCCCGGGTAGCGGCACCACGCGGATCCCCGCTTCACGGc                                                 <  5:112652/148‑1 (MQ=255)
                                                                                                cTTTTAGCGCATCACGGCGGCCTTTTGATTTGGCAGGTAAAAAGCCTTCGTAACAGAAACGGTCAGAGGGTAAACCCGCTGCGCTTAACGCAGTGATAGCAGCACACGGCCCGGGTAGCGGCACCACGCGGATCCCCGCTTCACGGCAg                                               <  4:192725/149‑1 (MQ=255)
                                                                                                cTTTTAGCGCATCACGGCGGCCTTTTGATTTGGCAGGTAAAAAGCCTTCGTAACAGAAACGGTCAGAGGGTAAACCCGCTGCGCTTAACGCAGTGATAGCAGCACACGGCCCGGGTAGCGGCACCACGCGGATCCCCGCTTCACGGCAg                                               <  3:197241/149‑1 (MQ=255)
                                                                                                      gcgcATCACGGCGGCCTTTTGATTTGGCAGGTAAAAAGCCTTCGTAACAGAAACGGTCAGAGGGTAAACCCGCTGCGCTTAACGCAGTGATAGCAGCACACGGCCCGGGTAGCGGCACCACGCGGATCCCCGCTTCACGGCAGGTACGc                                         >  2:365/1‑149 (MQ=255)
                                                                                                          aTCACGGCGGCCTTTTGATTTGGCAGGTAAAAAGCCTTCGTAACAGAAACGGTCAGAGGGTAAACCCGCTGCGCTTAACGCAGTGATAGCAGCACACGGCCCGGGTAGCGGCACCACGCGGATCCCCGCTTCACGGCAGGTACGCAc                                       >  7:370300/1‑147 (MQ=255)
                                                                                                                 cggcCTTTTGATTTGGCAGGTAAAAAGCCTTCGTAACAGAAACGGTCAGAGGGTAAACCCGCTGCGCTTAACGCAGTGATAGCAGCACACGGCCCGGGTAGCGGCACCACGCGGATCCCCGCTTCACGGCAGGTACGCACCAGATGGTa                              <  5:53387/149‑1 (MQ=255)
                                                                                                                    ccTTTTGATTTGGCAGGTAAAAAGCCTTCGTAACAGAAACGGTCAGAGGGTAAACCCGCTGCGCTTAACGCAGTGATAGCAGCACACGGCCCGGGTAGCGGCACCACGCGGATCCCCGCTTCACGGCAGGTACGCACCAGATGGTAGcc                           <  5:346532/149‑1 (MQ=255)
                                                                                                                                   ggTAAAAAGCCTTCGTAACAGAAACGGTCAGAGGGTAAACCCGCTGCGCTTAACGCAGTGATAGCAGCACACGGCCCGGGTAGCGGCACCACGCGGATCCCCGCTTCACGGCAGGTACGCACCAGATGGTAGCCAGGATCGTTAATTAg            >  4:226166/1‑149 (MQ=255)
                                                                                                                                     tAAAAAGCCTTCGTAACAGAAACGGTCAGAGGGTAAACCCGCTGCGCTTAACGCAGTGATAGCAGCACACGGCCCGGGTAGCGGCACCACGCGGATCCCCGCTTCACGGCAGGTACGCACCAGATGGTAGCCAGGATCGTTAATTAgcg          >  2:161361/1‑149 (MQ=255)
                                                                                                                                       aaaaGCCTTCGTAACAGAAACGGTCAGAGGGTAAACCCGCTGCGCTTAACGCAGTGATAGCAGCACACGGCCCGGGTAGCGGCACCACGCGGATCCCCGCTTCACGGCAGGTACGCACCAGATGGTAGCCAGGATCGTTAATTAgcggc        <  6:18995/149‑1 (MQ=255)
                                                                                                                                             cTTCGTAACAGAAACGGTCAGAGGGTAAACCCGCTGCGCTTAACGCAGTGATAGCAGCACACGGCCCGGGTAGCGGCACCACGCGGATCCCCGCTTCACGGCAGGTACGCACCAGATGGTAGCCAGGATCGTTAATTAGCGGCGTTCCg  >  4:159192/1‑149 (MQ=255)
                                                                                                                                             cTTCGTAACAGAAACGGTCAGAGGGTAAACCCGCTGCGCTTAACGCAGTGATAGCAGCACACGGCCCGGGTAGCGGCACCACGCGGATCCCCGCTTCACGGCAGGTACGCACCAGATGGTAGCCAGGATCGTTAATTAGCGGCGTTCCg  <  8:370300/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
GGGATTCGCCTAATACCGCAACGATATCTTCCAGGCTATCTAACAGACGGTGGGTAGATTCATAAAAAATCAGCGTGCGCGGCTCCGCTTCAATGGCTTTTAGCGCATCACGGCGGCCTTTTGATTTGGCAGGTAAAAAGCCTTCGTAACAGAAACGGTCAGAGGGTAAACCCGCTGCGCTTAACGCAGTGATAGCAGCACACGGCCCGGGTAGCGGCACCACGCGGATCCCCGCTTCACGGCAGGTACGCACCAGATGGTAGCCAGGATCGTTAATTAGCGGCGTTCCG  >  NC_000913/3292768‑3293057

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: ‑