Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,571,096	A→C	63.4%	D334A (GAT→GCT)	yibH →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,571,096	0	A	C	63.4%	62.3 / 53.5	71	D334A (GAT→GCT)	yibH	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (10/16);  new base C (21/24);  total (31/40)
Fisher's exact test for biased strand distribution p-value = 6.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.64e-01

TTGCAATCATTAACGGTCGTGCCCGGCACGGACGGTGTGCTGGGAACCATTGAACTGGACCCTAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATT  >  minE/2571029‑2571155                                                                    |                                                            ttGCAATCATTAACGGTCGTGCCCGGCACGGACGGTGTGCTGGGAACCATTGAACTGGACCCTAACgat                                                            <  1:1073325/69‑1 (MQ=255) ttGCAATCATTAACGGTCGTGCCCGGCACGGACGGTGTGCTGGGAACCATTGAACTGGACCCTAACgat                                                            <  1:108121/69‑1 (MQ=255) ttGCAATCATTAACGGTCGTGCCCGGCACGGACGGTGTGCTGGGAACCATTGAACTGGACCCTAACgat                                                            <  1:1049136/69‑1 (MQ=255) ttGCAATCATTAACGGTCGTGCCCGGCACGGACGGTGTGCTGGGAACCATTGAACTGGACCCTAACgat                                                            <  1:1728355/69‑1 (MQ=255) ttGCAATCATTAACGGTCGTGCCCGGCACGGACGGTGTGCTGGGAACCATTGAACTGGACCCTAACgat                                                            <  1:1696955/69‑1 (MQ=255) ttGCAATCATTAACGGTCGTGCCCGGCACGGACGGTGTGCTGGGAACCATTGAACTGGACCCTAACgat                                                            <  1:2110145/69‑1 (MQ=255) ttGCAATCATTAACGGTCGTGCCCGGCACGGACGGTGTGCTGGGAACCATTGAACTGGACCCTAACgat                                                            <  1:453926/69‑1 (MQ=255)  tGCAATCATTAACGGTCGTGCCCGGCACGGACGGTGTGCTGGGAACCATTGAACTGGACCCTAACgat                                                            <  1:3224180/68‑1 (MQ=255)  tGCAATCATTAACGGTCGTGCACGGCACGGACGGTGTGCTGGGAACCATTGAACTGGACCCTAACgat                                                            <  1:2859614/68‑1 (MQ=255)        cATTAACGGTCGTGCCCGGCACGGACGGTGTGCTGGGAACCATTGAACTGGACCCTAACGATGATATCg                                                     >  1:1052176/1‑69 (MQ=255)        cATTAACGGTCGTGCCCGGCACGGACGGTGTGCTGGGAACCATTGAACTGGACCCTAACGATGATATCg                                                     >  1:178507/1‑69 (MQ=255)        cATTAACGGTCGTGCCCGGCACGGACGGTGTGCTGGGAACCATTGAACTGGACCCTAACGATGATATCg                                                     >  1:1587736/1‑69 (MQ=255)        cATTAACGGTCGTGCCCGGCACGGACGGTGTGCTGGGAACCATTGAACTGGACCCTAACGATGATATCg                                                     >  1:3026694/1‑69 (MQ=255)        cATTAACGGTCGTGCCCGGCACGGACGGTGTGCTGGGAACCATTGAACTGGACCCTAACGATGATATCg                                                     >  1:1372038/1‑69 (MQ=255)        cATTAACGGTCGTGCCCGGCACGGACGGTGTGCTGGGAACCATTGAACTGGACCCTAACGATGATATCg                                                     >  1:1320647/1‑69 (MQ=255)        cATTAACGGTCGTGCCCGGCACGGACGGTGTGCTGGGAACCATTGAACTGGACCCTAACGATGATATCg                                                     >  1:754897/1‑69 (MQ=255)        cATTAACGGTCGTGCCCGGCACGGACGGTGTGCTGGGAACCATTGAACTGGACCCTAACGATGATATCg                                                     >  1:1223308/1‑69 (MQ=255)            aaCGGTCGTGCCCGGCACGGACGGTGTGCTGGGAACCATTGAACTGGACCCTAACgat                                                            >  1:381164/1‑58 (MQ=255)            aaCGGTCGTGCCCGGCACGGACGGTGTGCTGGGAACCATTGAACTGGACCCTAACgat                                                            >  1:2889564/1‑58 (MQ=255)            aaCGGTCGTGCCCGGCACGGACGGTGTGCTGGGAACCATTGAACTGGACCCTAACgat                                                            <  1:2887302/58‑1 (MQ=255)            aaCGGTCGTGCCCGGCACGGACGGTGTGCTGGGAACCATTGAACTGGACCCTAACgat                                                            <  1:2605332/58‑1 (MQ=255)            aaCGGTCGTGCCCGGCACGGACGGTGTGCTGGGAACCATTGAACTGGACCCTAACgat                                                            <  1:55831/58‑1 (MQ=255)            aaCGGTCGTGCCCGGCACGGACGGTGTGCTGGGAACCATTGAACTGGACCCTAACgat                                                            >  1:2044224/1‑58 (MQ=255)            aaCGGTCGTGCCCGGCACGGACGGTGTGCTGGGAACCATTGAACTGGACCCTAACgat                                                            <  1:960360/58‑1 (MQ=255)            aaCGGTCGTGCCCGGCACGGACGGTGTGCTGGGAACCATTGAACTGGACCCTAACgat                                                            >  1:1746245/1‑58 (MQ=255)            aaCGGTCGTGCCCGGCACGGACGGTGTGCTGGGAACCATTGAACTGGACCCTAACgat                                                            <  1:1503639/58‑1 (MQ=255)            aaCGGTCGTGCCCGGCACGGACGGTGTGCTGGGAACCATTGAACTGGACCCTAACgat                                                            >  1:1556120/1‑58 (MQ=255)            aaCGGTCGTGCCCGGCACGGACGGTGTGCTGGGAACCATTGAACTGGACCCTAACGATGATATCGATGc                                                 <  1:1511825/69‑1 (MQ=255)            aaCGGTCGTGCCCGGCACGGACGGTGTGCTGGGAACCATTGAACTGGACCCTAACGATGATATCGATGc                                                 <  1:1963684/69‑1 (MQ=255)             aCGGTCGTGCCCGGCACGGACGGTGTGCTGGGAACCATTGAACTGGACCCTAACgat                                                            <  1:1593165/57‑1 (MQ=255)              cGGTCGTGCGCGGCTCGGACGGTGTGCTGGGAACCATTGAACTGGACCCTAACgat                                                            <  1:509119/56‑1 (MQ=255)                     gCCCGGCACGGACGGTGTGCTGGGAACCATTGAACTGGACCCTAACgat                                                            <  1:1921457/49‑1 (MQ=255)                         ggCACGGACGGTGTGCTGGGAACCATTGAACTGGACCCTAACGATGATATCGATGCCTTACCCGACGGc                                    <  1:1438398/69‑1 (MQ=255)                         ggCACGGACGGTGTGCTGGGAACCATTGAACTGGACCCTAACGATGATATCGATGCCTTACCCGACGGc                                    <  1:3282736/69‑1 (MQ=255)                         ggCACGGACGGTGTGCTGGGAACCATTGAACTGGACCCTAACGATGATATCGATGCCTTACCCGACGGc                                    <  1:502888/69‑1 (MQ=255)                          gCACGGACGGTGTGCTGGGAACCATTGAACTGGACCCTAACGATGATATCGATGCCTTACCCGACGGc                                    <  1:1423605/68‑1 (MQ=255)                                          gggAACCATTGAACTGGACCCTAACGCTGATATCGATGCCTTAc                                            >  1:3006530/1‑44 (MQ=255)                                          gggAACCATTGAACTGGACCCTAACGCTGATATCGATGCCTTAc                                            >  1:3022487/1‑44 (MQ=255)                                          gggAACCATTGAACTGGACCCTAACGCTGATATCGATGCCTTAc                                            >  1:288414/1‑44 (MQ=255)                                          gggAACCATTGAACTGGACCCTAACGCTGATATCGATGCCTTAc                                            >  1:2726689/1‑44 (MQ=255)                                          gggAACCATTGAACTGGACCCTAACGCTGATATCGATGCCTTAc                                            >  1:3177884/1‑44 (MQ=255)                                          gggAACCATTGAACTGGACCCTAACGCTGATATCGATGCCTTAc                                            >  1:3385941/1‑44 (MQ=255)                                          gggAACCATTGAACTGGACCCTAACGCTGATATCGATGCCTTAc                                            >  1:3218622/1‑44 (MQ=255)                                          gggAACCATTGAACTGGACCCTAACGCTGATATCGATGCCTTAc                                            >  1:3261248/1‑44 (MQ=255)                                          gggAACCATTGAACTGGACCCTAACGCTGATATCGATGCCTTAc                                            >  1:1126145/1‑44 (MQ=255)                                          gggAACCATTGAACTGGACCCTAACGCTGATATCGATGCCTTAc                                            >  1:384410/1‑44 (MQ=255)                                          gggAACCATTGAACTGGACCCTAACGCTGATATCGATGCCTTAc                                            >  1:511566/1‑44 (MQ=255)                                          gggAACCATTGAACTGGACCCTAACGCTGATATCGATGCCTTAc                                            >  1:574046/1‑44 (MQ=255)                                          gggAACCATTGAACTGGACCCTAACGCTGATATCGATGCCTTAc                                            >  1:585273/1‑44 (MQ=255)                                          gggAACCATTGAACTGGACCCTAACGCTGATATCGATGCCTTAc                                            >  1:75131/1‑44 (MQ=255)                                          gggAACCATTGAACTGGACCCTAACGCTGATATCGATGCCTTAc                                            >  1:1305413/1‑44 (MQ=255)                                          gggAACCATTGAACTGGACCCTAACGCTGATATCGATGCCTTAc                                            >  1:803823/1‑44 (MQ=255)                                          gggAACCATTGAACTGGACCCTAACGCTGATATCGATGCCTTAc                                            >  1:2008791/1‑44 (MQ=255)                                          gggAACCATTGAACTGGACCCTAACGCTGATATCGATGCCTTAc                                            >  1:1732737/1‑44 (MQ=255)                                          gggAACCATTGAACTGGACCCTAACGCTGATATCGATGCCTTAc                                            >  1:1912056/1‑44 (MQ=255)                                          gggAACCATTGAACTGGACCCTAACGCTGATATCGATGCCTTAc                                            >  1:1904705/1‑44 (MQ=255)                                          gggAACCATTGAACTGGACCCTAACGCTGATATCGATGCCTTAc                                            >  1:2356523/1‑44 (MQ=255)                                             aaCCATTGAACTGGACCCTAACGATGATATCGATGCCTTACCCGACGGc                                    <  1:1972849/49‑1 (MQ=255)                                             aaCCATTGAACTGGACCCTAACGATGATATCGATGCCTTACCCGACGGc                                    <  1:1917789/49‑1 (MQ=255)                                              aCCATTGAACTGGACCCTAACGATGATATCGATGCCTTACCCGACGGc                                    <  1:3012561/48‑1 (MQ=255)                                              aCCATTGAACTGGACCCTAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGtt               <  1:2177660/69‑1 (MQ=255)                                                    gAACTGGACCCTAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTcc         >  1:2609844/1‑69 (MQ=255)                                                    gAACTGGACCCTAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTcc         >  1:2386703/1‑69 (MQ=255)                                                       cTGGACCCTAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGaa      <  1:2662970/69‑2 (MQ=255)                                                       cTGGACCCTAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGAc      <  1:3044734/69‑1 (MQ=255)                                                       cTGGACCCTAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGAc      <  1:2420071/69‑1 (MQ=255)                                                       cTGGACCCTAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGAc      <  1:3201244/69‑1 (MQ=255)                                                       cTGGACCCTAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGAc      <  1:1753786/69‑1 (MQ=255)                                                        tGGACCCTAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGAc      <  1:1625449/68‑1 (MQ=255)                                                            cccTAACGCTGATATCGATGCCTTACCCGACGGCATATACGCCCAGGTGGCAGTTTACTCCGACCAtt  <  1:1314753/68‑1 (MQ=255)                                                            cccTAACGCTGATATCGATGCCTTACCCGACGGCATATACGCCCAGGTGGCAGTTTACTCCGACCAtt  <  1:1870297/68‑1 (MQ=255)                                                            cccTAACGCTGATATCGATGCCTTACCCGACGGCATATACGCCCAGGTGGCAGTTTACTCCGACCAtt  <  1:2675242/68‑1 (MQ=255)                                                            cccTAACGCTGATATCGATGCCTTACCCGACGGCATATACGCCCAGGTGGCAGTTTACTCCGACCAtt  <  1:174407/68‑1 (MQ=255)                                                            cccTAACGCTGATATCGATGCCTTACCCGACGGCATATACGCCCAGGTGGCAGTTTACTCCGACCAtt  <  1:1405274/68‑1 (MQ=255)                                                            cccTAACGCTGATATCGATGCCTTACCCGACGGCATATACGCCCAGGTGGCAGTTTACTCCGACCAtt  <  1:170362/68‑1 (MQ=255)                                                            cccTAACGCTGATATCGATGCCTTACCCGACGGCATATACGCCCAGGTGGCAGTTTACTCCGACCAtt  <  1:2054277/68‑1 (MQ=255)                                                            cccTAACGCTGATATCGATGCCTTACCCGACGGCATATACGCCCAGGTGGCAGTTTACTCCGACCAtt  <  1:375940/68‑1 (MQ=255)                                                            cccTAACGCTGATATCGATGCCTTACCCGACGGCATATACGCCCAGGTGGCAGTTTACTCCGACCAtt  <  1:3393745/68‑1 (MQ=255)                                                            cccTAACGCTGATATCGATGCCTTACCCGACGGCATATACGCCCAGGTGGCAGTTTACTCCGACCAtt  <  1:3332440/68‑1 (MQ=255)                                                            cccTAACGCTGATATCGATGCCTTACCCGACGGCATATACGCCCAGGTGGCAGTTTACTCCGACCAtt  <  1:3282956/68‑1 (MQ=255)                                                            cccTAACGCTGATATCGATGCCTTACCCGACGGCATATACGCCCAGGTGGCAGTTTACTCCGACCAtt  <  1:2120203/68‑1 (MQ=255)                                                            cccTAACGCTGATATCGATGCCTTACCCGACGGCATATACGCCCAGGTGGCAGTTTACTCCGACCAtt  <  1:2146007/68‑1 (MQ=255)                                                            cccTAACGCTGATATCGATGCCTTACCCGACGGCATATACGCCCAGGTGGCAGTTTACTCCGACCAtt  <  1:3235239/68‑1 (MQ=255)                                                            cccTAACGCTGATATCGATGCCTTACCCGACGGCATATACGCCCAGGTGGCAGTTTACTCCGACCAtt  <  1:2232817/68‑1 (MQ=255)                                                            cccTAACGCTGATATCGATGCCTTACCCGACGGCATATACGCCCAGGTGGCAGTTTACTCCGACCAtt  <  1:3213657/68‑1 (MQ=255)                                                            cccTAACGCTGATATCGATGCCTTACCCGACGGCATATACGCCCAGGTGGCAGTTTACTCCGACCAtt  <  1:2282167/68‑1 (MQ=255)                                                            cccTAACGCTGATATCGATGCCTTACCCGACGGCATATACGCCCAGGTGGCAGTTTACTCCGACCAtt  <  1:3166402/68‑1 (MQ=255)                                                            cccTAACGCTGATATCGATGCCTTACCCGACGGCATATACGCCCAGGTGGCAGTTTACTCCGACCAtt  <  1:3155664/68‑1 (MQ=255)                                                            cccTAACGCTGATATCGATGCCTTACCCGACGGCATATACGCCCAGGTGGCAGTTTACTCCGACCAtt  <  1:2372135/68‑1 (MQ=255)                                                            cccTAACGCTGATATCGATGCCTTACCCGACGGCATATACGCCCAGGTGGCAGTTTACTCCGACCAtt  <  1:2568138/68‑1 (MQ=255)                                                            cccTAACGCTGATATCGATGCCTTACCCGACGGCATATACGCCCAGGTGGCAGTTTACTCCGACCAtt  <  1:1616423/68‑1 (MQ=255)                                                            cccTAACGCTGATATCGATGCCTTACCCGACGGCATATACGCCCAGGTGGCAGTTTACTCCGACCAtt  <  1:2806855/68‑1 (MQ=255)                                                            cccTAACGCTGATATCGATGCCTTACCCGACGGCATATACGCCCAGGTGGCAGTTTACTCCGACCAtt  <  1:2788928/68‑1 (MQ=255)                                                                    |                                                            TTGCAATCATTAACGGTCGTGCCCGGCACGGACGGTGTGCTGGGAACCATTGAACTGGACCCTAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATT  >  minE/2571029‑2571155

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