Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	125,165	A→G	41.5%	G267G (GGT→GGC)	gcd ←	glucose dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	125,165	0	A	G	41.5%	27.2 / 57.7	81	G267G (GGT→GGC)	gcd	glucose dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (11/36);  new base G (18/15);  total (30/51)
Fisher's exact test for biased strand distribution p-value = 8.73e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.45e-01

ACGCGGGCAATCCGCCATCACTTCCGGCGAAGCGGTTTCTGCTTTGGCTTCATGATAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCATTTCTCTTT  >  minE/125102‑125233                                                                |                                                                     acgCGGGCAATCCGCCATCACTTCCGGCGAAGCGGTTTCTGCTTTGGCTTCATGATAAGAGACACCACg                                                                 <  1:3251571/69‑1 (MQ=255)   gcgGGAAATCCGCCATCACTTCCGGCGAAGCGGTTTCTGCTTTGGCTTCATGATAAGAGACACCACGGc                                                               <  1:2645657/69‑1 (MQ=255)   gcgGGAAATCCGCCATCACTTCCGGCGAAGCGGTTTCTGCTTTGGCTTCATGATAAGAGACACCACGGc                                                               <  1:3426939/69‑1 (MQ=255)                   cACTTCCGGCGAAGCGGTTTCTGCTTTGGCTTCATGATAAGAGACACCACGGCAGGTTa                                                         <  1:2257107/59‑1 (MQ=255)                    aCTTCCGGCGAAGCGGTTTCTGCTTTGGCTTCATGATAAGAGACACCACGGCAGGTTACGTGCTGGaaa                                              >  1:1923513/1‑69 (MQ=255)                    aCTTCCGGCGAAGCGGTTTCTGCTTTGGCTTCATGATAAGAGACACCACGGCAGGTTACGTGCTGGaaa                                              >  1:2828240/1‑69 (MQ=255)                    aCTTCCGGCGAAGCGGTTTCTGCTTTGGCTTCATGATAAGAGACACCACGGCAGGTTACGTGCTGGaaa                                              >  1:1441309/1‑69 (MQ=255)                    aCTTCCGGCGAAGCGGTTTCTGCTTTGGCTTCATGATAAGAGACACCACGGCAGGTTACGTGCTGGaaa                                              >  1:690617/1‑69 (MQ=255)                    aCTTCCGGCGAAGCGGTTTCTGCTTTGGCTTCATGATAAGAGACACCACGGCAGGTTACGTGCTGGaaa                                              >  1:602357/1‑69 (MQ=255)                    aCTTCCGGCGAAGCGGTTTCTGCTTTGGCTTCATGATAAGAGACACCACGGCAGGTTACGTGCTGGaaa                                              >  1:1713968/1‑69 (MQ=255)                    aCTTCCGGCGAAGCGGTTTCTGCTTTGGCTTCATGATAAGAGACACCACGGCAGGTTACGTGCTGGaaa                                              >  1:3279672/1‑69 (MQ=255)                    aCTTCCGGCGAAGCGGTTTCTGCTTTGGCTTCATGATAAGAGACACCACGGCAGGTTACGTGCTGGaaa                                              >  1:1994441/1‑69 (MQ=255)                    aCTTCCGGCGAAGCGGTTTCTGCTTTGGCTTCATGATAAGAGACACCACGGCAGGTTACGTGCTGGaa                                               >  1:2847457/1‑68 (MQ=255)                    aCTTCCGGCGAAGCAGTTTCTGCTTTGGCTTCATGATAAGAGACGCCACGGCAGGTTACATGCTGGaag                                              >  1:2896890/1‑68 (MQ=255)                    aCTTCCGGCGAAGCAGTTTCTGCTTTGGCTTCATGATAAGAGACGCCACGGCAGGTTACATGCTGGaag                                              >  1:3083329/1‑68 (MQ=255)                    aCTTCCGGCGAAGCAGTTTCTGCTTTGGCTTCATGATAAGAGACGCCACGGCAGGTTACATGCTGGaag                                              >  1:1946691/1‑68 (MQ=255)                    aCTTCCGGCGAAGCAGTTTCTGCTTTGGCTTCATGATAAGAGACGCCACGGCAGGTTACATGCTGGaag                                              >  1:2992560/1‑68 (MQ=255)                    aCTTCCGGCGAAGCAGTTTCTGCTTTGGCTTCATGATAAGAGACGCCACGGCAGGTTACATGCTGGaag                                              >  1:2066344/1‑68 (MQ=255)                    aCTTCCGGCGAAGCAGTTTCTGCTTTGGCTTCATGATAAGAGACGCCACGGCAGGTTACATGCTGGaag                                              >  1:1058307/1‑68 (MQ=255)                    aCTTCCGGCGAAGCAGTTTCTGCTTTGGCTTCATGATAAGAGACGCCACGGCAGGTTACATGCTGGaag                                              >  1:2624873/1‑68 (MQ=255)                    aCTTCCGGCGAAGCAGTTTCTGCTTTGGCTTCATGATAAGAGACGCCACGGCAGGTTACATGCTGGaag                                              >  1:2541583/1‑68 (MQ=255)                    aCTTCCGGCGAAGCAGTTTCTGCTTTGGCTTCATGATAAGAGACGCCACGGCAGGTTACATGCTGGaag                                              >  1:2463468/1‑68 (MQ=255)                    aCTTCCGGCGAAGCAGTTTCTGCTTTGGCTTCATGATAAGAGACGCCACGGCAGGTTACATGCTGGaag                                              >  1:1910709/1‑68 (MQ=255)                    aCTTCCGGCGAAGCAGTTTCTGCTTTGGCTTCATGATAAGAGACGCCACGGCAGGTTACATGCTGGaag                                              >  1:324447/1‑68 (MQ=255)                    aCTTCCGGCGAAGCAGTTTCTGCTTTGGCTTCATGATAAGAGACGCCACGGCAGGTTACATGCTGGaag                                              >  1:3350012/1‑68 (MQ=255)                    aCTTCCGGCGAAGCAGTTTCTGCTTTGGCTTCATGATAAGAGACGCCACGGCAGGTTACATGCTGGaag                                              >  1:337114/1‑68 (MQ=255)                    aCTTCCGGCGAAGCAGTTTCTGCTTTGGCTTCATGATAAGAGACGCCACGGCAGGTTACATGCTGGaag                                              >  1:588566/1‑68 (MQ=255)                    aCTTCCGGCGAAGCAGTTTCTGCTTTGGCTTCATGATAAGAGACGCCACGGCAGGTTACATGCTGGaag                                              >  1:134501/1‑68 (MQ=255)                    aCTTCCGGCGAAGCAGTTTCTGCTTTGGCTTCATGATAAGAGACGCCACGGCAGGTTACATGCTGGaag                                              >  1:774835/1‑68 (MQ=255)                    aCTTCCGGCGAAGCAGTTTCTGCTTTGGCTTCATGATAAGAGACGCCACGGCAGGTTACATGCTGGaag                                              >  1:1263375/1‑68 (MQ=255)                    aCTTCCGGCGAAGCAGTTTCTGCTTTGGCTTCATGATAAGAGACGCCACGGCAGGTTACATGCTGGaa                                               >  1:2860879/1‑68 (MQ=255)                                gCGGTTTCTGCTTTGGCTTCATGATAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGtt                                  >  1:3145991/1‑69 (MQ=255)                                     ttCTGCTTTGGCTTCATGATAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCa                              >  1:164464/1‑68 (MQ=255)                                     ttCTGCTTTGGCTTCATGATAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAg                             <  1:2585977/69‑1 (MQ=255)                                     ttCTGCTTTGGCTTCATGATAAGAGACACCAAGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAg                             <  1:2699784/69‑1 (MQ=255)                                     ttCTGCTTTGGCTTCATGATAAGAGA‑‑‑CACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTc                          >  1:822324/1‑69 (MQ=255)                                         gCTTTGGCTTCATGATAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTc                          <  1:3144600/68‑1 (MQ=255)                                         gCTTTGGCTTCATGATAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTc                          <  1:1362334/68‑1 (MQ=255)                                         gCTTTGGCTTCATGATAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTc                          <  1:1452121/68‑1 (MQ=255)                                          cTTTGGCTTCATGATAAGAGACACCACGGCAGGTTACGTGCTGGa                                                <  1:119935/45‑1 (MQ=255)                                            ttGGCTTCATGATAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAgg                       <  1:3251234/68‑1 (MQ=255)                                            ttGGCTTCATGATAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAgg                       <  1:2360497/68‑1 (MQ=255)                                            ttGGCTTCATGATAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAgg                       <  1:319048/68‑1 (MQ=255)                                            ttGGCTTCATGATAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAgg                       <  1:2743821/68‑1 (MQ=255)                                            ttGGCTTCATGATAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAgg                       <  1:3393325/68‑1 (MQ=255)                                            ttGGCTTCATGATAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAgg                       <  1:3048684/68‑1 (MQ=255)                                            ttGGCTTCATGATAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAgg                       <  1:2064846/68‑1 (MQ=255)                                            ttGGCTTCATGATAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAgg                       <  1:212890/68‑1 (MQ=255)                                                 ttCATGATAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGt                  <  1:1800757/68‑1 (MQ=255)                                                 ttCATGATAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGt                  <  1:1852745/68‑1 (MQ=255)                                                 ttCATGATAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGt                  <  1:1797281/68‑1 (MQ=255)                                                 ttCATGATAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGt                  <  1:3168844/68‑1 (MQ=255)                                                 ttCATGATAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGt                  <  1:1753021/68‑1 (MQ=255)                                                 ttCATGATAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGt                  <  1:1505843/68‑1 (MQ=255)                                                 ttCATGATAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGt                  <  1:736960/68‑1 (MQ=255)                                                 ttCATGATAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGt                  <  1:2217404/68‑1 (MQ=255)                                                 ttCATGATAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGt                  <  1:2788455/68‑1 (MQ=255)                                                 ttCATGATAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGt                  <  1:563938/68‑1 (MQ=255)                                                 ttCATGATAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGt                  <  1:2586276/68‑1 (MQ=255)                                                 ttCATGATAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGt                  <  1:1734844/68‑1 (MQ=255)                                                       aTAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCa           <  1:2970460/69‑1 (MQ=255)                                                       aTAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCa           <  1:1746104/69‑1 (MQ=255)                                                       aTAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCa           <  1:519770/69‑1 (MQ=255)                                                        tAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGt                  <  1:2239916/61‑1 (MQ=255)                                                          agagACGCCACGGCAGGTTACATGCTGGAAGGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCAt          <  1:3433508/67‑1 (MQ=255)                                                          agagACGCCACGGCAGGTTACATGCTGGAAGGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCAt          <  1:411243/67‑1 (MQ=255)                                                          agagACGCCACGGCAGGTTACATGCTGGAAGGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCAt          <  1:424362/67‑1 (MQ=255)                                                          agagACGCCACGGCAGGTTACATGCTGGAAGGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCAt          <  1:1526929/67‑1 (MQ=255)                                                          agagACGCCACGGCAGGTTACATGCTGGAAGGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCAt          <  1:1881451/67‑1 (MQ=255)                                                          agagACGCCACGGCAGGTTACATGCTGGAAGGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCAt          <  1:3176260/67‑1 (MQ=255)                                                          agagACGCCACGGCAGGTTACATGCTGGAAGGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCAt          <  1:2986548/67‑1 (MQ=255)                                                          agagACGCCACGGCAGGTTACATGCTGGAAGGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCAt          <  1:719336/67‑1 (MQ=255)                                                          agagACGCCACGGCAGGTTACATGCTGGAAGGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCAt          <  1:1322401/67‑1 (MQ=255)                                                          agagACGCCACGGCAGGTTACATGCTGGAAGGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCAt          <  1:2704519/67‑1 (MQ=255)                                                          agagACGCCACGGCAGGTTACATGCTGGAAGGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCAt          <  1:2266631/67‑1 (MQ=255)                                                          agagACGCCACGGCAGGTTACATGCTGGAAGGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCAt          <  1:2272983/67‑1 (MQ=255)                                                          agagACGCCACGGCAGGTTACATGCTGGAAGGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCAt          <  1:2431273/67‑1 (MQ=255)                                                          agagACGCCACGGCAGGTTACATGCTGGAAGGACTCGTTGCTTTTCAGCTCAGGATCGTAATGCCAt          <  1:1840574/67‑1 (MQ=255)                                                           gagaCGCCACGGCAGGTTACATGCTGGAAGGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCAt          <  1:186527/66‑1 (MQ=255)                                                             gaCACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCATTtctc     <  1:669712/69‑1 (MQ=255)                                                             gaCACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCATTtctc     <  1:828197/69‑1 (MQ=255)                                                              acacCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGt                  <  1:2165924/55‑1 (MQ=255)                                                              acacCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCa           <  1:746231/62‑1 (MQ=255)                                                              acacCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCATTtc       <  1:2232742/66‑1 (MQ=255)                                                              acacCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCATTTctct    <  1:2416226/69‑1 (MQ=255)                                                               caccacGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGt                  <  1:1258666/54‑1 (MQ=255)                                                                accacGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCATTTCTCttt  >  1:3320879/1‑69 (MQ=255)                                                                accacGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCATTTCTCttt  >  1:3273713/1‑69 (MQ=255)                                                                accacGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCATTTCTCttt  >  1:1921597/1‑69 (MQ=255)                                                                accacGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCATTTCTCttt  >  1:2888365/1‑69 (MQ=255)                                                                |                                                                     ACGCGGGCAATCCGCCATCACTTCCGGCGAAGCGGTTTCTGCTTTGGCTTCATGATAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCATTTCTCTTT  >  minE/125102‑125233

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