Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,006,546	G→A	23.2%	A143A (GCG→GCA)	ydgH →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,006,546	0	G	A	23.2%	103.5 / 22.3	77	A143A (GCG→GCA)	ydgH	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (27/32);  new base A (13/5);  total (40/37)
Fisher's exact test for biased strand distribution p-value = 6.20e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.99e-01

CACGGTCACCGTCCAGGGCTTCTATCGTAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACATCGTTCGTCAAATCGATGCCAACCAGGGCGGCAAC  >  minE/1006483‑1006612                                                                |                                                                   cacGGTCACCGTCCATGGCTTCTATCGTAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCaa                                                                <  1:1906550/68‑1 (MQ=255) cacGGTCACCGTCCAGGGCTTCTATCGTAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCaa                                                                <  1:1677444/68‑1 (MQ=255) cacGGTCACCGTCCAGGGCTTCTATCGTAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCaa                                                                <  1:1537558/68‑1 (MQ=255) cacGGTCACCGTCCAGGGCTTCTATCGTAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCaa                                                                <  1:127747/68‑1 (MQ=255) cacGGTCACCGTCCAGGGCTTCTATCGTAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCaa                                                                <  1:3285404/68‑1 (MQ=255) cacGGTCACCGTCCAGGGCTTCTATCGTAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCaa                                                                <  1:216618/68‑1 (MQ=255)  acGGTCACCGTCCAGGGCTTCTATCGTAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCaaa                                                               <  1:2195944/68‑1 (MQ=255)  acGGTCACCGTCCAGGGCTTCTATCGTAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCaaa                                                               <  1:600724/68‑1 (MQ=255)  acGGTCACCGTCCAGGGCTTCTATCGTAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCaaa                                                               <  1:555331/68‑1 (MQ=255)  acGGTCACCGTCCAGGGCTTCTATCGTAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCaaa                                                               <  1:3381754/68‑1 (MQ=255)  acGGTCACCGTCCAGGGCTTCTATCGTAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCaaa                                                               <  1:2828369/68‑1 (MQ=255)  acGGTCACCGTCCAGGGCTTCTATCGTAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCaaa                                                               <  1:2265946/68‑1 (MQ=255)  acGGTCACCGTCCAGGGCTTCTATCGTAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCaaa                                                               <  1:3271870/68‑1 (MQ=255)  acGGTCACCGTCCAGGGCTTCTATCGTAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCaaa                                                               <  1:3085573/68‑1 (MQ=255)  acGGTCACCGTCCAGGGCTTCTATCGTAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCaaa                                                               <  1:3041672/68‑1 (MQ=255)  acGGTCACCGTCCAGGGCTTCTATCGTAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCaaa                                                               <  1:2643622/68‑1 (MQ=255)  acGGTCACCGTCCAGGGCTTCTATCGTAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCaaa                                                               <  1:2771380/68‑1 (MQ=255)  acGGTCACCGTCCAGGGCTTCTATCGTAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCaaa                                                               <  1:2806756/68‑1 (MQ=255)  acGGTCACCGTCCAGGGCTTCTATCGTAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCaaa                                                               <  1:1663216/68‑1 (MQ=255)  acGGTCACCGTCCAGGGCTTCTATCGTAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCaaa                                                               <  1:1358182/68‑1 (MQ=255)  acGGTCACCGTCCAGGGCTTCTATCGTAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCaaa                                                               <  1:1755807/68‑1 (MQ=255)  acGGTCACCGTCCAGGGCTTCTAACGTAGCCAGCCAGAAGCCAATGATGCCATCACCAAAGCGGCaaa                                                               <  1:3150841/68‑1 (MQ=255)       cgccGTCCAGGGCTTCTATCGTAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGaa                                                         <  1:138254/67‑1 (MQ=255)       cACCGTCCAGGGCTTCTATCGTAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGaa                                                         <  1:119237/69‑1 (MQ=255)       cACCGTCCAGGGCTTCTATCGTAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGaa                                                         <  1:216147/69‑1 (MQ=255)       cACCGTCCAGGGCTTCTATCGTAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGaa                                                         <  1:870710/69‑1 (MQ=255)             ccAGGGCTTCTATCGTAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAGGTGc                                                   <  1:1231775/69‑1 (MQ=255)              cAGGGCTTCTATCGTAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAGGTGcc                                                  <  1:1435646/69‑1 (MQ=255)              cAGGGCTTCTATCGTAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAGGTGc                                                   <  1:1475187/68‑1 (MQ=255)                      cTATCGTAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCAGCaaaa                                                              >  1:2355869/1‑49 (MQ=255)                      cTATCGTAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCAGCAAAAGCGAAAGGTGCCTActct                                            >  1:3031753/1‑67 (MQ=255)                      cTATCGTAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCAGCAAAAGCGAAAGGTGCCTACTCtt                                           >  1:1104879/1‑68 (MQ=255)                      cTATCGTAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCAGCAAAAGCGAAAGGTGCCTACTCtt                                           >  1:3024311/1‑68 (MQ=255)                      cTATCGTAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCAGCAAAAGCGAAAGGTGCCTACTCtt                                           >  1:2726339/1‑68 (MQ=255)                      cTATCGTAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCAGCAAAAGCGAAAGGTGCCTACTCtt                                           >  1:2543771/1‑68 (MQ=255)                      cTATCGTAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCAGCAAAAGCGAAAGGTGCCTACTCtt                                           >  1:2245072/1‑68 (MQ=255)                      cTATCGTAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCAGCAAAAGCGAAAGGTGCCTACTCtt                                           >  1:63316/1‑68 (MQ=255)                      cTATCGTAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCAGCAAAAGCGAAAGGTGCCTACTCtt                                           >  1:214933/1‑68 (MQ=255)                      cTATCGTAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCAGCAAAAGCGAAAGGTGCCTACTCtt                                           >  1:3382216/1‑68 (MQ=255)                      cTATCGTAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCAGCAAAAGCGAAAGGTGCCTACTCtt                                           >  1:546470/1‑68 (MQ=255)                      cTATCGTAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCAGCAAAAGCGAAAGGTGCCTACTCtt                                           >  1:1848231/1‑68 (MQ=255)                      cTATCGTAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCAGCAAAAGCGAAAGGTGCCTACTCtt                                           >  1:1814440/1‑68 (MQ=255)                         tCGTAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTa                                       >  1:2660450/1‑69 (MQ=255)                         tCGTAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTa                                       >  1:1154075/1‑69 (MQ=255)                            tAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCaa                                                                <  1:1494994/41‑1 (MQ=255)                                           cAATGATGCCATCACCAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACATCGTTCGTCAAATCGa                     >  1:909654/1‑69 (MQ=255)                                                  gCCATCACCAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACATCGTTCGTCAAATCGATGCCa                >  1:616259/1‑67 (MQ=255)                                                       cACCAAAGGGGCAAAAGCGAAAGGTGCCTACTCTTTCTACATCGTTCGTCAAATCGa                     >  1:955674/1‑57 (MQ=255)                                                       cACCAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACATCTTTCGTCAAATCGa                     >  1:608649/1‑57 (MQ=255)                                                       cACCAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACATCGTTCGTCAAATCGa                     >  1:2763832/1‑57 (MQ=255)                                                       cACCAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACATCGTTCGTCAAATCGa                     >  1:754755/1‑57 (MQ=255)                                                       cACCAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACATCGTTCGTCAAATCGa                     >  1:1200018/1‑57 (MQ=255)                                                       cACCAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACATCGTTCGTCAAATCGa                     >  1:2065181/1‑57 (MQ=255)                                                       cACCAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACATCGTTCGTCAAATCGa                     >  1:226609/1‑57 (MQ=255)                                                       cACCAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACATCGTTCGTCAAATCGa                     >  1:3095214/1‑57 (MQ=255)                                                       cACCAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACATCGTTCGTCAAATCGa                     >  1:2560186/1‑57 (MQ=255)                                                       cACCAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACATCGTTCGTCAAATCGa                     >  1:2580400/1‑57 (MQ=255)                                                       cACCAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACATCGTTCGTCAAATCGa                     >  1:2865785/1‑57 (MQ=255)                                                       cACCAAAGCGGCAAAAGCGAAAGGTGCATACTCTTTCTACATCGTTCGTCAAATCGa                     >  1:716682/1‑57 (MQ=255)                                                        aCCAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACATCGTTCGTCAAATCGATGCCAACCAGggc        >  1:743224/1‑69 (MQ=255)                                                        aCCAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACATCGTTCGTCAAATCGATGCCAACCAGggc        >  1:195396/1‑69 (MQ=255)                                                        aCCAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACATCGTTCGTCAAATCGATGCCAACCAGggc        >  1:1855517/1‑69 (MQ=255)                                                        aCCAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACATCGTTCGTCAAATCGATGCCAACCAGggc        >  1:1764565/1‑69 (MQ=255)                                                          cAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACAt                                    >  1:1435345/1‑39 (MQ=255)                                                          cAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACATCGTTCGTCAAATCGATGCCAAc              >  1:1286130/1‑61 (MQ=255)                                                          cAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACATCGTTCGTCAAATCGATGCCAAc              >  1:3366599/1‑61 (MQ=255)                                                          cAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACATCGTTCGTCAAATCGATGCCAAc              >  1:2562985/1‑61 (MQ=255)                                                          cAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACATCGTTCGTCAAATCGATGCCAAc              >  1:3338812/1‑61 (MQ=255)                                                          cAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACATCGTTCGTCAAATCGATGCCAACCAGggcg       >  1:768315/1‑68 (MQ=255)                                                          cAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACATCGTTCGTCAAATCGATGCCAACCAGggcg       >  1:243640/1‑68 (MQ=255)                                                            aaGCAGCAAAAGCGAAAGGTGCCTACTCTTTCTACATCGTTCGTCAAATCGATGCCAACCAGGGCGGCa    <  1:10547/69‑1 (MQ=255)                                                            aaGCAGCAAAAGCGAAAGGTGCCTACTCTTTCTACATCGTTCGTCAAATCGATGCCAACCAGGGCGGCa    <  1:1836254/69‑1 (MQ=255)                                                            aaGCAGCAAAAGCGAAAGGTGCCTACTCTTTCTACATCGTTCGTCAAATCGATGCCAACCAGGGCGGCa    <  1:1936032/69‑1 (MQ=255)                                                            aaGCAGCAAAAGCGAAAGGTGCCTACTCTTTCTACATCGTTCGTCAAATCGATGCCAACCAGGGCGGCa    <  1:3359805/69‑1 (MQ=255)                                                            aaGCAGCAAAAGCGAAAGGTGCCTACTCTTTCTACATCGTTCGTCAAATCGATGCCAACCAGGGCGGCa    <  1:3120839/69‑1 (MQ=255)                                                              gCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACATCGTTCGTCAAATCGATGCCAACCAGGGCGGCAAc  <  1:1633477/69‑1 (MQ=255)                                                              gCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACATCGTTCGTCAAATCGATGCCAACCAGGGCGGCAAc  <  1:809225/69‑1 (MQ=255)                                                                |                                                                   CACGGTCACCGTCCAGGGCTTCTATCGTAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACATCGTTCGTCAAATCGATGCCAACCAGGGCGGCAAC  >  minE/1006483‑1006612

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