Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,924,237	G→A	18.6%	P329S (CCT→TCT)	valS ←	valyl‑tRNA synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,924,237	0	G	A	18.6%	129.2 / 22.0	59	P329S (CCT→TCT)	valS	valyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base G (26/22);  new base A (6/5);  total (32/27)
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.83e-01

AGCGCGTCAACTGCGGCAACGACTGCTTTACGTGCAGCAAAACGCTCCAGTTTCTGGAACTCTGCAGGGATTTCGCTGGAATAAACGTCAGATTCGTTACCTTTGGTATCGAACACCTGGGCGCTTTCACGGATATCG  >  minE/2924170‑2924307                                                                    |                                                                       aGCGCGTCAACTGCGGCAACGACTGCTTTACGTGCAGCAAAACGCTCCAGTTTCTGGAACTCTGCAGGGAt                                                                     >  1:1097608/1‑71 (MQ=255) aGCGCGTCAACTGCGGCAACGACTGCTTTACGTGCAGCAAAACGCTCCAGTTTCTGGAACTCTGCAGGGAt                                                                     >  1:326600/1‑71 (MQ=255)     cgTCAACTGCGGCAACGACTGCTTTACGTGCAGCAAAACGCTCCAGTTTCTGGAACTCTGCAGGGATTTCg                                                                 <  1:1031641/71‑1 (MQ=255)     cgTCAACTGCGGCAACGACTGCTTTACGTGCAGCAAAACGCTCCAGTTTCTGGAACTCTGCAGAGATTTCg                                                                 <  1:3793258/71‑1 (MQ=255)      gTCAACTGCGGCAACGACTGCTTTACGTGCAGCAAAACGCTCCAGTTTCTGGAACTCTGCAGGGATTTCGc                                                                <  1:3424835/71‑1 (MQ=255)        cAACTGCGGCAACGACTGCTTTACGTGCAGCAAAACGCTCCAGTTTCTGGAACTCTGCAGGGATTTCGCTg                                                              >  1:2926241/1‑71 (MQ=255)             gCGGCAACGACTGCTTTACGTGCAGCAAAACGCTCCAGTTTCTGGAACTCTGCAGGGATTTCGCTGGAAt                                                          >  1:196367/1‑70 (MQ=255)             gCGGCAACGACTGCTTTACGTGCAGCAAAACGCTCCAGTTTCTGGAACTCTGCAGAGATTTCGCTGGAAt                                                          >  1:149101/1‑70 (MQ=255)              cGGCAACGACTGCTTTACGTGCAGCAAAACGCTCCAGTTTCTGGAACTCTGCAGGGATTTCGCTGGAATa                                                         <  1:2294561/70‑1 (MQ=255)                gCAACGACTGCTTTACGTGCAGCAAAACGCTCCAGTTTCTGGAACTCTGCAGAGAtt                                                                    <  1:3653818/57‑1 (MQ=255)                 cAACGACTGCTTTACGTGCAGCAAAACGCTCCAGTTTCTGGAACTCTGCAGGGATTTCg                                                                 <  1:1923119/59‑1 (MQ=255)                     gACTGCTTTACGTGCAGCAAAACGCTCCAGTTTCTGGAACTCTGCAGGGATTTCGCTGGAATAAACGt                                                    >  1:2851970/1‑68 (MQ=255)                         gCTTTACGTGCAGCAAAACGCTCCAGTTTCTGGAACTCTGCAGGGATTTCGCTGGAATAAACGTCAGATTc                                             >  1:3561048/1‑71 (MQ=255)                          cTTTACGTGCAGCAAAACGCTCCAGTTTCTGGAACTCTGCAGGGATTTCGCTGGAATAAACGTCAGATTCg                                            >  1:142647/1‑71 (MQ=255)                             tACGTGCAGCAAAACGCTCCAGTTTCTGGAACTCTGCAGAGATTTCGCTGGAATAAACGTCAGATTCGTTa                                         >  1:2471548/1‑71 (MQ=255)                              aCGTGCAGCAAAACGCTCCAGTTTCTGGAACTCTGCAGGGATTTCGCTGGAATAAACGTCAGATTCGTTAc                                        >  1:2131386/1‑71 (MQ=255)                               cGTGCAGCAAAACGCTCCAGTTTCTGGAACTCTGCAGGGATTTCGCTGGAATAAACGTCAGATTc                                             >  1:457189/1‑65 (MQ=255)                                gTGCAGCAAAACGCTCCAGTTTCTGGAACTCTGCAGGGATTTCg                                                                 >  1:3236681/1‑44 (MQ=255)                                 tGCAGCAAAACGCTCCAGTTTCTGGAACTCTGCAGGGATTTCGCTGGAGTAAACGTCAGATTCGTTACCtt                                     <  1:3596839/71‑1 (MQ=255)                                 tGCAGCAAAACGCTCCAGTTTCTGGAACTCTGCAGGGATTTCGCTGGAATAAACGTCAGATTCGTTACCtt                                     <  1:2253778/71‑1 (MQ=255)                                   cagcaAAACGCTCCAGTTTCTGGAACTCTGCAGGGATTTCGCTGGAATAAACGTCAGATTCGTTACCttt                                    <  1:2065891/70‑1 (MQ=255)                                    agcaAAACGCTCCAGTTTCTGGAACTCTGCAGGGATTTCGCTGGAATAAACGTCAGATTCGTTACCTTTgg                                  <  1:420238/71‑1 (MQ=255)                                      caAAACGCTCCAGTTTCTGGAACTCTGCAGGGAttt                                                                   <  1:906569/36‑1 (MQ=255)                                      caAAACGCTCCAGTTTCTGGAACTCTGCAGGGATTTCGCTGGAATAAACGTCAGATTCGTTACCTTTGGTa                                >  1:3412989/1‑71 (MQ=255)                                      caAAACGCTCCAGTTTCTGGAACTCTGCAGGGATTTCGCTGGAATAAACGTCAGATTCGTTACCTTTGGTa                                >  1:3408395/1‑71 (MQ=255)                                      caAAACGCTCCAGTTTCTGGAACTCTGCAGAGATTTCGCTGGAATAAACGTCAGATTCGTTACCTTTGGTa                                <  1:3495328/71‑1 (MQ=255)                                        aaaCGCTCCAGTTTCTGGAACTCTGCAGGGATTTCGCTGGAATAAACGTCAGAt                                               >  1:3719030/1‑54 (MQ=255)                                        aaaCGCTCCAGTTTCTGGAACTCTGCAGGGATTTCGCTGGAATAAACGTCAGATTCGTTACCTTTg                                   >  1:3442886/1‑66 (MQ=255)                                        aaaCGCTCCAGTTTCTGGAACTCTGCAGGGATTTCGCTGGAATAAACGTCAGATTCGTTACCTTTGGTATc                              <  1:1513546/71‑1 (MQ=255)                                         aaCGCTCCAGTTTCTGGAACTCTGCAGGGATTTCGCTGGAATAAACGTCAGATTCGTTACCTTTGGTATCg                             <  1:109023/71‑1 (MQ=255)                                         aaCGCTCCAGTTTCTGGAACTCTGCAGGGATTTCGCTGGAATAAACGTCAGATTCGTTACCTTTGGTATCg                             <  1:3061/71‑1 (MQ=255)                                          aCGCTCCAGTTTCTGGAACTCTGCAGGGATTTCGCTGGAATAAACGTCAGATTCGTTACCTTTGGTATc                              >  1:3095923/1‑69 (MQ=255)                                            gCTCCAGTTTCTGGAACTCTGCAGGGATTTCGCTGGAATAAACGTCAGATTCGTTACCTTTGGTATCGaa                           <  1:2702647/70‑1 (MQ=255)                                              tCCAGTATCTGGAACTCTGCAGGGATTTCGCTGGAATAAACGTCAGATTCGTTACCTTTGGTATCGAacac                        >  1:510722/1‑71 (MQ=255)                                                cAGTTTCTGGAACTCTGCAGGGATTTCGCTGGAATAAACGTCAGATTCGTTACCTTTGGTATCGAACACCt                      <  1:2233277/71‑1 (MQ=255)                                                cAGTTTCTGGAACTCTGCAGGGATTTCGCTGGAATAAACGTCAGATTCGTTACCTTTGGTATCGAACACCt                      <  1:1853945/71‑1 (MQ=255)                                                 aGTTTCTGGAACTCTGCAGAGATTTCGCTGGAATAAACGTCAGATTCGTTACCttt                                    >  1:3898755/1‑56 (MQ=255)                                                 aGTTTCTGGAACTCTGCAGAGATTTCGCTGGAATAAACGTCAGATTCGTTACCTTTGGTATCGAACACCTg                     >  1:2175277/1‑71 (MQ=255)                                                  gTTTCTGGAACTCTGCAGGGATTTCGCTGGAATAAACGTCAGATTCGTTACCt                                      >  1:1135702/1‑53 (MQ=255)                                                  gTTTCTGGAACTCTGCAGGGATTTCGCTGGAATAAACGTCAGATTCGTTACCTTTGGTATCGAACACCTgg                    >  1:2021009/1‑71 (MQ=255)                                                    ttCTGGAACTCTGCAGGGATTTCGCTGGAATAAACGTCAGATTCGTTACCTTTGGTATCGAACACCTGGgc                  >  1:2950928/1‑71 (MQ=255)                                                    ttCTGGAACTCTGCAGAGATTTCGCTGGAATAAACGTCAGATTCGTTACCTTTGGTATCGAACACCTGGgc                  >  1:395186/1‑71 (MQ=255)                                                     tCTGGAACTCTGCAGAGATTTCGCTGGAATAAACGTCAGATTCGTTACCTTTGGTATCGAACACCTGGgcg                 >  1:319112/1‑71 (MQ=255)                                                      cTGGAACTCTGCAGAGATTTCGCTGGAATAAACGTCAGATTCGTTACCTTTGGTATCGAACACCTGGgcg                 <  1:2997702/70‑1 (MQ=255)                                                       tGGAACTCTGCAGGGATTTCGCTGGAATAAACGTCAGATTCGTTACCTTTGGTa                                <  1:2430159/54‑1 (MQ=255)                                                        ggAACTCTGCAGAGATTTCGCTGGAATAAACGTCAGATTCGTTACCTTTGGTa                                <  1:1116636/53‑1 (MQ=255)                                                         gAACTCTGCAGGGATTTCGCTGGAATAAACGTCAGATTCGTTACCTTTGGTATCGAACACCTGGGCGCttt             >  1:1285058/1‑71 (MQ=255)                                                          aaCTCTGCAGGGATTTCGCTGGAATAAACGTCAGATTCGTTACCTTTGGTATCGAACACCTGGGCGCttt             >  1:1381433/1‑70 (MQ=255)                                                           aCTCTGCAGGGATTTCGCTGGAATAAACGTCAGATTCGTTACCTTTGGTATCGAACACCTGGGCGCTTTCa           >  1:3932799/1‑71 (MQ=255)                                                           aCTCTGCAGGGATTTCGCTGGAATAAACGTCAGATTCGTTACCTTTGGTATCGAACACCTGGGCGCTTTCa           <  1:1134554/71‑1 (MQ=255)                                                           aCTCTGCAGGGATTTCGCTGGAATAAACGTCAGATTCGTTACCTTTGGTATCGAACACCTGGGCGCTTTCa           >  1:829549/1‑71 (MQ=255)                                                            ctctGCAGGGATTTCGCTGGAATAAACGTCAGATTCGTTACCttt                                    >  1:366371/1‑45 (MQ=255)                                                               tGCAGGGATTTCGCTGGAATAAACGTCAGATTCGTTACCTTTGGTATc                              <  1:1992276/48‑1 (MQ=255)                                                               tGCAGGGATTTCGCTGGAATAAACGTCAGATTCGTTACCTTTGGTATCGAACACCt                      <  1:679125/56‑1 (MQ=255)                                                                gCAGGGATTTCGCTGGAATAAACGTCAGATTCGTTACCTTTGGTATc                              <  1:2156960/47‑1 (MQ=255)                                                                 cAGGGATTTCGCTGGAATAAACGTCAGATTCGTTACCTTTGGTATCGAACACCTGGGCGCTTTCACGGata     <  1:2346652/71‑1 (MQ=255)                                                                 cAGGGATTTCGCTGGAATAAACGTCAGATTCGTTACCTTTGGTATCGAACACCTGGGCGCTTTCACGGata     >  1:3472480/1‑71 (MQ=255)                                                                 cAGGGATTTCGCTGGAATAAACGTCAGATTCGTTACCTTTGGTATCGAACACCTGGGCGCTTTCACGGata     >  1:276337/1‑71 (MQ=255)                                                                  aGGGATTTCGCTGGAATAAACGTCAGATTCGTTACCTTTGGTAt                               <  1:3853789/44‑1 (MQ=255)                                                                    ggATTTCGCTGGAATAAACGTCAGATTCGTTACCTTTGGTATCGAACACCTGGGCGCTTTCACGGATATCg  >  1:1711821/1‑71 (MQ=255)                                                                    |                                                                       AGCGCGTCAACTGCGGCAACGACTGCTTTACGTGCAGCAAAACGCTCCAGTTTCTGGAACTCTGCAGGGATTTCGCTGGAATAAACGTCAGATTCGTTACCTTTGGTATCGAACACCTGGGCGCTTTCACGGATATCG  >  minE/2924170‑2924307

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