Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,478,879 T→C V136A (GTT→GCT)  rraB → protein inhibitor of RNase E

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,478,8790TC91.9% 231.6 / 5.4 74V136A (GTT→GCT) rraBprotein inhibitor of RNase E
Reads supporting (aligned to +/- strand):  ref base T (3/3);  new base C (34/34);  total (37/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAA  >  NC_000913/4478785‑4478976
                                                                                              |                                                                                                 
cttccgatctACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGcc                                                                                              <  1:1812017/91‑1 (MQ=255)
    ttGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTa                                                                                          >  2:1632640/1‑100 (MQ=255)
       aGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAAtt                                                                                       <  2:1503736/100‑1 (MQ=255)
         tACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTaa                                                                                     >  1:775799/1‑100 (MQ=255)
          acgacgGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCagatcgaaga                                                                                    >  2:1812028/1‑91 (MQ=255)
          acgacgGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAg                                                                                    >  2:1630423/1‑100 (MQ=255)
           cgacgGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGc                                                                                   >  1:841686/1‑100 (MQ=255)
            gacgGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCg                                                                                  >  1:1428627/1‑100 (MQ=255)
            gacgGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCg                                                                                  <  1:1604303/100‑1 (MQ=255)
               ggATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAg                                                                               >  1:1552392/1‑100 (MQ=255)
                 aTGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAgcg                                                                             >  2:1239964/1‑100 (MQ=255)
                 aTGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAgcg                                                                             <  1:1680741/100‑1 (MQ=255)
                  tGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAgcgc                                                                            >  1:1915004/1‑100 (MQ=255)
                   ggggCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCa                                                                           <  2:945600/100‑1 (MQ=255)
                     ggCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATc                                                                         >  2:864434/1‑100 (MQ=255)
                     ggCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATc                                                                         <  1:1937035/100‑1 (MQ=255)
                       cACTTACTTTGAAGATCCTAACGGCGAAGATCGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAg                                                                       <  1:499889/100‑1 (MQ=255)
                        acttacttTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAgg                                                                      <  1:1790922/100‑1 (MQ=255)
                         cttacttTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGc                                                                     >  1:422022/1‑100 (MQ=255)
                         cttacttTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGc                                                                     >  2:429233/1‑100 (MQ=255)
                         cttacttTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGAAGCTCGCCACTAATTAAGCGTAGCGCATCAGGc                                                                     >  2:647416/1‑100 (MQ=255)
                           tacttTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGccc                                                                   >  2:533058/1‑100 (MQ=255)
                           tacttTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGccc                                                                   >  2:1656627/1‑100 (MQ=255)
                            acttTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCa                                                                  >  1:1091666/1‑100 (MQ=255)
                             cttTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAg                                                                 <  2:376568/100‑1 (MQ=255)
                             cttTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAg                                                                 <  2:557328/100‑1 (MQ=255)
                             cttTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAg                                                                 <  1:1181845/100‑1 (MQ=255)
                              tttGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAgg                                                                >  2:1398681/1‑100 (MQ=255)
                               ttGAAGATCCTAACGGCGAAGATGGCGAGGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGt                                                               <  1:429772/100‑1 (MQ=255)
                                   aGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGtt                                                           >  2:1942711/1‑100 (MQ=255)
                                   aGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGtt                                                           >  1:74885/1‑100 (MQ=255)
                                     aTCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTcc                                                         <  1:947805/100‑1 (MQ=255)
                                      tCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCt                                                        >  1:1414656/1‑100 (MQ=255)
                                       ccTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTg                                                       <  1:361004/100‑1 (MQ=255)
                                        cTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGc                                                      >  1:47183/1‑100 (MQ=255)
                                          aaCGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCaa                                                    >  1:270366/1‑100 (MQ=255)
                                            cGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAAcc                                                  <  2:1414664/100‑1 (MQ=255)
                                             ggCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCt                                                 <  1:1013741/100‑1 (MQ=255)
                                              gCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTg                                                >  2:1929492/1‑100 (MQ=255)
                                               cGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTgg                                               >  2:312092/1‑100 (MQ=255)
                                               cGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTgg                                               <  2:989609/100‑1 (MQ=255)
                                                   atctgCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGaa                                           <  2:1041302/96‑1 (MQ=255)
                                                    aTGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAAt                                          >  1:1789736/1‑100 (MQ=255)
                                                     tGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAAtt                                         <  1:652535/100‑1 (MQ=255)
                                                       gCGACGATGAAGATTTTGTCGATGAAGACGATGACGGCGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGaaagat                                       >  1:1041301/1‑96 (MQ=255)
                                                       gCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTg                                       >  2:1639117/1‑100 (MQ=255)
                                                        cgacgaTGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGt                                      >  1:1120093/1‑100 (MQ=255)
                                                         gacgaTGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTa                                     <  1:312090/100‑1 (MQ=255)
                                                            gaTGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTTATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGcc                                  <  1:1239960/100‑1 (MQ=255)
                                                             aTGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGccc                                 >  1:369466/1‑100 (MQ=255)
                                                             aTGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGccc                                 <  2:469948/100‑1 (MQ=255)
                                                             aTGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGccc                                 >  2:405859/1‑100 (MQ=255)
                                                              tGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCg                                <  2:422024/100‑1 (MQ=255)
                                                               gAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCgg                               >  1:1723981/1‑100 (MQ=255)
                                                                aaGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGa                              <  1:1632630/100‑1 (MQ=255)
                                                                 aGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGAt                             <  2:956675/100‑1 (MQ=255)
                                                                  gATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATa                            <  2:878089/100‑1 (MQ=255)
                                                                    ttttGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCTCAGGTAGTTCATGCAACCTGGTGAATTTGTAGCCCGGATatg                          <  2:1313377/100‑3 (MQ=255)
                                                                     tttGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAgg                         <  1:429231/100‑1 (MQ=255)
                                                                          cGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTc                    >  2:1602915/1‑100 (MQ=255)
                                                                           gATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCa                   >  1:1311167/1‑100 (MQ=255)
                                                                           gATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCa                   >  2:1853074/1‑100 (MQ=255)
                                                                             tGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATg                 <  2:775800/100‑1 (MQ=255)
                                                                               aaGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGTCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTc               <  1:1650986/100‑1 (MQ=255)
                                                                               aaGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTc               >  2:447958/1‑100 (MQ=255)
                                                                               aaGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTc               <  1:1219925/100‑1 (MQ=255)
                                                                                aGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTca              <  2:74885/100‑1 (MQ=255)
                                                                                aGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTca              >  2:22562/1‑100 (MQ=255)
                                                                                   cGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACAt           <  1:237035/100‑1 (MQ=255)
                                                                                     aTGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCATATcc         <  1:447956/100‑1 (MQ=255)
                                                                                      tGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCg        >  2:230929/1‑100 (MQ=255)
                                                                                         cGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCa     <  2:1264439/100‑1 (MQ=255)
                                                                                          ggAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCAt    <  1:1639107/100‑1 (MQ=255)
                                                                                            aGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATaa  <  2:884825/100‑1 (MQ=255)
                                                                                              |                                                                                                 
GACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAA  >  NC_000913/4478785‑4478976

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