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,8670CA75.7% 209.8 / 77.4 107A378D (GCT→GAT) gabT4‑aminobutyrate aminotransferase GabT
Reads supporting (aligned to +/- strand):  ref base C (20/6);  new base A (53/28);  total (73/34)
Fisher's exact test for biased strand distribution p-value = 3.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.21e-01

Only 100 of 108 total aligned reads displayed.

GACGGCGATCACAACAAGCCGGACGCCAAACTCACCGCCGAGATCGTGGCTCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCGGCCCGTAT  >  NC_000913/2793818‑2793913
                                                 |                                              
gACGGCGATCACAACAAGCCGGACGCCAAACTCACCGCCGAGATCGTGGc                                                >  2:550356/1‑50 (MQ=255)
  cGGCGATCACAACAAGCCGGACGCCAAACTCACCGCCGAGATCGTGGATc                                              <  2:3153690/50‑1 (MQ=255)
   ggCGATCACAACAAGCCGGACGCCAAACTCACCGCCGAGATCGTGGCTcg                                             >  1:560537/1‑50 (MQ=255)
   ggCGATCACAACAAGCCGGACGCCAAACTCACCGCCGAGATCGTGGATcg                                             >  2:2568579/1‑50 (MQ=255)
     cGATCACAACAAGCCGGACGCCAAACTCACCGCCGAGATCGTGGCTcgcg                                           >  1:2536092/1‑50 (MQ=255)
      gATCACAACAAGCCGGACGCCAAACTCACCGCCGAGATCGTGGCTCgcgc                                          >  1:5078942/1‑50 (MQ=255)
      gATCACAACAAGCCGGACGCCAAACTCACCGCCGAGATCGTGGCTCgcgc                                          >  1:3449335/1‑50 (MQ=255)
      gATCACAACAAGCCGGACGCCAAACTCACCGCCGAGATCGTGGCTCgcgc                                          >  1:2921154/1‑50 (MQ=255)
      gATCACAACAAGCCGGACGCCAAACTCACCGCCGAGATCGTGGCTCgagc                                          >  2:3631220/1‑50 (MQ=255)
      gATCACAACAAGCCGGACGCCAAACTCACCGCCGAGATCGTGGATCgcgc                                          >  1:569827/1‑50 (MQ=255)
        tCACAACAAGCCGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGccc                                        <  1:3808735/50‑1 (MQ=255)
        tCACAACAAGCCGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGccc                                        >  1:1972827/1‑50 (MQ=255)
         cacaACAAGCCGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCcg                                       >  2:2746787/1‑50 (MQ=255)
          ccaacaaCCCGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCcgc                                      <  2:143970/49‑1 (MQ=255)
          acaacaAGCCGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCcgc                                      <  2:474534/50‑1 (MQ=255)
          acaacaAGCCGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCcgc                                      >  2:3123958/1‑50 (MQ=255)
          acaacaAGCCGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCcgc                                      >  2:4753766/1‑50 (MQ=255)
           caacaaGCCGGACGCCAAACTCACTGCCGAGATCGTGGATCGCGCCcgc                                      >  2:4699555/1‑49 (MQ=38)
           caacaaGCCGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCcgc                                      >  2:1128086/1‑49 (MQ=255)
            aacaaGCCGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGa                                    >  1:2725014/1‑50 (MQ=255)
            aacaaGCCGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGa                                    >  1:4834372/1‑50 (MQ=255)
              caaGCCGGGCGACAAACCCCCCGCCGAGATCGTGGATCGCGCCCGCGATa                                  <  2:5321070/50‑1 (MQ=255)
              caaGCCGGACGCCAAACTCACCGCCGAGATCGTGGCTCGCGCCCGCGAt                                   >  2:1773281/1‑49 (MQ=255)
              caaGCCGGACGCCAAACTCACCGCCGAGATCGTGGCTCGCGCCCGCGATa                                  <  1:3931370/50‑1 (MQ=255)
              caaGCCGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGAt                                   >  2:2581268/1‑49 (MQ=255)
              caaGCCGAACGCCAAACTCACCGCCGAGATCGTGGCTCGCGCCCGCGATa                                  <  1:1028884/50‑1 (MQ=255)
              caaGACGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATa                                  <  1:3143887/50‑1 (MQ=255)
               aaGCCGGACGCCAAACTCGCCGCCGAGATCGTGGATCGCGCCCGCGATaa                                 <  2:663155/50‑1 (MQ=255)
               aaGCCGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATa                                  >  2:1909399/1‑49 (MQ=255)
                aGCCGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATaaa                                <  1:1401715/50‑1 (MQ=255)
                 gCCGGGCGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAg                               <  2:1444098/50‑1 (MQ=255)
                 gCCGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAg                               >  1:4787892/1‑50 (MQ=255)
                 gCCGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAg                               >  1:4840120/1‑50 (MQ=255)
                  ccGGACGCCACACTCACCGCCGAGATCGTGGCTCGCGCCCCCGATAAAgg                              >  2:2029319/1‑50 (MQ=255)
                  ccGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAgg                              >  1:3690070/1‑50 (MQ=255)
                   cGGACGCCAAACTCACCGCCGAGATCGTGGCTCGCGCCCGCGATAAAGGc                             >  2:3850671/1‑50 (MQ=255)
                   cGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGc                             >  1:2188917/1‑50 (MQ=255)
                   cGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGc                             <  1:2050068/50‑1 (MQ=255)
                   cGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGc                             >  2:3016370/1‑50 (MQ=255)
                   cGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGc                             >  1:5324649/1‑50 (MQ=255)
                    ggACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGcc                            <  1:4251728/50‑1 (MQ=255)
                    ggACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGcc                            >  1:3233153/1‑50 (MQ=255)
                     ggcGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCt                           <  2:281181/48‑1 (MQ=255)
                     gACGCCAAACTCTCCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCt                           <  2:889203/50‑1 (MQ=39)
                     gACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAACGGCCt                           >  1:3261165/1‑50 (MQ=255)
                     gACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCt                           >  1:4949902/1‑50 (MQ=255)
                     gACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCt                           <  1:3749159/50‑1 (MQ=255)
                      aCGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCt                           >  1:1818943/1‑49 (MQ=255)
                       cGCCAAACTCACCGCCGAGATCGTGGCTCGCGCCCGCGATAAAGGCCTGa                         <  1:5145350/50‑1 (MQ=255)
                       cGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGa                         <  2:2580493/50‑1 (MQ=255)
                        gCCAATCTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGAt                        <  1:4906783/50‑1 (MQ=255)
                        gCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGAt                        >  2:3870081/1‑50 (MQ=255)
                        gCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGAt                        >  1:1997779/1‑50 (MQ=255)
                        gCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGAt                        <  2:5093674/50‑1 (MQ=255)
                        gCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGAt                        >  2:3067156/1‑50 (MQ=255)
                          cAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGAttc                      >  1:2606561/1‑50 (MQ=255)
                           aaaCTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGAttct                     <  2:2896122/50‑1 (MQ=255)
                           aaaCTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGAttct                     >  2:4091585/1‑50 (MQ=255)
                            aaCTCACCGCCGAGATCGTGGCTCGCGCCCGCGATAAAGGCCTGAttctt                    >  2:1571454/1‑50 (MQ=255)
                            aaCTCACCGCCGAGATCGTGGCTCGCGCCCGCGATAAAGGCCTGAttctt                    >  1:2517109/1‑50 (MQ=255)
                            aaCTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGAttctt                    >  2:589305/1‑50 (MQ=255)
                            aaCTAACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGAttctt                    >  2:2686373/1‑50 (MQ=255)
                             aCTCACCGCCGAGATCGTGGCTCGCGCCCGCGATAAAGGCCTGAttcttc                   >  1:4805836/1‑50 (MQ=255)
                             aCTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGAttcttc                   >  1:2042185/1‑50 (MQ=255)
                              cTCACCGCCGAGATCGTGGCTCGCGCCCGCGATAAAGGCCTGATtcttct                  >  1:5246028/1‑50 (MQ=255)
                               tCACCGCCGAGATCGTGGCTCGCGCCCGCGATAAAGGCCTGATtcttctt                 >  2:3268357/1‑49 (MQ=255)
                               tCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTtctc                 >  2:762395/1‑50 (MQ=255)
                               tCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTtctc                 >  1:2454925/1‑50 (MQ=255)
                                cACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTtctct                >  2:4906396/1‑50 (MQ=255)
                                 aCCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTtctctc               <  1:1850921/50‑1 (MQ=255)
                                 aCCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTtctct                <  2:1961568/49‑1 (MQ=255)
                                 aCCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTtctct                <  1:1328086/49‑1 (MQ=255)
                                  ccgccgAGATCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTcc              >  1:4652088/1‑50 (MQ=255)
                                  ccgccgAGATCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTcc              <  1:4699556/50‑1 (MQ=255)
                                   cgccgAGATCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCt             >  2:3614605/1‑50 (MQ=255)
                                    gccgAGATCGTGGCTCGCGCCCGCGATAAAGGCCTGATTTTTCTCTCCTg            >  2:349047/1‑50 (MQ=255)
                                    gccgAGATCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTg            <  2:2541870/50‑1 (MQ=255)
                                    gccgAGATCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTg            >  2:5332555/1‑50 (MQ=255)
                                    gccgAGATCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTg            >  2:3994872/1‑50 (MQ=255)
                                     ccgAGATCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGc           >  2:1016523/1‑50 (MQ=255)
                                     ccgAGATCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGc           >  1:174961/1‑50 (MQ=255)
                                      cgAGATCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCg          <  2:4878685/50‑1 (MQ=255)
                                       gagaTCGTGGCTCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCgg         >  2:3931999/1‑50 (MQ=255)
                                       gagaTCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCgg         >  1:3243136/1‑50 (MQ=255)
                                        agaTCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCGGc        >  2:2811446/1‑50 (MQ=255)
                                        agaTCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCGGc        >  1:2128803/1‑50 (MQ=255)
                                         gaTCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCGGcc       <  1:4155798/50‑1 (MQ=255)
                                         gaTCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCGGcc       >  2:5142624/1‑50 (MQ=255)
                                          aTCGTGGCTCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCGGccc      >  2:5323450/1‑50 (MQ=255)
                                          aTCGTGGCTCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCGGccc      >  1:4129682/1‑50 (MQ=255)
                                           tCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCGGCCCg     >  1:1230571/1‑50 (MQ=255)
                                           tCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCGGCCCg     <  1:2236794/50‑1 (MQ=255)
                                           tCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCGGCCCg     >  2:5078866/1‑50 (MQ=255)
                                            agTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCGGCCCg     >  2:4576285/2‑49 (MQ=255)
                                             gTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCGGCCCGt    <  1:2537367/49‑1 (MQ=255)
                                             gTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCGGCCCGTa   >  1:4354385/1‑50 (MQ=255)
                                             gTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCGGCCCGTa   >  2:785449/1‑50 (MQ=255)
                                              tGGCTCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCGGCCCGTAt  >  2:3212802/1‑50 (MQ=255)
                                              tGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCGGCCCGTAt  >  1:450053/1‑50 (MQ=255)
                                              tGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCGGCCCGTAt  >  1:56572/1‑50 (MQ=255)
                                                 |                                              
GACGGCGATCACAACAAGCCGGACGCCAAACTCACCGCCGAGATCGTGGCTCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCGGCCCGTAT  >  NC_000913/2793818‑2793913

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