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
*	minE	2,923,265	0	G	T	12.5%	212.9 / 11.3	64	L653M (CTG→ATG)	valS	valyl‑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: ‑