Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 729,087 C→T 7.1% H191H (CAC→CAT)  rluD → 23S rRNA pseudouridine(1911/1915/1917) synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947729,0870CT 7.1% 161.3 / 3.1 56H191H (CAC→CAT) rluD23S rRNA pseudouridine(1911/1915/1917) synthase
Reads supporting (aligned to +/- strand):  ref base C (26/26);  new base T (2/2);  total (28/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 = 6.29e-01

CGGCCTGATGGTGGTGGCCAAGACCCTGCAGGCGCAGACCAAGCTGGTCGAGCAGATGCAGGCCCGCAAGGTCAGCCGCATCTACGAGTGCATCGTGATTGGCGTGGTGACTGCCGGTGGCAAGATCGACGCGCCCATCGGCCGCCACGGCGGCATGCGCCAGCGCATGGCTGTAACCGACGGCGGCAAGCCGGCGGTCAGTCACTACCGCGTACTCAAGCGCTTCCGCTCGCACACCCATGTGCGGGTCAAGCTGGAAACCGGCCGTACCCACCAGATTCGCGTGCACATGGCCCA  >  NC_002947/728940‑729236
                                                                                                                                                   |                                                                                                                                                     
cGGCCTGATGGTGGTGGCCAAGACCCTGCAGGCGCAGACCAAGCTGGTCGAGCAGATGCAGGCCCGCAAGGTCAGCCGCATCTACGAGTGCATCGTGATTGGCGTGGTGACTGCCGGTGGCAAGATCGACGCGCCCATCGGCCGCCAcgg                                                                                                                                                     <  1:350312/150‑1 (MQ=255)
        tggtggtggCCAAGACCCTGCAGGCGCAGACCAAGCTGGTCGAGCAGATGCAGGCCCGCAAGGTCAGCCGCATCTACGAGTGCATCGTGATTGGCGTGGTGACTGCCGGTGGCAAGATCGACGCGCCCATCGGCCGCCACGGCGGCATgc                                                                                                                                             <  1:722594/150‑1 (MQ=255)
        tggtggtggCCAAGACCCTGCAGGCGCAGACCAAGCTGGTCGAGCAGATGCAGGCCCGCAAGGTCAGCCGCATCTACGAGTGCATCGTGATTGGCGTGGTGACTGCCGGTGGCAAGATCGACGCGCCCATCGGCCGCCACGGCGGCATgc                                                                                                                                             <  2:913219/150‑1 (MQ=255)
                  cAAGACCCTGCAGGCGCAGACCAAGCTGGTCGAGCAGATGCAGGCCCGCAAGGTCAGCCGCATCTACGAGTGCATCGTGATTGGCGTGGTGACTGCCGGTGGCAAGATCGACGCGCCCATCGGCCGCCACGGCGGCATGCGCCAGCGCAt                                                                                                                                   >  1:890722/1‑150 (MQ=255)
                    aGACCCTGCAGGCGCAGACCAAGCTGGTCGAGCAGATGCAGGCCCGCAAGGTCAGCCGCATCTACGAGTGCATCGTGATTGGCGTGGTGACTGCCGGTGGCAAGATCGACGCGCCCATCGGCCGCCACGGCGGCATGCGCCAGCGCATgg                                                                                                                                 <  1:924858/150‑1 (MQ=255)
                     gACCCTGCAGGCGCAGACCAAGCTGGTCGAGCAGATGCAGGCCCGCAAGGTCAGCCGCATCTACGAGTGCATCGTGATTGGCGTGGTGACTGCCGGTGGCAAGATCGACGCGCCCATCGGCCGCCACGGCGGCATGCGCCAGCGCATgg                                                                                                                                 >  1:828138/1‑149 (MQ=255)
                     gACCCTGCAGGCGCAGACCAAGCTGGTCGAGCAGATGCAGGCCCGCAAGGTCAGCCGCATCTACGAGTGCATCGTGATTGGCGTGGTGACTGCCGGTGGCAAGATCGACGCGCCCATCGGCCGCCACGGCGGCATGCGCCAGCGCATgg                                                                                                                                 <  2:828138/149‑1 (MQ=255)
                     gACCCTGCAGGCGCAGACCAAGCTGGTCGAGCAGATGCAGGCCCGCAAGGTCAGCCGCATCTACGAGTGCATCGTGATTGGCGTGGTGACTGCCGGTGGCAAGATCGACGCGCCCATCGGCCGCCACGGCGGCATGCGCCAGCGCATGGc                                                                                                                                <  1:1190838/150‑1 (MQ=255)
                      aCCCTGCAGGCGCAGACCAAGCTGGTCGAGCAGATGCAGGCCCGCAAGGTCAGCCGCATCTACGAGTGCATCGTGATTGGCGTGGTGACTGCCGGTGGCAAGATCGACGCGCCCATCGGCCGCCACGGCGGCATGCGCCAGCGCATGGCt                                                                                                                               >  2:624948/1‑150 (MQ=255)
                         cTGCAGGCGCAGACCAAGCTGGTCGAGCAGATGCAGGCCCGCAAGGTCAGCCGCATCTACGAGTGCATCGTGATTGGCGTGGTGACTGCCGGTGGCAAGATCGACGCGCCCATCGGCCGCCACGGCGGCATGCGCCAGCGCATGGCt                                                                                                                               >  2:266291/1‑147 (MQ=255)
                         cTGCAGGCGCAGACCAAGCTGGTCGAGCAGATGCAGGCCCGCAAGGTCAGCCGCATCTACGAGTGCATCGTGATTGGCGTGGTGACTGCCGGTGGCAAGATCGACGCGCCCATCGGCCGCCACGGCGGCATGCGCCAGCGCATGGCt                                                                                                                               <  1:266291/147‑1 (MQ=255)
                           gCAGGCGCAGACCAAGCTGGTCGAGCAGATGCAGGCCCGCAAGGTCAGCCGCATCTACGAGTGCATCGTGATTGGCGTGGTGACTGCCGGTGGCAAGATCGACGCTCCCATCGGCCGCCACGGCTGCATGCTCCAGCGCATGGCTGTAAc                                                                                                                          <  2:184239/150‑1 (MQ=255)
                             agccgcAGACCAAGCTGTTCGAGCAGATGCAGGCCCGCAAGGTCAGCCGCATCTACGAGTGCATCGTGATTGGCGTGGTGACTGCCGGTGGCAAGATCGACGCGCCCATCGGCCGCCACGGCGGCATGCGCCAGCGCATGGCTGTAACCg                                                                                                                        >  2:866580/4‑150 (MQ=255)
                                   aGACCAAGCTGGTCGAGCAGATGCAGGCCCGCAAGGTCAGCCGCATCTACGAGTGCATCGTGATTGGCGTGGTGACTGCCGGTGGCAAGATCGACGCGCCCATCGGCCGCCACGGCGGCATGCGCCAGCGCATGGCTGTAACCGAcggcg                                                                                                                  <  1:793151/150‑1 (MQ=255)
                                       cAAGCTGGTCGAGCAGATGCAGGCCCGCAAGGTCAGCCGCATCTACGAGTGCATCGTGATTGGCGTGGTGACTGCCGGTGGCAAGATCGACGCGCCCATCGGCCGCCACGGCGGCAt                                                                                                                                               >  1:136613/1‑117 (MQ=255)
                                       cAAGCTGGTCGAGCAGATGCAGGCCCGCAAGGTCAGCCGCATCTACGAGTGCATCGTGATTGGCGTGGTGACTGCCGGTGGCAAGATCGACGCGCCCATCGGCCGCCACGGCGGCAt                                                                                                                                               <  2:136613/117‑1 (MQ=255)
                                       cAAGCTGGTCGAGCAGATGCAGGCCCGCAAGGTCAGCCGCATCTACGAGTGCATCGTGATTGGCGTGGTGACTGCCGGTGGCAAGATCGACGCGCCCATCGGCCGCCACGGCGGCATGCGCCa                                                                                                                                         >  2:672402/1‑123 (MQ=255)
                                       cAAGCTGGTCGAGCAGATGCAGGCCCGCAAGGTCAGCCGCATCTACGAGTGCATCGTGATTGGCGTGGTGACTGCCGGTGGCAAGATCGACGCGCCCATCGGCCGCCACGGCGGCATGCGCCa                                                                                                                                         <  1:672402/123‑1 (MQ=255)
                                       cAAGCTGGTCGAGCAGATGCAGGCCCGCAAGGTCAGCCGCATCTACGAGTGCATCGTGATTGGCGTGGTGACTGCCGGTGGCAAGATCGACGCGCCCATCGGCCGCCACGGCGGCATGCGCCAGCGCATGGCTGTAACCGACGGCGGCaa                                                                                                              <  2:890722/150‑1 (MQ=255)
                                          gCTGGTCGAGCAGATGCAGGCCCGCAAGGTCAGCCGCATCTACGAGTGCATCGTGATTGGCGTGGTGACTGCCGGTGGCAAGATCGACGCGCCCATCGGCCGCCACGGCGGCATGCGCCAGCGCATGGCTGTAACCGACGGCGGCAAGcc                                                                                                           >  2:64358/1‑150 (MQ=255)
                                                  aGCAGATGCAGGCCCGCAAGGTCAGCCGCATCTACGAGTGCATCGTGATTGGCGTGGTGACTGCCGGTGGCAAGATCGACGCGCCCATCGGCCGCCACGGCGGCATGCGCCAGCGCATGGCTGTAACCGACGGCGGCAAGCCGGCGgtca                                                                                                   <  1:788271/150‑1 (MQ=255)
                                                     aGATGCAGGCCCGCAAGGTCAGCCGCATCTACGAGTGCATCGCGATTGGCGTGGTGACTGCCGGCGGCAAGATCGACGCGCCCATCGGCCGCCACGGCGGCATGCGCCAGCGCATGGCTGTAACCGACGGCGGCAAGCCGGCGgtcagtc                                                                                                <  2:1055092/150‑1 (MQ=255)
                                                         gCAGGCCCGCAAGGTCAGCCGCATCTACGAGTGCATCGTGATTGGCGTGGTGACTGCCGGTGGCAAGATCGATGCGCCCATCGGCCGCCACGGCGGCATGCGCCAGCGCATGGCTGTAACCGACGGCGGCAAGCCGGCGGTCAGTCACTa                                                                                            <  2:700120/150‑1 (MQ=255)
                                                           aGGCCCGCAAGGTCAGCCGCATCTACGAGTGCATCGTGATTCGCGTGGAGACAAAAAGTGGCAGGAACGACAACCACATAGGACGCCACGGCGGCATGCGCCAGCGAATGGCTGTACCCTCCGGCGGCAAGCCGGCGGTTAGTCGCTAcc                                                                                          >  2:42379/1‑150 (MQ=255)
                                                             gCCCGCAAGGTCAGCCGCATCTACGAGTGCATCGTGATTGGCGTGGTGACTGCCGGTGGCAAGATCGACGCGCCCATCGGCCGCCACGGCGGCATGCGCCAgc                                                                                                                                       <  1:44011/103‑1 (MQ=255)
                                                             gCCCGCAAGGTCAGCCGCATCTACGAGTGCATCGTGATTGGCGTGGTGACTGCCGGTGGCAAGATCGACGCGCCCATCGGCCGCCACGGCGGCATGCGCCAgc                                                                                                                                       >  2:44011/1‑103 (MQ=255)
                                                             gCCCGCAAGGTCAGCCGCATCTACGAGTGCATCGTGATTGGCGTGGTGACTGCCGGTGGCAAGATCGACGCGCCCATCGGCCGCCACGGCGGCATGCGATAGCGCATGGCTGTAACCGACGGCGGCAAGCCGGCGGTCAGTCACTACcgc                                                                                        <  1:220524/150‑1 (MQ=255)
                                                                cGCAAGGTCAGCCGCATCTACGAGTGCATCGTGATTGGCGTGGTGACTGCCGGTGGCAAGATCGACGCGCCCATCGGCCGCCACGGCGGCATGCGCCAGCGCATGGCTGTAACCGACGGCGGCAAGCCGGCGGTCAGTCACTACCGCGTa                                                                                     >  2:758936/1‑150 (MQ=255)
                                                                  cAAGGTCAGCCGCATCTACGAGTGCATCGTGATTGGCGTGGTGACTGCCGGTGGCAAGATCGACGCGCCCATCGGCCGCCACGGCGGCATGCGCCAGCGCATGGCTGTAACCGACGGCGGCAAGCCGGCGGTCAGTCACTACCGCGTACt                                                                                   <  1:611366/150‑1 (MQ=255)
                                                                    aGGTCAGCCGCATCTACGAGTGCATCGTGATTGGCGTGGTGACTGCCGGTGGCAAGATCGACGCGCCCATCGGCCGCCACGGCGGCATGCGCCAGCGCATGGCTGTAACCGACGGCGGCAAGCCGGCGGTCAGTCACTACCGCGTACTCa                                                                                 <  1:649622/150‑1 (MQ=255)
                                                                           ccGCATCTACGAGTGCATCGTGATTGGCGTGGTGACTGCCGGTGGCAAGATCGACGCGCCCATCGGCCGCCACGGCGGCATGCGCCAGCGCATGGCTGTAACCGACGGCGGCAAGCCGGCGGTCAGTCACTACCGCGTACTCAAGCGCtt                                                                          >  1:1023714/1‑150 (MQ=255)
                                                                                  tACGAGTGCATCGTGATTGGCGTGGTGACTGCCGGTGGCAGGATCGACGCGCCCATCGGCCGCCACGGCGGCATGCGCCAGCGCATGGCTGTAACCGACGGCGGCAAGCCGGCGGTCAGTCACTACCGCGTACTCAAGCGCTTCCGCTCg                                                                   >  1:1154091/1‑150 (MQ=255)
                                                                                       gTGCATCGTGATTGGCGTGGTGACTGCCGGTGGCAAGATCGACGCGCCCATCGGCCGCCATGGCGGCATGCGCCAGCGCATGGCTGTAACCGACGGCGGCAAGCCGGCGGTCAGTCACTACCGCGTACTCAAGCGCTTCCGCTCGCacac                                                              >  1:1073667/1‑150 (MQ=255)
                                                                                       gTGCATCGTGATTGGCGTGGTGACTGCCGGTGGCAAGATCGACGCGCCCATCGGCCGCCATGGCGGCATGCGCCAGCGCATGGCTGTAACCGACGGCGGCAAGCCGGCGGTCAGTCACTACCGCGTACTCAAGCGCTTCCGCTCGCacac                                                              >  1:1072507/1‑150 (MQ=255)
                                                                                       gTGCATCGTGATTGGCGTGGTGACTGCCGGTGGCAAGATCGACGCGCCCATCGGCCGCCACGGCGGCATGCGCCAGCGCATGGCTGTAACCGACGGCGGCAAGCCGGCGGTCAGTCACTACCGCGTACTCAAGCGCTTCCGCTCGCacac                                                              <  1:758936/150‑1 (MQ=255)
                                                                                           aTCGTGATTGGCGTGGTGACTGCCGGTGGCAAGATCGACGCGCCCATCGGCCGCCACGGCGGCATGCGCCAGCGCATgg                                                                                                                                 >  1:419454/1‑79 (MQ=255)
                                                                                           aTCGTGATTGGCGTGGTGACTGCCGGTGGCAAGATCGACGCGCCCATCGGCCGCCACGGCGGCATGCGCCAGCGCATgg                                                                                                                                 <  2:419454/79‑1 (MQ=255)
                                                                                                 aTTGGCGTGGTGACTGCCGGTGGCAAGATCGACGCGCCCATCGGCCGCCACGGCGGCATGCGCCAGCGCATGGCTGTAACCGACGGCGGCAAGCCGGCGGTCAGTCACTACCGCGTACTCAAGCGCTTCCGCTCGCACACCCATGTGCgg                                                    >  1:804898/1‑150 (MQ=255)
                                                                                                           tgACTGCCGGTGGCAAGATCGACGCGCCCATCGGCCGCCACGGCGGCATGCGCCAGCGCATGGCTGTAACCGACGGCGGCAAGCCGGCGGTCAGTCACTACCGCGTACTCAAGCGCTTCCGCTCGCACACCCATGTGCGGGTCAAGCTgg                                          >  1:944205/1‑150 (MQ=255)
                                                                                                             aCTGCCGGTGGCAAGATCGACGCGCCCATCGGCCGCCACGGCGGCATGCGCCAGCGCATGGCTGTAACCGACGGCGGCAAGCCGGCGGTCAGTCACTACCGCGTa                                                                                     <  1:420135/105‑1 (MQ=255)
                                                                                                             aCTGCCGGTGGCAAGATCGACGCGCCCATCGGCCGCCACGGCGGCATGCGCCAGCGCATGGCTGTAACCGACGGCGGCAAGCCGGCGGTCAGTCACTACCGCGTa                                                                                     >  2:420135/1‑105 (MQ=255)
                                                                                                              cTGCCGGTGGCAAGATCGACGCGCCCATCGGCCGCCACGGCGGCATGCGCCAGCGCATGGCTGTAACCGACGGCGGCAAGCCGGCGGTCAGTCACTACCGCGTACTCAAGCGCTTCCGCt                                                                     >  1:1049489/1‑120 (MQ=255)
                                                                                                              cTGCCGGTGGCAAGATCGACGCGCCCATCGGCCGCCACGGCGGCATGCGCCAGCGCATGGCTGTAACCGACGGCGGCAAGCCGGCGGTCAGTCACTACCGCGTACTCAAGCGCTTCCGCt                                                                     <  2:1049489/120‑1 (MQ=255)
                                                                                                               tGCCGGTGGCAAGATCGACGCGCCCATCGGCCGCCACGGCGGCATGCGCCAGCGCATGGCTGTAACCGACGGCGGCAAGCCGGCGGTCAGTCACTACCGCGTACTCAAGCGCTTCCGCTCGCACACCCATGTGCGGGt                                                  >  2:1196483/1‑138 (MQ=255)
                                                                                                               tGCCGGTGGCAAGATCGACGCGCCCATCGGCCGCCACGGCGGCATGCGCCAGCGCATGGCTGTAACCGACGGCGGCAAGCCGGCGGTCAGTCACTACCGCGTACTCAAGCGCTTCCGCTCGCACACCCATGTGCGGGt                                                  <  1:1196483/138‑1 (MQ=255)
                                                                                                                 ccGGTGGCAAGATCGACGCGCCCATCGGCCGCCATGGCGGCATGCGCCAGCGCATGGCTGTAACCGACGGCGGCAAGCCGGCGGTCAGTCACTACCGCGTACTCAAGCGCTTCCGCTCGCACACCCATGTGCGGGTCAAGCTGGAAACCg                                    <  2:1073667/150‑1 (MQ=255)
                                                                                                                 ccGGTGGCAAGATCGACGCGCCCATCGGCCGCCATGGCGGCATGCGCCAGCGCATGGCTGTAACCGACGGCGGCAAGCCGGCGGTCAGTCACTACCGCGTACTCAAGCGCTTCCGCTCGCACACCCATGTGCGGGTCAAGCTGGAAACCg                                    <  2:1072507/150‑1 (MQ=255)
                                                                                                                     tGGCAAGATCGACGCGCCCATCGGCCGCCACGGCGGCATGCGCCAGCGCATGGCTGTAACCGACGGCGGCAAGCCGGCGGTCAGTCACTACCGCGTACTCAAGCGCTTCCGCTCGCACACCCATGTGCGGGTCAAGCTGGAAACCGGCCg                                >  1:400516/1‑150 (MQ=255)
                                                                                                                       gCAAGATCGACGCGCCCATCGGCCGCCACGGCGGCATGCGCCAGCGCATGGCTGTAACCGACGGCGGCAAGCcggcgg                                                                                                      <  1:994550/78‑1 (MQ=255)
                                                                                                                       gCAAGATCGACGCGCCCATCGGCCGCCACGGCGGCATGCGCCAGCGCATGGCTGTAACCGACGGCGGCAAGCcggcgg                                                                                                      >  2:994550/1‑78 (MQ=255)
                                                                                                                        cAAGATCGACGCGCCCATCGGCCGCCACGGCGGCATGCGCCAGCGCATGGCTGTAACCGACGGCGGCAAGCCGGCGGTCAGTCACTACCGCGTACTCAAGCGCTTCCGCTCGCACACCCATGTGCGGGTCAAGCTGGAAACCGGCCGTAc                             >  1:14597/1‑150 (MQ=255)
                                                                                                                          aGATCGACGCGCCCATCGGCGGCCACGGCGGCATGCGCCAGCGCATGGCTGTAACCGACGGCGGCAAGCCGGCGGTCAGTCACTACCGCGTGCTCAAGCGCTTCCGCGCGCACACCCATGTGCGGGTCGAGCTGGAAACCGGCCGTAccc                           >  1:198573/1‑150 (MQ=255)
                                                                                                                                aCGCGCCCATCGGCCGCCACGGCGGCATGCGCCAGCGCATGGCTGTAACCGACGGCGGCAAGCCGGCGGTCAGTCACTACCGCGTACTCAAGCGCTTCCGCTCGCACACCCATGTGCGGGTCAAGCTGGAAACCGGCCGTACCCACCAGa                     <  2:1023714/150‑1 (MQ=255)
                                                                                                                                  gcacccATCGGCCGCCACGGCGGCATGCGCCAGCGCATGGCTGTAACCGACGGCGGCAAGCCGGCGGTCAGTCACTACCGCGTACTCAAGCGCTTCCGCTCGCACACCCATGTGCGGGTCAAGCTGGAAACCGGCCGTACCCACCAGAtt                   >  2:1141811/4‑150 (MQ=255)
                                                                                                                                           ggCCGCCACGGCGGCATGCGCCAGCGCATGGCTGTAACCGACGGCGGCAAGCCGGCGGTCAGTCACTACCGCGTACTCAAGCGCTTCCGCTCGCACACCCATGTGCGGGTCAAGCTGGAAACCGGCCGTACCCACCAGATTCGCGTGcac          >  2:768592/1‑150 (MQ=255)
                                                                                                                                                ccACGGCGGCATGCGCCAGCGCATGGCTGTAACCGACGGCGGCAAGCCGGCGGTCAGTCACTACCGCGTACTCAAGCGCTTCCGCTCGCACACCCATGTGCGGGTCAAGCTGGAAACCGGCCGTACCCACCAGATTcgc                <  2:708946/139‑1 (MQ=255)
                                                                                                                                                ccACGGCGGCATGCGCCAGCGCATGGCTGTAACCGACGGCGGCAAGCCGGCGGTCAGTCACTACCGCGTACTCAAGCGCTTCCGCTCGCACACCCATGTGCGGGTCAAGCTGGAAACCGGCCGTACCCACCAGATTcgc                >  1:708946/1‑139 (MQ=255)
                                                                                                                                                   cggcggCATGCGCCAGCGCATGGCTGTAACCGACGGCGGCAAGCCGGCGGTCAGTCACTACCGCGTACTCAAACGCTCCCGCTCGCACACCCATGTGCGGGTCAAGCTGGAAACCGGCCGTACCCACCAGATTCGCGTGCACATGGCCCa  >  2:3935/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
CGGCCTGATGGTGGTGGCCAAGACCCTGCAGGCGCAGACCAAGCTGGTCGAGCAGATGCAGGCCCGCAAGGTCAGCCGCATCTACGAGTGCATCGTGATTGGCGTGGTGACTGCCGGTGGCAAGATCGACGCGCCCATCGGCCGCCACGGCGGCATGCGCCAGCGCATGGCTGTAACCGACGGCGGCAAGCCGGCGGTCAGTCACTACCGCGTACTCAAGCGCTTCCGCTCGCACACCCATGTGCGGGTCAAGCTGGAAACCGGCCGTACCCACCAGATTCGCGTGCACATGGCCCA  >  NC_002947/728940‑729236

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