Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,886,319	G→A	12.3%	G191S (GGT→AGT)	ytfI →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,886,319	0	G	A	12.3%	184.3 / 15.0	65	G191S (GGT→AGT)	ytfI	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (24/33);  new base A (3/5);  total (27/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.37e-01

GGTTATGCCTTATTGGAGGATATAGTAAGCAAATTAAAACGGGGTAAATTTGAACTCCAAATACCAGGAGGTGGAATTAAAGAATGCGATGGATGGTATATTTACGAAAAAATAATTGATGATAATTTTGCCATAGTCA  >  minE/2886250‑2886388                                                                      |                                                                      ggTTATGCCTTATTGGAGGATATAGTAAGCAAATTAAAACGGGGTAAATTTGAACTCCAAATACCaggag                                                                       >  1:2470351/1‑70 (MQ=255) ggTTATGCCTTATTGGAGGATATAGTAAGCAAATTAAAACGGGGTAAATTTGAACTCCAAATACCaggag                                                                       >  1:3297673/1‑70 (MQ=255) ggTTATGCCTTATTGGAGGATATAGTAAGCAAATTAAAACGGGGTAAATTTGAACTCCAAATACCaggag                                                                       >  1:1716433/1‑70 (MQ=255)    tATGCCTTATTGGAGGATATAGTAAGCAAATTAAAACGGGGTAAATTTGAACTCCAAATACCAGGAGGTg                                                                    >  1:942348/1‑70 (MQ=255)    tATGCCTTATTGGAGGATATAGTAAGCAAATTAAAACGGGGTAAATTTGAACTCCAAATACCAGGAAGTg                                                                    >  1:1730986/1‑70 (MQ=255)        ccTTATTGGAGGATATAGTAAGCAAATTAAAACGGGGTAAATTTGAACTCCAAATACCAGGAGGTGGAAtt                                                               <  1:3426434/71‑1 (MQ=255)        ccTTATTGGAGGATATAGTAAGCAAATTAAAACGGGGTAAATTTGAACTCCAAATACCAGGAGGTGGAAtt                                                               <  1:1116427/71‑1 (MQ=255)        ccTTATTGGAGGATATAGTAAGCAAATTAAAACGGGGTAAATTTGAACTCCAAATACCAGGAGGTGGAAt                                                                >  1:2345021/1‑70 (MQ=255)         cTTATTGGAGGATATAGTAAGCAAATTAAAACGGGGTAAATTTGAACTCCAAATACCAGGAGGTGGAAtt                                                               <  1:1446391/70‑1 (MQ=255)         cTTATTGGAGGATATAGTAAGCAAATTAAAACGGGGTAAATTTGAACTCCAAATACCAGGAGGTGGAATTa                                                              <  1:2416635/71‑1 (MQ=255)            aTTGGAGGATATAGTAAGCAAATTAAAACGGGGTAAATTTGAACTCCAAATACCAGGAGGTGGAATTa                                                              >  1:2729624/1‑68 (MQ=255)            aTTGGAGGATATAGTAAGCAAATTAAAACGGGGTAAATTTGAACTCCAAATACCAGGAGGTGGAATTa                                                              >  1:703335/1‑68 (MQ=255)            aTTGGAGGATATAGTAAGCAAATTAAAACGGGGTAAATTTGAACTCCAAATACCAGGAGGTGGAATTa                                                              >  1:947748/1‑68 (MQ=255)            aTTGGAGGATATAGTAAGCAAATTAAAACGGGGTAAATTTGAACTCCAAATACCAGGAAGTGGAATTa                                                              >  1:1651484/1‑68 (MQ=255)             ttGGAGGATATAGTAAGCAAATTAAAACGGGGTAAATTTGAACTCCAAATACCAGGAGGTGGAATTAAAg                                                           <  1:3583869/70‑1 (MQ=255)             ttGGAGGATATAGTAAGCAAATTAAAACGGGGTAAATTTGAACTCCAAATACCAGGAGGTGGAATTAAAg                                                           <  1:2618500/70‑1 (MQ=255)               ggaggaTATAGTAAGCAAATTAAAACGGGGTAAATTTGAACTCCAAATACCAGGAGGTGGAATTAAAGAAt                                                        >  1:1303478/1‑71 (MQ=255)                     tataGTAAGCAAATTAAAACGGGGTAAATTTGAACTCCAAATACCAGGAGGTGGAATTAAAGAATGCgatg                                                  <  1:3273829/71‑1 (MQ=255)                       taGTAAGCAAATTAAAACGGGGTAAATTTGAACTCCAAATACCAGGAGGTGGAATTAAAGAATGCgatgga                                                >  1:1544992/1‑71 (MQ=255)                       taGTAAGCAAATTAAAACGGGGTAAATTTGAACTCCAAATACCAGGAGGTGGAATTAAAGAATGCgatgga                                                >  1:161736/1‑71 (MQ=255)                       taGTAAGCAAATTAAAACGGGGTAAATTTGAACTCCAAATACCAGGAGGTGGAATTAAAGAATGCgatgga                                                >  1:717672/1‑71 (MQ=255)                         gTAAGAAAATTAAAACGGGGTAAATTTGAACTCCAAATACCAGGAAGTGGAATTAAAGaa                                                         >  1:1204976/1‑60 (MQ=255)                           aaGCAAATTAAAACGGGGTAAATTTGAACTCCAAATACCAGGAGGTGGAATTAAAGAATGCGatggatgg                                             >  1:3188469/1‑70 (MQ=255)                           aaGCAAATTAAAACGGGGTAAATTTGAACTCCAAATACCAGGAGGTGGAATTAAAGAATGCGatggatgg                                             >  1:3533859/1‑70 (MQ=255)                            aGCAAATTAAAACGGGGTAAATTTGAACTCCAAATACCAGGAGGTGGAATTAAAGAATGCGATGGATGGta                                           <  1:3332319/71‑1 (MQ=255)                            aGCAAATTAAAACGGGGTAAATTTGAACTCCAAATACCAGGAGGTGGAATTAAAGAATGCGATGGATGGta                                           <  1:198801/71‑1 (MQ=255)                            aGCAAATTAAAACGGGGTAAATTTGAACTCCAAATACCAGGAGGTGGAATTAAAGAATGCGATGGATGGta                                           <  1:2822150/71‑1 (MQ=255)                             gCAAATTAAAACGGGGTAAATTTGAACTCCAAATACCAGGAGGTGGAATTAAAGAATGCGATGGATGGta                                           <  1:2070958/70‑1 (MQ=255)                               ataTTAAAACGGGGTAAATTTGAACTCCAAATACCAGGAGGTGGAATTAAAGAATGCGATGGATGGta                                           <  1:3170935/66‑1 (MQ=255)                               aaaTTAAAACGGGGTAAATTTGAACTCCAAATACCAGGAGGTGGAATTAAAGAATGCGATGGATGGTatat                                        >  1:1434257/1‑71 (MQ=255)                                aaTTAAAACGGGGTAAATTTGAACTCCAAATACCAGGAGGTGGAATTAAAGAATGCgatgga                                                >  1:1146853/1‑62 (MQ=255)                                aaTTAAAACGGGGTAAATTTGAACTCCAAATACCAGGAGGTGGAATTAAAGAATGCGATGGATGGTATAtt                                       <  1:2744121/71‑1 (MQ=255)                                aaTTAAAACGGGGTAAATTTGAACTCCAAATACCAGGAAGTGTAAtt                                                               <  1:714158/47‑1 (MQ=255)                                aaTTAAAACGGGGTAAATTTGAACTCCAAATACCAGGAAGTGGAAtt                                                               <  1:1304048/47‑1 (MQ=255)                                aaTTAAAACGGGGTAAATTTGAACTCCAAATACCAGGAAGTGGAAtt                                                               <  1:2395035/47‑1 (MQ=255)                                aaTTAAAACGGGGTAAATTTGAACTCCAAATACCAGGAAGTGGAAtt                                                               <  1:2151206/47‑1 (MQ=255)                                aaTTAAAACGGGGTAAATTTGAACTCCAAATACCAGGAAGTGGAAtt                                                               <  1:2071199/47‑1 (MQ=255)                                 aTTAAAACGGGGTAAATTTGAACTCCAAATACCaggag                                                                       >  1:1693476/1‑38 (MQ=255)                                 aTTAAAACGGGGTAAATTTGAACTCCAAATACCAGGAGGTGGAATTAAAGAATGCGATGGATGGTATAtt                                       >  1:182573/1‑70 (MQ=255)                                  ttAAAACGGGGTAAATTTGAACTCCAAATACCAGGAGGTGGAATTa                                                              <  1:3567926/46‑1 (MQ=255)                                  ttAAAACGGGGTAAATTTGAACTCCAAATACCAGGAGGTGGAATTAAAGAATGCGATGGATGGTATATTTa                                     <  1:703949/71‑1 (MQ=255)                                           ggTAAATTTGAACTCCAAATACCAGGAGGTGGAATTAAAGAATGCGATGGATGGTATATTTACGAAAAaat                            <  1:962577/71‑1 (MQ=255)                                            gTAAATTTGAACTCCAAATACCAGGAGGTGGAATTAAAGAATGCGATGGATGGTATATTTACGAAAAaata                           <  1:1324773/71‑1 (MQ=255)                                             tAAATTTGAACTCCAAATACCAGGAGGTGGAATTAAAGAATGCGATGGATGGTATATTTACGaaaa                               <  1:3121765/66‑1 (MQ=255)                                             tAAATTTGAACTCCAAATACCAGGAGGTGGAATTAAAGAATGCGATGGATGGTATATTTACGaaaa                               <  1:74136/66‑1 (MQ=255)                                               aaTTTGAACTCCAAATACCAGGAGGTGGAATTAAAGAATGCGGTGGATGGTATATTTACGAAAAAATAAtt                        >  1:3497185/1‑71 (MQ=255)                                               aaTTTGAACTCCAAATACCAGGAGGTGGAATTAAAGAATGCGATGGATGGtat                                          <  1:2705493/53‑1 (MQ=255)                                               aaTTTGAACTCCAAATACCAGGAGGTGGAATTAAAGAATGCGATGGATGGTATATTTACGAAAAAATAAtt                        >  1:1062404/1‑71 (MQ=255)                                                   tGAACTCCAAATACCAGGAGGTGGAATTAAAGAATGCGATGGATGGt                                            >  1:895614/1‑47 (MQ=255)                                                      aCTCCAAATACCAGGAGGTGGAATTAAAGAATGCGATGGATGGTATATTTACGAAAAAATAATTGATGATa                 <  1:1072101/71‑1 (MQ=255)                                                      aCTCCAAATACCAGGAGGTGGAATTAAAGAATGCGATGGATGGTATATTTACGAAAAAATAATTGATGATa                 <  1:402353/71‑1 (MQ=255)                                                       cTCCAAATACCAGGAGGTGGAATTAAAGAATGCGATGGATGGTATATTTACGAAAAAATAATTGATGATaa                <  1:3607305/71‑1 (MQ=255)                                                       cTCCAAATACCAGGAGGTGGAATTAAAGAATGCGATGGATGGTATATTTACGAAAAAATAATTGATGATaa                <  1:696734/71‑1 (MQ=255)                                                       cTCCAAATACCAGGAGGTGGAATTAAAGAATGCGATGGATGGTATATTTACGAAAAAATAATTGATGATaa                >  1:3609602/1‑71 (MQ=255)                                                           aaaTACCAGGAGGTGGAATTAAAGAATGCGATGGATGGTATATTTACGAAAAAATAATTGATGATAAt               <  1:1058714/68‑1 (MQ=255)                                                            aaTACCAGGAGGTGGAATTAAAGAATGCGATGGATGGTATATTTACGaaaaa                              >  1:2459668/1‑52 (MQ=255)                                                            aaTACCAGGAGGTGGAATTAAAGAATGCGATGGATGGTATATTTACGaaaaa                              >  1:3336132/1‑52 (MQ=255)                                                            aaTACCAGGAGGTGGAATTAAAGAATGCGATGGATGGTATATTTACGaaaaa                              >  1:337484/1‑52 (MQ=255)                                                            aaTACCAGGAGGTGGAATTAAAGAATGCGATGGATGGTATATTTACGAAAAaat                            <  1:3467906/54‑1 (MQ=255)                                                            aaTACCAGGAGGTGGAATTAAAGAATGCGATGGATGGTATATTTACGAAAAaat                            <  1:1365481/54‑1 (MQ=255)                                                            aaTACCAGGAGGTGGAATTAAAGAATGCGATGGATGGTATATTTACGAAAAaat                            <  1:2064038/54‑1 (MQ=255)                                                            aaTACCAGGAGGTGGAATTAAAGAATGCGATGGATGGTATATTTACGAAAAAATAATTGATGATAATTTTg           <  1:2321484/71‑1 (MQ=255)                                                              tACCAGGAGGTGGAATTAAAGAATGCGATGGATGGTATATTTACGaaaaa                              <  1:3001077/50‑1 (MQ=255)                                                              tACCAGGAGGTGGAATTAAAGAATGCGATGGATGGTATATTTACGaaaaa                              <  1:2874894/50‑1 (MQ=255)                                                                ccAGGAGGTGGAATTAAAGAATGCGATGGATGGTATATTTACGAAAAAATAATtga                      >  1:3598347/1‑56 (MQ=255)                                                                ccAGGAGGTGGAATTAAAGAATGCGATGGATGGTATATTTACGAAAAAATAATtga                      >  1:1105034/1‑56 (MQ=255)                                                                ccAGGAGGTGGAAATAAAGAATGCGATGGATGGTATATTTACGAAAAAATAAttc                       >  1:1146568/1‑54 (MQ=255)                                                                 cAGGAGGTGGAATTAAAGAATGCGATGGATGGTATATTTACGAAAAAATAATTGATGATAATTTTGCCATa      <  1:3823299/71‑1 (MQ=255)                                                                 cAGGAGGTGGAATTAAAGAATGCGATGGATGGTATATTTACGAAAAAATAATTGATGATAATTTTGCCATa      <  1:2714514/71‑1 (MQ=255)                                                                  aggaggTGGAATTAAAGAATGCGATGGATGGTATAtt                                       <  1:1114005/37‑1 (MQ=255)                                                                    gaggTGGAATTAAAGAATGCGATGGATGGTATATTTa                                     >  1:2547426/1‑37 (MQ=255)                                                                    gaggTGGAATTAAAGAATGCGATGGATGGTATATTTa                                     >  1:256734/1‑37 (MQ=255)                                                                     aggTGGAATTAAAGAATGCGATGGATGGTATATTTACGAAAAAATAATTGATGATa                 <  1:3245663/56‑1 (MQ=255)                                                                     aggTGGAATTAAAGAATGCGATGGATGGTATATTTACGAAAAAATAATTGATGATAATTTTGCCATAGTCa  <  1:2639098/71‑1 (MQ=255)                                                                      ggTGGAATTAAAGAATGCGATGGATGGTATATTTACGAAAAAATAATTGATGATAATTTTGCCATAGTCa  >  1:3538236/1‑70 (MQ=255)                                                                      |                                                                      GGTTATGCCTTATTGGAGGATATAGTAAGCAAATTAAAACGGGGTAAATTTGAACTCCAAATACCAGGAGGTGGAATTAAAGAATGCGATGGATGGTATATTTACGAAAAAATAATTGATGATAATTTTGCCATAGTCA  >  minE/2886250‑2886388

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