Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 125,477 T→A 24.9% V163V (GTA→GTT gcd ← glucose dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE125,4770TA24.9% 44.9 / 11.8 32V163V (GTA→GTTgcdglucose dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (20/4);  new base A (6/2);  total (26/6)
Fisher's exact test for biased strand distribution p-value = 6.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTGTTTCAGCGGCGAAAAGCGTTGACCTTCCTGATTACGACCATAGGCAGGCCAGTCCTGATCGGCTACGGGGGAGATAGCTTCAGCAGGTGTGGCATCGGCGCTTAAGGTGCCGTTGATCTCCTGCGGAT  >  minE/125410‑125541
                                                                   |                                                                
tttGTTTCAGCGGCGAAAAGCGTTGACCTTCCTGATTACGACCATAGGCAGGCCAGTCCTGATCGGCTa                                                                 >  1:1073343/1‑69 (MQ=255)
tttGTTTCAGCGGCGAAAAGCGTTGACCTTCCTGATTACGACCATAGGCAGGCCAGTCCTGATCGGCTa                                                                 >  1:397131/1‑69 (MQ=255)
tttGTTTCAGCGGCGAAAAGCGTTGACCTTCCTGATTACGACCATAGGCAGGCCAGTCCTGATCGGCTa                                                                 >  1:151714/1‑69 (MQ=255)
                       tGACCTTCCTGATTGCGACCATAGGCAGGCCAGTCCTGATCGGCAACGGGGGAGATGGCTTCTGCTggc                                          >  1:2802582/1‑68 (MQ=255)
                       tGACCTTCCTGATTGCGACCATAGGCAGGCCAGTCCTGATCGGCAACGGGGGAGATGGCTTCTGCTggc                                          >  1:776173/1‑68 (MQ=255)
                       tGACCTTCCTGATTGCGACCATAGGCAGGCCAGTCCTGATCGGCAACGGGGGAGATGGCTTCTGCTggc                                          >  1:2395529/1‑68 (MQ=255)
                       tGACCTTCCTGATTGCGACCATAGGCAGGCCAGTCCTGATCGGCAACGGGGGAGATGGCTTCTGCTggc                                          >  1:1891082/1‑68 (MQ=255)
                       tGACCTTCCTGATTGCGACCATAGGCAGGCCAGTCCTGATCGGCAACGGGGGAGATGGCTTCTGCTggc                                          >  1:1884412/1‑68 (MQ=255)
                       tGACCTTCCTGATTGCGACCATAGGCAGGCCAGTCCTGATCGGCAACGGGGGAGATGGCTTCTGCTgg                                           >  1:1444255/1‑68 (MQ=255)
                        gACCTTCCTGATTGCGACCATAGGCAGGCCAGTCCTGATCGGCAACGGGGGAGATGGCTTCTGCTggcg                                         <  1:3258584/69‑3 (MQ=255)
                        gACCTTCCTGATTGCGACCATAGGCAGGCCAGTCCTGATCGGCAACGGGGGAGATGGCTTCTGCTggcg                                         <  1:949902/69‑3 (MQ=255)
                                    tACGACCATAGGCAGGCCAGTCCTGATCGGCTACGGGGGAGATAGCTTCAGCAGGTGTgg                                      >  1:1693892/1‑60 (MQ=255)
                                    tACGACCATAGGCAGGCCAGTCCTGATCGGCTACGGGGGAGATAGCTTCAGCAGGTGTGGCATCGgcg                              >  1:703967/1‑68 (MQ=255)
                                    tACGACCATAGGCAGGCCAGTCCTGATCGGCTACGGGGGAGATAGCTTCAGCAGGTGTGGCATCGgcg                              >  1:1532830/1‑68 (MQ=255)
                                    tACGACCATAGGCAGGCCAGTCCTGATCGGCTACGGGGGAGATAGCTTCAGCAGGTGTGGCATCGgcg                              >  1:749228/1‑68 (MQ=255)
                                    tACGACCATAGGCAGGCCAGTCCTGATCGGCTACGGGGGAGATAGCTTCAGCAGGTGTGGCATCGgcg                              >  1:631306/1‑68 (MQ=255)
                                    tACGACCATAGGCAGGCCAGTCCTGATCGGCTACGGGGGAGATAGCTTCAGCAGGTGTGGCATCGgcg                              >  1:1917902/1‑68 (MQ=255)
                                    tACGACCATAGGCAGGCCAGTCCTGATCGGCTACGGGGGAGATAGCTTCAGCAGGTGTGGCATCGgcg                              >  1:1333462/1‑68 (MQ=255)
                                    tACGACCATAGGCAGGCCAGTCCTGATCGGCTACGGGGGAGATAGCTTCAGCAGGTGTGGCATCGgcg                              >  1:2446655/1‑68 (MQ=255)
                                    tACGACCATAGGCAGGCCAGTCCTGATCGGCTACGGGGGAGATAGCTTCAGCAGGTGTGGCATCGgcg                              >  1:1146608/1‑68 (MQ=255)
                                     aCGACCATAGGCAGGCCAGTCCTGATCGGCTACGGGGGAGATAGCTTCAGCAGGTGTGGCATCGgcgc                             >  1:1745596/1‑68 (MQ=255)
                                       gACCATAGGCAGGCCAGTCCTGATCGGCTACGGGGGAGATAGCTTCAGCAgg                                           <  1:1627714/52‑1 (MQ=255)
                                                     cAGTCCTGATCGGCTACGGGGGAGATAGCTTCAGCAGGTGTGGCATCGGCGCTTAAGGTGCCGTTGAtc            <  1:623050/69‑1 (MQ=255)
                                                     cAGTCCTGATCGGCTACGGGGGAGATAGCTTCAGCAGGTGTGGCATCGGCGCTTAAGGTGCCGTTGAtc            <  1:1099748/69‑1 (MQ=255)
                                                      aGTCCTGATCGGCTACGGGGGAGATAGCTTCAGCAGGTGTGGCATCGGCGCTTAAGGTGCCGTTGAtc            <  1:1144812/68‑1 (MQ=255)
                                                         ccTGATCGGCTACGGGGGAGATAGCTTCAGCAGGTGTGGCATCGGCGCTTAAGGTGCCGTTGATCTCCt        >  1:1915954/1‑69 (MQ=255)
                                                          cTGATCGGCTACGGGGGAGATAGCTTCAGCAGGTGTGGCATCGGCGCTTAAGGTGCCGTTGATCTCCTg       >  1:122985/1‑69 (MQ=255)
                                                             aTCGGCTACGGGGGAGATAGCTTCAGCAGGTGTGGCATCGGCGCTTaa                         >  1:2948590/1‑48 (MQ=255)
                                                             aTCGGCTACGGGGGAGATAGCTTCAGCAGGTGTGGCATCGGCGCTTaa                         >  1:401159/1‑48 (MQ=255)
                                                             aTCGGCTACGGGGGAGATAGCTTCAGCAGGTGTGGCATCGGCGCTTaa                         >  1:1870494/1‑48 (MQ=255)
                                                             aTCGGCTACGGGGGAGATAGCTTCAGCAGGTGTGGCATCGGCGCTTaa                         >  1:1500874/1‑48 (MQ=255)
                                                               cGGCTACGGGGGAGATAGCTTCAGCAGGTGTGGCATCGGCGCTTAAGGTGCCGTTGATCTCCTGCGGAt  >  1:2094718/1‑69 (MQ=255)
                                                                   |                                                                
TTTGTTTCAGCGGCGAAAAGCGTTGACCTTCCTGATTACGACCATAGGCAGGCCAGTCCTGATCGGCTACGGGGGAGATAGCTTCAGCAGGTGTGGCATCGGCGCTTAAGGTGCCGTTGATCTCCTGCGGAT  >  minE/125410‑125541

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: