Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,923,276	C→G	12.1%	S649T (AGC→ACC)	valS ←	valyl‑tRNA synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,923,276	0	C	G	12.1%	190.8 / 16.1	66	S649T (AGC→ACC)	valS	valyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base C (29/29);  new base G (3/5);  total (32/34)
Fisher's exact test for biased strand distribution p-value = 7.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.35e-01

GACAGCGTCATTTCGCCGCCGTTGAAGCCGCAATCCTGACCTTCTGTGTTCATCAGCACAAAGCGGCTGGCGTTCCACAGCTTGTTACAGAAGTTACGGTAACCTTCCAGACGCTTCATATCCCAGTTGATG  >  minE/2923210‑2923341                                                                   |                                                                  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:1846282/71‑1 (MQ=255)   cAGCGTCATTTCGCCGCCGTTGAAGCCGCAATCCTGACCTTCTGTGTTCATCAGCACAAAGCGGCTGGCGt                                                             >  1:1802274/1‑71 (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:2014190/71‑1 (MQ=255)     gCGTCATTTCGCCGCCGTTGAAGCCGCAATCCTGACCTTCTGTGTTCATCAGCACAAAGCGGCTGGCGTTc                                                           <  1:3266870/71‑1 (MQ=255)       gTCATTTCGCCGCCGTTGAAGCCGCAATCCTGACCTTCTGTGTTCATCAGCACAAAGCGGCTGGCGTTCca                                                         <  1:2976266/71‑1 (MQ=255)       gTCATTTCGCCGCCGTTGAAGCCGCAATCCTGACCTTCTGTGTTCATCAGCACAAAGCGGCTGGCGTTCca                                                         <  1:1374307/71‑1 (MQ=255)          aTTTCGCCGCCGTTGAAGCCGCAATCCTGACCTTCTGTGTTCATCAGCACAAAGCGGCTGGCGTTCCACAg                                                      >  1:1908120/1‑71 (MQ=255)           tttCGCCGCCGTTGAAGCCGCAATCCTGACCTTCTGTGTTCATCAGCACAAAGCGGCTGGCGTTCCACAGc                                                     >  1:1802991/1‑71 (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)                       tGAAGCCGCAATCCTGACCTTCTGTGTTCATCAGCACAAAGCGGc                                                                   >  1:1505174/1‑45 (MQ=255)                       tGAAGCCGCAATCCTGACCTTCTGTGTTCATCAGCACAAAGCGGc                                                                   >  1:1433483/1‑45 (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:2506201/71‑1 (MQ=255)                                       aCCTTCTGTGTTCATCAGCACAAAGCGGCTGGCGTTCCACAGCTTGTTACAGAAGTTACGGTAACCTTCCa                         <  1:3081431/71‑1 (MQ=255)                                       aCCTTCTGTGTTCATCAGCACAAAGCGGCTGGCGTTCCACAGCTTGTTACAGAAGTTACGGTAACCTTCCa                         <  1:1476123/71‑1 (MQ=255)                                        ccTTCTGTGTTCATCAGCACAAAGCGGCTGGCGTTCCACAGCTTGTTACAGAAGTTACGGTAACCTTCCAg                        <  1:1039188/71‑1 (MQ=255)                                          ttCGGTGTTCATCATCACGAAGCTGGTCGCATTCCACAGCTTGTTGCAGa                                           <  1:1666247/50‑1 (MQ=255)                                          ttCGGTGTTCATCATCACGAAGCGGGTCGCATTCCACAGCTTGTTGCAGa                                           <  1:1915694/50‑1 (MQ=255)                                          ttCGGTGTTCATCATCACGAAGCGGGTCGCATTCCACAGCTTGTTGCAGa                                           <  1:638743/50‑1 (MQ=255)                                          ttCGGTGTTCATCATCACGAAGCGGGTCGCATTCCACAGCTTGTTGCAGa                                           <  1:1853491/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:2921775/1‑71 (MQ=255)                                            cTGTGTTCATCAGCACAAAGCGGCTGGCGTTCCACAGCTTGTTACAGAAGTTACGGTAACCTTCCAGACGc                    >  1:1226892/1‑71 (MQ=255)                                             tgtgTTCATCAGCACAAAGCGGCTGGCGTTCCACAGCTTGTTACAGAAGTTACGGTAACCTTCCa                         >  1:522003/1‑65 (MQ=255)                                             tgtgTTCATCAGCACAAAGCGGCTGGCGTTCCACAGCTTGTTACAGAAGTTACGGTAACCTTCCa                         >  1:2185184/1‑65 (MQ=255)                                             tgtgTTCATCAGCACAAAGCGGCTGGCGTTCCACAGCTTGTTACAGAAGTTACGGTAACCTTCCa                         >  1:433400/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)                                                             aaGCGGCTGGCGTTCCACAGCTTGTTACAGAAGTTACGGTAACCTTCCAGACGCTTCatat             >  1:3187046/1‑61 (MQ=255)                                                             aaGCGGCTGGCGTTCCACAGCTTGTTACAGAAGTTACGGTAACCTTCCAGACGCTTCATATCCCAGTTGa    >  1:2678563/1‑70 (MQ=255)                                                             aaGCGGCTGGCGTTCCACAGCTTGTTACAGAAGTTACGGTAACCTTCCAGACGCTTCATATCCCAGTTGa    >  1:1464942/1‑70 (MQ=255)                                                             aaGCGGCTGGCGTTCCACAGCTTGTTACAGAAGTTACGGTAACCTTCCAGACGCTTCATATCCCAGTTGAt   >  1:1549479/1‑71 (MQ=255)                                                              aGCGGCTGGCGTTCCACAGCTTGTTACAGAAGTTACGGTAACCTTCCAGACGCTTCATATCCCAGTTGATg  <  1:3418372/71‑1 (MQ=255)                                                                   |                                                                  GACAGCGTCATTTCGCCGCCGTTGAAGCCGCAATCCTGACCTTCTGTGTTCATCAGCACAAAGCGGCTGGCGTTCCACAGCTTGTTACAGAAGTTACGGTAACCTTCCAGACGCTTCATATCCCAGTTGATG  >  minE/2923210‑2923341

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