Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	869,449	T→A	33.3%	T58T (ACA→ACT)	yciN ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	869,449	0	T	A	33.3%	88.2 / 70.5	93	T58T (ACA→ACT)	yciN	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (27/35);  new base A (15/16);  total (42/51)
Fisher's exact test for biased strand distribution p-value = 6.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.72e-01

AGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTAAAAAGTAGTCGCCGGTAAATACCAGCACGC  >  minE/869381‑869501                                                                     |                                                     atatgatATTTTTCTTAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCt                                                      <  1:2609807/67‑1 (MQ=255) aGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAATATATTAAACACCGCGGTGCTTTTGGCt                                                      <  1:1600314/69‑1 (MQ=255) aGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCt                                                      <  1:2665216/69‑1 (MQ=255)   gtgatATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCAGt                                                    <  1:2281396/68‑1 (MQ=255)   gtgatATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCAGt                                                    <  1:2421942/68‑1 (MQ=255)   gtgatATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCAGt                                                    <  1:1368089/68‑1 (MQ=255)   gtgatATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCAGt                                                    <  1:135173/68‑1 (MQ=255)   gtgatATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCAGt                                                    <  1:2943197/68‑1 (MQ=255)   gtgatATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCAGt                                                    <  1:3151016/68‑1 (MQ=255)   gtgatATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCAGt                                                    <  1:2170479/68‑1 (MQ=255)   gtgatATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCAGt                                                    <  1:3420861/68‑1 (MQ=255)   gtgatATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCAGt                                                    <  1:440171/68‑1 (MQ=255)   gtgatATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCAGt                                                    <  1:456493/68‑1 (MQ=255)   gtgatATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCAGt                                                    <  1:977200/68‑1 (MQ=255)   gtgatATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCAGt                                                    <  1:2153259/68‑1 (MQ=255)   gtgatATTTTTCTAAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCAGt                                                    <  1:1577174/68‑1 (MQ=255)   gtgatATTTATCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCAGt                                                    <  1:80878/68‑1 (MQ=255)    tgatATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCAGt                                                    <  1:1777644/68‑1 (MQ=255)     gatATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCAGt                                                    <  1:1663102/67‑1 (MQ=255)                            gCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCt                                                      <  1:2942652/42‑1 (MQ=255)                                 gtgtTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTAAAAAGTAGTc                      >  1:1043630/1‑69 (MQ=255)                                  tgtTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTAAAAAGTAGTCg                     >  1:251224/1‑69 (MQ=255)                                  tgtTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTAAAAAGTAGTCg                     >  1:2956517/1‑69 (MQ=255)                                  tgtTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTAAAAAGTAGTCg                     >  1:3287090/1‑69 (MQ=255)                                  tgtTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTAAAAAGTAGTCg                     >  1:2081097/1‑69 (MQ=255)                                  tgtTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTAAAAAGTAGTCg                     >  1:1078983/1‑69 (MQ=255)                                  tgtTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTAAAAAGTAGTCg                     >  1:1304604/1‑69 (MQ=255)                                      tGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTaaaaatt                          <  1:1724163/60‑3 (MQ=255)                                      tGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTa                                <  1:2695767/54‑1 (MQ=255)                                      tGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTa                                >  1:3431370/1‑54 (MQ=255)                                      tGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTa                                >  1:2486252/1‑54 (MQ=255)                                      tGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTa                                <  1:2792024/54‑1 (MQ=255)                                      tGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTa                                <  1:2346674/54‑1 (MQ=255)                                      tGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTa                                >  1:3380750/1‑54 (MQ=255)                                      tGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTa                                >  1:527434/1‑54 (MQ=255)                                      tGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTa                                <  1:1996800/54‑1 (MQ=255)                                      tGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTa                                >  1:1937560/1‑54 (MQ=255)                                      tGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTa                                <  1:1216259/54‑1 (MQ=255)                                      tGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTa                                >  1:3001904/1‑54 (MQ=255)                                      tGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTa                                >  1:1834157/1‑54 (MQ=255)                                      tGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTAAAAagt                          <  1:2019434/60‑1 (MQ=255)                                      tGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTAAAAagt                          <  1:2772463/60‑1 (MQ=255)                                      tGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTAAAAagt                          <  1:1763922/60‑1 (MQ=255)                                      tGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTAAAAagt                          <  1:61402/60‑1 (MQ=255)                                      tGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTAAAAagt                          <  1:525185/60‑1 (MQ=255)                                      tGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTAAAAagt                          <  1:2533742/60‑1 (MQ=255)                                      tGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTAAAAagt                          <  1:2219155/60‑1 (MQ=255)                                      tGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTAAAAagt                          <  1:1436396/60‑1 (MQ=255)                                      tGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTAAAAagt                          <  1:2784841/60‑1 (MQ=255)                                      tGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTAAAAAGTAGTCGCCg                  >  1:1904153/1‑68 (MQ=255)                                      tGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTAAAAAGTAGTCGCCg                  >  1:3263192/1‑68 (MQ=255)                                      tGAACATATTAAACACCGCGGTGCTTTTGGCAGTTTGTAATCCTTGTTCATCTAAAAagta                         >  1:2985217/1‑61 (MQ=255)                                      tGAACATATTAAACACCGCGGTGCTTTTGGCAGTTGGTAATCCTTGTTCATCTAAAAagta                         >  1:3118441/1‑61 (MQ=255)                                      tGAACATATTAAACACCGCGGTGCTTTTGGCAGTTGGTAATCCTTGTTCATCTAAAAagta                         >  1:2935950/1‑61 (MQ=255)                                      tGAACATATTAAACACCGCGGTGCTTTTGGCAGTTGGTAATCCTTGTTCATCTAAAAagta                         >  1:1310409/1‑61 (MQ=255)                                      tGAACATATTAAACACCGCGGTGCTTTTGGCAGTTGGTAATCCTTGTTCATCTAAAAagta                         >  1:1270133/1‑61 (MQ=255)                                      tGAACATATTAAACACCGCGGTGCTTTTGGCAGTTGGTAATCCTTGTTCATCTAAAAagta                         >  1:3191846/1‑61 (MQ=255)                                      tGAACATATTAAACACCGCGGTGCTTTTGGCAGTTGGTAATCCTTGTTCATCTAAAAagta                         >  1:353891/1‑61 (MQ=255)                                      tGAACATATTAAACACCGCGGTGCTTTTGGCAGTTGGTAATCCTTGTTCATCTAAAAagta                         >  1:1554745/1‑61 (MQ=255)                                      tGAACATATTAAACACCGCGGTGCTTTTGGCAGTTGGTAATCCTTGTTCATCTAAAAagta                         >  1:124759/1‑61 (MQ=255)                                      tGAACATATTAAACACCGCGGTGCTTTTGGCAGTTGGTAATCCTTGTTCATCTAAAAagta                         >  1:1242459/1‑61 (MQ=255)                                      tGAACATATTAAACACCGCGGTGCTTTTGGCAGTTGGTAATCCTTGTTCATCTAAAAagta                         >  1:2187092/1‑61 (MQ=255)                                      tGAACATATTAAACACCGCGGTGCTTTTGGCAGTTGGTAATCCTTGTTCATCTAAAAagta                         >  1:2097284/1‑61 (MQ=255)                                      tGAACATATTAAACACCGCGGTGCTTTTGGCAGTTGGTAATCCTTGTTCATCTAAAAagta                         >  1:2071971/1‑61 (MQ=255)                                      tGAACATATTAAACACCGCGGTGCTTTTGGCAGTTGGTAATCCTTGTTCATCTAAAAagta                         >  1:720381/1‑61 (MQ=255)                                      tGAACATATTAAACACCGCGGTGCTTTTGGCAGTTGGTAATCCTTGTTCATCTAAAAagta                         >  1:1842757/1‑61 (MQ=255)                                      tGAACATATTAAACACCGAGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTa                                >  1:1264201/1‑54 (MQ=255)                                         aCATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTAAAAAGTAGTCTCCGGTa               <  1:3272/68‑1 (MQ=255)                                         aCATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTAAAAAGTAGTCGCCGGTa               <  1:3389822/68‑1 (MQ=255)                                         aCATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTAAAAAGTAGTCGCCGGTa               <  1:574899/68‑1 (MQ=255)                                         aCATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTAAAAAGTAGTCGCCGGTa               <  1:63669/68‑1 (MQ=255)                                         aCATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTAAAAAGTAGTCGCCGGTa               <  1:3264472/68‑1 (MQ=255)                                         aCATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTAAAAAGTAGTCGCCGGTa               <  1:1796072/68‑1 (MQ=255)                                         aCATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTAAAAAGTAGTCGCCGGTa               <  1:3152350/68‑1 (MQ=255)                                         aCATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTAAAAAGTAGTCGCCGGTa               <  1:1286120/68‑1 (MQ=255)                                         aCATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTAAAAAGTAGTCGCCGGTa               <  1:3064270/68‑1 (MQ=255)                                         aCATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTAAAAAGTAGTCGCCGGTa               <  1:1884327/68‑1 (MQ=255)                                         aCATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTAAAAAGTAGTCGCCGGTa               <  1:135252/68‑1 (MQ=255)                                         aCATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTAAAAAGTAGTCGCCGGTa               <  1:2495774/68‑1 (MQ=255)                                         aCATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTAAAAAGTAGTCGCCGGTa               <  1:1653344/68‑1 (MQ=255)                                          cATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTAAAAAGTAGTCGCCGGTa               <  1:1848398/67‑1 (MQ=255)                                            tatTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTAAAAAGTAGTCGCCGGTa               <  1:1219429/65‑1 (MQ=255)                                             atTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTAAAAAGTAGTCGCCGGTAGAt            >  1:943436/1‑67 (MQ=255)                                             atTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTAAAAAGTAGTCGCCGGTAAAt            >  1:3333754/1‑67 (MQ=255)                                             atTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTAAAAAGTAGTCGCCGGTAAAt            >  1:2885144/1‑67 (MQ=255)                                             atTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTAAAAAGTAGTCGCCGGTAAAt            >  1:1675085/1‑67 (MQ=255)                                             atTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTAAAAAGTAGTCGCCGGTAAAt            >  1:220917/1‑67 (MQ=255)                                                aaaCACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTAAAAAGTAGTCGCCGGTAAATACCAg       <  1:1908038/69‑1 (MQ=255)                                                aaaCACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTAAAAAGTAGTCGCCGGTAAATACCAg       <  1:1332116/69‑1 (MQ=255)                                                 aaCACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTa                                <  1:2000827/43‑1 (MQ=255)                                                   cacCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTAAAAAGTAGTCGCCGGTa               >  1:1795159/1‑58 (MQ=255)                                                   cacCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTAAAAAGTAGTCGCCGGTa               >  1:1933677/1‑58 (MQ=255)                                                   cacCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTAAAAAGTAGTCGCCGGTa               >  1:1300728/1‑58 (MQ=255)                                                   cacCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTAAAAAGTAGTCGCCGGTa               >  1:3112032/1‑58 (MQ=255)                                                   cacCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTAAAAAGTAGTCGCCGGTa               >  1:3184723/1‑58 (MQ=255)                                                     ccGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTAAAAAGTAGTCGCCGGTAAATACCAGCACGc  <  1:3080970/69‑1 (MQ=255)                                                     ccGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTAAAAAGTAGGCGCCGGTAAATACCAGCACGc  <  1:2649947/69‑1 (MQ=255)                                                                     |                                                     AGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTAAAAAGTAGTCGCCGGTAAATACCAGCACGC  >  minE/869381‑869501

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