Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,771,354 C→T 57.3% A308A (GCC→GCT pheA → fused chorismate mutase P and prephenate dehydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,771,3540CT57.3% 20.9 / 62.6 82A308A (GCC→GCTpheAfused chorismate mutase P and prephenate dehydratase
Reads supporting (aligned to +/- strand):  ref base C (20/15);  new base T (16/31);  total (36/46)
Fisher's exact test for biased strand distribution p-value = 4.51e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAGCCATTAACGTGTCTGATCAGGTTCCGGCGAAAACCACGTTGTTAATGGCGACCGGGCAACAAGCCGGTGCGCTGGTTGAAGCGTTGCTGGTACTGCGCAACCACAATCTGATTATG  >  minE/1771287‑1771405
                                                                   |                                                   
aaGCCATTAACGTGTCTGATCAGGTTCCGGCGAAAACCACGTTGTTAATGGCGACCGGGCAACAAGCCg                                                    >  1:3372630/1‑69 (MQ=255)
aaGCCATTAACGTGTCTGATCAGGTTCCGGCGAAAACCACGTTGTTAATGGCGACCGGGCAACAAGCCg                                                    >  1:2224302/1‑69 (MQ=255)
aaGCCATTAACGTGTCTGATCAGGTTCCGGCGAAAACCACGTTGTTAATGGCGACCGGGCAACAAGCCg                                                    >  1:1188807/1‑69 (MQ=255)
aaGCCATTAACGTGTCTGATCAGGTTCCGGCGAAAACCACGTTGTTAATGGCGACCGGGCAACAAGCCg                                                    >  1:1117048/1‑69 (MQ=255)
        aaCGTGTCTGATCAGGTTCCGGCGAAAACCACGTTGTTAATGGCGACCGGGCAACAAGCCGGTGCGCTg                                            >  1:1968513/1‑69 (MQ=255)
        aaCGTGTCTGATCAGGTTCCGGCGAAAACCACGTTGTTAATGGCGACCGGGCAACAAGCCGGTGCGCTg                                            >  1:734919/1‑69 (MQ=255)
        aaCGTGTCTGATCAGGTTCCGGCGAAAACCACGTTGTTAATGGCGACCGGGCAACAAGCCGGTGCGCTg                                            >  1:1095703/1‑69 (MQ=255)
               cTGACCAGGTTCCGGCGAAAACGACGTTGTTAATGGCGACCGGACAACAAGCTGGTgcac                                              <  1:269201/60‑3 (MQ=255)
               cTGACCAGGTTCCGGCGAAAACGACGTTGTTAATGGCGACCGGACAACAAGCTGGTgcac                                              <  1:277542/60‑3 (MQ=255)
               cTGACCAGGTTCCGGCGAAAACGACGTTGTTAATGGCGACCGGACAACAAGCTGGTgcac                                              <  1:1828267/60‑3 (MQ=255)
               cTGACCAGGTTCCGGCGAAAACGACGTTGTTAATGGCGACCGGACAACAAGCTGGTgcac                                              <  1:1741580/60‑3 (MQ=255)
               cTGACCAGGTTCCGGCGAAAACGACGTTGTTAATGGCGACCGGACAACAAGCTGGTgcac                                              <  1:2900506/60‑3 (MQ=255)
               cTGACCAGGTTCCGGCGAAAACGACGTTGTTAATGGCGACCGGACAACAAGCTGGTgcac                                              <  1:1668904/60‑3 (MQ=255)
               cTGACCAGGTTCCGGCGAAAACGACGTTGTTAATGGCGACCGGACAACAAGCTGGTgcac                                              <  1:3178099/60‑3 (MQ=255)
               cTGACCAGGTTCCGGCGAAAACGACGTTGTTAATGGCGACCGGACAACAAGCTGGTgcac                                              <  1:1517701/60‑3 (MQ=255)
               cTGACCAGGTTCCGGCGAAAACGACGTTGTTAATGGCGACCGGACAACAAGCTGGTgcac                                              <  1:2190951/60‑3 (MQ=255)
               cTGACCAGGTTCCGGCGAAAACGACGTTGTTAATGGCGACCGGACAACAAGCTGGTgcac                                              <  1:1967681/60‑3 (MQ=255)
               cTGACCAGGTTCCGGCGAAAACGACGTTGTTAATGGCGACCGGACAACAAGCTGGTgcac                                              <  1:10778/60‑3 (MQ=255)
               cTGACCAGGTTCCGGCGAAAACGACGTTGTTAATGGCGACCGGACAACAAGCTGGTgcac                                              <  1:747767/60‑3 (MQ=255)
               cTGACCAGGTTCCGGCGAAAACGACGTTGTTAATGGCGACCGGACAACAAGCTGGTgcac                                              <  1:836072/60‑3 (MQ=255)
               cTGACCAGGTTCCGGCGAAAACGACGTTGTTAATGGCGACCGGACAACAAGCTGGTgcac                                              <  1:1328209/60‑3 (MQ=255)
               cTGACCAGGTTCCGGCGAAAACGACGTTGTTAATGGCGACCGGACAACAAGCTGGTgcac                                              <  1:1313670/60‑3 (MQ=255)
               cTGACCAGGTTCCGGCGAAAACGACGTTGTTAATGGCGACCGGACAACAAGCTGGTgcac                                              <  1:1251039/60‑3 (MQ=255)
                tGACCAGGTTCCGGCTAAAACGACGTTGTTAATGGCGACCGGACAACAAGCTGGTgcac                                              <  1:1517943/59‑3 (MQ=255)
                tGACCAGGTTCCGGCGAAAACGACGTTGTTAATGGCGACCGGACAACAAGCTGGTgcac                                              <  1:113709/59‑3 (MQ=255)
                    cAGGTTCCGGCGAAAACCACGTTGTTAATGGCGACCGGGCAACAAGCCGGTGCGCTGGTTGAAGCGTTg                                >  1:1457783/1‑69 (MQ=255)
                    cAGGTTCCGGCGAAAACCACGTTGTTAATGGCGACCGGGCAACAAGCCGGTGCGCTGGTTGAAGCGTTg                                >  1:1001385/1‑69 (MQ=255)
                    cAGGTTCCGGCGAAAACCACGTTGTTAATGGCGACCGGGCAACAAGCCGGTGCGCTGGTTGAAGCGTTg                                >  1:883445/1‑69 (MQ=255)
                    cAGGTTCCGGCGAAAACCACGTTGTTAATGGCGACCGGGCAACAAGCCGGTGCGCTGGTTGAAGCGTTg                                >  1:2727331/1‑69 (MQ=255)
                    cAGGTTCCGGCGAAAACCACGTTGTTAATGGCGACCGGGCAACAAGCCGGTGCGCTGGTTGAAGCGTTg                                >  1:1787649/1‑69 (MQ=255)
                    cAGGTTCCGGCGAAAACCACGTTGTTAATGGCGACCGGGCAACAAGCCGGTGCGCTGGTTGAAGCGTTg                                >  1:14571/1‑69 (MQ=255)
                    cAGGTTCCGGCGAAAACCACGTTGTTAATGGCGACCGGGCAACAAGCCGGTGCGCTGGTTGAAGCGTTg                                >  1:432234/1‑69 (MQ=255)
                        ttCCGGCGAAAACCACGTTGTTAATGGCGACCGGGCAACAAGCCGGTGCGCTGGTTGAAGCGTTGCTgg                            <  1:1576004/69‑1 (MQ=255)
                        ttCCGGAGAAAACCACGTTGTTAATGGCGACCGGGCAACAAGCCGGTGCGCTGGTTGAAGCGTTGCTgg                            <  1:1025498/69‑1 (MQ=255)
                          ccGGCGAAAACGACGTTGTTAATGGCGACCGGACAACAAGCTGGTgcac                                              <  1:3250107/49‑3 (MQ=255)
                            ggCGAAAACCACGTTGTTAATGGCGACCGGGCAACAAGCCGGTGCGCTGGTTGAAGCGTTGCTGGTAc                         <  1:2690177/68‑1 (MQ=255)
                               gAAACCGACGTTGTTAATGGCGACCGGACAACAAGCTGGTGCACTGGTTGAAGCGTTGCTGGTACTgc                      >  1:852876/1‑68 (MQ=255)
                               gAAAACGCCGTTGTTAATGGCGACCGGACAACAAGCTGGTGCACTGGTTGAAGCGTTGCTGGTACTgc                      >  1:2597290/1‑68 (MQ=255)
                               gAAAACGACGTTGTTAATGGCGACCGGACAACAAGCTGGTGCACTGGTTGAGGCGTTGCTGGTACTgc                      >  1:2978361/1‑68 (MQ=255)
                               gAAAACGACGTTGTTAATGGCGACCGGACAACAAGCTGGTGCACTGGTTGAAGCGTTGCTGGTACTgc                      >  1:2040256/1‑68 (MQ=255)
                               gAAAACGACGTTGTTAATGGCGACCGGACAACAAGCTGGTGCACTGGTTGAAGCGTTGCTGGTACTgc                      >  1:918222/1‑68 (MQ=255)
                               gAAAACGACGTTGTTAATGGCGACCGGACAACAAGCTGGTGCACTGGTTGAAGCGTTGCTGGTACTgc                      >  1:1478192/1‑68 (MQ=255)
                               gAAAACGACGTTGTTAATGGCGACCGGACAACAAGCTGGTGCACTGGTTGAAGCGTTGCTGGTACTgc                      >  1:1497305/1‑68 (MQ=255)
                               gAAAACGACGTTGTTAATGGCGACCGGACAACAAGCTGGTGCACTGGTTGAAGCGTTGCTGGTACTgc                      >  1:374823/1‑68 (MQ=255)
                               gAAAACGACGTTGTTAATGGCGACCGGACAACAAGCTGGTGCACTGGTTGAAGCGTTGCTGGTACTgc                      >  1:3318248/1‑68 (MQ=255)
                               gAAAACGACGTTGTTAATGGCGACCGGACAACAAGCTGGTGCACTGGTTGAAGCGTTGCTGGTACTgc                      >  1:3187071/1‑68 (MQ=255)
                               gAAAACGACGTTGTTAATGGCGACCGGACAACAAGCTGGTGCACTGGTTGAAGCGTTGCTGGTACTgc                      >  1:2999208/1‑68 (MQ=255)
                               gAAAACGACGTTGTTAATGGCGACCGGACAACAAGCTGGTGCACTGGTTGAAGCGTTGCTGGTACTgc                      >  1:2868517/1‑68 (MQ=255)
                               gAAAACGACGTTGTTAATGGCGACCGGACAACAAGCTGGTGCACTGGTTGAAGCGTTGCTGGTACTgc                      >  1:1850802/1‑68 (MQ=255)
                               gAAAACGACGTTGTTAATGGCGACCGGACAACAAGCTGGTGCACTGGTTGAAGCGTTGCTGGTACTgc                      >  1:2255611/1‑68 (MQ=255)
                               gAAAACGACGTTGTTAATGGCGACCGGACAACAAGCTGGTGCACTGGTTGAAGCGTTGCTGGTACTgc                      >  1:1924985/1‑68 (MQ=255)
                               gAAAACGACGTTGTTAATGGCGACCGGACAACAAGCTGGTGCACTGGTTGAAGCGTTGCTGGTACTg                       >  1:1920307/1‑67 (MQ=255)
                                   aCCACGTTGTTAATGGCGACCGGGCAACAAGCCGGTGCGCTGGTTGAAGCGTTGCTGGTACTGCGCAAc                 <  1:1409380/69‑1 (MQ=255)
                                   aCCACGTTGTTAATGGCGACCGGGCAACAAGCCGGTGCGCTGGTTGAAGCGTTGCTGGTACTGCGCAAc                 <  1:756077/69‑1 (MQ=255)
                                   aCCACGTTGTTAATGGCGACCGGGCAACAAGCCGGTGCGCTGGTTGAAGCGTTGCTGGTACTGCGCAAc                 <  1:1413439/69‑1 (MQ=255)
                                   aCCACGTTGTTAATGGCGACCGGGCAACAAGCCGGTGCGCTGGTTGAAGCGTTGCTGGTACTGCGCAAc                 <  1:2307794/69‑1 (MQ=255)
                                   aCCACGTTGTTAATGGCGACCGGGCAACAAGCCGGTGCGCTGGTTGAAGCGTTGCTGGTACTGCGCAAc                 <  1:1681917/69‑1 (MQ=255)
                                                 ggCGACCGGGCAACAAGCCGGTGCGCTGGTTGAAGCGTTGCTGGTACTGCGCAACCACAATCTGATTAt   >  1:3223971/1‑69 (MQ=255)
                                                 ggCGACCGGGCAACAAGCCGGTGCGCTGGTTGAAGCGTTGCTGGTACTGCGCAACCACAATCTGATTAt   >  1:3280695/1‑69 (MQ=255)
                                                 ggCGACCGGGCAACAAGCCGGTGCGCTGGTTGAAGCGTTGCTGGTACTGCGCAACCACAATCTGATTAt   >  1:2189281/1‑69 (MQ=255)
                                                 ggCGACCGGGCAACAAGCCGGTGCGCTGGTTGAAGCGTTGCTGGTACTGCGCAACCACAATCTGATTAt   >  1:2695919/1‑69 (MQ=255)
                                                 ggCGACCGGGCAACAAGCCGGTGCGCTGGTTGAAGCGTTGCTGGTACTGCGCAACCACAATCTGATTAt   >  1:1055888/1‑69 (MQ=255)
                                                 ggCGACCGGACAACAAGCTGGTGCACTGGTTGAAGCGTTGCTGGTACTGCGCAATCACAGTCTGATTAt   <  1:2158596/69‑1 (MQ=255)
                                                 ggCGACCGGACAACAAGCTGGTGCACTGGTTGAAGCGTTGCTGGTACTGCGCAATCACAGTCTGATTAt   <  1:809228/69‑1 (MQ=255)
                                                 ggCGACCGGACAACAAGCTGGTGCACTGGTTGAAGCGTTGCTGGTACTGCGCAATCACAGTCTGATTAt   <  1:328031/69‑1 (MQ=255)
                                                 ggCGACCGGACAACAAGCTGGTGCACTGGTTGAAGCGTTGCTGGTACTGCGCAATCACAGTCTGATTAt   <  1:1946027/69‑1 (MQ=255)
                                                 ggCGACCGGACAACAAGCTGGTGCACTGGTTGAAGCGTTGCTGGTACTGCGCAATCACAGTCTGATTAt   <  1:2228325/69‑1 (MQ=255)
                                                 ggCGACCGGACAACAAGCTGGTGCACTGGTTGAAGCGTTGCTGGTACTGCGCAATCACAGTCTGATTAt   <  1:2591925/69‑1 (MQ=255)
                                                 ggCGACCGGACAACAAGCTGGTGCACTGGTTGAAGCGTTGCTGGTACTGCGCAATCACAGTCTGATTAt   <  1:2283863/69‑1 (MQ=255)
                                                 ggCGACCGGACAACAAGCTGGTGCACTGGTTGAAGCGTTGCTGGTACTGCGCAATCACAGTCTGATTAt   <  1:2916082/69‑1 (MQ=255)
                                                 ggCGACCGGACAACAAGCTGGTGCACTGGTTGAAGCGTTGCTGGTACTGCGCAATCACAGTCTGATTAt   <  1:29021/69‑1 (MQ=255)
                                                 ggCGACCGGACAACAAGCTGGTGCACTGGTTGAAGCGTTGCTGGTACTGCGCAATCACAGTCTGATTAt   <  1:2837372/69‑1 (MQ=255)
                                                 ggCGACCGGACAACAAGCTGGTGCACTGGTTGAAGCGTTGCTGGTACTGCGCAATCACAGTCTGATTAt   <  1:2804956/69‑1 (MQ=255)
                                                 ggCGACCGGACAACAAGCTGGTGCACTGGTTGAAGCGTTGCTGGTACTGCGCAATCACAGTCTGATTAt   <  1:1841868/69‑1 (MQ=255)
                                                  gCGACCGGGCAACAAGCCGGTGCGCTGGTTGAAGCGTTGCTGGTACTGCGCAACCACAATCTGATTATg  <  1:1520223/69‑1 (MQ=255)
                                                  gCGACCGGGCAACAAGCCGGTGCGCTGGTTGAAGCGTTGCTGGTACTGCGCAACCACAATCTGATTATg  <  1:3163087/69‑1 (MQ=255)
                                                  gCGACCGGGCAACAAGCCGGTGCGCTGGTTGAAGCGTTGCTGGTACTGCGCAACCACAATCTGATTATg  <  1:1505286/69‑1 (MQ=255)
                                                  gCGACCGGGCAACAAGCCGGTGCGCTGGTTGAAGCGTTGCTGGTACTGCGCAACCACAATCTGATTATg  <  1:2153870/69‑1 (MQ=255)
                                                  gCGACCGGGCAACAAGCCGGTGCGCTGGTTGAAGCGTTGCTGGTACTGCGCAACCACAATCTGATTATg  <  1:1829095/69‑1 (MQ=255)
                                                  gCGACCGGGCAACAAGCCGGTGCGCTGGTTGAAGCGTTGCTGGTACTGCGCAACCACAATCTGATTATg  <  1:2103613/69‑1 (MQ=255)
                                                  gCGACCGGGCAACAAGCCGGTGCGCTGGTTGAAGCGTTGCTGGTACTGCGCAACCACAATCTGATTATg  >  1:2359305/1‑69 (MQ=255)
                                                  gCGACCGGGCAACAAGCCGGTGCGCTGGTTGAAGCGTTGCTGGTACTGCGCAACCACAATCTGATTATg  <  1:2576352/69‑1 (MQ=255)
                                                                   |                                                   
AAGCCATTAACGTGTCTGATCAGGTTCCGGCGAAAACCACGTTGTTAATGGCGACCGGGCAACAAGCCGGTGCGCTGGTTGAAGCGTTGCTGGTACTGCGCAACCACAATCTGATTATG  >  minE/1771287‑1771405

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: