Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	779,653	G→A	48.1%	L12L (CTG→TTG)	ycgK ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	779,653	0	G	A	48.1%	3.7 / 70.7	81	L12L (CTG→TTG)	ycgK	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (21/21);  new base A (21/18);  total (42/39)
Fisher's exact test for biased strand distribution p-value = 8.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.39e-01

CTTTTCTGAACTCAACATTGACATTCTTACCCGCAGCAAAGCTCGTGGAGGTTAGCAGTGCAAGCAGCAATACAGCCTTCCTGATGCTCTTGATTTTCATTGGTAGTCTCCTGTTCAGAATGAGTACGATG  >  minE/779587‑779717                                                                   |                                                                 cTTTTCTGAACTCAACATTGACATTCTTACCCGCAGCAAAGCTCGTGGAGGTTAGCATTGCAagcagc                                                                 <  1:614384/68‑1 (MQ=255) cTTTTCTGAACTCAACATTGACATTCTTACCCGCAGCAAAGCTCGTGGAGGTTAGCAGTGCAagcagc                                                                 <  1:2663227/68‑1 (MQ=255) cTTTTCTGAACTCAACATTGACATTCTTACCCGCAGCAAAGCTCGTGGAGGTTAGCAGTGCAagcagc                                                                 <  1:1473159/68‑1 (MQ=255) cTTTTCTGAACTCAACATTGACATTCTTACCCGCAGCAAAGCTCGTGGAGGTTAGCAGTGCAagcagc                                                                 <  1:2630194/68‑1 (MQ=255) cTTTTCTGAACTCAACATTGACATTCTTACCCGCAGCAAAGCTCGTGGAGGTTAGCAGTGCAagcagc                                                                 <  1:433218/68‑1 (MQ=255)  ttttCTGAACTCAACATTGACATTCTTACCCGCAGCAAAGCTCGTGGAGGTTAGCAGTGCAAGCAGCaa                                                               <  1:549081/69‑1 (MQ=255)  ttttCTGAACTCAACATTGACATTCTTACCCGCAGCAAAGCTCGTGGAGGTTAGCAGTGCAAGCAGCaa                                                               <  1:1632587/69‑1 (MQ=255)  ttttCTGAACTCAACATTGACATTCTTACCCGCAGCAAAGCTCGTGGAGGTTAGCAGTGCAAGCAGCaa                                                               <  1:1485170/69‑1 (MQ=255)  ttttCTGAACTCAACATTGACATTCTTACCCGCAGCAAAGCTCGTGGAGGTTAGCAGTGCAAGCAGCaa                                                               <  1:3218945/69‑1 (MQ=255)  ttttCTGAACTCAACATTGACATTCTTACCCGCAGCAAAGCTCGTGGAGGTTAGCAGTGCAAGCAGCaa                                                               <  1:2370948/69‑1 (MQ=255)  ttttCTGAACTCAACATTGACATTCTTACCCGCAGCAAAGCTCGTGGAGGTTAGCAGTGCAAGCAGCaa                                                               <  1:961721/69‑1 (MQ=255)          aCTCAACATTGACATTCTTACCCGCAGCAAAGCTCGTGGAGGTTAGCAGTGCAAGCAACAATACaaga                                                        >  1:2605852/1‑65 (MQ=255)          aCTCAACATTGACATTCTTACCCGCAGCAAAGCTCGTGGAGGTTAGCAGTGCAAGCAACAATACaaga                                                        >  1:2790690/1‑65 (MQ=255)          aCTCAACATTGACATTCTTACCCGCAGCAAAGCTCGTGGAGGTTAGCAGTGCAAGCAACAATACaaga                                                        >  1:2225847/1‑65 (MQ=255)          aCTCAACATTGACATTCTTACCCGCAGCAAAGCTCGTGGAGGTTAGCAGTGCAAGCAACAATACaaga                                                        >  1:2211269/1‑65 (MQ=255)          aCTCAACATTGACATTCTTACCCGCAGCAAAGCTCGTGGAGGTTAGCAGTGCAAGCAACAATACaaga                                                        >  1:2640696/1‑65 (MQ=255)          aCTCAACATTGACATTCTTACCCGCAGCAAAGCTCGTGGAGGTTAGCAGTGCAAGCAACAATACaaga                                                        >  1:1698896/1‑65 (MQ=255)          aCTCAACATTGACATTCTTACCCGCAGCAAAGCTCGTGGAGGTTAGCAGTGCAAGCAACAATACaaga                                                        >  1:1697660/1‑65 (MQ=255)          aCTCAACATTGACATTCTTACCCGCAGCAAAGCTCGTGGAGGTTAGCAGTGCAAGCAACAATACaaga                                                        >  1:2198225/1‑65 (MQ=255)          aCTCAACATTGACATTCTTACCCGCAGCAAAGCTCGTGGAGGTTAGCAGTGCAAGCAACAATACaaga                                                        >  1:214030/1‑65 (MQ=255)          aCTCAACATTGACATTCTTACCCGCAGCAAAGCTCGTGGAGGTTAGCAGTGCAAGCAACAATACaaga                                                        >  1:2111947/1‑65 (MQ=255)          aCTCAACATTGACATTCTTACCCGCAGCAAAGCTCGTGGAGGTTAGCAGTGCAAGCAACAATACaaga                                                        >  1:2738992/1‑65 (MQ=255)          aCTCAACATTGACATTCTTACCCGCAGCAAAGCTCGTGGAGGTTAGCAGTGCAAGCAACAATACaaga                                                        >  1:790187/1‑65 (MQ=255)          aCTCAACATTGACATTCTTACCCGCAGCAAAGCTCGTGGAGGTTAGCAGTGCAAGCAACAATACaaga                                                        >  1:1209933/1‑65 (MQ=255)          aCTCAACATTGACATTCTTACCCGCAGCAAAGCTCGTGGAGGTTAGCAGTGCAAGCAACAATACaaga                                                        >  1:1272730/1‑65 (MQ=255)          aCTCAACATTGACATTCTTACCCGCAGCAAAGCTCGTGGAGGTTAGCAGTGCAAGCAACAATACaaga                                                        >  1:3086588/1‑65 (MQ=255)          aCTCAACATTGACATTCTTACCCGCAGCAAAGCTCGTGGAGGTTAGCAGTGCAAGCAACAATACaaga                                                        >  1:2923179/1‑65 (MQ=255)          aCTCAACATTGACATTCTTACCCGCAGCAAAGCTCGTGGAGGTTAGCAGTGCAAGCAACAATACaaga                                                        >  1:2919292/1‑65 (MQ=255)          aCTCAACATTGACATTCTTACCCGCAGCAAAGCTCGTGGAGGTTAGCAGTGCAAGCAACAATACaaga                                                        >  1:2797716/1‑65 (MQ=255)          aCTCAACATTGACATTCTTACCCGCAGCAAAGCTCGTGGAGGTTAGCAGTGCAAGCAACAATACaag                                                         >  1:2360055/1‑65 (MQ=255)          aCTCAACATTGACATTCTTACCCGCAGCAAAGCTCGTGGAGGTTAGCAGTGCAAGCAACAATACaag                                                         >  1:1092122/1‑65 (MQ=255)          aCTCAACATTGACATTCTTACCCGCAGCAAAGCTCGTGGAAGTTAGCAGTGCAAGCAACAATACaaga                                                        >  1:3078320/1‑65 (MQ=255)           cTCAACATTGACATTCTTACCCGCAGCAAAGCTCGTGGAGGTTAGCAGTGCAAGCAGCAATACAGcc                                                        <  1:199670/67‑1 (MQ=255)           cTCAACATTGACATTCTTACCCGCAGCAAAGCTCGTGGAGGTTAGCAGTGCAAGCAGCAATACAGcc                                                        <  1:2113172/67‑1 (MQ=255)           cTCAACATTGACATTCTTACCCGCAGCAAAGCTCGTGGAGGTTAGCAGTGCAAGCAGCAATACAGcc                                                        <  1:714845/67‑1 (MQ=255)           cTCAACATTGACATTCTTACCCGCAGCAAACCTCGTGGAGGTTAGCAGTGCAAGCAGCAATACAGcc                                                        <  1:1094431/67‑1 (MQ=255)             cAACATTGACATTCTTACCCGCAGCAAAGCTCGTGGAGGTTAGCAGTGCAAGCAGCAATACAGCCtt                                                      >  1:1949039/1‑67 (MQ=255)             cAACATTGACATTCTTACCCGCAGCAAAGCTCGTGGAGGTTAGCAGTGCAAGCAGCAATACAGCCTTc                                                     >  1:1002717/1‑68 (MQ=255)             cAACATTGACATTCTTACCCGCAGCAAAGCTCGTGGAGGTTAGCAGTGCAAGCAGCAATACAGCCTTc                                                     >  1:2488388/1‑68 (MQ=255)             cAACATTGACATTCTTACCCGCAGCAAAGCTCGTGGAGGTTAGCAGTGCAAGCAGCAATACAGCCTTc                                                     >  1:966163/1‑68 (MQ=255)             cAACATTGACATTCTTACCCGCAGCAAAGCTCGTGGAGGTTAGCAGTGCAAGCAGCAATACAGCCTTc                                                     >  1:27185/1‑68 (MQ=255)             cAACATTGACATTCTTACCCGCAGCAAAGCTCGTGGAGGTTAGCAGTGCAAGCAGCAATACAGCCTTc                                                     >  1:1290664/1‑68 (MQ=255)             cAACATTGACATTCTTACCCGCAGCAAAGCTCGTGGAGGTTAGCAGTGCAAGCAGCAATACAGCCTTc                                                     >  1:1312499/1‑68 (MQ=255)             cAACATTGACATTCTTACCCGCAGCAAAGCTCGTGGAGGTTAGCAGTGCAAGCAGCAATACAGCCTTc                                                     >  1:1904274/1‑68 (MQ=255)             cAACATTGACATTCTAACCCGCAGCAAAGCTCGTGGAGGTTAGCAGTGCAAGCAGCAATACAGCCTTc                                                     >  1:3134286/1‑68 (MQ=255)                cATTGACATTCTTACCCGCAGCAAAGCTCGTGGAGGTTAGCAGTGCAAGCAGCAATACAGCCTTCCt                                                   <  1:2792654/67‑1 (MQ=255)                    gACATTCTTACCCGCAGCAAAGCTCGTGGAGGTTAGCAGTGCAAGCAGCAATACAGCCt                                                       >  1:1811147/1‑59 (MQ=255)                    gACATTCTTACCCGCAGCAAAGCTCGTGGAGGTTAGCAGTGCAAGCAGCAATACAGCCt                                                       <  1:1839151/59‑1 (MQ=255)                    gACATTCTTACCCGCAGCAAAGCTCGTGGAGGTTAGCAGTGCAAGCAGCAATACAGCCt                                                       <  1:81680/59‑1 (MQ=255)                        ttCTTACCCGCAGCAAAGCTCGTGGAGGTTAGCAGTGCAAGCAGCAATACAGCCTTCCTGATGCTCTTg                                         <  1:1331733/69‑1 (MQ=255)                          cTTACCCGCAGCAAAGCTCGTGGAGGTTAGCAGTGCAAGCAGCAATa                                                             >  1:1384333/1‑47 (MQ=255)                          cTTACCCGCAGCAAAGCTCGTGGAGGTTAGCAGTGCAAGCAGCAATACAGCCTTCCTTATGCTCTTGAt                                       >  1:628452/1‑69 (MQ=255)                          cTTACCCGCAGCAAAGCTCGTGGAGGTTAGCAGTGCAAGCAGCAATACAGCCTTCCTGATGCTCTTGAt                                       >  1:490405/1‑69 (MQ=255)                          cTTACCCGCAGCAAAGCTCGTGGAGGTTAGCAGTGCAAGCAGCAATACAGCCTTCCTGATGCTCTTGAt                                       >  1:1793894/1‑69 (MQ=255)                          cTTACCCGCAGCAAAGCTCGTGGAGGTTAGCAGTGCAAGCAGCAATACAGCCTTCCTGATGCTCTTGAt                                       >  1:346853/1‑69 (MQ=255)                          cTTACCCGCAGCAAAGCTCGTGGAGGTTAGCAGTGCAAGCAGCAATACAGCCTTCCTGATGCTCTTGAt                                       >  1:1348060/1‑69 (MQ=255)                          cTTACCCGCAGCAAAGCTCGTGGAGGTTAGCAGTGCAAGCAGCAATACAGCCTTCCTGATGCTCTTGAt                                       >  1:1916591/1‑69 (MQ=255)                          cTTACCCGCAGCAAAGCTCGTGGAGGTTAGCAGTGCAAGCAGCAATACAGCCTTCCTGATGCTCTTGAt                                       >  1:717443/1‑69 (MQ=255)                          cTTACCCGCAGCAAAGCTCGTGGAGGTTAGCAGTGCAAGCAGCAATACAGCCTTCCTGATGCTCTTGAt                                       >  1:1017637/1‑69 (MQ=255)                               ccGCAGCAAAGCTCGTGGAGGTTAGCAGTGCAAGCAGCAATACAGCCTTCCTGATGCTCTTGAtttt                                    <  1:2432168/67‑1 (MQ=255)                               ccGCAGCAAAGCTCGTGGAGGTTAGCAGTGCAAGCAGCAATACAGCCTTCCTGATGCTCTTGAtttt                                    <  1:652102/67‑1 (MQ=255)                                   agcaAAGCTCGTGGAGGTTAGCGGTGCAAGCAGCAATACAGCCTTCCt                                                   <  1:1375535/48‑1 (MQ=255)                                      aaaGCTCGTGGAGGTTAGCAGTGCAAGCAGCAATACAGCCt                                                       <  1:2987252/41‑1 (MQ=255)                                          cTCGTGGAGGTTAGCAGTGCAAGCAGCAATACAGCCTTCCTGATGCTCTTGATTTTCATTGGTAGtctc                       <  1:1360784/69‑1 (MQ=255)                                          cTCGTGGAGGTTAGCAGTGCAAGCAGCAATACAGCCTTCCTGATGCTCTTGATTTTCATTGGTAGtctc                       <  1:1932656/69‑1 (MQ=255)                                             gTGGAGGTTAGCAGTGCAAGCAACAATACAGCCTTACTGATGCCCTTAATTTTCATTg                               <  1:387501/58‑1 (MQ=255)                                             gTGGAGGTTAGCAGTGCAAGCAACAATACAGCCTTACTGATGCCCTTAATTTTCATTg                               <  1:51326/58‑1 (MQ=255)                                             gTGGAGGTTAGCAGTGCAAGCAACAATACAGCCTTACTGATGCCCTTAATTTTCATTg                               <  1:276425/58‑1 (MQ=255)                                             gTGGAGGTTAGCAGTGCAAGCAACAATACAGCCTTACTGATGCCCTTAATTTTCATTg                               <  1:2888983/58‑1 (MQ=255)                                             gTGGAGGTTAGCAGTGCAAGCAACAATACAGCCTTACTGATGCCCTTAATTTTCATTg                               <  1:2854916/58‑1 (MQ=255)                                             gTGGAGGTTAGCAGTGCAAGCAACAATACAGCCTTACTGATGCCCTTAATTTTCATTg                               <  1:680700/58‑1 (MQ=255)                                             gTGGAGGTTAGCAGTGCAAGCAACAATACAGCCTTACTGATGCCCTTAATTTTCATTg                               <  1:1199588/58‑1 (MQ=255)                                             gTGGAGGTTAGCAGTGCAAGCAACAATACAGCCTTACTGATGCCCTTAATTTTCATTg                               <  1:1620679/58‑1 (MQ=255)                                             gTGGAGGTTAGCAGTGCAAGCAACAATACAGCCTTACTGATGCCCTTAATTTTCATTg                               <  1:1628608/58‑1 (MQ=255)                                             gTGGAGGTTAGCAGTGCAAGCAACAATACAGCCTTACTGATGCCCTTAATTTTCATTg                               <  1:1711410/58‑1 (MQ=255)                                             gTGGAGGTTAGCAGTGCAAGCAACAATACAGCCTTACTGATGCCCTTAATTTTCATTg                               <  1:2523353/58‑1 (MQ=255)                                             gTGGAGGTTAGCAGTGCAAGCAACAATACAGCCTTACTGATGCCCTTAATTTTCATTg                               <  1:1851492/58‑1 (MQ=255)                                             gTGGAGGTTAGCAGTGCAAGCAACAATACAGCCTTACTGATGCCCTTAATTTTCATTg                               <  1:89592/58‑1 (MQ=255)                                             gTGGAGGTTAGCAGTGCAAGCAACAATACAGCCTTACTGATGCCCTTAATTTTCATTg                               <  1:2278281/58‑1 (MQ=255)                                             gTGGAGGTTAGCAGTGCAAGCAACAATACAGCCTTACTGATGCCCTTAATTTTCATTg                               <  1:197509/58‑1 (MQ=255)                                             gTGGAGGTTAGCAGTGCAAGCAACAATACAGCCTTACTGATGCCCTTAATTTTCATTg                               <  1:2058086/58‑1 (MQ=255)                                             gTGGAGGTTAGCAGTGCAAGCAACAAAACAGCCTTACTGATGCCCTTAATTTTCATTg                               <  1:2015591/58‑1 (MQ=255)                                                    ttAGCAGTGCAAGCAACAATACAGCCTTACTGATGCCCTTAATTTTCATTg                               <  1:2821493/51‑1 (MQ=255)                                                            gCAAGCAGCAATACAGCCTTCCTGATGCTCTTGATTTTCATTGGTAGtctc                       <  1:863381/51‑1 (MQ=255)                                                             cAAGCAGCAATACAGCCTTCCTGATGCTCTTGATTTTCATTGGTAGTCTCCTGTTCAGAATGAGTACg     >  1:242909/1‑68 (MQ=255)                                                             cAAGCAGCAATACAGCCTTCCTGATGCTCTTGATTTTCATTGGTAGTCTCCTGTTCAGAATGAGTACg     >  1:2200677/1‑68 (MQ=255)                                                               agcagcAATACAGCCTTCCTGATGCTCTTGATTTTCATTGGTAGTCTCCTGTTCAGAATGAGTACGATg  <  1:625156/69‑1 (MQ=255)                                                               agcagcAATACAGCCTTCCTGATGCTCTTGATTTTCATTGGTAGTCTCCTGTTCAGAATGAGTACGATg  <  1:1422417/69‑1 (MQ=255)                                                               agcagcAATACAGCCTTCCTGATGCTCTTGATTTTCATTGGTAGTCTCCTGTTCAGAATGAGTACGATg  <  1:1517512/69‑1 (MQ=255)                                                               agcagcAATACAGCCTTCCTGATGCTCTTGATTTTCATTGGTAGTCTCCTGTTCAGAATGAGTACGATg  <  1:1519965/69‑1 (MQ=255)                                                               agcagcAATACAGCCTTCCTGATGCTCTTGATTTTCATTGGTAGTCTCCTGTTCAGAATGAGTACGATg  <  1:1570830/69‑1 (MQ=255)                                                               agcagcAATACAGCCTTCCTGATGCTCTTGATTTTCATTGGTAGTCTCCTGTTCAGAATGAGTACGATg  <  1:1185749/69‑1 (MQ=255)                                                               agcagcAATACAGCCTTCCTGATGCTCTTGATTTTCATTGGTAGTCTCCTGTTCAGAATGAGTACGATg  <  1:888946/69‑1 (MQ=255)                                                               agcagcAATACAGCCTTCCTGATGCTCTTGATTTTCATTGGTAGTCTCCTGTTCAGAATGAGTACGATg  <  1:190003/69‑1 (MQ=255)                                                               agcagcAATACAGCCTTCCTGATGCTCTTGATTTTCATTGGTAGTCTCCTGTTCAGAATGAGTACGATg  <  1:1055839/69‑1 (MQ=255)                                                                   |                                                                 CTTTTCTGAACTCAACATTGACATTCTTACCCGCAGCAAAGCTCGTGGAGGTTAGCAGTGCAAGCAGCAATACAGCCTTCCTGATGCTCTTGATTTTCATTGGTAGTCTCCTGTTCAGAATGAGTACGATG  >  minE/779587‑779717

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