Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 78,003 A→C E128A (GAG→GCG)  yabM → putative transport protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091378,0030AC99.2% 356.0 / NA 131E128A (GAG→GCG) yabMputative transport protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base C (67/63);  minor base G (0/1);  total (67/64)
Fisher's exact test for biased strand distribution p-value = 4.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 134 total aligned reads displayed.

TGACGCTTATCACCTGTGGTGTGCTTCTGGCATCTCTGGCCAATACGGCAATGCCACAGTTATTTGCTCTGGAGCGGGAATATGCGGATAACTCGGCGCGAGAAGTGGTGATGTTTAGCTCGGTGATGCGTGCGCAGCTTTCT  >  NC_000913/77931‑78073
                                                                        |                                                                      
tGACGCTTATCACCTGTGGTGTGCTTCTGGCCTCTCTGGCCACTCCGGCACTGCCACCGATCTTTGCTCCGGcgcg                                                                     >  1:4871366/1‑76 (MQ=255)
tGACGCTTATCACCTGTGGTGTGCTTCTGGCATCTCTGGCCAATACGGCAATGCCACAGTTATTTGCTCTGGcgcg                                                                     >  1:750824/1‑76 (MQ=255)
  aCGCTTATCACCTGTGGTGTGCTTCTGGCATCTCTGGCCAATACGGCAATGCCACAGTTATTTGCTCTGGCGCggg                                                                   >  1:1184920/1‑76 (MQ=255)
   cGCTTATCACCTGTGGTGTGCTTCTGGCATCTCTGGCCAATACGGCAATGCCACAGTTATTTGCTCTGGCGCGGGa                                                                  <  2:3202602/76‑1 (MQ=255)
   cGCTTATCACCTGTGGTGTGCTTCTGGCATCTCTGGCCAATACGGCAAGGCCACAGTTATTTGCCCTcgcgcgcca                                                                  >  2:57035/1‑73 (MQ=255)
    gCTTATCACCTGTGGTGTGCTTCTGGCATCTCTGGCCAATACGGCAATGCCACAGTTATTTGCTCTGGCGCGGGaa                                                                 >  1:1956780/1‑76 (MQ=255)
    gCTTATCACCTGTGGTGTGCTTCTGGCATCTCTGGCCAATACGGCAATGCCACAGTTATTTGCTCTGGCGCGGGaa                                                                 >  1:1607870/1‑76 (MQ=255)
     cTTATCACCTGTGGTGTGCTTCTGGCATCTCTGGCCAATACGGCAATGCCACAGTTATTTGCTCTGGCGCGGGAat                                                                >  2:4043401/1‑76 (MQ=255)
     cTTATCACCTGTGGTGTGCTTCTGGCATCTCTGGCCAATACGGCAATGCCACAGTTATTTGCTCTGGCGCGGGAat                                                                <  2:4158359/76‑1 (MQ=255)
      ttATCACCTGTTGTGTGCTGCTGGCATCTCTGGCCAATACGGCAATGCCACAGTTATTTGCTCTGGCGCGGGAata                                                               <  1:4083863/76‑1 (MQ=255)
      ttATCACCTGTGGTGTGCTTCTGGCATCTCTGGCCAATACGGCAATGCCACAGTTATTTGCTCTGGCGCGGGActa                                                               >  2:4157689/1‑76 (MQ=255)
      ttATCACCTGTGGTGTGCTTCTGGCATCTCTGGCCAATACGGCAATGCCACAGTTATTTGCTCTGGCGCGGGAata                                                               >  1:4262622/1‑76 (MQ=255)
      ttATCACCTGTGGTGTGCTTCTGGCATCTCTGGCCAATACGGCAATGCCACAGTTATTTGCTCTGGCGCGGGAata                                                               <  1:269672/76‑1 (MQ=255)
        aTCACCTGTGGTGTGCTTCTGGCATCTCTGGCCAATACGGCAATGCCACAGTTATTTGCTCTGGCGCGGGAATATg                                                             >  2:3890639/1‑76 (MQ=255)
         tCACCTGTGGTGTGCTTCTGGCATCTCTGGCCAATACGGCAATGCCACAGTTATTTGCTCTGGCGCGGGAATATGc                                                            <  2:4741228/76‑1 (MQ=255)
          cACCTGTGGTGTGCTTCTGGCATCTCTGGCCAATACGGCAATGCCACAGTTATTTGCTCTGGCGCGGGAATATGCg                                                           >  2:4467767/1‑76 (MQ=255)
          cACCTGTGGTGTGCTTCTGGCATCTCTGGCCAATACGGCAATGCCACAGTTATTTGCTCTGGCGCGGGAATATGCg                                                           <  1:918418/76‑1 (MQ=255)
           aCCTGTGGTGTGCTTCTGGCATCTCTGGCCAATACGGCAATGCCACAGTTATTTGCTCTGGCGCGGGAATATGCgg                                                          >  2:2742868/1‑76 (MQ=255)
           aCCTGTGGTGTGCTTCTGGCATCTCTGGCCAATACGGCAATGCCACAGTTATTTGCTCTGGCGCGGGAATATGCgg                                                          <  1:1011354/76‑1 (MQ=255)
            ccTGTGGTGTGCTTCTGGCATCTCTGGCCAATACGGCAATGCCACAGTTATTTGCTCTGGCGCGGGAATATGCGGa                                                         >  1:3971713/1‑76 (MQ=255)
              tgtgGTGTGCTTCTGGCATCTCTGGCCAATACGGCAATGCCACAGTTATTTGCTCTGGCGCGGGAATATGCGGATa                                                       <  2:3438569/76‑1 (MQ=255)
               gtggtgTGTTGCGGGCATCTCTGGCCAATACGGCAATGCCACAGTTATTTGCTCTGGCGCGGGAATATGCGGATaa                                                      <  2:9613/76‑1 (MQ=255)
               gtggtgTGCTTCTGGCATCTCTGGCCAATACGGCAATGCCACAGTTATTTGCTCTGGCGCGGGAATATGCGGATaa                                                      <  1:2544115/76‑1 (MQ=255)
                  gtgtGCTTCTGGCATCTCTGGCCAATACGGCAATGCCACAGTTATTTGCTCTGGCGCGGGAATATGCGGATAACTc                                                   >  2:1260024/1‑76 (MQ=255)
                  gtgtGCTTCTGGCATCTCTGGCCAATACGGCAATGCCACAGTTATTTGCTCTGGCGCGGGAATATGCGGATAACTc                                                   <  1:1073583/76‑1 (MQ=255)
                  gtgtGCTTCTGGCATCTCTGGCCAATACGGCAATGCCACAGTTATTTGCTCTGGCGCGGGAATATGCGGATAACTc                                                   <  1:4040387/76‑1 (MQ=255)
                  gtgtGCTTCTGGCATCTCTGGCCAATACGGCAATGCCACAGTTATTTGCTCTGGCGCGGGAATATGCGGATAACTc                                                   <  1:3306950/76‑1 (MQ=255)
                  gtgtGCGGCTGGCATCTCTGGCCAAGACGGCAATGCCACAGTTAGTTGCTCTGGCGCGGGAATATGCGGATAAcgc                                                   <  1:3343764/76‑3 (MQ=255)
                   tgtgCTTCTGGCATCTCTGGCCAATACGGCAATGCCACAGTTATTTGCTCTGGCGCGGGAATATGCGGATAACTCg                                                  >  1:2935474/1‑76 (MQ=255)
                   tgtgCTTCTGGCATCTCTGGCCAATACGGCAATGCCACAGTTATTTGCTCTGGCGCGGGAATATGCGGATAACTCg                                                  <  2:2597403/76‑1 (MQ=255)
                   ggggCTTCTGGCATCTCTGGCCAATACGGCAATGCCACAGTTATTTGCTCTGGCGCGGGAATATGCGGATAACTCg                                                  <  2:283669/73‑1 (MQ=255)
                    gtgCTTCTGGCATCTCTGGCCAATACGGCAATGCCACAGTTATTTGCTCTGGCGCGGTAATATGCGGATCACTCgg                                                 >  2:3198755/1‑76 (MQ=255)
                    gtgCTTCTGGCATCTCTGGCCAATACGGCAATGCCACAGTTATTTGCTCTGGCGCGGGAATATGCGGATAACTCgg                                                 >  1:2613822/1‑76 (MQ=255)
                    gtgCTTCTGGCATCTCTGGCCAATACGGCAATGCCACAGTTATTTGCTCTGGCGCGGGAATATGCGGATAACTCgg                                                 <  1:3235250/76‑1 (MQ=255)
                    gtgCTTCTGGCATCTCTGGCCAATACGGCAATGCCACAGTTATTTGCTCTGGCGCGGGAATATGCGGATAACTCgg                                                 >  2:3195877/1‑76 (MQ=255)
                     tgCTTCTGGCATCTCTGGCCAATACGGCAATGCCACAGTTATTTGCTCTGGCGCGGGAATATGCGGATAGCTCGgc                                                <  2:111548/76‑1 (MQ=255)
                     tgCTTCTGGCATCTCTGGCCAATACGGCAATGCCACAGTTATTTGCTCTGGCGCGGGAATATGCGGATAACTCGgc                                                <  1:3546374/76‑1 (MQ=255)
                     tgCTTCTGGCATCTCTGGCCAATACGGCAATGCCACAGTTATTTGCTCTGGCGCGGGAATATGCGGATAACTCGgc                                                <  2:3790164/76‑1 (MQ=255)
                       cTTCTGGCATCTCTGGCCAATACGGCAATGCCACAGTTATTTGCTCTGGCGCGGGAATATGCGGATAACTCGgcgc                                              >  2:1497178/1‑76 (MQ=255)
                        ttCTGGCATCTCTGGCCAATACGGCAATGCCACAGTTATTTGCTCTGGCGCGGGAATATGCGGATAACTCGGcgcg                                             >  2:965026/1‑76 (MQ=255)
                        ttCTGGCATCTCTGGCCAATACGGCAATGCCACAGTTATTTGCTCTGGCGCGGGAATATGCGGATAACTCGGcgcg                                             >  2:2509453/1‑76 (MQ=255)
                        ttCTGGCATCTCTGGCCAATACGGCAATGCCACAGTTATTTGCTCTGGCGCGGGAATATGCGGATAACTCGGcgcg                                             >  2:3406803/1‑76 (MQ=255)
                         tCTGGCATCTCTGGCCAATACGGCAATGCCACAGTTATTTGCTCTGGCGCGGGAATATGCGGATAACTCGGcgcgg                                            >  1:421535/1‑75 (MQ=255)
                         tCTGGCATCTCTGGCCAATACGGCAATGCCACAGTTATTTGCTCTGGCGCGGGAATATGCGGATAACTCGGcgcgg                                            >  2:2114840/1‑75 (MQ=255)
                         tCTGGCATCTCTGGCCAATACGGCAATGCCACAGTTATTTGCTCTGGCGCGGGAATATGCGGATAACTCGGcgcgg                                            >  2:1779503/1‑75 (MQ=255)
                         tCTGGCATCTCTGGCCAATACGGCAATGCCACAGTTATTTGCTCTGGCGCGGGAATATGCGGATAACTCGGCGCga                                            <  1:2742868/76‑1 (MQ=255)
                         tCTGGCATCTCTGGCCAATACGGCAATGCCACAGTTATTTGCTCTGGCGCGGGAATATGCGGATAACTCGGCGCga                                            <  1:1668641/76‑1 (MQ=255)
                         tCTGGCATCTCTGGCCAATACGGCAATGCCACAGTTATTTGCTCTGGCGCGGGAATATGCGGATAACTCGGCGCga                                            >  2:1425890/1‑76 (MQ=255)
                          cTGGCATCTCTGGCCAATACGGCAATGCCACAGTTATTTGCTCTGGCGCGGGAATATGCGGATAACTCGGCGCgag                                           >  2:4281117/1‑76 (MQ=255)
                           tGGCATCTCTGGCCAATACGGCAATGCCACAGTTATTTGCTCTGGCGCGGGAATATGCGGATAACTCGGCGCgaga                                          >  1:1674531/1‑76 (MQ=255)
                            ggCATCTCTGGCCAATACGGCAATGCCACAGTTATTTGCTCTGGCGCGGGAATATGCGGATAACTCGGCGCGAGaa                                         <  2:850649/76‑1 (MQ=255)
                              gaTCTCTGGCCAATACGGCAATGCCACAGTTATTTGCTCTGGCGCGGGAATATGCGGATAACTCGGCGCGAGAAgt                                       <  2:2199827/75‑1 (MQ=255)
                               aTCTCTGGCCAATACGGCAATGCCACAGTTATTTGCTCTGGCGCGGGAATATGCGGATAACTCGGCGCGAGAAgtg                                      <  2:1094976/76‑1 (MQ=255)
                               aTCTCTGGCCAATACGGCAATGCCACAGTTATTTGCTCTGGCGCGGGAATATGCGGATAACTCGGCGCGAGAAgtg                                      >  1:3411596/1‑76 (MQ=255)
                               aTCTCTGGCCAATACGGCAATGCCACAGTTATTTGCTCTGGCGCGGGAATATGCGGATAACTCGGCGCGAGAAgtg                                      <  1:4765751/76‑1 (MQ=255)
                                  tctGGCCAATACGGCAATGCCACAGTTATTTGCTCTGGCGCGGGAATATGCGGATAACTCGGCGCGAGAAgtggtg                                   >  1:1672309/1‑76 (MQ=255)
                                  tctGGCCAATACGGCAATGCCACAGTTATTTGCTCTGGCGCGGGAATATGCGGATAACTCGGCGCGAGAAgtggtg                                   >  1:4563865/1‑76 (MQ=255)
                                  tctGGCCAATACGGCAATGCCACAGTTATTTGCTCTGGCGCGGGAATATGCGGATAACTCGGCGCGAGAAgtggtg                                   >  2:2218617/1‑76 (MQ=255)
                                   ctGGCCAATACGGCAATGCCACAGTTATTTGCTCTGGCGCGGGAATATGCGGATAACTCGGCGCGAGAAGTGGTGa                                  >  1:4765289/1‑76 (MQ=255)
                                   ctGGCCAATACGGCAATGCCACAGTTATTTGCTCTGGCGCGGGAATATGCGGATAACTCGGCGCGAGAAGTGGTGa                                  <  1:3918450/76‑1 (MQ=255)
                                   ctGGCCAATACGGCAATGCCACAGTTATTTGCTCTGGCGCGGGAATATGCGGATAACTCGGCGCGAGAAGTGGTGa                                  >  2:2159088/1‑76 (MQ=255)
                                   ctGGCCAATACGGCAATGCCACAGTTATTTGCTCTGGCGCGGGAATATGCGGATAACTCGGCGCGAGAAGTGGTGa                                  >  1:1261387/1‑76 (MQ=255)
                                       ccAATACGGCAATGCCACAGTTATTTGCTCTGGCGCGGGAATATGCGGATAACTCGGCGCGAGAAGTGGTGATGtt                              <  2:888224/76‑1 (MQ=255)
                                       ccAATACGGCAATGCCACAGTTATTTGCTCTGGCGCGGGAATATGCGGATAACTCGGCGCGAGAAGTGGTGATGtt                              >  1:968469/1‑76 (MQ=255)
                                         aaTACGGCAATGCCACAGTTATTTGCTCTGGCGCGGGAATATGCGGATAACTCGGCGCGAGAAGTGGTGATGTTTa                            >  2:715543/1‑76 (MQ=255)
                                         aaTACGGCAATGCCACAGTTATTTGCTCTGGCGCGGGAATATGCGGATAACTCGGCGCGAGAAGTGGTGATGTTTa                            <  1:2223812/76‑1 (MQ=255)
                                         aaTACGGCAATGCCACAGTTATTTGCTCTGGCGCGGGAATATGCGGATAACTCGGCGCGAGAAGTGGTGATGTTTa                            <  1:2199365/76‑1 (MQ=255)
                                              ggCAATGCCACAGTTATTTGCTCTGGCGCGGGAATATGCGGATAACTCGGCGCGAGAAGTGGTGATGTTTAGCTCg                       <  1:3675405/76‑1 (MQ=255)
                                              ggCAATGCCACAGTTATTTGCTCTGGCGCGGGAATATGCGGATAACTCGGCGCGAGAAGTGGTGATGTTTAGCTCg                       <  2:2407086/76‑1 (MQ=255)
                                               gCAATGCCACAGTTATTTGCTCTGGCGCGGGAATATGCGGATAACTCGGCGCGAGAAGTGGTGATGTTTAGCTCgg                      >  2:3502104/1‑76 (MQ=255)
                                               gCAATGCCACAGTTATTTGCTCTGGCGCGGGAATATGCGGATAACTCGGCGCGAGAAGTGGTGATGTTTAGCTCgg                      <  1:4080493/76‑1 (MQ=255)
                                               gCAATGCCACAGTTATTTGCTCTGGCGCGGGAATATGCGGATAACTCGGCGCGAGAAGTGGTGATGTTTAGCTCgg                      >  1:2323647/1‑76 (MQ=255)
                                               gCAATGCCACAGTTATTTGCTCTGGCGCGGGAATATGCGGATAACTCGGCGCGAGAAGTGGTGATGTTTAGCTCgg                      >  1:736728/1‑76 (MQ=255)
                                                cAATGCCACAGTTATTTGCTCTGGCGCGGGAATATGCGGATAACTCGGCGCGAGAAGTGGTGATGTTTAGCTCGGt                     >  1:1634735/1‑76 (MQ=255)
                                                 aaTGCCACAGTTATTTGCTCTGGCGCGGGAATATGCGGATAACTCGGCGCGAGAAGTGGTGATGTTTAGCTCGGTg                    <  1:2846874/76‑1 (MQ=255)
                                                 aaTGCCACAGTTATTTGCTCTGGCGCGGGAATATGCGGATAACTCGGCGCGAGAAGTGGTGATGTTTAGCTCGGTg                    >  2:2969043/1‑76 (MQ=255)
                                                  aTGCCACAGATATTTGCTCTCGCGCGGGAATAGGCCGATAACTCGGCGCGAGACGAGGCGATGTTTAGCTCGGTgg                   >  2:934892/1‑75 (MQ=255)
                                                   tGCCACAGTTATTTGCTCTGGCGCGGGAATATGCGGATAACTCGGCGCGAGAAGTGGTGATGTTTAGCTCGGTGAt                  >  1:1206269/1‑76 (MQ=255)
                                                    gCCACAGTTATTTGCTCTGGCGCGGGAATATGCGGATAACTCGGCGCGAGAAGTGGTGATGTTTAGCTCGGTGAtg                 <  1:2023947/76‑1 (MQ=255)
                                                     ccgcaGTTATTTGCTCTGGCGCGGGAATATGCGGATAACTCGGCGCGAGAAGTGGTGATGTTTAGCTCGGTGAtgc                <  1:4731126/73‑1 (MQ=255)
                                                     ccACAGTTATTTGCTCTGGCGCGGGAATATGCGGATAACTCGGCGCGAGAAGTGGTGATGTTTAGCTCGGTGAtgc                <  2:3673743/76‑1 (MQ=255)
                                                     ccACAGTTATTTGCTCTGGCGCGGGAATATGCGGATAACTCGGCGCGAGAAGTGGTGATGTTTAGCTCGGTGAtgc                <  1:1425891/76‑1 (MQ=255)
                                                      cacaGTTATTTGCTCTGGCGCGGGAATATGCGGATAACTCGGCGCGAGACGCGGCGATGTTTAGCTCGGTGAtgcg               >  2:4717767/1‑76 (MQ=255)
                                                      cacaGTTATTTGCTCTGGCGCGGGAATATGCGGATAACTCGGCGCGAGAAGTGGTGATGTTTAGCTCGGTGAtgcg               <  1:4378037/76‑1 (MQ=255)
                                                      cacaGTTATTTGCTCTGGCGCGGGAATATGCGGATAACTCGGCGCGAGAAGTGGTGATGTTTAGCTCGGTGAtgcg               <  2:4086145/76‑1 (MQ=255)
                                                      cacaGCTATTTGCTCTGGCGCGGGAATATGCGGATAACTCGGCGCGAGAAGTGGTGATGTTTAGCTCGGTGAtgcg               >  1:2453050/1‑76 (MQ=255)
                                                       acaGTTATTTGCTCTGGCGCGGGAATATGCGGATAACTCGGCGCGAGAAGTGGTGATGTTTAGCTCGGTGAtgcgt              >  2:1531328/1‑76 (MQ=255)
                                                         aGTTATTTGCTCTGGCGCGGGAATATGCGGATAACTCGGCGCGAGAAGTGGTGATGTTTAGCTCGGTGAtgcgtgc            <  2:4608803/76‑1 (MQ=255)
                                                          gTTATTTGCTCTGGCGCGGGAATATGCGGATAACTCGGCGCGAGAAGTGGTGATGTTTAGCTCGGTGAtgcgtgcg           <  1:4467767/76‑1 (MQ=255)
                                                          gTTATTTGCTCTGGCGCGGGAATATGCGGATAACTCGGCGCGAGAAGTGGTGATGTTTAGCTCGGTGAtgcgtgcg           <  1:1354268/76‑1 (MQ=255)
                                                           ttATTTGCTCTGGCGCGGGAATATGCGGATAACTCGGCGCGAGAAGTGGTGATGTTTAGCTCGGTGATGCGTgcgc          <  1:2975780/76‑1 (MQ=255)
                                                           gtATTTGCTCTGGCGCGGGAATATGCGGATAACTCGGCGCGAGAAGTGGTGATGTTTAGCTCGGTGATGCGTgcgc          <  1:295070/75‑1 (MQ=255)
                                                            tATTTGCTCTGGCGCGGGAATATGCGGATAACTCGGCGCGAGAAGTGGTGATGTTTAGCTCGGTGATGCGTGCGCa         <  1:4512842/76‑1 (MQ=255)
                                                             aTTTGCTCTGGCGCGGGAATATGCGGATAACTCGGCGCGAGAAGTGGTGATGTTTAGCTCGGTGATGCGTGCGCAg        >  2:3632995/1‑76 (MQ=255)
                                                                tGCTCTGGGGCGGGAATATGCGGATAACTCGGCGCGAGAAGTGGTGATGTTTAGCTCGGTGATGCGTGCGCAGCtt     <  2:71911/76‑1 (MQ=255)
                                                                 gCTCTGGCGCGGGAATATGCGGATAACTCGGCGCGAGAAGTGGTGATGTTTAGCTCGGTGATGCGTGCGCAGCttt    <  1:1998718/76‑1 (MQ=255)
                                                                   tctGGCGCGGGAATATGCGGATAACTCGGCGCGAGAAGTGGTGATGTTTAGCTCGGTGAtgc                <  2:2099011/62‑1 (MQ=255)
                                                                   tctGGCGCGGGAATATGCGGATAACTCGGCGCGAGAAGTGGTGATGTTTAGCTCGGTGATGCGTGCGCAGCTTtct  >  2:4746628/1‑76 (MQ=255)
                                                                   tctGGCGCGGGAATATGCGGATAACTCGGCGCGAGAAGTGGTGATGTTTAGCTCGGTGATGCGTGCGCAGCTTtct  >  2:4400095/1‑76 (MQ=255)
                                                                   tctGGCGCGGGAATATGCGGATAACTCGGCGCGAGAAGTGGTGATGTTTAGCTCGGTGATGCGTGCGCAGCTTtct  <  1:4502361/76‑1 (MQ=255)
                                                                        |                                                                      
TGACGCTTATCACCTGTGGTGTGCTTCTGGCATCTCTGGCCAATACGGCAATGCCACAGTTATTTGCTCTGGAGCGGGAATATGCGGATAACTCGGCGCGAGAAGTGGTGATGTTTAGCTCGGTGATGCGTGCGCAGCTTTCT  >  NC_000913/77931‑78073

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