Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,143,125 C→T D415N (GAC→AAC)  rne ← endoribonuclease; RNA‑binding protein;RNA degradosome binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,143,1250CT98.2% 190.4 / ‑5.0 56D415N (GAC→AAC) rneendoribonuclease; RNA‑binding protein;RNA degradosome binding protein
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (26/29);  total (27/29)
Fisher's exact test for biased strand distribution p-value = 4.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

GGCACAGGAACAATGGCGTGAACTTCCTGGGTGTTCTCTTTCAGCGCTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTCACGCACGGTGCCAGTACCAGAACAACGCGGACAAACGTGATGACTGGATTCACCCAGTGATGGGCTCAGGCGCTGACGGGACATTTCCAGCAGG  >  NC_000913/1143031‑1143219
                                                                                              |                                                                                              
ggCACAGGAACAATGGCGTGAACTTCCTGGGTGTTCTCTTTCAGCGCTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTtacgca                                                                                           >  2:2013131/1‑100 (MQ=255)
ggCACAGGAACAATGGCGTGAACTTCCTGGGTGTTCTCTTTCAGCGCTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTtacgca                                                                                           >  2:1180025/1‑100 (MQ=255)
    caGGAACAATGGCGTGAACTTCCTGGGTGTTCTCTTTCAGCGCTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTTACGCACGGt                                                                                       >  2:231668/1‑100 (MQ=255)
     aGGAACAATGGCGTGAACTTCCTGGGTGTTCTCTTTCAGCGCTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTTACGCACGGTg                                                                                      >  2:68998/1‑100 (MQ=255)
      ggAACAATGGCGTGAACTTCCTGGGTGTTCTCTTTCAGCGCTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTTACGCACGGTGc                                                                                     >  1:140479/1‑100 (MQ=255)
       gAACAATGGCGTGAACTTCCTGGGTGTTCTCTTTCAGCGCTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTTACGCACGGTGcc                                                                                    <  1:2004802/100‑1 (MQ=255)
       gAACAATGGCGTGAACTTCCTGGGTGTTCTCTTTCAGCGCTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTTACGCACGGTGcc                                                                                    <  1:717256/100‑1 (MQ=255)
        tcttccgatctGAACTTCCTGGGTGTTCTCTTTCAGCGCTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTTACGCACGGTGCCa                                                                                   <  1:1016048/90‑1 (MQ=255)
         aCAATGGCGTGAACTTCCTGGGTGTTCTCTTTCAGCGCTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTTACGCACGGTGCCAg                                                                                  <  1:429425/100‑1 (MQ=255)
          cAATGGCGTGAACTTCCTGGGTGTTCTCTTTCAGCGCTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTTACGCACGGTGCCAGt                                                                                 <  1:231376/100‑1 (MQ=255)
            aTGGCGTGAACTTCCTGGGTGTTCTCTTTCAGCGCTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTTACGCACGGTGCCAGTAc                                                                               <  1:1047006/100‑1 (MQ=255)
            aTGGCGTGAACTTCCCGGGTGTTCTCTTTCAGCGCTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTTACGCACGGTGCCAGTAc                                                                               <  1:2013149/100‑1 (MQ=255)
             tGGCGTGAACTTCCTGGGTGTTCTCTTTCAGCGCTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTTACGCACGGTGCCAGTAcc                                                                              <  2:140478/100‑1 (MQ=255)
                cGTGAACTTCCTGGGTGTTCTCTTTCAGCGCTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTTACGCACGGTGCCAGTACCAGa                                                                           <  1:128382/100‑1 (MQ=255)
                cGTGAACTTCCTGGGTGTTCTCTTTCAGCGCTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTTACGCACGGTGCCAGTACCAGa                                                                           >  2:810069/1‑100 (MQ=255)
                 gTGAACTTCCTGGGTGTTCTCTTTCAGCGCTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTTACGCACGGTGCCAGTACCAGaa                                                                          <  2:676264/100‑1 (MQ=255)
                 gTGAACTTCCTGGGTGTTCTCTTTCAGCGCTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTTACGCACGGTGCCAGTACCAGaa                                                                          >  2:1149662/1‑100 (MQ=255)
                 gTGAACTTCCTGGGTGTTCTCTTTCAGCGCTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTTACGCACGGTGCCAGTACCAGaa                                                                          >  1:1361805/1‑100 (MQ=255)
                  tGAACTTCCTGGGTGTTCTCTTTCAGCGCTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTTACGCACGGTGCCAGTACCAGaac                                                                         <  1:1308870/100‑1 (MQ=255)
                     aCTTCCTGGGTGTTCTCTTTCAGCGCTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTTACGCACGGTGCCAGTACCAGaacaac                                                                      >  2:1560189/1‑100 (MQ=255)
                       ttCCTGGGTGTTCTCTTTCAGCGCTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTTACGCACGGTGCCAGTACCAGAACAAcgc                                                                    >  2:1469134/1‑100 (MQ=255)
                                  tctcTTTCAGCGCTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTTACGCACGGTGCCAGTACCAGAACAACGCGGACAAACGtg                                                         >  2:218286/1‑100 (MQ=255)
                                  tctcTTTCAGCGCTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTTACGCACGGTGCCAGTACCAGAACAACGCGGACAAACGtg                                                         >  2:1975554/1‑100 (MQ=255)
                                  tctcTTTCAGCGCTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTTACGCACGGTGCCAGTACCAGAACAACGCGGACAAACGtg                                                         >  2:1045260/1‑100 (MQ=255)
                                     ctTTCAGCGCTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTCACGCACGGTGCCAGTACCAGAACAACGCGGACAAACGtgatg                                                      >  2:861307/1‑100 (MQ=255)
                                      tttCAGCGCTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTTACGCACGGTGCCAGTACCAGAACAACGCGGACAAACGtgatga                                                     >  2:326760/1‑100 (MQ=255)
                                       ttCAGCGCTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTTACGCACGGTGCCAGTACCAGAACAACGCGGACAAACGTGATGAc                                                    <  1:164973/100‑1 (MQ=255)
                                        tCAGCGCTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTTACGCACGGTGCCAGTACCAGAACAACGCGGACAAACGTGATGACt                                                   >  1:1881490/1‑100 (MQ=255)
                                         cAGCGCTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTTCCGCACGGTGCCAGTACCAGAACAACGCAGACAAACGTGATGACTg                                                  <  1:1180031/100‑1 (MQ=255)
                                         cAGCGCTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTTACGCACGGTGCCAGTACCAGAACAACGCGGACAAACGTGATGACTg                                                  <  1:1149668/100‑1 (MQ=255)
                                          aGCGCTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTTACGCACGGTGCCAGTACCAGAACAACGCGGACAAACGTGATGACTgg                                                 <  2:904555/100‑1 (MQ=255)
                                          aGCGCTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTTACGCACGGTGCCAGTACCAGAACAACGCGGACAAACGTGATGACTgg                                                 >  1:862580/1‑100 (MQ=255)
                                              cttcttcttCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTTACGCACGGTGCCAGTACCAGAACAACGCGGACAAACGTGATGACTGGATTc                                             >  2:1727527/1‑100 (MQ=255)
                                                tcttcttcGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTTACGCACGGTGCCAGTACCAGAACAACGCGGACAAACGTGATGACTGGATTCAc                                           <  2:643774/100‑1 (MQ=255)
                                                tcttcttcGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTTACGCACGGTGCCAGTACCAGAACAACGCGGACAAACGTGATGACTGGATTCAc                                           <  1:1324253/100‑1 (MQ=255)
                                                 cttcttcGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTTACGCACGGTGCCAGTACCAGAACAACGCGGACAAACGTGCTGACTGCATTCAcc                                          >  2:1385695/1‑100 (MQ=255)
                                                  ttcttcGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTTACGCACGGTGCCAGTACCAGAACAACGCGGACAAACGTGATGACTGGATTCAccc                                         <  2:1879947/100‑1 (MQ=255)
                                                  ttcttcGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTTACGCACGGTGCCAGTACCAGAACAACGCGGACAAACGTGATGACTGGATTCAccc                                         <  1:1681683/100‑1 (MQ=255)
                                                   tcttcGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTTACGCACGGTGCCAGTACCAGAACAACGCGGACAAACGTGATGACTGGATTCACCCa                                        >  1:1319835/1‑100 (MQ=255)
                                                      tcGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTTACGCACGGTGCCAGTACCAGAACAACGCGGACAAACGTGATGACTGGATTCACCCAGTg                                     >  2:1865885/1‑100 (MQ=255)
                                                         aTCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTTACGCACGGTGCCAGTACCAGAACAACGCGGACAAACGTGATGACTGGATTCACCCAGTGATg                                  <  2:733590/100‑1 (MQ=255)
                                                          tCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTTACGCACGGTGCCAGTACCAGAACAACGCGGACAAACGTGATGACTGGATTCACCCAGTGATgg                                 <  1:816814/100‑1 (MQ=255)
                                                               cGCAGAATAGAGAGCGACAGCGATTCGTTGTTACGCACGGTGCCAGTACCAGAACAACGCGGACAAACGTGATGACTGGATTCACCCAGTGATGGGCTCa                            <  2:229256/100‑1 (MQ=255)
                                                                 cAGAATAGAGAGCGACAGCGATTCGTTGTTACGCACGGTGCCAGTACCAGAACAACGCGGACAAACGTGATGACTGGATTCACCCAGTGATGGGCTCAgg                          >  2:502629/1‑100 (MQ=255)
                                                                   gAATAGAGAGCGACAGCGATTCGTTGTTACGCACGGTGCCAGTACCAGAACAACGCGGACAAACGTGATGACTGGATTCACCCAGTGATGGGCTCAGgcg                        <  2:1149261/100‑1 (MQ=255)
                                                                        gagagCGACAGCGATTCGTTGTTACGCACGGTGCCAGTACCAGAACAACGCGGACAAACGTGATGACTGGATTCACCCAGTGATGGGCTCAGGCGCTGa                    >  1:89150/1‑99 (MQ=255)
                                                                        gagagCGACAGCGATTCGTTGTTACGCACGGTGCCAGTACCAGAACAACGCGGACAAACGTGATGACTGGATTCACCCAGTGATGGGCTCAGGCGCTGAc                   >  2:99397/1‑100 (MQ=255)
                                                                        gagagCGACAGCGATTCGTTGTTACGCACGGTGCCAGTACCAGAACAACGCGGACAAACGTGATGACTGGATTCACCCAGTGATGGGCTCAGGCGCTGAc                   >  1:699002/1‑100 (MQ=255)
                                                                          gagCGACAGCGATTCGTTGTTACGCACGGTGCCAGTACCAGAACAACGCGGACAAACGTGATGACTGGATTCACCCAGTGATGGGCTCAGGCGCTGACgg                 >  2:1351719/1‑100 (MQ=255)
                                                                            gCGACAGCGATTCGTTGTTACGCACGGTGCCAGTACCAGAACAACGCGGACAAACGTGATGACTGGATTCACCCAGTGATGGGCTCAGGCGCTGACGGGa               <  2:1306681/100‑1 (MQ=255)
                                                                             ctaCAGCGATGCGTTGTTACGTACGGTGCCAGTACCAGAACAACGCGGACAAACGTGATGACTGGATTCACCCAGTGATGGGCTCAGGCGCTGACGGGAc              >  2:1322134/3‑100 (MQ=255)
                                                                             cGACAGCGATTCGTTGTTACGCACGGTGCCAGTACCAGAACAACGCGGACAAACGTGATGACTGGATTCACCCAGTGATGGGCTCAGGCGCTGACggggc              >  1:1058499/1‑98 (MQ=255)
                                                                             cGACAGCGATTCGTTGTTACGCACGGTGCCAGTACCAGAACAACGCGGACAAACGTGATGACTGGATTCACCCAGTGATGGGCTCAGGCGCTGACGGGAc              <  2:512787/100‑1 (MQ=255)
                                                                             cGACAGCGATTCGTTGTTACGCACGGTGCCAGTACCAGAACAACGCGGACAAACGTGATGACTGGATTCACCCAGTGATGGGCTCAGGCGCTGACGGGAc              <  1:670157/100‑1 (MQ=255)
                                                                                   cGATTCGTTGTTACGCACGGTGCCAGTACCAGAACAACGCGGACAAACGTGATGACTGGATTCACCCAGTGATGGGCTCAGGCGCTGACGGGACATTTcc        <  1:810072/100‑1 (MQ=255)
                                                                                     aTTCGTTGTTACGCACGGTGCCAGTACCAGAACAACGCGGACAAACGTGATGACTGGATTCACCCAGTGATGGGCTCAGGCGCTGACGGGACATTTCcag      <  2:466864/100‑1 (MQ=255)
                                                                                     aTTCGTTGTTACGCACGGTGCCAGTACCAGAACAACGCGGACAAACGTGATGACTGGATTCACCCAGTGATGGGCTCAGGCGCTGACGGGACATTTCcag      >  1:266031/1‑100 (MQ=255)
                                                                                      ttCGTTGTTACGCACGGTGCCAGTACCAGAACAACGCGGACAAACGTGATGACTGGATTCACCCAGTGATGGGCTCAGGCGCTGACGGGACATTTCcagc     <  1:326761/100‑1 (MQ=255)
                                                                                         gttgttACGCACGGTGCCAGTACCAGAACAACGCGGACAAACGTGATGACTGGATTCACCCAGTGATGGGCTCAGGCGCTGACGGGACATTTCCAGCAgg  >  2:1593205/1‑100 (MQ=255)
                                                                                              |                                                                                              
GGCACAGGAACAATGGCGTGAACTTCCTGGGTGTTCTCTTTCAGCGCTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTCACGCACGGTGCCAGTACCAGAACAACGCGGACAAACGTGATGACTGGATTCACCCAGTGATGGGCTCAGGCGCTGACGGGACATTTCCAGCAGG  >  NC_000913/1143031‑1143219

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