Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 125,420 C→T 51.0% P182P (CCG→CCA)  gcd ← glucose dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE125,4200CT51.0% ‑4.0 / 92.7 100P182P (CCG→CCA) gcdglucose dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (15/34);  new base T (36/15);  total (51/49)
Fisher's exact test for biased strand distribution p-value = 1.20e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.37e-01

GTACGGAACACCCAGGCTTCTTTCAGATTATGGACGTTATCGGCGTTAATTTGTTTCAGCGGCGAAAAGCGTTGACCTTCCTGATTACGACCATAGGCAGGCCAGTCCTGATCGGCTA  >  minE/125361‑125478
                                                           |                                                          
gTACGGAACACCCAGGCTTCTTTCAGATTATTGACGTTATCGGCGTTAATTTGTTTCAGCGGCGAAAAg                                                   <  1:2765523/69‑1 (MQ=255)
gTACGGAACACCCAGGCTTCTTTCAGATTATGGACGTTATCGGCGTTAATTTGTTTCAGCGGCGAAAAg                                                   <  1:2487930/69‑1 (MQ=255)
gTACGGAACACCCAGGCTTCTTTCAGATTATGGACGTTATCGGCGTTAATTTGTTTCAGCGGCGAAAAg                                                   <  1:1970540/69‑1 (MQ=255)
gTACGGAACACCCAGGCTTCTTTCAGATTATGGACGTTATCGGCGTTAATTTGTTTCAGCGGCGAAAAg                                                   <  1:358125/69‑1 (MQ=255)
gTACGGAACACCCAGGCTTCTTTCAGATTATGGACGTTATCGGCGTTAATTTGTTTCAGCGGCGAAAAg                                                   <  1:1615944/69‑1 (MQ=255)
gTACGGAACACCCAGGCTTCTTTCAGATTATGGACGTTATCGGCGTTAATTTGTTTCAGCGGCGAAAAg                                                   <  1:237929/69‑1 (MQ=255)
gTACGGAACACCCAGGCTTCTTTCAGATTATGGACGTTATCGGCGTTAATTTGTTTCAGCGGCGAAAAg                                                   <  1:1448472/69‑1 (MQ=255)
gTACGGAACACCCAGGCTTCTTTCAGATTATGGACGTTATCGGCGTTAATTTGTTTCAGCGGCGAAAAg                                                   <  1:2740856/69‑1 (MQ=255)
gTACGGAACACCCAGGCTTCTTTCAGATTATGGACGTTATCGGCGTTAATTTGTTTCAGCGGCGAAAAg                                                   <  1:3212579/69‑1 (MQ=255)
gTACGGAACACCCAGGCTTCTTTCAGATTATGGACGTTATCGGCGTTAATTTGTTTCAGCGGCGAAAAg                                                   <  1:536265/69‑1 (MQ=255)
  aCGGAACCCCCAGGCTTCTTTCAGCTTGTGGACGTTATCAGCGTTAATTTGCTTCAGTGGTGAAAAGCg                                                 >  1:882678/1‑69 (MQ=255)
  aCGGAACACCCAGGCTTCTTTCAGCTTGTGGACGTTATCAGCGTTAATTTGCTTCGGTGGTGAAAAGCg                                                 >  1:2630930/1‑69 (MQ=255)
  aCGGAACACCCAGGCTTCTTTCAGCTTGTGGACGTTATCAGCGTTAATTTGCTTCAGTGGTGAAAAGc                                                  >  1:859148/1‑68 (MQ=255)
  aCGGAACACCCAGGCTTCTTTCAGCTTGTGGACGTTATCAGCGTTAATTTGCTTCAGTGGTGAAAAGCg                                                 >  1:2494485/1‑69 (MQ=255)
  aCGGAACACCCAGGCTTCTTTCAGCTTGTGGACGTTATCAGCGTTAATTTGCTTCAGTGGTGAAAAGCg                                                 >  1:1056445/1‑69 (MQ=255)
  aCGGAACACCCAGGCTTCTTTCAGCTTGTGGACGTTATCAGCGTTAATTTGCTTCAGTGGTGAAAAGCg                                                 >  1:1045577/1‑69 (MQ=255)
  aCGGAACACCCAGGCTTCTTTCAGCTTGTGGACGTTATCAGCGTTAATTTGCTTCAGTGGTGAAAAGCg                                                 >  1:2576589/1‑69 (MQ=255)
  aCGGAACACCCAGGCTTCTTTCAGCTTGTGGACGTTATCAGCGTTAATTTGCTTCAGTGGTGAAAAGCg                                                 >  1:460851/1‑69 (MQ=255)
  aCGGAACACCCAGGCTTCTTTCAGCTTGTGGACGTTATCAGCGTTAATTTGCTTCAGTGGTGAAAAGCg                                                 >  1:2377523/1‑69 (MQ=255)
  aCGGAACACCCAGGCTTCTTTCAGCTTGTGGACGTTATCAGCGTTAATTTGCTTCAGTGGTGAAAAGCg                                                 >  1:2826107/1‑69 (MQ=255)
  aCGGAACACCCAGGCTTCTTTCAGCTTGTGGACGTTATCAGCGTTAATTTGCTTCAGTGGTGAAAAGCg                                                 >  1:28759/1‑69 (MQ=255)
  aCGGAACACCCAGGCTTCTTTCAGCTTGTGGACGTTATCAGCGTTAATTTGCTTCAGTGGTGAAAAGCg                                                 >  1:1119219/1‑69 (MQ=255)
  aCGGAACACCCAGGCTTCTTTCAGCTTGTGGACGTTATCAGCGTTAATTTGCTTCAGTGGTGAAAAGCg                                                 >  1:3065593/1‑69 (MQ=255)
  aCGGAACACCCAGGCTTCTTTCAGCTTGTGGACGTTATCAGCGTTAATTTGCTTCAGTGGTGAAAAGCg                                                 >  1:3380948/1‑69 (MQ=255)
  aCGGAACACCCAGGCTTCTTTCAGATTATGGACGTTATCGGCGTTAATTTGTTTCAGCGGCGAAAAGCg                                                 >  1:792122/1‑69 (MQ=255)
  aCGGAACACCCAGGCTTCTTTCAGATTATGGACGTTATCGGCGTTAATTTGTTTCAGCGGCGAAAAGCg                                                 >  1:2426253/1‑69 (MQ=255)
   cGGAACACCCAGGCTTCTTTCAGATTATGGACGTTATCGGCGTTAATTTGTTTCAGCGGCGAAAAGCg                                                 >  1:2983539/1‑68 (MQ=255)
   cGGAACACCCAGGCTTCTTTCAGATTATGGACGTTATCGGCGTTAATTTGTTTCAGCGGCGAAAAGCg                                                 >  1:1904489/1‑68 (MQ=255)
   cGGAACACCCAGGCTTCTTTCAGATTATGGACGTTATCGGCGTTAATTTGTTTCAGCGGCGAAAAGCg                                                 >  1:1197319/1‑68 (MQ=255)
   cGGAACACCCAGGCTTCTTTCAGATTATGGACGTTATCGGCGTTAATTTGTTTCAGCGGCGAAAAGCg                                                 >  1:228711/1‑68 (MQ=255)
   cGGAACACCCAGGCTTCTTTCAGATTATGGACGTTATCGGCGTTAATTTGTTTCAGCGGCGAAAAGCg                                                 >  1:2305292/1‑68 (MQ=255)
   cGGAACACCCAGGCTTCTTTCAGATTATGGACGTTATCGGCGTTAATTTGTTTCAGCGGCGAAAAGCg                                                 >  1:1387996/1‑68 (MQ=255)
   cGGAACACCCAGGCTTCTTTCAGATTATGGACGTTATCGGCGTTAATTTGTTTCAGCGGCGAAAAGCg                                                 >  1:1499356/1‑68 (MQ=255)
      aaCACCCAGGCTTCTTTCAGATTATGGACGTTATCGGCGTTAATTTGTTTCAGCGGCGAAAAGCGTTGa                                             <  1:2625045/69‑1 (MQ=255)
      aaCACCCAGGCTTCTTTCAGATTATGGACGTTATCGGCGTTAATTTGTTTCAGCGGCGAAAAGCGTTGa                                             <  1:2680212/69‑1 (MQ=255)
      aaCACCCAGGCTTCTTTCAGATTATGGACGTTATCGGCGTTAATTTGTTTCAGCGGCGAAAAGCGTTGa                                             <  1:810253/69‑1 (MQ=255)
      aaCACCCAGGCTTCTTTCAGATTATGGACGTTATCGGCGTTAATTTGTTTCAGCGGCGAAAAGCGTTGa                                             <  1:3146023/69‑1 (MQ=255)
      aaCACCCAGGCTTCTTTCAGATTATGGACGTTATCGGCGTTAATTTGTTTCAGCGGCGAAAAGCGTTGa                                             <  1:2358052/69‑1 (MQ=255)
      aaCACCCAGGCTTCTTTCAGATTATGGACGTTATCGGCGTTAATTTGTTTCAGCGGCGAAAAGCGTTGa                                             <  1:2271238/69‑1 (MQ=255)
      aaCACCCAGGCTTCTTTCAGATTATGGACGTTATCGGCGTTAATTTGTTTCAGCGGCGAAAAGCGTTGa                                             <  1:66677/69‑1 (MQ=255)
      aaCACCCAGGCTTCTTTCAGATTATGGACGTTATCGGCGTTAATTTGTTTCAGCGGCGAAAAGCGTTGa                                             <  1:684776/69‑1 (MQ=255)
      aaCACCCAGGCTTCTTTCAGATTATGGACGTTATCGGCGTTAATTTGTTTCAGCGGCGAAAAGCGTTGa                                             <  1:1971037/69‑1 (MQ=255)
                  tcttTCAGCTTGTGGACGTTATCAGCGTTAATTTGCTTCAGTGGTGAAAAGCGTTGACCTTCCTGAttg                                 <  1:2758796/69‑2 (MQ=255)
                  tcttTCAGCTTGTGGACGTTATCAGCGTTAATTTGCTTCAGTGGTGAAAAGCGTTGACCTTCCTGAttg                                 <  1:299409/69‑2 (MQ=255)
                  tcttTCAGCTTGTGGACGTTATCAGCGTTAATTTGCTTCAGTGGTGAAAAGCGTTGACCTTCCTGAttg                                 <  1:3221688/69‑2 (MQ=255)
                  tcttTCAGCTTGTGGACGTTATCAGCGTTAATTTGCTTCAGTGGTGAAAAGCGTTGACCTTCCTGAttg                                 <  1:325195/69‑2 (MQ=255)
                  tcttTCAGCTTGTGGACGTTATCAGCGTTAATTTGCTTCAGTGGTGAAAAGCGTTGACCTTCCTGAttg                                 <  1:2434310/69‑2 (MQ=255)
                  tcttTCAGCTTGTGGACGTTATCAGCGTTAATTTGCTTCAGTGGTGAAAAGCGTTGACCTTCCTGAttg                                 <  1:2592203/69‑2 (MQ=255)
                  tcttTCAGCTTGTGGACGTTATCAGCGTTAATTTGCTTCAGTGGTGAAAAGCGTTGACCTTCCTGAttg                                 <  1:1402113/69‑2 (MQ=255)
                  tcttTCAGCTTGTGGACGTTATCAGCGTTAATTTGCTTCAGTGGTGAAAAGCGTTGACCTTCCTGAttg                                 <  1:1521373/69‑2 (MQ=255)
                  tcttTCAGCTTGTGGACGTTATCAGCGTTAATTTGCTTCAGTGGTGAAAAGCGTTGACCTTCCTGAttg                                 <  1:1942959/69‑2 (MQ=255)
                  tcttTCAGCTTGTGGACGTTATCAGCGTTAATTTGCTTCAGTGGTGAAAAGCGTTGACCTTCCTGAttg                                 <  1:689254/69‑2 (MQ=255)
                  tcttTCAGCTTGTGGACGTTATCAGCGTTAATTTGCTTCAGTGGTGAAAAGCGTTGACCTTCCTGAttg                                 <  1:547953/69‑2 (MQ=255)
                  tcttTCAGCTTGTGGACGTTATCAGCGTTAATTTGCTTCAGTGGTGAAAAGCGTTGACCTTCCTGAttg                                 <  1:2499633/69‑2 (MQ=255)
                  tcttTCAGCTTGTGGACGTTATCAGCGTTAATTTGCTTCAGTGGTGAAAAGCGTTGACCTTCCTGAttg                                 <  1:2588433/69‑2 (MQ=255)
                   cttTCAGCTTGTGGACGTTATCAGCGTTAATTTGCTTCAGTGGTGAAAAGCGTTGACCTTCCTGAttg                                 <  1:3214108/68‑2 (MQ=255)
                     ttCAGCTTGTGGACGTTATCAGCGTTAATTTGCTTCAGTGGTGAAAAGCGTTGACCTTCCTGAttg                                 <  1:1353458/66‑2 (MQ=255)
                     ttCAGATTATGGACGTTATCGGCGTTAATTTGTTTCAGCGGCGAAAAGCGTTGACCTTCCTGATTAc                                <  1:2330705/67‑1 (MQ=255)
                     ttCAGATTATGGACGTTATCGGCGTTAATTTGTTTCAGCGGCGAAAAGCGTTGACCTTCCTGATTAc                                <  1:3369575/67‑1 (MQ=255)
                     ttCAGATTATGGACGTTATCGGCGTTAATTTGTTTCAGCGGCGAAAAGCGTTGACCTTCCTGATTAc                                <  1:295906/67‑1 (MQ=255)
                     ttCAGATTATGGACGTTATCGGCGTTAATTTGTTTCAGCGGCGAAAAGCGTTGACCTTCCTGATTAc                                <  1:1650719/67‑1 (MQ=255)
                     ttCAGATTATGGACGTTATCGGCGTTAATTTGTTTCAGCGGCGAAAAGCGTTGACCTTCCTGATTAc                                <  1:2620457/67‑1 (MQ=255)
                                   gTTATCAGCGTTGATTTGCTTCAGTGGTGAAAAGCGTTGACCTTCCTGATTGCGACCATAGGCAGGCCa                >  1:2418676/1‑69 (MQ=255)
                                   gTTATCAGCGTTAATTTGCTTCAGTGGTGAAAAGCGTTGACCTTCCTGATTGCGACCATaggcaggcga                >  1:1377747/1‑67 (MQ=255)
                                   gTTATCAGCGTTAATTTGCTTCAGTGGTGAAAAGCGTTGACCTTCCTGATTGCGACCATag                        >  1:3277944/1‑61 (MQ=255)
                                   gTTATCAGCGTTAATTTGCTTCAGTGGTGAAAAGCGTTGACCTTCCTGATTGCGACCATAGGCAGGCCa                >  1:3384448/1‑69 (MQ=255)
                                   gTTATCAGCGTTAATTTGCTTCAGTGGTGAAAAGCGTTGACCTTCCTGATTGCGACCATAGGCAGGCCa                >  1:3315230/1‑69 (MQ=255)
                                   gTTATCAGCGTTAATTTGCTTCAGTGGTGAAAAGCGTTGACCTTCCTGATTGCGACCATAGGCAGGCCa                >  1:880417/1‑69 (MQ=255)
                                   gTTATCAGCGTTAATTTGCTTCAGTGGTGAAAAGCGTTGACCTTCCTGATTGCGACCATAGGCAGGCCa                >  1:79244/1‑69 (MQ=255)
                                   gTTATCAGCGTTAATTTGCTTCAGTGGTGAAAAGCGTTGACCTTCCTGATTGCGACCATAGGCAGGCCa                >  1:37225/1‑69 (MQ=255)
                                   gTTATCAGCGTTAATTTGCTTCAGTGGTGAAAAGCGTTGACCTTCCTGATTGCGACCATAGGCAGGCCa                >  1:675752/1‑69 (MQ=255)
                                   gTTATCAGCGTTAATTTGCTTCAGTGGTGAAAAGCGTTGACCTTCCTGATTGCGACCATAGGCAGGCCa                >  1:371806/1‑69 (MQ=255)
                                   gTTATCAGCGTTAATTTGCTTCAGTGGTGAAAAGCGTTGACCTTCCTGATTGCGACCATAGGCAGGCCa                >  1:2695279/1‑69 (MQ=255)
                                   gTTATCAGCGTTAATTTGCTTCAGTGGTGAAAAGCGTTGACCTTCCTGATTGCGACCATAGGCAGGCCa                >  1:1347683/1‑69 (MQ=255)
                                   gTTATCAGCGTTAATTTGCTTCAGTGGTGAAAAGCGTTGACCTTCCTGATTGCGACCATAGGCAGGCCa                >  1:200541/1‑69 (MQ=255)
                                   gTTATCAGCGTTAATTTGCTTCAGTGGTGAAAAGCGTTGACCTTCCTGATTGCGACCATAGGCAGGCCa                >  1:3174308/1‑69 (MQ=255)
                                   gTTATCAGCGTTAATTTGCTTCAGTGGTGAAAAGCGTTGACCTTCCTGATTGCGACCATAGGCAGGCCa                >  1:2219600/1‑69 (MQ=255)
                                   gTTATCAGCGTTAATTTGCTTCAGTGGTGAAAAGCGTTGACCTTCCTGATTGCGACCATAGGCAGGCCa                >  1:3023737/1‑69 (MQ=255)
                                   gTTATCAGCGTTAATTTGCTTCAGTGGTGAAAAGCGTTGACCTTCCTGATTGCGACCATAGGCAGGCCa                >  1:2905413/1‑69 (MQ=255)
                                   gTTATCAGCGTTAATTTGCTTCAGTGGTGAAAAGCGTTGACCTTCCTGATTGCGACCATAGGCAGGCCa                >  1:2783464/1‑69 (MQ=255)
                                   gTTATCAGCGTTAATTTGCTTCAGTGGTGAAAAGCGTTGACCTTCCTGATTGCGACCATAGGCAGGCCa                >  1:2472213/1‑69 (MQ=255)
                                   gTTATCAGCGTTAATTTGCTTCAGTGGTGAAAAGCGTTGACCTTCCTGATTGCGACCATAGGCAGGCCa                >  1:261269/1‑69 (MQ=255)
                                   gTTATCAGCGTTAATTTGCTTCAGTGGTGAAAAGCGTTGACCTTCCTGATTGCGACCATAGGCAGGCCa                >  1:2658298/1‑69 (MQ=255)
                                   gTTATCAGCGTTAATTTGCTTCAGTGGTGAAAAGCGTTGACCTTCCTGATTGAGACCATAGGCAGGCCa                >  1:499470/1‑69 (MQ=255)
                                      aTCGGCGTTAATTTGTTTCAGCGGCGAAAAGCGTTGACCTTCCTGATTACGACCATa                         >  1:3383527/1‑57 (MQ=255)
                                      aTCGGCGTTAATTTGTTTCAGCGGCGAAAAGCGTTGACCTTCCTGATTACGACCATAGGCAGGCCAGTc             >  1:2107416/1‑69 (MQ=255)
                                      aTCGGCGTTAATTTGTTTCAGCGGCGAAAAGCGTTGACCTTCCTGATTACGACCATAGGCAGGCCAGTc             >  1:3268000/1‑69 (MQ=255)
                                         ggCGTTAATTTGTTTCAGCGGCGAAAAGCGTTGACCTTCCTGATTACGACCATAGGCAGGCCAGTCCTg          <  1:2741277/69‑1 (MQ=255)
                                         ggCGTTAATTTGTTTCAGCGGCGAAAAGCGTTGACCTTCCTGATTACGACCATAGGCAGGCCAGTCCTg          <  1:2118272/69‑1 (MQ=255)
                                          gCGTTAATTTGTTTCAGCGGCGAAAAGCGTTGACCTTCCTGATTACGACCATAGGCAGGCCAGTCCTg          <  1:1379588/68‑1 (MQ=255)
                                               aaTTTGTTTCAGCGGCGAAAAGCGTTGACCTTCCTGATTACGACCATAGGCAGGCCAGTCCTGATCGGc    <  1:12607/69‑1 (MQ=255)
                                                 tttGTTTCAGCGGCGAAAAGCGTTGACCTTCCTGATTACGACCATAGGCAGGCCAGTCCTGATCGGCTa  >  1:397131/1‑69 (MQ=255)
                                                 tttGTTTCAGCGGCGAAAAGCGTTGACCTTCCTGATTACGACCATAGGCAGGCCAGTCCTGATCGGCTa  >  1:1073343/1‑69 (MQ=255)
                                                 tttGTTTCAGCGGCGAAAAGCGTTGACCTTCCTGATTACGACCATAGGCAGGCCAGTCCTGATCGGCTa  >  1:151714/1‑69 (MQ=255)
                                                   tGTTTCAGCGGCGAAAAGCGTTGACCTTCCTGATTACGACCATaggcag                    <  1:918427/49‑1 (MQ=255)
                                                   tGTTTCAGCGGCGAAAAGCGTTGACCTTCCTGATTACGACCATaggcag                    <  1:3310131/49‑1 (MQ=255)
                                                   tGTTTCAGCGGCGAAAAGCGTTGACCTTCCTGATTACGACCATaggcag                    <  1:863758/49‑1 (MQ=255)
                                                   tGTTTCAGCGGCGAAAAGCGTTGACCTTCCTGATTACGACCATaggcag                    <  1:1764743/49‑1 (MQ=255)
                                                   tGTTTCAGCGGCGAAAAGCGTTGACCTTCCTGATTACGACCATaggcag                    <  1:791661/49‑1 (MQ=255)
                                                   tGTTTCAGCGGCGAAAAGCGTTGAACTTCCTGATTACGACCATagggag                    <  1:3292707/49‑1 (MQ=38)
                                                           |                                                          
GTACGGAACACCCAGGCTTCTTTCAGATTATGGACGTTATCGGCGTTAATTTGTTTCAGCGGCGAAAAGCGTTGACCTTCCTGATTACGACCATAGGCAGGCCAGTCCTGATCGGCTA  >  minE/125361‑125478

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