Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,125,754 C→T 8.2% A294T (GCC→ACC)  pheS ← phenylalanine tRNA synthetase, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,125,7540CT 8.2% 229.8 / 8.9 73A294T (GCC→ACC) pheSphenylalanine tRNA synthetase, alpha subunit
Reads supporting (aligned to +/- strand):  ref base C (34/33);  new base T (2/4);  total (36/37)
Fisher's exact test for biased strand distribution p-value = 6.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.15e-01

CGAAGAATGAACGCAGGTCGGTGACGCCGTAACGCAACATAGTCAGACGCTCCATCCCCATCCCGAAGGCGAAACCAGAGTAAACTTCCGGGTCGATGCCAACGTTACGCAACACGTTCGGATGCACCATCCCG  >  minE/1125685‑1125818
                                                                     |                                                                
cGAAGAATGAACGCAGGTCGGTGACGCCGTAACGCAACATAGTCAGACGCTCCATCCCCATCCCGAAGGc                                                                  <  1:237916/70‑1 (MQ=255)
      aTGAACGCAGGTCGGTGACGCCGTAACGCAACATAGTCAGACGCTCCATCCCCATCCCGAAGGTGAAACCa                                                           <  1:3456110/71‑1 (MQ=255)
      aTGAACGCAGGTCGGTGACGCCGTAACGCAACATAGTCAGACGCTCCATCCCCATCCCGAAGGCGAAACCa                                                           <  1:3769093/71‑1 (MQ=255)
      aTGAACGCAGGTCGGTGACGCCGTAACGCAACATAGTCAGACGCTCCATCCCCATCCCGAAGGCGAAACCa                                                           <  1:3561657/71‑1 (MQ=255)
       tGAACGCAGGTCGGTGACGCCGTAACGCAACATAGTCAGACGCTCCATCCCCATCCCGAAGGCGAAACCag                                                          <  1:1023552/71‑1 (MQ=255)
       tGAACGCAGGTCGGTGACGCCGTAACGCAACATAGTCAGACGCTCCATCCCCATCCCGAAGGCGAAACCag                                                          <  1:796821/71‑1 (MQ=255)
          aCGCAGGTCGGTGACGCCGTAACGCAACATAGTCAGACGCTCCATCCCCATCCCGAAGGCTAAACCagag                                                        >  1:3263031/1‑70 (MQ=255)
          aCGCAGGTCGGTGACGCCGTAACGCAACATAGTCAGACGCTCCATCCCCATCCCGAAGGCGAAAGCAGAGt                                                       <  1:727/71‑1 (MQ=255)
          aCGCAGGTCGGTGACGCCGTAACGCAACATAGTCAGACGCTCCATCCCCATCCCGAAGGCGAAACCagag                                                        >  1:757022/1‑70 (MQ=255)
           cGCAGGTCGGTGACGCCGTAACGCAACATAGTCAGACGCTCCATCCCCATCCCGAAGGCGAAACCAGACTa                                                      <  1:2350776/71‑1 (MQ=255)
              aGGTCGGTGACGCCGTAACGCAACATAGTCAGACGCTCCATCCCCATCCCGAAGGCGAAACCAGAGTaa                                                     >  1:781505/1‑69 (MQ=255)
                 tCGGTGACGCCGTAACGCAACATAGTCAGACGCTCCATCCCCATCCCGAAGGCGAAACCAGAGTAAACTTc                                                <  1:3333122/71‑1 (MQ=255)
                  cGGTCACGCCGTAACGCAACATAGTCAGACGCTCCATCCCCATCCCGAAGGCGaaa                                                              >  1:3507934/1‑56 (MQ=255)
                   ggTGACGCCGTAACGCAACATAGTCAGACGCTCCATCCCCTTCCCGAAGGCGAAACCAGAGTAAACTTcc                                               <  1:647437/70‑1 (MQ=255)
                   ggTGACGCCGTAACGCAACATAGTCAGACGCTCCATCCCCATCCCGAAGGCGAAACCAGAGTAAACTTcc                                               <  1:3048488/70‑1 (MQ=255)
                   ggTGACGCCGTAACGCAACATAGTCAGACGCTCCATCCCCATCCCGAAGGCGAAACCAGAGTAAACTTCCg                                              >  1:1547924/1‑71 (MQ=255)
                   ggTGACGCCGTAACGCAACATAGTCAGACGCTCCATCCCCATCCCGAAGGCGAAACCAGAGTAAACTTCCg                                              >  1:1356695/1‑71 (MQ=255)
                    gTGACGCCGTAACGCAACATAGTCAGACGCTCCATCCCCATCCCGAAGGCGAAACCAGAGTAAACTTCCgg                                             <  1:3039758/71‑1 (MQ=255)
                       aCGCCGTAACGCAACATAGTCAGACGCTCCATCCCCATCCCGAAGGTGAAACCAGAGTAAACTTCCGGGTc                                          <  1:2386924/71‑1 (MQ=255)
                       aCGCCGTAACGCAACATAGTCAGACGCTCCATCCCCATCCCGAAGGCGAAACCAGAGTAAACTTCCGGGTc                                          <  1:1996570/71‑1 (MQ=255)
                         gCCGTAACGCAACATAGTCAGACGCTCCATCCCCATCCCGAAGGCGAAACCAGAGTAAACTTCCGGGTCGa                                        <  1:1000120/71‑1 (MQ=255)
                         gCCGTAACGCAACATAGTCAGACGCTCCATCCCCATCCCGAAGGCGAAACCAGAGTAAACTTCCGGGTCGa                                        <  1:1070330/71‑1 (MQ=255)
                         gCCGTAACGCAACATAGTCAGACGCTCCATCCCCATCCCGAAGGCGAAACCAGAGTAAACTTCCGGGTCGa                                        <  1:2227245/71‑1 (MQ=255)
                          ccGTAACGCAACATAGTCAGACGCTCCATCCCCATCCCGAAGGCGAAACCAGAGTa                                                      <  1:2017155/56‑1 (MQ=255)
                           cGTAACGCAACATAGTCAGACGCTCCATCCCCATCCCGAAGGCGAAACCa                                                           <  1:127502/50‑1 (MQ=255)
                           cGTAACGCAACATAGTCAGACGCTCCATCCCCATCCCGAAGGCGAAACCAGAGTAAACTTCCGGGTCGATg                                      >  1:2359973/1‑71 (MQ=255)
                                 gCAACATAGTCAGACGCTCCATCCCCATCCCGAAGGCGAAACCAGAGTAAACTTCCGGGTCGATGCCAACg                                >  1:3588549/1‑71 (MQ=255)
                                  cAACATAGTCAGACGCTCCATCCCCATCCCGAAGTCGAAACCAGAGTAAACTTCCGGGTCGATGCCAACGt                               <  1:2924199/71‑1 (MQ=255)
                                  cAACATAGTCAGACGCTCCATCCCCATCCCGAAGGCGAAACCAGAGTAAACTTCCGGGTCGATGCCAACGt                               <  1:2497240/71‑1 (MQ=255)
                                      aTAGTCAGACGCTCCATCCCCATCCCGAAGGCGAAACCAGAGTAAACTTCCGGGTCGATGCCAACGTTACg                           >  1:2525025/1‑71 (MQ=255)
                                       tAGTCAGACGCTCCATCCCCATCCCGAAGGCGAAACCa                                                           <  1:3408405/38‑1 (MQ=255)
                                         gTCAGACGCTCCATCCCCATCCCGACGGCGAAACCAGAGTAAACTTCCGGGTCGATGCCAACGTTACGCa                         <  1:2688737/70‑1 (MQ=255)
                                         gTCAGACGCTCCATCCCCATCCCGAAGGCGAAACCAGAGTAAACTTCCGGGTCGATGCCAACGTTACGCa                         <  1:2328360/70‑1 (MQ=255)
                                            aGACGCTCCATCCCCATCCCGAAGGTGAAACCAGAGTAAACTTCCggg                                            >  1:2307966/1‑48 (MQ=255)
                                            aGACGCTCCATCCCCATCCCGAAGGCGAAACCAGAGTAAACTTCCGGGTCGATGCCAACGTTACGCAacac                     >  1:3876163/1‑71 (MQ=255)
                                            aGACGCTCCATCCCCATCCCGAAGGCGAAACCAGAGTAAACTTCCGGGTCGATGCCAACGTTACGCAacac                     >  1:2939741/1‑71 (MQ=255)
                                             gACGCTCCATCCCCATCCCGAAGGCGAAACCAGAGTAAACTTCCGGGTCGATGCCAACGTTACGCAACACg                    <  1:1042582/71‑1 (MQ=255)
                                               cGCTCCATCCCCATCCCGAAGGCGAAACCAGAGTAAACTTCCGGGTCGATGCCAACGTTACGCa                         >  1:2540988/1‑64 (MQ=255)
                                                 cTCCATCCCCATCCCGAAGGCGAAACCAGAGTAAACTTCCGGGTCGATGCCAACGTTACGCAACACGTTCg                <  1:724606/71‑1 (MQ=255)
                                                   ccatccccatccCGAAGGCGAAACCAGAGTAAACTTCCGGGTCGATGCCaa                                  >  1:3127384/1‑51 (MQ=255)
                                                    catccccatcccGAAGGCGAAACCAGAGTAAACTTCCGGGTCTATGCCAACGTTACGCAACACGTTCGGAt             <  1:1794574/71‑1 (MQ=255)
                                                    catccccatcccGAAGGCGAAACCAGAGTAAACTTCCGGGTCGATGCCAACGTTACGCAACACGTTCGGa              <  1:3488152/70‑1 (MQ=255)
                                                    catccccatcccGAAGGCGAAACCAGAGTAAACTTCCGGGTCGATGCCAACGTTACGCAACACGTTCGGa              <  1:760300/70‑1 (MQ=255)
                                                    catccccatcccGAAGGCGAAACCAGAGTAAACTTCCGGGTCGATGCCAACGTTACGCAACACGTTCGGAt             >  1:3664721/1‑71 (MQ=255)
                                                      tccccatcccGAAGGCGAAACCAGAGTAAACTTCCGGGTCGATGCCAACGTTACGc                          <  1:1514064/56‑1 (MQ=255)
                                                       ccccatcccGAAGGCGAAACCAGAGTAAACTTCCGGGTCGATGCCAACGTTACGCAACACGTTCGGATGCa          >  1:2401950/1‑71 (MQ=255)
                                                       ccccatcccGAAGGCGAAACCAGAGTAAACTTCCGGGTCGATGCCAACGTTACGCAACACGTTCGGATGCa          >  1:3990082/1‑71 (MQ=255)
                                                       ccccatcccGAAGGCGAAACCAGAGTAAACTTCCGGGTCGATGCCAACGTTACGCAACACGTTCGGATGCa          >  1:491151/1‑71 (MQ=255)
                                                       ccccatcccGAAGGCGAAACCAGAGTAAACTTCCGGGTCGATGCCAACGTTACGCAACACGTTCGGATGCa          >  1:2570021/1‑71 (MQ=255)
                                                        cccatcccGAAGGTGAAACCAGAGTAAACTTCCGGGTCGATGCCAACGTTACGCAACACGTTAGGATGCAc         >  1:3465058/1‑71 (MQ=255)
                                                        cccatcccGAAGGCGAAACCAGAGTAAACTTCCGGGTCGATGCCAACGTTACGCAACACGTTGGGAt             >  1:3651304/1‑67 (MQ=255)
                                                        cccatcccGAAGGCGAAACCAGAGTAAACTTCCGGGTCGATGCCAACGTTACGCAACACGTTCGGATGCa          <  1:1950097/70‑1 (MQ=255)
                                                        cccatcccGAAGGCGAAACCAGAGTAAACTTCCGGGTCGATGCCAACGTTACGCAACACGTTCGGATGCa          <  1:3880197/70‑1 (MQ=255)
                                                           atcccGAAGGCGAAACCAGAGTAAACTTCCGGGTCGATGCCAACGTTACGCAACACGTTCGGATGCACCAt      >  1:656985/1‑71 (MQ=255)
                                                           atcccGAAGGCGAAACCAGAGTAAACTTCCGGGTCGATGCCAACGTTACGCAACACGTTCGGATGCACCAt      >  1:3166429/1‑71 (MQ=255)
                                                           atcccGAAGGCGAAACCAGAGTAAACTTCCGGGTCGATGCCAACGTTACGCAACACGTTCGGATGCACCAt      >  1:2171650/1‑71 (MQ=255)
                                                           ataccGAAGGCGAAACCAGAGTAAACTTCCGGGTCGATGCCAACGTTACGCAACACGTTCGGATGCACCAt      >  1:3463903/4‑71 (MQ=255)
                                                            tcccGAAGGTGAAACCAGAGTAAACTTCCGGGTCGATGCCAAc                                 <  1:3286892/43‑1 (MQ=255)
                                                            tcccGAAGGCGAAACCAGAGTAAACTTCCGGGTCGATga                                     >  1:787099/1‑38 (MQ=255)
                                                            tcccGAAGGCGAAACCAGAGTAAACTTCCGGGTCGATGCCAACGTTACGCAACACGTTCGGATGCACCAt      <  1:2002800/70‑1 (MQ=255)
                                                            tcccGAAGGCGAAACCAGAGTAAACTTCCGGGTCGATGCCAACGTTACGCAACACGTTCGGATGCACCAt      <  1:2007551/70‑1 (MQ=255)
                                                            tcccGAAGGCGAAACCAGAGTAAACTTCCGGGTCGATGCCAACGTTACGCAACACGTTCGGATGCACCATc     >  1:2240198/1‑71 (MQ=255)
                                                            tcccGAAGGCGAAACCAGAGTAAACTTCCGGGTCGATGCCAACGTTACGCAACACGTTCGGATGCACCATc     >  1:567011/1‑71 (MQ=255)
                                                              ccGAAGGTGAAACCAGAGTAAACTTCCGGGTCGATGCCAACGTTACGCAACACGTTCGGATGCACCATccc   <  1:2681396/71‑1 (MQ=255)
                                                              ccGAAGGCGAAACCAGAGTAAACTTCCGGGTCGATGCCAACGTTACGCAACACGtt                  >  1:2665135/1‑56 (MQ=255)
                                                              ccGAAGGCGAAACCAGAGTAAACTTCCGGGTCGATGCCAACGTTACGCAACACGTTCGGATGCACCATcc    >  1:3578573/1‑70 (MQ=255)
                                                              ccGAAGGCGAAACCAGAGTAAACTTCCGGGTCGATGCCAACGTTACGCAACACGTTCGGATGCACCATcc    >  1:2905231/1‑70 (MQ=255)
                                                              ccGAAGGCGAAACCAGAGTAAACTTCCGGGTCGATGCCAACGTTACGCAACACGTTCGGATGCACCATc     >  1:2324052/1‑69 (MQ=255)
                                                              acGAAGGCGAAACCAGAGTAAACTTCCGGGTCGATGCCAACGTTACGc                          >  1:1474624/2‑48 (MQ=255)
                                                                gAAGGCGAAACCAGAGTAAACTTCCGGGTCGATGCCAACGTTACGCaa                        <  1:237353/48‑1 (MQ=255)
                                                                 aaGGCGAAACCAGAGTAAACTTCCGGGTCGATGCCAACGTTACGCaa                        >  1:1756843/1‑47 (MQ=255)
                                                                 aaGGCGAAACCAGAGTAAACTTCCGGGTCGATGCCAACGTTACGCaa                        >  1:2951798/1‑47 (MQ=255)
                                                                 aaGGCGAAACCAGAGTAAACTTCCGGGTCGATGCCAACGTTACGCAACACGTTCg                <  1:1973378/55‑1 (MQ=255)
                                                                 aaGGCGAAACCAGAGTAAACTTCCGGGTCGATGCCAACGTTACGCAACACGTTCGGATGCACCATCCCg  >  1:2152374/1‑69 (MQ=255)
                                                                     |                                                                
CGAAGAATGAACGCAGGTCGGTGACGCCGTAACGCAACATAGTCAGACGCTCCATCCCCATCCCGAAGGCGAAACCAGAGTAAACTTCCGGGTCGATGCCAACGTTACGCAACACGTTCGGATGCACCATCCCG  >  minE/1125685‑1125818

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