Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,280,025 G→A 8.4% Q153* (CAG→TAG)  rsuA ← 16S rRNA pseudouridine(516) synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,280,0250GA 8.4% 124.2 / 2.0 48Q153* (CAG→TAG) rsuA16S rRNA pseudouridine(516) synthase
Reads supporting (aligned to +/- strand):  ref base G (20/24);  new base A (2/2);  total (22/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.73e-01

GATGCAGCTCAACCACGTGGTTACCCACGGCGGCGAACATGCGTTTCACCTGATGATAACGCCCTTCGCTGATGGTCAGACGAACCTGCGTTGGGGTAATCACTTCCAGCACCGCAGGCTTAGTGAGATCTTTTTCGTTATGCAGCTGCACGCCTTTAGCAAATTGCTCTGCCGTATCGTCAGCTACAGGTGATTCCAGTGTCACCAGATAGGTCTTCTCGCAATGATGGCGCGGAGAAGTAATGCGGTGCGACCACTGACCATCATCAGTCATCAGCACCAGAC  >  NC_000913/2279878‑2280162
                                                                                                                                                   |                                                                                                                                         
gATGCAGCTCAACCACGTGGTTACCCACGGCGGCGAACATGCGTTTCACCTGATGATAACGCCCTTCGCTGGTGGTCAGACGAACCTGCGTTGGGGTAATCACTTCCAGCACCGCCGGCTAAGTGAGATCTTTTTCGTTATTCAGCTGc                                                                                                                                          >  1:66151/1‑149 (MQ=255)
 aTGCAGCTCAACCACGTGGTTACCCACGGTGGCGAACATGCGTTTCACCTGATGATAACGCCCTTCGCTGATGGTCAGACGAACCTGCGTTGGGGTACTCCCTTCCAGCACCGCAGGCCACGTGAGACCTTTTTCGTTATGCAGCTgac                                                                                                                                         >  2:146481/1‑147 (MQ=255)
               cGTGGTTACCCACGGCGGCGAACATGCGTTTCACCTGATGATAACGCCCTTCGCTTATGGTCAGACGAACCTGCGTTGGGGTAATCACTTCCAGCACCGCAGGCTTAGTGAGATCTTTTTCGTTATGCAGCTGCACGCCTTTAGCAAAt                                                                                                                           <  1:787245/149‑1 (MQ=255)
               cGTGGTTACCCACGGCGGCGAACATGCGTTTCACCTGATGATAACGCCCTTCGCTGATGGTCAGACGAACCTGCGTTGGGGTAATCACTTCCAGCACCGCAGGCTTAGTGAGATCTTTTTCGTTATGCAGCTGCACGCCTTTAGCAAAt                                                                                                                           <  1:721325/149‑1 (MQ=255)
                     tACCCACGGCGGCGAACATGCGTTTCACCTGATGATAACGCCCTTCGCTGATGGTCAGACGAACCTGCGTTGGGGTAATCACTTCCAGCACCGCAGGCTTAGTGAGATCTTTTTCGTTATGCAGCTGCACGCCTTTAGCAAATTGctct                                                                                                                     >  1:894242/1‑149 (MQ=255)
                        ccACGGCGGCGAACATGCGTTTCACCTGATGATAACGCCCTTCGCTGATGGTCAGACGAACCTGCGTTGGGGTAATCACTTCCAGCACCGCAGGCTTAGTGCGATATTTTTCGTTATGACGCTGCACGCCTTGAGCACATTGCTCTTcc                                                                                                                  >  2:456276/1‑149 (MQ=255)
                        ccACGGCGGCGAACATGCGTTTCACCTGATGATAACGCCCTTCGCTGATGGTCAGACGAACCTGCGTTGGGGTAATCACTTCCAGCACCGCAGGCTTAGTGAGATCTTTTTCGTTATGCAGCTGCACGCCTTTAGCAAATTGCTCTGcc                                                                                                                  >  1:813872/1‑149 (MQ=255)
                          aCGGCGGCGAACATGCGTTTCACCTGATGATAACGCCCTTCGCTGATGGTCAGACGAACCTGCGTTGGGGTAATCACTTCCAGCACCGCAGGCTTAGTGAGATCTTTTTCGTTATGCAGCTGCACGCCTTTAGCAAATTGCTCTGCCGt                                                                                                                <  1:323922/149‑1 (MQ=255)
                          aCGGCGGCGAACATGCGTTTCACCTGATGATAACGCCCTTCGCTGATGGTCAGACGAACCTGCGTTGGGGTAATCACTTCCAGCACCGCAGGCTTAGTGAGATCTTTTTCGTTATGCAGCTGCACGCCTTTAGCAAATTGCTCTGCCGt                                                                                                                >  2:841820/1‑149 (MQ=255)
                                gcgAACATGCGTTTCACCTGATGATAACGCCCTTCGCTGATGGTCAGACGAACCTGCGTTGGGGTAATCACTTCCAGCACCGCAGGCTTAGTGAGATCTTTTTCGTTATGCAGCTGCACGCCTTTAGCAAATTGCTCTGCCGTATCGTc                                                                                                          >  1:693040/1‑149 (MQ=255)
                                 cgAACATGCGTTTCACCTGATGATAACGCCCTTCGCTGATGGTCAGACGAACCTGCGTTGGGGTAATCACTTCCAGCACCGCAGGCTTAGTGAGATCTTTTTCGTTATGCAGCTGCACGCCTTTAGCAAATTGCTCTGCCGTATCGTCa                                                                                                         <  1:824517/149‑1 (MQ=255)
                                  gAACATGCGTTTCACCTGATGATAACGCCCTTCGCTGATGGTCAGACGAACCTGCGTTGGGGTAATCACTTCCAGCACCGCAGGCTTAGTGAGATCTTTTTCGTTATGCAGCTGCACGCCTTTAGCAAATTGCTCTGCCGTATCGTCAg                                                                                                        <  2:197274/149‑1 (MQ=255)
                                           tttCACCTGATGATAACGCCCTTCGCTGATGGTCAGACGAACCTGCGTTGGGGTAATCACTTCCAGCACCGCAGGCTTAGTGAGATCTTTTTCGTTATGCAGCTGCACGCCTTTAGCAAATTGCTCTGCCGTATCGTCAGCTACAAGTg                                                                                               >  2:39584/1‑149 (MQ=255)
                                            ttCACCTTATCATAACGCCCTTCGCTGATGGTCAGACGAACCTTCGTTGGGGTAATCACTTCCAGCACCGCAGGCTTAGTGAGATCTTTTTCGTTATGCAGCTGCACGCCTTTAGCAAATTGCTCTGCCGTATCGTCAGCTACAGGTGa                                                                                              <  2:749377/149‑1 (MQ=255)
                                                ccTGATGATAACGCCCTTCGCTGATGGTCAGACGAACCTGCGTTGGGGTAATCACTTCCAGCACCGCAGGCTTAGTGAGATCTTTTGCGTTCTGCAGCTGCACGCCTGTAGCAAATTGCTCTGCGGTATCGTCAGCTACAGGTGATTcc                                                                                          >  2:389328/1‑149 (MQ=255)
                                                   gatgatAACGCCCTTCGCTGATGGTCAGACGAACCTGCGTTGGGGTAATCACTTCCAGCACCGCAGGCTTAGTGAGATCTTTTTCGTTATGCAGCTGCACGCCTTTAGCAAATTGCTCTGCCGTATCGTCAGCTACAGGTGATTCCAgt                                                                                       <  2:121700/149‑1 (MQ=255)
                                                    atgatAACGCCCTTCGCTGATGGTCAGACGAACCTGCGTTGGGGTAATCACTTCCAGCACCGCAGGCTTAGTGAGATCTTTTTCGTTATGCAGCTGCACGCCTGTAGGAAATTGGTCTGCCGACCCGTCCGCAACAGGTGATTccggga                                                                                      >  2:596856/1‑145 (MQ=255)
                                                        tAACGCCCTTCGCTGATGGTCAGACGAACCTGCGTTGGGGTAATCACTTCCAGCACCGCAGGCTTAGTGAGATCTTTTTCGTTATGCAGCTGCACGCCTTTAGCAAATTGCTCTGCCGTATCGTGAGCTACAGGTGATTCCAGTGTCAc                                                                                  <  1:676651/149‑1 (MQ=255)
                                                         aaCGCCCTTCGCTGATGGTCAGACGAACCTGCGTTGGGGTAATCACTTCCAGCACCGCAGGCTTAGTGAGATCTTTTTCGTTATGCAGCTGCACGCCTTTAGCAAATTGCTCTGCCGTATCGTCAGCTACAGGTGATTCCAGTGTCAcc                                                                                 >  1:785117/1‑149 (MQ=255)
                                                           cgcgCTTCGCTGATGGTCAGAGGAACCTGGGTTGGGGTAATCACTTGCAGCACCGCAGGCTTAGTGGGATCTTTTTCGTTATGCAGCTGCACGCCTTTAGCAAATTGCTCTGCGGTGTCGGCGGGTACAGGTGCTTCGCGTGTGACCAg                                                                               >  2:45768/1‑149 (MQ=255)
                                                               cTTCGCTGATGGTCAGACGAACCGGCGTTGGGGTAATCACTTCCAGCACCGCAGGCTTAGTGAGATCTTTTTCGTTATGCAGCTGCACGCCTTTAGCAAATTGCTCTGCCGTATCGTCAGCTACAGGTGATTCCAGTGTCACCAGATAg                                                                           <  1:259217/149‑1 (MQ=255)
                                                                ttCGCTGATGGTCAGACGAACCTGCGTTGGGGTAATCACTTCCAGCACCGCAGGCTTAGTGAGATCTTTTTCGTTATGCAGCTACACGCCTTTAGCAAATTGCTCTGCCGTATCGTCAGCTACAGGTGATTCCAGTGTCACCAGATAgg                                                                          >  2:833038/1‑149 (MQ=255)
                                                                ttCGCTGATGGTCAGACGAACCTGCGTTGGGGTAATCACTTCCAGCACCGCAGGCTTAGTGAGATCTTTTTCGTTATGCAGCTACACGCCTTTAGCAAATTGCTCTGCCGTATCGTCAGCTACAGGTGATTCCAGTGTCACCAGATAgg                                                                          >  2:832698/1‑149 (MQ=255)
                                                                 tCGCTGATGGTCAGACGAACCTGCGTTGGGGTAATCACTTCCAGCACCGCAGGCTTAGTGAGATCTTTTTCGTTATGCAGCTGCACGCCTTTAGCAAATTGCTCTGTCGTATCGTCAGCTACAGGTGATTCCAGTGTCACCAGATAGGt                                                                         >  1:704180/1‑149 (MQ=255)
                                                                 tCGCTGATGGTCAGACGAACCTGCGTTGGGGTAATCACTTCCAGCACCGCAGGCTTAGTGAGATCTTTTTCGTTATGCAGCTGCACGCCTTTAGCAAATTGCTCTGCCGTATCGTCAGCTACAGGTGATTCCAGTGTCACCAGATAGGt                                                                         <  2:728073/149‑1 (MQ=255)
                                                                        tGGTCAGACGAACCTGCGTTGGGGTAATCACTTCCAGCACCGCAGGCTTAGTGAGATCTTTTTCGTTATGCAGCGGCACGCCTTTAGCAAATTGCTCTGCCGTATCGTCAGCTACAGGTGATTCCCGTGTCACCAGATAGGACTTCTCg                                                                  >  1:409306/1‑149 (MQ=255)
                                                                              gACGAACCTGTGTTGGGGTAATCACTTCCAGCCTCGCAGGCTTAGTGATATCTTTTTCGTTATGCAGCTGCACGCCTTTAGCAAATTGCTCTGCCGTATCGTCAGCTACAGGTGATTCCAGTGTCACCAGATAGGTCTTCTCGCAatga                                                            <  2:419290/149‑1 (MQ=255)
                                                                                 gagccTGGGTTGGGGTAAGCACTTCCAGCACGGCAGGCTTAGTGAGATCTTTTTCGTTATGCAGCTGGACGCCTTTAGCAAATTGCTCGGCCGTATCGTCAGCTACAGGTGATTCCAGTGTCACCAGATAGGTCTTCTCGCAATGATgg                                                         <  2:145263/146‑1 (MQ=255)
                                                                                       gCGTTGGGGTAAGCACTGCCGGCACCGCAGGCTTAGTGAGATCTTTTTCGTTATGCAGCTGCACGCCTTTAGCAAATTGCTCTGCCGTATCGTCAGCTACAGGTGATTCCAGTGTCACCAGATAGGTCTTCTCGCAATGATGGCGCGga                                                   <  2:893702/149‑1 (MQ=255)
                                                                                          tgggggTAATCACTTCCAGCACCGCAGGCTTAGTGAGATCTTTTTCGTTATGCAGCTGCACGCCTTTAGCAAATTGCTCTGCCGTATCGTCAGCTACAGGTGATTCCAGTGTCACCAGATAGGTCTTCTCGCAATGATGGCGCGGAGaa                                                <  2:92112/147‑1 (MQ=255)
                                                                                            ggggTAACCACTTCCAGCACCGCAGGCTTATTTAGATCTTTTTCGTTATGCAGCTGCACGCCTTTATCAAATTGCTCGGCTGTATCGTCAGCTACAGGTGATTCCAGTGTCACCAGATAGGTCTTCTCGCAATGATGGCGCGGAGAAGt                                              <  2:224515/149‑1 (MQ=255)
                                                                                              ggTAATCACTTCCAGCACCGCAGGCTTAGTGAGATCTTTTTCGTTATGCAGCTGCACGCCTTTAGCAAATTGCTCTGCCGTATCGTCAGCTACAGGTGATTCCAGTGTCACCAGATAGGTCTTCTCGCAATGATGGCGCGGAGAAGTaa                                            >  2:747185/1‑149 (MQ=255)
                                                                                              ggTAATCACTTCCAGCACCGCAGGCTTAGTGAGATCTTTTGCGTTATGCAGCTGCACGCCTTTAGCAAATTGCTCTGCCGTATCGTCAGCTACAGGTGATTCCAGTGTCACCAGATCGGTCTTCTCGCGATGATGGATCGGAGAAGtga                                            >  2:607199/1‑147 (MQ=255)
                                                                                                  aTCACGTCCAGCACCGCAGGCTTAGTGAGATCTTTTTCGTTATGCAGCTGCACGCCTTTAGCAAATTGCTCTGCCGTATCGTCAGCTACAGGTGATTCCAGTGTCACCAGATAGGTCTTCTCGCAATGATGGCGCGGAGAAGTAAGGCg                                        >  2:581002/1‑149 (MQ=255)
                                                                                                   tCACTTCCAGCACCGCAGGCTTAGTGAGATCTTTTTCGTTATGCAGCTGCACGCCTTTAGCAAATTGCTCTGCCGTATCGTCAGCTACAGGTGATTCCAGTGTCACCAGATAGGTCTTCTCGCAATGATGGCGCGGAGAAGTAATGCgg                                       >  1:595574/1‑149 (MQ=255)
                                                                                                    ccattCCAGCACCGCAGGCTTAGTGAGATCTTTTTCGTTATGAAGCTGCACGCCTTTAGCAAATTGCTCTGCCGTATCGTCAGCTACAGGTGATTCCAGTGTCACCAGATAGGTCTTCTCGCAATGATGGCGCGGAGAAGTAATGCGGt                                      <  2:642361/146‑1 (MQ=255)
                                                                                                           aGCACCGCAGGCTTATGGAGATCTTTTTCGTTATGCAGCTGCACGCCTTTAGCAAATTGCTCTGCCGTATCGTCAGCTACAGGTGATTCCAGTGTCACCAGATAGGTCTTCTCGCAATGATGGCGCGGAGAAGTAATGCGGTGCGACCa                               <  2:391241/149‑1 (MQ=255)
                                                                                                               ccGCAGGCTTAGTGAGATCTTTTTCGTTATGCAGCTGCACGCCTTTAGCAAATTGCTCTGCCGTATCGTCAGCTACAGGTGATTCCAGTGTCACCAGATAGGTCTTCTCGCAATGATGGCGCGGAGAAGTAATGCGGTGCGACCACTGa                           <  1:707926/149‑1 (MQ=255)
                                                                                                               ccGCAGGCTGAGTGAGATCTTTTTCGTTATGCAGCTGCACGCCTTTAGCAAATTGCTCTGCCGTATCGTCAGCTACAGGTGATTCCAGTGTCACCAGATAGGTCTTCTCGCAATGATGGCGCGGAGAAGTAATGCGGTGCGACCACTGa                           <  2:186613/149‑1 (MQ=255)
                                                                                                                cGCACGCTTAGTGAGCTATTGTTCCGTATGTAGCTGCACGCCTTTAGCAAATTGCTCTGCCGTATCGTCAGCTACAGGTGATTCCAGTGTCACCAGATATGTCTTCTCGCAATGATGGCGCGGAGAAGTAATGCGGTGCGACCACTGAc                          <  1:276596/149‑1 (MQ=255)
                                                                                                                      cTTAGTGAGATCCTTTTCGTGATTCAGCTGCACGCCTTTAGCAAAGTGCTCTGCCGTATCGTCAGCTACAGGTGATTCCAGTGTCACCAGATAGGTCTTCTCGCAATGATGGCGCGGAGAAGTAATGCGGTGCGACCACTGACcatcat                    <  2:909801/149‑1 (MQ=255)
                                                                                                                        tAGTGAGATCTTTTTCGTTATGCAGCTGCACGCCTTTAGCAAATTGCTCTGCCGTATCGTCAGCTACAGGTGATTCCAGTGTCACCAGATAGGTCTTCTCGCAATGATGGCGCGGAGAAGTAATGCGGTGCGACCACTGACCAtcatca                  <  1:709620/149‑1 (MQ=255)
                                                                                                                          gTGAGATCTTTTTCGTTATGCAGCTGCACGCCTTTAGCAAATTGCTCTGCCGTATCGTCAGCTACAGGTGATTCCAGTGTCACCAGATAGGTCTTCTCGCAATGATGGCGCGGAGAAGTAATGCGGTGCGACCACTGACCAtcatcagt                >  1:394887/1‑149 (MQ=255)
                                                                                                                           tGAGATCTTTTTCGTTATGCAGCTGCACGCCTTTAGCAAATTGCTCTGCCGTATCGTCAGCTACAGGTGATTCCAGTGTCACCAGATAGGTCTTCTCGCAATGATGGCGCGGAGAAGTAATGCGGTGCGACCACTGAc                          >  2:306440/1‑138 (MQ=255)
                                                                                                                           tGAGATCTTTTTCGTTATGCAGCTGCACGCCTTTAGCAAATTGCTCTGCCGTATCGTCAGCTACAGGTGATTCCAGTGTCACCAGATAGGTCTTCTCGCAATGATGGCGCGGAGAAGTAATGCGGTGCGACCACTGAc                          <  1:306440/138‑1 (MQ=255)
                                                                                                                           tGAGATCTTTTTCGTTATGCAGCTACACGCCTTTAGCAAATTGCTCTGCCGTATCGTCAGCTACAGGTGATTCCAGTGTCACCAGATAGGTCTTCTCGCAATGATGGCGCGGAGAAGTAATGCGGTGCGACCACTGACCAtcatcagtc               <  1:833038/149‑1 (MQ=255)
                                                                                                                           tGAGATCTTTTTCGTTATGCAGCTACACGCCTTTAGCAAATTGCTCTGCCGTATCGTCAGCTACAGGTGATTCCAGTGTCACCAGATAGGTCTTCTCGCAATGATGGCGCGGAGAAGTAATGCGGTGCGACCACTGACCAtcatcagtc               <  1:832698/149‑1 (MQ=255)
                                                                                                                                     ttCGTTATGCAGCTGCACGCCTTTAGCAAATTGCTCTGCCGTATCGTCAGCTACAGGTGATTCCAGTGTCACCAGATAGGTCTTCTCGCAATGATGGCGCGGAGAAGTAATGCGGTGCGACCACTGACCATCATCAGTCATCAGCACCa     <  2:305051/149‑1 (MQ=255)
                                                                                                                                        gTTATGCAGCTGCACGCCTTTAGCAAATTGCTCTGCCGTATCGTCAGCTACAGGTGATTCCAGTGTCACCAGATAGGTCTTCTCGCAATGATGGCGCGGAGAAGTAATGCGGTGCGACCACTGACCATCCTCAGTCATCCGCACCAGAc  >  2:65466/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                         
GATGCAGCTCAACCACGTGGTTACCCACGGCGGCGAACATGCGTTTCACCTGATGATAACGCCCTTCGCTGATGGTCAGACGAACCTGCGTTGGGGTAATCACTTCCAGCACCGCAGGCTTAGTGAGATCTTTTTCGTTATGCAGCTGCACGCCTTTAGCAAATTGCTCTGCCGTATCGTCAGCTACAGGTGATTCCAGTGTCACCAGATAGGTCTTCTCGCAATGATGGCGCGGAGAAGTAATGCGGTGCGACCACTGACCATCATCAGTCATCAGCACCAGAC  >  NC_000913/2279878‑2280162

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