Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 125,180 G→A 33.2% H262H (CAC→CAT gcd ← glucose dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE125,1800GA33.2% 84.8 / 53.8 102H262H (CAC→CATgcdglucose dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (16/52);  new base A (18/16);  total (34/68)
Fisher's exact test for biased strand distribution p-value = 3.99e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.62e-01

Only 100 of 102 total aligned reads displayed.

ACTTCCGGCGAAGCGGTTTCTGCTTTGGCTTCATGATAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCATTTCTCTTT  >  minE/125121‑125233
                                                           |                                                     
aCTTCCGGCGAAGCGGTTTCTGCTTTGGCTTCATGATAAGAGACACCACGGCAGGTTACGTGCTGGaaa                                              >  1:690617/1‑69 (MQ=255)
aCTTCCGGCGAAGCGGTTTCTGCTTTGGCTTCATGATAAGAGACACCACGGCAGGTTACGTGCTGGaaa                                              >  1:602357/1‑69 (MQ=255)
aCTTCCGGCGAAGCGGTTTCTGCTTTGGCTTCATGATAAGAGACACCACGGCAGGTTACGTGCTGGaaa                                              >  1:1713968/1‑69 (MQ=255)
aCTTCCGGCGAAGCGGTTTCTGCTTTGGCTTCATGATAAGAGACACCACGGCAGGTTACGTGCTGGaaa                                              >  1:1994441/1‑69 (MQ=255)
aCTTCCGGCGAAGCGGTTTCTGCTTTGGCTTCATGATAAGAGACACCACGGCAGGTTACGTGCTGGaaa                                              >  1:1441309/1‑69 (MQ=255)
aCTTCCGGCGAAGCGGTTTCTGCTTTGGCTTCATGATAAGAGACACCACGGCAGGTTACGTGCTGGaaa                                              >  1:1923513/1‑69 (MQ=255)
aCTTCCGGCGAAGCGGTTTCTGCTTTGGCTTCATGATAAGAGACACCACGGCAGGTTACGTGCTGGaaa                                              >  1:2828240/1‑69 (MQ=255)
aCTTCCGGCGAAGCGGTTTCTGCTTTGGCTTCATGATAAGAGACACCACGGCAGGTTACGTGCTGGaaa                                              >  1:3279672/1‑69 (MQ=255)
aCTTCCGGCGAAGCGGTTTCTGCTTTGGCTTCATGATAAGAGACACCACGGCAGGTTACGTGCTGGaa                                               >  1:2847457/1‑68 (MQ=255)
aCTTCCGGCGAAGCAGTTTCTGCTTTGGCTTCATGATAAGAGACGCCACGGCAGGTTACATGCTGGaag                                              >  1:3083329/1‑68 (MQ=255)
aCTTCCGGCGAAGCAGTTTCTGCTTTGGCTTCATGATAAGAGACGCCACGGCAGGTTACATGCTGGaag                                              >  1:588566/1‑68 (MQ=255)
aCTTCCGGCGAAGCAGTTTCTGCTTTGGCTTCATGATAAGAGACGCCACGGCAGGTTACATGCTGGaag                                              >  1:2992560/1‑68 (MQ=255)
aCTTCCGGCGAAGCAGTTTCTGCTTTGGCTTCATGATAAGAGACGCCACGGCAGGTTACATGCTGGaag                                              >  1:2896890/1‑68 (MQ=255)
aCTTCCGGCGAAGCAGTTTCTGCTTTGGCTTCATGATAAGAGACGCCACGGCAGGTTACATGCTGGaag                                              >  1:1910709/1‑68 (MQ=255)
aCTTCCGGCGAAGCAGTTTCTGCTTTGGCTTCATGATAAGAGACGCCACGGCAGGTTACATGCTGGaag                                              >  1:1946691/1‑68 (MQ=255)
aCTTCCGGCGAAGCAGTTTCTGCTTTGGCTTCATGATAAGAGACGCCACGGCAGGTTACATGCTGGaag                                              >  1:2066344/1‑68 (MQ=255)
aCTTCCGGCGAAGCAGTTTCTGCTTTGGCTTCATGATAAGAGACGCCACGGCAGGTTACATGCTGGaag                                              >  1:2624873/1‑68 (MQ=255)
aCTTCCGGCGAAGCAGTTTCTGCTTTGGCTTCATGATAAGAGACGCCACGGCAGGTTACATGCTGGaag                                              >  1:1058307/1‑68 (MQ=255)
aCTTCCGGCGAAGCAGTTTCTGCTTTGGCTTCATGATAAGAGACGCCACGGCAGGTTACATGCTGGaag                                              >  1:2463468/1‑68 (MQ=255)
aCTTCCGGCGAAGCAGTTTCTGCTTTGGCTTCATGATAAGAGACGCCACGGCAGGTTACATGCTGGaag                                              >  1:1263375/1‑68 (MQ=255)
aCTTCCGGCGAAGCAGTTTCTGCTTTGGCTTCATGATAAGAGACGCCACGGCAGGTTACATGCTGGaag                                              >  1:337114/1‑68 (MQ=255)
aCTTCCGGCGAAGCAGTTTCTGCTTTGGCTTCATGATAAGAGACGCCACGGCAGGTTACATGCTGGaag                                              >  1:774835/1‑68 (MQ=255)
aCTTCCGGCGAAGCAGTTTCTGCTTTGGCTTCATGATAAGAGACGCCACGGCAGGTTACATGCTGGaag                                              >  1:134501/1‑68 (MQ=255)
aCTTCCGGCGAAGCAGTTTCTGCTTTGGCTTCATGATAAGAGACGCCACGGCAGGTTACATGCTGGaag                                              >  1:3350012/1‑68 (MQ=255)
aCTTCCGGCGAAGCAGTTTCTGCTTTGGCTTCATGATAAGAGACGCCACGGCAGGTTACATGCTGGaag                                              >  1:324447/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:1452121/68‑1 (MQ=255)
                     gCTTTGGCTTCATGATAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTc                          <  1:1362334/68‑1 (MQ=255)
                     gCTTTGGCTTCATGATAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTc                          <  1:3144600/68‑1 (MQ=255)
                      cTTTGGCTTCATGATAAGAGACACCACGGCAGGTTACGTGCTGGa                                                <  1:119935/45‑1 (MQ=255)
                        ttGGCTTCATGATAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAgg                       <  1:3393325/68‑1 (MQ=255)
                        ttGGCTTCATGATAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAgg                       <  1:2360497/68‑1 (MQ=255)
                        ttGGCTTCATGATAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAgg                       <  1:2743821/68‑1 (MQ=255)
                        ttGGCTTCATGATAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAgg                       <  1:3251234/68‑1 (MQ=255)
                        ttGGCTTCATGATAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAgg                       <  1:3048684/68‑1 (MQ=255)
                        ttGGCTTCATGATAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAgg                       <  1:319048/68‑1 (MQ=255)
                        ttGGCTTCATGATAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAgg                       <  1:2064846/68‑1 (MQ=255)
                        ttGGCTTCATGATAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAgg                       <  1:212890/68‑1 (MQ=255)
                             ttCATGATAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGt                  <  1:1852745/68‑1 (MQ=255)
                             ttCATGATAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGt                  <  1:1800757/68‑1 (MQ=255)
                             ttCATGATAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGt                  <  1:1797281/68‑1 (MQ=255)
                             ttCATGATAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGt                  <  1:1753021/68‑1 (MQ=255)
                             ttCATGATAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGt                  <  1:3168844/68‑1 (MQ=255)
                             ttCATGATAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGt                  <  1:2788455/68‑1 (MQ=255)
                             ttCATGATAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGt                  <  1:1734844/68‑1 (MQ=255)
                             ttCATGATAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGt                  <  1:2586276/68‑1 (MQ=255)
                             ttCATGATAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGt                  <  1:736960/68‑1 (MQ=255)
                             ttCATGATAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGt                  <  1:1505843/68‑1 (MQ=255)
                             ttCATGATAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGt                  <  1:2217404/68‑1 (MQ=255)
                             ttCATGATAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGt                  <  1:563938/68‑1 (MQ=255)
                                   aTAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCa           <  1:1746104/69‑1 (MQ=255)
                                   aTAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCa           <  1:519770/69‑1 (MQ=255)
                                   aTAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCa           <  1:2970460/69‑1 (MQ=255)
                                    tAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGt                  <  1:2239916/61‑1 (MQ=255)
                                      agagACGCCACGGCAGGTTACATGCTGGAAGGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCAt          <  1:1322401/67‑1 (MQ=255)
                                      agagACGCCACGGCAGGTTACATGCTGGAAGGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCAt          <  1:1526929/67‑1 (MQ=255)
                                      agagACGCCACGGCAGGTTACATGCTGGAAGGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCAt          <  1:3433508/67‑1 (MQ=255)
                                      agagACGCCACGGCAGGTTACATGCTGGAAGGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCAt          <  1:411243/67‑1 (MQ=255)
                                      agagACGCCACGGCAGGTTACATGCTGGAAGGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCAt          <  1:3176260/67‑1 (MQ=255)
                                      agagACGCCACGGCAGGTTACATGCTGGAAGGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCAt          <  1:424362/67‑1 (MQ=255)
                                      agagACGCCACGGCAGGTTACATGCTGGAAGGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCAt          <  1:2986548/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)
                                      agagACGCCACGGCAGGTTACATGCTGGAAGGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCAt          <  1:719336/67‑1 (MQ=255)
                                      agagACGCCACGGCAGGTTACATGCTGGAAGGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCAt          <  1:1881451/67‑1 (MQ=255)
                                      agagACGCCACGGCAGGTTACATGCTGGAAGGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCAt          <  1:2704519/67‑1 (MQ=255)
                                      agagACGCCACGGCAGGTTACATGCTGGAAGGACTCGTTGCTTTTCAGCTCAGGATCGTAATGCCAt          <  1:1840574/67‑1 (MQ=255)
                                       gagaCGCCACGGCAGGTTACATGCTGGAAGGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCAt          <  1:186527/66‑1 (MQ=255)
                                         gaCACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCATTtctc     <  1:828197/69‑1 (MQ=255)
                                         gaCACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCATTtctc     <  1:669712/69‑1 (MQ=255)
                                          acgccacGGCAGGTTACATGCTGGAAGGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCAt          <  1:1550887/60‑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:2888365/1‑69 (MQ=255)
                                            accacGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCATTTCTCttt  >  1:3273713/1‑69 (MQ=255)
                                            accacGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCATTTCTCttt  >  1:1921597/1‑69 (MQ=255)
                                                cGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCATTtctc     <  1:1913353/62‑1 (MQ=255)
                                                       ttACGTGCTTGAAAGACTCGTTGGTTTTCAGCTCAGGATCGt                  <  1:397986/42‑1 (MQ=255)
                                                       ttACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGt                  <  1:257548/42‑1 (MQ=255)
                                                       ttACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGt                  <  1:811431/42‑1 (MQ=255)
                                                       ttACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGt                  <  1:324414/42‑1 (MQ=255)
                                                       ttACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGt                  <  1:2961013/42‑1 (MQ=255)
                                                       ttACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGt                  <  1:2464512/42‑1 (MQ=255)
                                                       ttACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGt                  <  1:2146761/42‑1 (MQ=255)
                                                       ttACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGt                  <  1:2143799/42‑1 (MQ=255)
                                                       ttACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGt                  <  1:1916184/42‑1 (MQ=255)
                                                       ttACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGt                  <  1:1644124/42‑1 (MQ=255)
                                                       ttACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGt                  <  1:1641187/42‑1 (MQ=255)
                                                       ttACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGt                  <  1:1158550/42‑1 (MQ=255)
                                                       ttACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCa           <  1:1747331/49‑1 (MQ=255)
                                                           |                                                     
ACTTCCGGCGAAGCGGTTTCTGCTTTGGCTTCATGATAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCATTTCTCTTT  >  minE/125121‑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: