Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,923,265 G→T 12.5% L653M (CTG→ATG)  valS ← valyl‑tRNA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,923,2650GT12.5% 212.9 / 11.3 64L653M (CTG→ATG) valSvalyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base G (30/26);  new base T (3/5);  total (33/31)
Fisher's exact test for biased strand distribution p-value = 4.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

CAGCGGTCCGCCAGCGACAGCGTCATTTCGCCGCCGTTGAAGCCGCAATCCTGACCTTCTGTGTTCATCAGCACAAAGCGGCTGGCGTTCCACAGCTTGTTACAGAAGTTACGGTAACCTTCCAGACGCTTCATATCC  >  minE/2923195‑2923332
                                                                      |                                                                   
cAGCGGTCCGCCAGCGACAGCGTCATTTCGCCGCCGTTGAAGCCGCAATCCTGACCTTCTGTGTTCATCAg                                                                     <  1:2927101/71‑1 (MQ=255)
cAGCGGTCCGCCAGCGACAGCGTCATTTCGCCGCCGTTGAAGCCGCAATCCTGACCTTCTGTGTTCATCAg                                                                     <  1:2019174/71‑1 (MQ=255)
         gCCAGCGACAGCGTCATTTCGCCGCCGTTGAAGCCGCAATCCTGACCTTCTGTGTTCATCAGCACAAAGCg                                                            <  1:462274/71‑1 (MQ=255)
               gACAGCGTCATTTCGCCGCCGTTGAAGCCGCAATCCTGACCTTCTGTGTTCATCAGCACAAAGCGGCTGGc                                                      >  1:1113180/1‑71 (MQ=255)
                aCAGCGTCATTTCGCCGCCGTTGAAGCCGCAATCCTGACCTTCTGTGTTCATCAGCACAAAGCGGCTGGCg                                                     >  1:2412348/1‑71 (MQ=255)
                aCAGCGTCATTTCGCCGCCGTTGAAGCCGCAATCCTGACCTTCTGTGTTCATCAGCACAAAGCGGCTGGCg                                                     >  1:2950589/1‑71 (MQ=255)
                 cAGCGTCATTTCGCCGCCGTTGAAGCCGCAATCCTGACCTTCTGTGTTCATCAGCACAAAGCGGCTGGCGt                                                    <  1:494890/71‑1 (MQ=255)
                 cAGCGTCATTTCGCCGCCGTTGAAGCCGCAATCCTGACCTTCTGTGTTCATCAGCACAAAGCGGCTGGCGt                                                    >  1:1802274/1‑71 (MQ=255)
                 cAGCGTCATTTCGCCGCCGTTGAAGCCGCAATCCTGACCTTCTGTGTTCATCAGCACAAAGCGGCTGGCGt                                                    <  1:1846282/71‑1 (MQ=255)
                  aGCGTCATTTCGCCGCCGTTGAAGCCGCAATCCTGACCTTCTGTGTTCATCAGCACAAAGCGGCTGGCg                                                     >  1:722794/1‑69 (MQ=255)
                   gCGTCGTTTCGCCGCCGTTGAAGCCGCAATCCTGACCTTCTGTGTTCATCAGCACAAAGCGGCTGGCGTTc                                                  <  1:1652077/71‑1 (MQ=255)
                   gCGTCATTTCGCCGCCGTTGAAGCCGCAATCCTGACCTTCTGTGTTCATCAGCACAAAGCGGCTGGCGTTc                                                  <  1:492703/71‑1 (MQ=255)
                   gCGTCATTTCGCCGCCGTTGAAGCCGCAATCCTGACCTTCTGTGTTCATCAGCACAAAGCGGCTGGCGTTc                                                  <  1:3266870/71‑1 (MQ=255)
                   gCGTCATTTCGCCGCCGTTGAAGCCGCAATCCTGACCTTCTGTGTTCATCAGCACAAAGCGGCTGGCGTTc                                                  <  1:2014190/71‑1 (MQ=255)
                     gTCATTTCGCCGCCGTTGAAGCCGCAATCCTGACCTTCTGTGTTCATCAGCACAAAGCGGCTGGCGTTCca                                                <  1:1374307/71‑1 (MQ=255)
                     gTCATTTCGCCGCCGTTGAAGCCGCAATCCTGACCTTCTGTGTTCATCAGCACAAAGCGGCTGGCGTTCca                                                <  1:2976266/71‑1 (MQ=255)
                        aTTTCGCCGCCGTTGAAGCCGCAATCCTGACCTTCTGTGTTCATCAGCACAAAGCGGCTGGCGTTCCACAg                                             >  1:1908120/1‑71 (MQ=255)
                         tttCGCCGCCGTTGAAGCCGCAATCCTGACCTTCTGTGTTCATCAGCACAAAGCGGCTGGCGTTCCACAGc                                            >  1:1802991/1‑71 (MQ=255)
                          ttCGCCGCCGTTGAAGCCGCAATCCTGACCTTCTGTGTTCATCAGca                                                                   >  1:1796719/1‑47 (MQ=255)
                          ttCGCCGCCGTTGAAGCCGCAATCCTGACCTTCTGTGTTCATCAGca                                                                   >  1:2599439/1‑47 (MQ=255)
                          ttCGCCGCCGTTGAAGCCGCAATCCTGACCTTCTGTGTTCATCAGCACaaa                                                               >  1:2231413/1‑51 (MQ=255)
                          ttCGCCGCCGTTGAAGCCGCAATCCTGACCTTCTGTGTTCATCAGCACAAAGCGGCTGGc                                                      >  1:3016265/1‑60 (MQ=255)
                          ttCGCCGCCGTTGAAGCCGCAATCCTGACCTTCTGTGTTCATCAGCACAAAGCGGCTGGc                                                      >  1:2736982/1‑60 (MQ=255)
                             gccgccGTTGAAGCCGCAATCCTGACCTTCTGTGTTCATCAGCACAAAGCGGCTGGCGTTCCACAg                                             <  1:591889/66‑1 (MQ=255)
                               cgccgTTGAAGCCGCAATCCTGACCTTCTGTGTTCATCATCACAAAGCGGCTGGCGTTCCACAGCTTGTTa                                      <  1:3041354/71‑1 (MQ=255)
                               cgccgTTGAAGCCGCAATCCTGACCTTCTGTGTTCATCAGCACAAAGCGGCTGGCGTTCCACAGCt                                           >  1:27520/1‑66 (MQ=255)
                               cgccgTTGAAGCCGCAATCCTGACCTTCTGTGTTCATCAGCACAAAGCGGCTGGCGTTCCACAGCt                                           >  1:819965/1‑66 (MQ=255)
                               cgccgTTGAAGCCGCAATCCTGACCTTCTGTGTTCATCAGCACAAAGCGGCTGGCGTTCCACAGCTTGTTa                                      <  1:594717/71‑1 (MQ=255)
                               cgccgTTGAAGCCGCAATCCTGACCTTCTGTGTTCATCAGCACAAAGCGGCTGGCGTTCAACAGCt                                           >  1:2182499/1‑66 (MQ=255)
                                gccgTTGAAGCCGCAATCCTGACCTTCTGTGTTCATCAGCACAAAGCGGCTGGCGTTCCACAGCTTGTTa                                      >  1:2374124/1‑70 (MQ=255)
                                 ccgTTGAAGCCGCAATCCTGACCTTCTGTGTTCATCAGCACAAAGCGGCTGGCGTTCCACAGCTTGTTa                                      >  1:1063362/1‑69 (MQ=255)
                                 ccgTTGAAGCCGCAATCCTGACCTTCTGTGTTCATCAGCACAAAGCGGCTGGCGTTCCACAGCTTGTTa                                      >  1:1141252/1‑69 (MQ=255)
                                  cgTTGAAGCCGCAATCCTGACCTTCTGTGTTCATCAGCACAAAGCGGCTGGCGTTCCACAGCTTGTTa                                      >  1:605590/1‑68 (MQ=255)
                                    ttGAAGCCGCAATCCTGACCTTCTGTGTTCATCAGcaca                                                                 <  1:3143350/39‑1 (MQ=255)
                                     tGAAGCCGCAATCCTGACCTTCTGTGTTCATCAGCACAAAGCGGc                                                          >  1:1505174/1‑45 (MQ=255)
                                     tGAAGCCGCAATCCTGACCTTCTGTGTTCATCAGCACAAAGCGGc                                                          >  1:1433483/1‑45 (MQ=255)
                                      gAAGCCGCAATCCTGACCTTCTGTGTTCATCAGCACaaa                                                               >  1:138022/1‑39 (MQ=255)
                                         gCCGCAATCCTGACCTTCTGTGTTCATCAGCACAAAGCGGCTGGCGTTCCACAGCTTGt                                        <  1:3396312/59‑1 (MQ=255)
                                         gCCGCAATCCTGACCTTCTGTGTTCATCAGCACAAAGCGGCTGGCGTTCCACAGCTTGt                                        <  1:2360525/59‑1 (MQ=255)
                                           cGCAATCCTGACCTTCTGTGTTCATCAGCACAAAGCGGCTGGCGTTCCACAGCTTGTTACAGAAGTTACgg                          >  1:313003/1‑71 (MQ=255)
                                            gCAATCCTGACCTTCTGTGTTCATCAGCACAAAGCGGCTGGCGTTCCACAGCTTGTTACAg                                   >  1:2218269/1‑61 (MQ=255)
                                            gCAATCCTGACCTTCTGTGTTCATCAGCACAAAGCGGCTGGCGTTCCACAGCTTGTTACAGAAGTTACGGt                         <  1:876502/71‑1 (MQ=255)
                                            gCAATCCTGACCTTCTGTGTTCATCAGCACAAAGCGGCTGGCGTTCCACAGCTTGTTACAGAAGTTACGGt                         <  1:3054158/71‑1 (MQ=255)
                                             cAATCCTGACCTTCTGTGTTCATCAGCACAAAGCGGCTGGCGTTCCACAGCTTGTTACAGAAGTTACGGt                         <  1:2744344/70‑1 (MQ=255)
                                                     aCCTTCTGTGTTCATCAGCACAAAGCGGCTGGCGTTCCACAGCTTGTTACAGAAGTTACGGTAACCTTCCa                <  1:1476123/71‑1 (MQ=255)
                                                     aCCTTCTGTGTTCATCAGCACAAAGCGGCTGGCGTTCCACAGCTTGTTACAGAAGTTACGGTAACCTTCCa                <  1:2506201/71‑1 (MQ=255)
                                                     aCCTTCTGTGTTCATCAGCACAAAGCGGCTGGCGTTCCACAGCTTGTTACAGAAGTTACGGTAACCTTCCa                <  1:3081431/71‑1 (MQ=255)
                                                      ccTTCTGTGTTCATCAGCACAAAGCGGCTGGCGTTCCACAGCTTGTTACAGAAGTTACGGTAACCTTCCAg               <  1:1039188/71‑1 (MQ=255)
                                                        ttCGGTGTTCATCATCACGAAGCTGGTCGCATTCCACAGCTTGTTGCAGa                                  <  1:1666247/50‑1 (MQ=255)
                                                        ttCGGTGTTCATCATCACGAAGCGGGTCGCATTCCACAGCTTGTTGCAGa                                  <  1:1853491/50‑1 (MQ=255)
                                                        ttCGGTGTTCATCATCACGAAGCGGGTCGCATTCCACAGCTTGTTGCAGa                                  <  1:1915694/50‑1 (MQ=255)
                                                        ttCGGTGTTCATCATCACGAAGCGGGTCGCATTCCACAGCTTGTTGCAGa                                  <  1:638743/50‑1 (MQ=255)
                                                          cggtgtTCATCATCACGAAGCGGGTCGCATTCCACAGCTTGTTGCAGaa                                 >  1:1847163/3‑49 (MQ=255)
                                                          cggtgtTCATCATCACGAAGCGGGTCGCATTCCACAGCTTGTTGCAGAAg                                >  1:252185/3‑50 (MQ=255)
                                                          cggtgtTCATCATCACGAAGCGGGTCGCATTCCACAGCTTGTTGCAGAAg                                >  1:3149956/3‑50 (MQ=255)
                                                          cTGTGTTCATCAGCACAAAGCGGCTGGCGTTCCACAGCTTGTTACAGAAGTTACGGTAACCTTCCAGACg            <  1:2735248/70‑1 (MQ=255)
                                                          cTGTGTTCATCAGCACAAAGCGGCTGGCGTTCCACAGCTTGTTACAGAAGTTACGGTAACCTTCCAGACGc           >  1:1226892/1‑71 (MQ=255)
                                                          cTGTGTTCATCAGCACAAAGCGGCTGGCGTTCCACAGCTTGTTACAGAAGTTACGGTAACCTTCCAGACGc           >  1:2921775/1‑71 (MQ=255)
                                                           tgtgTTCATCAGCACAAAGCGGCTGGCGTTCCACAGCTTGTTACAGAAGTTACGGTAACCTTCCa                >  1:433400/1‑65 (MQ=255)
                                                           tgtgTTCATCAGCACAAAGCGGCTGGCGTTCCACAGCTTGTTACAGAAGTTACGGTAACCTTCCa                >  1:522003/1‑65 (MQ=255)
                                                           tgtgTTCATCAGCACAAAGCGGCTGGCGTTCCACAGCTTGTTACAGAAGTTACGGTAACCTTCCa                >  1:2185184/1‑65 (MQ=255)
                                                             tgtTCATCAGCACAAAGCGGCTGGCGTTCCACAGCTTGTTACAGAAGTTACGGTAACCTTCCAGACGc           <  1:1545896/68‑1 (MQ=255)
                                                             tgtTCATCAGCACAAAGCGGCTGGCGTTCCACAGCTTGTTACAGAAGTTACGGTAACCTTCCAGACGCtt         <  1:3341670/70‑1 (MQ=255)
                                                             tgtTCATCAGCACAAAGCGGCTGGCGTTCCACAGCTTGTTACAGAAGTTACGGTAACCTTCCAGACGCtt         <  1:2107639/70‑1 (MQ=255)
                                                             tgtTCATCAGCACAAAGCGGCTGGCGTTCCACAGCTTGTTACAGAAGTTACGGTAACCTTCCAGACGCt          <  1:3015916/69‑1 (MQ=255)
                                                                     aGCACAAAGCGGCTGGCGTTACACAGCTTGTTACAGAAGTTACGGTAACCTTCCAGACGCTTCATATcc  >  1:89855/1‑69 (MQ=255)
                                                                      |                                                                   
CAGCGGTCCGCCAGCGACAGCGTCATTTCGCCGCCGTTGAAGCCGCAATCCTGACCTTCTGTGTTCATCAGCACAAAGCGGCTGGCGTTCCACAGCTTGTTACAGAAGTTACGGTAACCTTCCAGACGCTTCATATCC  >  minE/2923195‑2923332

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