Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,835,651	G→A	22.7%	N82N (AAC→AAT)	yjeJ ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,835,651	0	G	A	22.7%	125.9 / 27.0	92	N82N (AAC→AAT)	yjeJ	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (31/40);  new base A (8/13);  total (39/53)
Fisher's exact test for biased strand distribution p-value = 8.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

TCCAGCGTGTTAACCCGACGATTAATATCGGCATTTTTAAGCTCATCAACCAGAATTTCTGGGATGTTTTCTGCCATTTTCTTAATGACTTTATCGCGAGCTTGCTCATACTGAAGGTG  >  minE/2835586‑2835704                                                                  |                                                      tCCAGCGTGTTAACCCGACGATTAATATCGGCATTTTTAAGCTCATCAACCAGAATTTCTGGGATGttt                                                    >  1:2672414/1‑69 (MQ=255)        tgtTAAGCCGACGATTAATATCGGCATTTTTAAGCTCATCAACCAGAATTTCTGGGATGTTTTCTGCCa                                             >  1:43160/1‑69 (MQ=255)        tgtTAACCCGACGATTAATATCGGCATTTTTAAGCTCATCAACCAGAATTTCTGGGATGTTTTCTGCCa                                             >  1:2421653/1‑69 (MQ=255)        tgtTAACCCGACGATTAATATCGGCATTTTTAAGCTCATCAACCAGAATTTCTGGGATGTTTTCTGCCa                                             >  1:357516/1‑69 (MQ=255)          ttAACCCGACGATTAATATCGGCATTTTTAAGCTCATCAACCAGAATTTCTGGGATGTTTTCTGCCAtt                                           >  1:2850976/1‑69 (MQ=255)          ttAACCCGACGATTAATATCGGCATTTTTAAGCTCATCAACCAGAATTTCTGGGATGTTTTCTGCCAtt                                           >  1:2403235/1‑69 (MQ=255)          ttAACCCGACGATTAATATCGGCATTTTTAAGCTCATCAACCAGAATTTCTGGGATGTTTTCTGCCAtt                                           <  1:1656586/69‑1 (MQ=255)          ttAACCCGACGATTAATATCGGCATTTTTAAGCTCATCAACCAGAATTTCTGGGATGTTTTCTGCCAtt                                           >  1:144532/1‑69 (MQ=255)          ttAACCCGACGATTAATATCGGCATTTTTAAGCTCATCAACCAGAATTTCTGGGATGTTTTCTGCCAtt                                           <  1:298252/69‑1 (MQ=255)          ttAACCCGACGATTAATATCGGCATTTTTAAGCTCATCAACCAGAATTTCTGGGATGTTTTCTGCCAt                                            >  1:385157/1‑68 (MQ=255)          ttAACCCGACGATTAATATCGGCATTTTTAAGCTCATCAACCAGAATTTCTGGGATGTTTTCTGCCAt                                            >  1:2972529/1‑68 (MQ=255)          ttAACCCGACGATTAATATCGGCATTTTTAAGCTCATCAACCAGAATTTCTGGGATGTTTTCTGCCAt                                            >  1:832826/1‑68 (MQ=255)               ccGACGATTAATATCGGCATTTTTAAGCTCATCAACCAGAATTTCTGGGATGTTTTCTGCCATTTTCtt                                      <  1:2321629/69‑1 (MQ=255)               ccGACGATTAATATCGGCATTTTTAAGCTCATCAACCAGAATTTCTGGGATGTTTTCTGCCATTTTCtt                                      <  1:142251/69‑1 (MQ=255)               ccGACGATTAATATCGGCATTTTTAAGCTCATCAACCAGAATTTCTGGGATGTTTTCTGCCATTTTCtt                                      <  1:1483120/69‑1 (MQ=255)               ccGACGATTAATATCGGCATTTTTAAGCTCATCAACCAGAATTTCTGGGATGTTTTCTGCCATTTTCtt                                      <  1:2341772/69‑1 (MQ=255)               ccGACGATTAATATCGGCATTTTTAAGCTCATCAACCAGAATTTCTGGGATGTTTTCTGCCATTTTCtt                                      <  1:3365575/69‑1 (MQ=255)               ccGACGATTAATATCGGCATTTTTAAGCTCATCAACCAGAATTTCTGGGATGTTTTCTGCCATTTTCtt                                      <  1:1208510/69‑1 (MQ=255)               ccGACGATTAATATCGGCATTTTTAAGCTCATCAACCAGAATTTCTGGGATGTTTTCTGCCATTTTCtt                                      <  1:2703790/69‑1 (MQ=255)               ccGACGATTAATATCGGCATTTTTAAGCTCATCAACCAGAATTTCTGGGATGTTTTCTGCCATTTTCtt                                      <  1:200606/69‑1 (MQ=255)                     aTTAATATCGGCATTTTTAAGCTCATCAACCAGACTTTCTGGGATGTTTTCTGCCATTTTCTTAATGAc                                <  1:3234700/69‑1 (MQ=255)                     aTTAATATCGGCATTTTTAAGCTCATCAACCAGAATTTCTGGGATGTTTTCTGCCATTTTCTTAATGAc                                <  1:2307319/69‑1 (MQ=255)                       tAATATCGGCATTTTTAAGCTCATCAACCAGAATTTCTGGGATGTTTTCTGCCATTTTCTTAATGACtt                              >  1:1789476/1‑69 (MQ=255)                       tAATATCGGCATTTTTAAGCTCATCAACCAGAATTTCTGGGATGTTTTCTGCCATTTTCTTAATGACtt                              >  1:228767/1‑69 (MQ=255)                       tAATATCGGCATTTTTAAGCTCATCAACCAGAATTTCTGGGATGTTTTCTGCCATTTTCTTAATGACtt                              >  1:2371510/1‑69 (MQ=255)                       tAATATCGGCATTTTTAAGCTCATCAACCAGAATTTCTGGGATGTTTTCTGCCATTTTCTTAATGACtt                              >  1:484164/1‑69 (MQ=255)                       tAATATCGGCATTTTTAAGCTCATCAACCAGAATTTCTGGGATGTTTTCTGCCATTTTCTTAATGACtt                              >  1:2890901/1‑69 (MQ=255)                            tCGGCATTTTTAAGCTCATCAACCAGAATTTCTGGGATGTTTTCTGCCATTTTCTTAATGACTTTATc                          >  1:2021458/1‑68 (MQ=255)                            tCGGCATTTTTAAGCTCATCAACCAGAATTTCTGGGATGTTTTCTGCCATTTTCTTAATGACTTTATc                          >  1:285813/1‑68 (MQ=255)                            tCGGCATTTTTAAGCTCATCAACCAGAATTTCTGGGATGTTTTCTGCCATTTTCTTAATGACTTTATc                          >  1:1545291/1‑68 (MQ=255)                                cATTTTTAAGCTCATCAACCAGAATTTCTGGGATGTTTTCTGCCATTTTCTTAATGACTTTATCGCGAg                     <  1:508822/69‑1 (MQ=255)                                 aTTTTTAAGCTCATGAACCAGAATTTCTGGGATGTTTTCTGCCAtt                                           <  1:3051006/46‑1 (MQ=255)                                 aTTTTTAAGCTCATCAACCAGAATTTCTGGGATATTTTCTg                                                >  1:1569559/1‑41 (MQ=255)                                 aTTTTTAAGCTCATCAACCAGAATTTCTGGGATATTTTCTGCCATTTTCTTAATGACTTTATCGCGAGc                    >  1:2365184/1‑69 (MQ=255)                                 aTTTTTAAGCTCATCAACCAGAATTTCTGGGATATTTTCTGCCATTTTCTTAATGACTTTATCGCGAGc                    >  1:3387066/1‑69 (MQ=255)                                 aTTTTTAAGCTCATCAACCAGAATTTCTGGGATATTTTCTGCCATTTTCTTAATGACTTTATCGCGAGc                    >  1:1901480/1‑69 (MQ=255)                                 aTTTTTAAGCTCATCAACCAGAATTTCTGGGATATTTTCTGCCATTTTCTTAATGACTTTATCGCGAGc                    >  1:341710/1‑69 (MQ=255)                                 aTTTTTAAGCTCATCAACCAGAATTTCTGGGATATTTTCTGCCATTTTCTTAATGACTTTATCGCGAGc                    >  1:1366767/1‑69 (MQ=255)                                 aTTTTTAAGCTCATCAACCAGAATTTCTGGGATATTTTCTGCCATTTTCTTAATGACTTTATCGCGAGc                    >  1:53663/1‑69 (MQ=255)                                 aTTTTTAAGCTCATCAACCAGAATTTCTGGGATATTTTCTGCCATTTTCTTAATGACTTTATCGCGAGc                    >  1:1182404/1‑69 (MQ=255)                                  tttttAAGCTCATCAACCAGAATTTCTGGGATGTTTTCTGCCATTTTCTTAATGACTTTATCGCGAGCt                   <  1:2748227/69‑1 (MQ=255)                                    tttAAGCTCATCAACCAGAATTTCTGGGATGTTTTCTGCCATTTTCTTAATGACTTTa                            <  1:3104024/58‑1 (MQ=255)                                    tttAAGCTCATCAACCAGAATTTCTGGGATGTTTTCTGCCATTTTCTTAATGACTTTa                            <  1:245318/58‑1 (MQ=255)                                    tttAAGCTCATCAACCAGAATTTCTGGGATGTTTTCTGCCATTTTCTTAATGACTTTa                            <  1:592017/58‑1 (MQ=255)                                    tttAAGCTCATCAACCAGAATTTCTGGGATGTTTTCTGCCATTTTCTTAATGACTTTa                            <  1:2938744/58‑1 (MQ=255)                                    tttAAGCTCATCAACCAGAATTTCTGGGATGTTTTCTGCCATTTTCTTAATGACTTTATcgcg                       >  1:218738/1‑63 (MQ=255)                                    tttAAGCTCATCAACCAGAATTTCTGGGATGTTTTCTGCCATTTTCTTAATGACTTTATCGCGAGCTTg                 >  1:336702/1‑69 (MQ=255)                                    tttAAGCTCATCAACCAGAATTTCTGGGATGTTTTCTGCCATTTTCTTAATGACTTTATCGCGAGCTTg                 >  1:1022593/1‑69 (MQ=255)                                    tttAAGCTCATCAACCAGAATTTCTGGGATGTTTTCTGCCATTTTCTTAATGACTTTATCGCGAGCTTg                 >  1:462575/1‑69 (MQ=255)                                    tttAAGCTCATCAACCAGAATTTCTGGGATGTTTTCTGCCATTTTCTTAATGACTTTATCGCGAGCTTg                 >  1:597492/1‑69 (MQ=255)                                    tttAAGCTCATCAACCAGAATTTCTGGGATGTTTTCTGCCATTTTCTTAATGACTTTATCGCGAGCTTg                 >  1:869678/1‑69 (MQ=255)                                    tttAAGCTCATCAACCAGAATTTCTGGGATGTTTTCTGCCATTTTCTTAATGACTTTATCGCGAGCTTg                 >  1:2154491/1‑69 (MQ=255)                                    tttAAGCTCATCAACCAGAATTTCTGGGATGTTTTCTGCCATTTTCTTAATGACTTTATCGCGAGCTTg                 >  1:1004129/1‑69 (MQ=255)                                    tttAAGCTCATCAACCAGAATTTCTGGGATGTTTTCTGCCATTTTCTTAATGACTTTATCGCGAGCTTg                 >  1:1127895/1‑69 (MQ=255)                                    tttAAGCTCATCAACCAGAATTTCTGGGATGTTTTCTGCCATTTTCTTAATGACTTTATCGCGAGCTTg                 >  1:2728221/1‑69 (MQ=255)                                    tttAAGCTCATCAACCAGAATTTCTGGGATGTTTTCTGCCATTTTCTTAATGACTTTATCGCGAGCTTg                 >  1:1331043/1‑69 (MQ=255)                                       aaGCTCATCAACCAGAATTTCTGGGATGTTTTCTGCCATTTTCTTAATGACTTTATCGCGAGCTTGCt               <  1:2990791/68‑1 (MQ=255)                                       aaGCTCATCAACCAGAATTTCTGGGATGTTTTCTGCCATTTTCTTAATGACTTTATCGCGAGCTTGCt               <  1:3348979/68‑1 (MQ=255)                                       aaGCTCATCAACCAGAATTTCTGGGATGTTTTCTGCCATTTTCTTAATGACTTTATCGCGAGCTTGCt               <  1:49112/68‑1 (MQ=255)                                       aaGCTCATCAACCAGAATTTCTGGGATGTTTTCTGCCATTTTCTTAATGACTTTATCGCGAGCTTGCt               <  1:2268788/68‑1 (MQ=255)                                       aaGCTCATCAACCAGAATTTCTGGGATGTTTTCTGCCATTTTCTTAATGACTTTATCGCGAGCTTGCt               <  1:1355350/68‑1 (MQ=255)                                       aaGCTCATCAACCAGAATTTCTGGGATGTTTTCTGCCATTTTCTTAATGACTTTATCGCGAGCTTGCt               <  1:432817/68‑1 (MQ=255)                                       aaGCTCATCAACCAGAATTTCTGGGATGTTTTCTGCCATTTTCTTAATGACTTTATCGCGAGCTTGCt               <  1:1556750/68‑1 (MQ=255)                                       aaGCTCATCAACCAGAATTTCTGGGATGTTTTCTGCCATTTTCTTAATGACTTTATCGCGAGCTTGCt               <  1:156291/68‑1 (MQ=255)                                         gCTCATCAACCAGAATTTCTGGGATGTTTTCTGCCTTTTTCTTAATGACTTTATCGCGAGCTTGCt               <  1:1416975/66‑1 (MQ=255)                                         gCTCATCAACCAGAATTTCTGGGATGTTTTCTGCCATTTTCTTAATGACTTTATCGCGAGCTTGCTCa             <  1:30890/68‑1 (MQ=255)                                         gCTCATCAACCAGAATTTCTGGGATGTTTTCTGCCATTTTCTTAATGACTTTATCGCGAGCTTGCTCa             <  1:1965484/68‑1 (MQ=255)                                         gCTCATCAACCAGAATTTCTGGGATGTTTTCTGCCATTTTCTTAATGACTTTATCGCGAGCTTGCTCa             <  1:210621/68‑1 (MQ=255)                                         gCTCATCAACCAGAATTTCTGGGATGTTTTCTGCCATTTTCTTAATGACTTTATCGCGAGCTTGCTCa             <  1:2856811/68‑1 (MQ=255)                                         gCTCATCAACCAGAATTTCTGGGATGTTTTCTGCCATTTTCTTAATGACTTTATCGCGAGCTTGCTCa             <  1:2359659/68‑1 (MQ=255)                                         gCTCATCAACCAGAATTTCTGGGATGTTTTCTGCCATTTTCTTAATGACTTTATCGCGAGCTTGCTCa             <  1:2624923/68‑1 (MQ=255)                                         gCTCATCAACCAGAATTTCTGGGATGTTTTCTGCCATTTTCTTAATGACTTTATCGCGAGCTTGCTCa             <  1:2861124/68‑1 (MQ=255)                                         gCTCATCAACCAGAATTTCTGGGATGTTTTCTGCCATTTTCTTAATGACTTTATCGCGAGCTTGCTCAt            >  1:2807949/1‑69 (MQ=255)                                         gCTCATCAACCAGAATTTCTGGGATGTTTTCTGCCATTTTCTTAATGACTTTATCGCGAGCTTGCTCAt            >  1:2811911/1‑69 (MQ=255)                                          cTCATCAACCAGAATTTCTGGGATGTTTTCTGCCATTTTCTTAATGACTTTATCGCGAGCTTGCTCa             <  1:3374591/67‑1 (MQ=255)                                                 aCCAGAATTTCTGGGATGTTTTCTGCCATTTTCTTAATGACTTTa                            <  1:2597152/45‑1 (MQ=255)                                                  ccAGAATTTCTGGGATATTTTCTGCCATTTTCTTAATGACTTTATCGCGAGCTTGTTCATACTGAAggc   <  1:963050/69‑2 (MQ=255)                                                  ccAGAATTTCTGGGATATTTTCTGCCATTTTCTTAATGACTTTATCGCGAGCTTGTTCATACTGAAggc   <  1:1683346/69‑2 (MQ=255)                                                  ccAGAATTTCTGGGATATTTTCTGCCATTTTCTTAATGACTTTATCGCGAGCTTGTTCATACTGAAggc   <  1:2928939/69‑2 (MQ=255)                                                  ccAGAATTTCTGGGATATTTTCTGCCATTTTCTTAATGACTTTATCGCGAGCTTGTTCATACTGAAggc   <  1:733946/69‑2 (MQ=255)                                                  ccAGAATTTCTGGGATATTTTCTGCCATTTTCTTAATGACTTTATCGCGAGCTTGTTCATACTGAAggc   <  1:72005/69‑2 (MQ=255)                                                  ccAGAATTTCTGGGATATTTTCTGCCATTTTCTTAATGACTTTATCGCGAGCTTGTTCATACTGAAggc   <  1:699341/69‑2 (MQ=255)                                                  ccAGAATTTCTGGGATATTTTCTGCCATTTTCTTAATGACTTTATCGCGAGCTTGTTCATACTGAAggc   <  1:1231824/69‑2 (MQ=255)                                                  ccAGAATTTCTGGGATATTTTCTGCCATTTTCTTAATGACTTTATCGCGAGCTTGTTCATACTGAAggc   <  1:1271975/69‑2 (MQ=255)                                                  ccAGAATTTCTGGGATATTTTCTGCCATTTTCTTAATGACTTTATCGCGAGCTTGTTCATACTGAAggc   <  1:3008960/69‑2 (MQ=255)                                                  ccAGAATTTCTGGGATATTTTCTGCCATTTTCTTAATGACTTTATCGCGAGCTTGTTCATACTGAAggc   <  1:3135102/69‑2 (MQ=255)                                                  ccAGAATTTCTGGGATATTTTCTGCCATTTTCTTAATGACTTTATCGCGAGCTTGTTCATACTGAAggc   <  1:1897617/69‑2 (MQ=255)                                                  ccAGAATTTCTGGGATATTTTCTGCCATTTTCTTAATGACTTTATCGCGAGCTTGTTCATACTGAAggc   <  1:405248/69‑2 (MQ=255)                                                   cAGAATTTCTGGGATGTTTTCTGCCATTTTCTTAATGACTTTATCGCGAGCTTGCTCATACTGAAGgtg  <  1:1382226/69‑1 (MQ=255)                                                   cAGAATTTCTGGGATATTTTCTGCCATTTTCTTAATGACTTTATCGCGAGCTTGTTCATACTGAAggc   <  1:2901062/68‑2 (MQ=255)                                                       aTTTCTGGGATGTTTTCTGCAATTTTCTTAATGACTTTATCGCGAGCTTGCt               <  1:3328470/52‑1 (MQ=255)                                                          tCTTGGATGTTTTCTGCCATTTTCTTAATGACTTTATCGCGAGCTTGCTCa             <  1:1424159/51‑1 (MQ=255)                                                                  |                                                      TCCAGCGTGTTAACCCGACGATTAATATCGGCATTTTTAAGCTCATCAACCAGAATTTCTGGGATGTTTTCTGCCATTTTCTTAATGACTTTATCGCGAGCTTGCTCATACTGAAGGTG  >  minE/2835586‑2835704

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