Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,793,867 C→A A378D (GCT→GAT)  gabT → 4‑aminobutyrate aminotransferase GabT

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,793,8670CA100.0% 507.3 / NA 129A378D (GCT→GAT) gabT4‑aminobutyrate aminotransferase GabT
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (86/43);  total (86/43)

Only 100 of 129 total aligned reads displayed.

CGGCGATCACAACAAGCCGGACGCCAAACTCACCGCCGAGATCGTGGCTCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCGGCCCGTAT  >  NC_000913/2793820‑2793913
                                               |                                              
cGGCGATCACAACAAGCCGGACGCCAAACTCACCGCCGAGATCGTGGATc                                              >  1:8219783/1‑50 (MQ=255)
cGGCGATCACAACAAGCCGGACGCCAAACTCACCGCCGAGATCGTGGATc                                              >  1:4735875/1‑50 (MQ=255)
 ggCGATCACAACAAGCCGGACGCCAAACTCACCGCCGAGATCGTGGATcg                                             >  1:7714440/1‑50 (MQ=255)
 ggCGATCACAACAAGCCGGACGCCAAACTCACCGCCGAGATCGTGGATcg                                             >  1:8281207/1‑50 (MQ=255)
  gCGATCGCAACAAGCCGGACGCCAAACTCACCGCCGAGATCGTGGATcgc                                            <  1:1769847/50‑1 (MQ=255)
  gCGATCACAACAAGCCGGACGCCAAACTCACCGCCGAGATCGTGGATcgc                                            >  1:7284449/1‑50 (MQ=255)
   cGATCACAACAAGCCGGACGCCAAACTCACCGCCGAGATCGTGGATcgcg                                           >  1:4315203/1‑50 (MQ=255)
   cGATCACAACAAGCCGGACGCCAAACTCACCGCCGAGATCGTGGATcgcg                                           >  1:5003771/1‑50 (MQ=255)
    gctCACAACAAGCCGGACGCCAAACTCACCGCCGAGATCGTGGATCgcgc                                          >  1:4724971/3‑50 (MQ=255)
    gATCACAACAAGCCGGACGCCAAACTCACCGCCGAGATCGTGGATCgcgc                                          >  1:3373603/1‑50 (MQ=255)
     aTCACAACAAGCCGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGcc                                         >  1:8429625/1‑50 (MQ=255)
      tCACACCAAGCCGCACGACAAACTCACCGCCGAGATCGTGGATCGCGccc                                        <  1:6542084/50‑1 (MQ=37)
      tCACAACAAGCCGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGccc                                        <  1:4013997/50‑1 (MQ=255)
       cacaACAAGCCGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCcg                                       >  1:6109767/1‑50 (MQ=255)
       cacaACAAGCCGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCcg                                       >  1:6426905/1‑50 (MQ=255)
        acaacaAGCCGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCcgc                                      <  1:201572/50‑1 (MQ=255)
        acaacaAGCCGGACGCAAAACTCACCGCCGAGATCGTGGATCGCGCCcgc                                      <  1:4429037/50‑1 (MQ=39)
        acaacaAGCCGGACGCAAAACTCACCGCCGAGATCGTGGATCGCGCCcg                                       <  1:2103430/49‑1 (MQ=38)
         caacaaGCCGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCcgc                                      >  1:8503942/1‑49 (MQ=255)
          aacaaGCCGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGa                                    >  1:7005778/1‑50 (MQ=255)
           acaaGCCGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGa                                    >  1:3298865/1‑49 (MQ=255)
           acaaGCCGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGa                                    >  1:6478612/1‑49 (MQ=255)
           acaaGCCGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGa                                    >  1:1544449/1‑49 (MQ=255)
            caaGCCGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATa                                  <  1:1378176/50‑1 (MQ=255)
             aaGCCGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATa                                  >  1:7860947/1‑49 (MQ=255)
              aGCCGGAGGCAAAACTATCCGCCGAGATCGTGGATCGCGCCCGCGATaaa                                <  1:2512257/50‑1 (MQ=255)
              aGCCGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATa                                  >  1:1993256/1‑48 (MQ=255)
              aGCCGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATa                                  >  1:4885869/1‑48 (MQ=255)
              aGCCGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATa                                  >  1:578858/1‑48 (MQ=255)
               gCCGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAg                               >  1:4954008/1‑50 (MQ=255)
                 cGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGc                             >  1:7222231/1‑50 (MQ=255)
                 cGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGc                             <  1:2445826/50‑1 (MQ=255)
                 cGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGc                             >  1:5600740/1‑50 (MQ=255)
                 cGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGc                             >  1:1608368/1‑50 (MQ=255)
                  ggACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGcc                            <  1:6228467/50‑1 (MQ=255)
                  ggACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGcc                            >  1:6468022/1‑50 (MQ=255)
                  ggACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGcc                            >  1:1402605/1‑50 (MQ=255)
                  ggACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGcc                            >  1:3500437/1‑50 (MQ=255)
                   gACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCt                           <  1:1547182/50‑1 (MQ=255)
                    aCGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCt                           >  1:8205149/1‑49 (MQ=255)
                    aCGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTg                          <  1:3041697/50‑1 (MQ=255)
                     cGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGa                         >  1:950953/1‑50 (MQ=255)
                      gCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGAt                        <  1:4879400/50‑1 (MQ=255)
                      gCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGAt                        <  1:3303555/50‑1 (MQ=255)
                      gCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGAt                        <  1:1964174/50‑1 (MQ=255)
                      gCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGAt                        >  1:4006953/1‑50 (MQ=255)
                       ccAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGAtt                       <  1:1250533/50‑1 (MQ=255)
                       ccAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGAtt                       >  1:2683390/1‑50 (MQ=255)
                        cAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGAttc                      <  1:1123383/50‑1 (MQ=255)
                        cAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGAttc                      <  1:1284515/50‑1 (MQ=255)
                        cAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGAttc                      <  1:2334282/50‑1 (MQ=255)
                         aaaCTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGAttct                     >  1:8587366/1‑50 (MQ=255)
                         aaaCTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGAttct                     >  1:7522917/1‑50 (MQ=255)
                         aaaCTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGAttct                     >  1:1886125/1‑50 (MQ=255)
                         aaaCTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGAttct                     >  1:2408850/1‑50 (MQ=255)
                          aaCTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGAttctt                    >  1:40361/1‑50 (MQ=255)
                          aaCTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGAttctt                    <  1:1723902/50‑1 (MQ=255)
                          aaCTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGAttctt                    <  1:3725042/50‑1 (MQ=255)
                          aaCTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGAttct                     >  1:6710206/1‑49 (MQ=255)
                           aCTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGAttcttc                   >  1:404000/1‑50 (MQ=255)
                           aCTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGAttcttc                   >  1:1361261/1‑50 (MQ=255)
                            cTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGAttcttca                  <  1:2420489/50‑2 (MQ=255)
                            cTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGATtcttct                  >  1:8221965/1‑50 (MQ=255)
                             tCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTtctc                 >  1:1856271/1‑50 (MQ=255)
                              cACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTtctct                >  1:5713568/1‑50 (MQ=255)
                              cACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTtctct                >  1:6980516/1‑50 (MQ=255)
                               aCCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTtctct                <  1:2900637/49‑1 (MQ=255)
                               aCCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTtctc                 <  1:7636059/48‑1 (MQ=255)
                               aCCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTtctc                 <  1:3548765/48‑1 (MQ=255)
                               aCCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTtctc                 <  1:6617242/48‑1 (MQ=255)
                                ccgccgAGATCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTcc              <  1:2652304/50‑1 (MQ=255)
                                ccgccgAGATCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTcc              >  1:3376636/1‑50 (MQ=255)
                                 cgccgAGATCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCt             >  1:4397414/1‑50 (MQ=255)
                                 cgccgAGATCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCt             >  1:8125349/1‑50 (MQ=255)
                                   ccgAGATCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGc           >  1:8171112/1‑50 (MQ=255)
                                   ccgAGATCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGc           >  1:6935676/1‑50 (MQ=255)
                                   ccgAGATCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGc           >  1:3050177/1‑50 (MQ=255)
                                   ccgAGATCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGc           >  1:5699455/1‑50 (MQ=255)
                                    cgAGATCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCg          >  1:3853722/1‑50 (MQ=255)
                                    cgAGATCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCg          <  1:2600198/50‑1 (MQ=255)
                                    cgAGATCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCg          <  1:3749678/50‑1 (MQ=255)
                                     gagaTCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCgg         >  1:2269699/1‑50 (MQ=255)
                                     gagaTCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCgg         <  1:328498/50‑1 (MQ=255)
                                     gagaTCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCgg         <  1:7928330/50‑1 (MQ=255)
                                      agaTCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCGGc        >  1:4880978/1‑50 (MQ=255)
                                      agaTCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCGGc        >  1:4262862/1‑50 (MQ=255)
                                      agaTCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCGGc        <  1:7873435/50‑1 (MQ=255)
                                        aTCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCGGccc      <  1:5684190/50‑1 (MQ=255)
                                         tCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCGGCCCg     >  1:8419620/1‑50 (MQ=255)
                                         tCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCGGCCCg     <  1:423279/50‑1 (MQ=255)
                                         tCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCGGCCCg     >  1:1464445/1‑50 (MQ=255)
                                         tCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCGGCCCg     >  1:2301092/1‑50 (MQ=255)
                                           gTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCGGCCCGt    <  1:7184940/49‑1 (MQ=255)
                                           gTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCGGCCCGt    <  1:6334251/49‑1 (MQ=255)
                                           gTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCGGCCCGTa   >  1:1720392/1‑50 (MQ=255)
                                           gTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCGGCCCGTa   >  1:2023797/1‑50 (MQ=255)
                                           gTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCGGCCCGTa   >  1:1292216/1‑50 (MQ=255)
                                            tGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCGGCCCGTAt  <  1:1949909/50‑1 (MQ=255)
                                            tGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCGGCCCGTAt  >  1:8370623/1‑50 (MQ=255)
                                            tGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCGGCCCGTAt  >  1:5993064/1‑50 (MQ=255)
                                               |                                              
CGGCGATCACAACAAGCCGGACGCCAAACTCACCGCCGAGATCGTGGCTCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCGGCCCGTAT  >  NC_000913/2793820‑2793913

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