Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,327,213 C→T 19.2% T121I (ACC→ATC)  rluB → 23S rRNA pseudouridine(2605) synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,327,2130CT19.2% 157.3 / 29.1 73T121I (ACC→ATC) rluB23S rRNA pseudouridine(2605) synthase
Reads supporting (aligned to +/- strand):  ref base C (31/28);  new base T (4/10);  total (35/38)
Fisher's exact test for biased strand distribution p-value = 1.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.43e-01

TCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGC  >  NC_000913/1327115‑1327311
                                                                                                  |                                                                                                  
tCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCaccacc                                                                                                   <  2:768732/100‑1 (MQ=255)
 ccAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCg                                                                                                  <  1:1089082/100‑1 (MQ=255)
       gtgtTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGt                                                                                             >  2:2386201/1‑99 (MQ=255)
             gACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCATCGGTGGTGAAcagc                                                                                      <  2:895515/100‑4 (MQ=255)
             gACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCATCGGTGGTGAAcagc                                                                                      <  1:78980/100‑4 (MQ=255)
             gACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCATCGGTGGTGAAcagc                                                                                      <  1:80764/100‑4 (MQ=255)
             gACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCATCGGTGGTGAAcagc                                                                                      <  2:208886/100‑4 (MQ=255)
             gACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCATCGGTGGTGAAcagc                                                                                      <  1:1601321/100‑4 (MQ=255)
             gACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTgg                                                                                      >  1:1079967/1‑100 (MQ=255)
             gACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTgg                                                                                      >  2:1080232/1‑100 (MQ=255)
              aCCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCATCGGTGGTGAAcagca                                                                                     <  1:845602/100‑5 (MQ=255)
              aCCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCATCGGTGGTGAAcagca                                                                                     <  2:1794337/100‑5 (MQ=255)
              aCCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGc                                                                                     <  1:1239907/100‑1 (MQ=255)
              aCCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTAGTCTGCTGCTGTTCACCACCGATGGTGAACTGGc                                                                                     <  1:1666696/100‑1 (MQ=255)
                cGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCaa                                                                                   >  2:363460/1‑100 (MQ=255)
                cGTCGGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGGCGCTTGGACGTTAATAGCTGTGGTCTGCTGCTGCTCACCATCGGTGGTGAAcagcagc                                                                                   <  2:102756/100‑7 (MQ=255)
                 gTCTGCCAAAACTGCGTGGCGCACGCTGGATTTCCGTGTGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCGCCATCGGTGGTGAAcagcag                                                                                    <  2:1536682/98‑6 (MQ=255)
                 gTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTTGGTCGGTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCaaa                                                                                  <  1:720780/100‑1 (MQ=255)
                  tCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCATCGGTGGTGAAcagcagcag                                                                                 <  2:316865/100‑9 (MQ=38)
                  tCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAAc                                                                                 <  1:2334034/100‑1 (MQ=255)
                     gCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGt                                                                              >  1:110073/1‑100 (MQ=255)
                     gCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGt                                                                              >  2:259052/1‑100 (MQ=255)
                     gCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGCGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGt                                                                              >  2:741096/1‑100 (MQ=255)
                      ccAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGtt                                                                             >  1:1733389/1‑100 (MQ=255)
                       cAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGttt                                                                            <  2:1591597/100‑1 (MQ=255)
                       cAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGttt                                                                            <  1:619241/100‑1 (MQ=255)
                        aaaaCTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTa                                                                           <  2:1126622/100‑1 (MQ=255)
                        aaaaCTGCGTGGCGCACGCTGGATCGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGTGGGTGAACTGGCAAACCGtt                                                                             <  2:1453132/98‑1 (MQ=255)
                           aCTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATg                                                                        <  2:1648569/100‑1 (MQ=255)
                               cGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCAcc                                                                    <  2:1709708/100‑1 (MQ=255)
                               cGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCAcc                                                                    >  1:322748/1‑100 (MQ=255)
                                  ggCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCaa                                                                 >  1:1187969/1‑100 (MQ=255)
                                    cgcACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGc                                                               >  1:23775/1‑100 (MQ=255)
                                        cGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTg                                                           <  1:1065171/100‑1 (MQ=255)
                                        cGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTg                                                           >  1:2035678/1‑100 (MQ=255)
                                         gCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGa                                                          >  1:1776696/1‑100 (MQ=255)
                                          cTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGaa                                                         >  2:2332938/1‑100 (MQ=255)
                                          cTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGaa                                                         >  2:2242541/1‑100 (MQ=255)
                                             gATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGtt                                                      >  1:1359902/1‑100 (MQ=255)
                                             gATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGtt                                                      <  2:2288676/100‑1 (MQ=255)
                                             gATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGtt                                                      <  1:819782/100‑1 (MQ=255)
                                                 gCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGaan                                                  <  1:393249/100‑2 (MQ=255)
                                                 gCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAAc                                                  <  1:1926435/100‑1 (MQ=255)
                                                  ccGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGCACTGGCAATACATTTATTTCACCCAAACAGTGAAGTTGAACg                                                 >  2:1673961/1‑100 (MQ=255)
                                                  ccGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACg                                                 <  2:2007207/100‑1 (MQ=255)
                                                  ccGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACg                                                 <  1:1869381/100‑1 (MQ=255)
                                                   cGTGTGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTGTAATGCACCCAAGCCGTGAAGTTGaa                                                   >  2:2268618/1‑97 (MQ=255)
                                                       ggTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGaa                                             <  2:2141763/99‑1 (MQ=255)
                                                        gTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCCAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAata                                           >  2:1223679/1‑100 (MQ=255)
                                                           gCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAACGCACCCAAGCCGTGAAGTTGAACGTGAATATGc                                        <  2:717004/100‑1 (MQ=255)
                                                            cttcgaggatAATACCTCTGGTCTCCTGCTGTTCGCCACCGATGGTGAAGTGGCAATCCGTTTAACGCTCCCAAGCCGGGAAGTTGATCGTGAATATGcc                                       <  2:1745559/91‑1 (MQ=255)
                                                            cTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGcc                                       >  1:990065/1‑100 (MQ=255)
                                                              tGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCcgt                                     <  2:455419/100‑1 (MQ=255)
                                                                  cGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCcgtgcg                                  >  1:1660015/1‑99 (MQ=255)
                                                                   gTTAATAGCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCcgtgcgtg                                <  2:1550627/100‑1 (MQ=255)
                                                                   gTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCcgtgcgtg                                <  2:520788/100‑1 (MQ=255)
                                                                    ttAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCgtgt                               >  1:1316699/1‑100 (MQ=255)
                                                                       aTACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTAtt                            >  1:1006419/1‑100 (MQ=255)
                                                                         aCCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTg                          >  2:1823223/1‑100 (MQ=255)
                                                                          ccTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTgg                         <  2:749631/100‑1 (MQ=255)
                                                                              tgGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAATCCGTTTAATGCACCCAAGTCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTc                       <  2:391275/98‑1 (MQ=255)
                                                                                     ctgctgTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTgacg              >  2:2291238/1‑100 (MQ=255)
                                                                                     ctgctgTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTgacg              >  1:2286020/1‑100 (MQ=255)
                                                                                       gctgTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTgacgac            >  1:1756012/1‑100 (MQ=255)
                                                                                       gctgTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTgacgac            >  2:1019745/1‑100 (MQ=255)
                                                                                       ccagTTCACCATCGGTGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTgacgac            >  1:774235/4‑100 (MQ=255)
                                                                                        ctgTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGacgacg           >  1:2013507/1‑100 (MQ=255)
                                                                                        cagTTCACCATCGGTGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGacgacg           >  1:551622/3‑100 (MQ=255)
                                                                                         tgTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCATGTATTTGGTCAGGTTGACGAcgc          <  2:322746/100‑1 (MQ=255)
                                                                                         agTTCACCATCGGTGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGAcgc          >  1:1785604/2‑100 (MQ=255)
                                                                                         agTTCACCATCGGTGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGAcgc          >  1:1308258/2‑100 (MQ=255)
                                                                                          gTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGAcgcg         <  1:2155643/100‑1 (MQ=255)
                                                                                          gTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGAcgcg         >  1:1077073/1‑100 (MQ=255)
                                                                                            tCACCACCGATGGTGAACCGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGaa       >  2:537256/1‑100 (MQ=255)
                                                                                               ccaccGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACt    >  2:1617896/1‑100 (MQ=255)
                                                                                                caccGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTg   <  1:2332959/100‑1 (MQ=255)
                                                                                                caccGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTg   >  2:1151924/1‑100 (MQ=255)
                                                                                                 accGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGc  >  1:386653/1‑100 (MQ=255)
                                                                                                  |                                                                                                  
TCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGC  >  NC_000913/1327115‑1327311

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