Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,923,269	A→G	11.9%	F651F (TTT→TTC)	valS ←	valyl‑tRNA synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,923,269	0	A	G	11.9%	165.6 / 12.6	59	F651F (TTT→TTC)	valS	valyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base A (26/26);  new base G (3/4);  total (29/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.49e-01

GCCAGCGACAGCGTCATTTCGCCGCCGTTGAAGCCGCAATCCTGACCTTCTGTGTTCATCAGCACAAAGCGGCTGGCGTTCCACAGCTTGTTACAGAAGTTACGGTAACCTTCCAGACGCTTCATATCCCAGTTGA  >  minE/2923204‑2923339                                                                  |                                                                       gCCAGCGACAGCGTCATTTCGCCGCCGTTGAAGCCGCAATCCTGACCTTCTGTGTTCATCAGCACAAAGCg                                                                   <  1:462274/71‑1 (MQ=255)       gACAGCGTCATTTCGCCGCCGTTGAAGCCGCAATCCTGACCTTCTGTGTTCATCAGCACAAAGCGGCTGGc                                                             >  1:1113180/1‑71 (MQ=255)        aCAGCGTCATTTCGCCGCCGTTGAAGCCGCAATCCTGACCTTCTGTGTTCATCAGCACAAAGCGGCTGGCg                                                            >  1:2950589/1‑71 (MQ=255)        aCAGCGTCATTTCGCCGCCGTTGAAGCCGCAATCCTGACCTTCTGTGTTCATCAGCACAAAGCGGCTGGCg                                                            >  1:2412348/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)                  ttCGCCGCCGTTGAAGCCGCAATCCTGACCTTCTGTGTTCATCAGCACaaa                                                                      >  1:2231413/1‑51 (MQ=255)                  ttCGCCGCCGTTGAAGCCGCAATCCTGACCTTCTGTGTTCATCAGCACAAAGCGGCTGGc                                                             >  1:2736982/1‑60 (MQ=255)                  ttCGCCGCCGTTGAAGCCGCAATCCTGACCTTCTGTGTTCATCAGCACAAAGCGGCTGGc                                                             >  1:3016265/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:1141252/1‑69 (MQ=255)                         ccgTTGAAGCCGCAATCCTGACCTTCTGTGTTCATCAGCACAAAGCGGCTGGCGTTCCACAGCTTGTTa                                             >  1:1063362/1‑69 (MQ=255)                          cgTTGAAGCCGCAATCCTGACCTTCTGTGTTCATCAGCACAAAGCGGCTGGCGTTCCACAGCTTGTTa                                             >  1:605590/1‑68 (MQ=255)                            ttGAAGCCGCAATCCTGACCTTCTGTGTTCATCAGcaca                                                                        <  1:3143350/39‑1 (MQ=255)                             tGAAGCCGCAATCCTGACCTTCTGTGTTCATCAGCACAAAGCGGc                                                                 >  1:1433483/1‑45 (MQ=255)                             tGAAGCCGCAATCCTGACCTTCTGTGTTCATCAGCACAAAGCGGc                                                                 >  1:1505174/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:3054158/71‑1 (MQ=255)                                    gCAATCCTGACCTTCTGTGTTCATCAGCACAAAGCGGCTGGCGTTCCACAGCTTGTTACAGAAGTTACGGt                                <  1:876502/71‑1 (MQ=255)                                     cAATCCTGACCTTCTGTGTTCATCAGCACAAAGCGGCTGGCGTTCCACAGCTTGTTACAGAAGTTACGGt                                <  1:2744344/70‑1 (MQ=255)                                             aCCTTCTGTGTTCATCAGCACAAAGCGGCTGGCGTTCCACAGCTTGTTACAGAAGTTACGGTAACCTTCCa                       <  1:1476123/71‑1 (MQ=255)                                             aCCTTCTGTGTTCATCAGCACAAAGCGGCTGGCGTTCCACAGCTTGTTACAGAAGTTACGGTAACCTTCCa                       <  1:3081431/71‑1 (MQ=255)                                             aCCTTCTGTGTTCATCAGCACAAAGCGGCTGGCGTTCCACAGCTTGTTACAGAAGTTACGGTAACCTTCCa                       <  1:2506201/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:638743/50‑1 (MQ=255)                                                ttCGGTGTTCATCATCACGAAGCGGGTCGCATTCCACAGCTTGTTGCAGa                                         <  1:1915694/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:522003/1‑65 (MQ=255)                                                   tgtgTTCATCAGCACAAAGCGGCTGGCGTTCCACAGCTTGTTACAGAAGTTACGGTAACCTTCCa                       >  1:433400/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)                                                               cacgAAGCGGGTCGCATTCCACAGCTTGTTGCAGa                                         <  1:3458943/35‑1 (MQ=16)                                                                acaAAGCGGCTGGCGTTCCACAGCTTGTTACAGAAGTTACGGTAACCTTCCAGACGCTTCATATCCCAGtt    <  1:2736336/71‑1 (MQ=255)                                                                acaAAGCGGCTGGCGTTCCACAGCTTGTTACAGAAGTTACGGTAACCTTCCAGACGCTTCATATCCCAGtt    <  1:1641966/71‑1 (MQ=255)                                                                 caAAGCGGCTGGCGTTCCACAGCTTGTTACAGAAGTTACGGTa                               <  1:217556/43‑1 (MQ=255)                                                                  aaaGCGGCTGGCGTTCCACAGCTTGTTACAGAAGTTACGGTAACCTTCCAGACGCTTCATATCCCAGTTGa  >  1:2555996/1‑71 (MQ=255)                                                                  |                                                                       GCCAGCGACAGCGTCATTTCGCCGCCGTTGAAGCCGCAATCCTGACCTTCTGTGTTCATCAGCACAAAGCGGCTGGCGTTCCACAGCTTGTTACAGAAGTTACGGTAACCTTCCAGACGCTTCATATCCCAGTTGA  >  minE/2923204‑2923339

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: ‑