Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,959,154 G→T 57.5% A306E (GCG→GAG) ‡ fbaA ← fructose‑bisphosphate aldolase, class II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,959,1540GT57.5% 19.5 / 83.0 80A306E (GCG→GAG) ‡fbaAfructose‑bisphosphate aldolase, class II
Reads supporting (aligned to +/- strand):  ref base G (14/20);  new base T (24/22);  total (38/42)
Fisher's exact test for biased strand distribution p-value = 3.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTCTTGTTCGGCTGATCTTCGCCTTTCGGGTTACCCAGCTGACCCTGCAGATAAGCTTCGTTCGCTTTGTAGTAGTTCAGAACGCCTTCCCAGGTTGCCCATTGGGTATCGGTATCGATGTTC  >  minE/1959090‑1959213
                                                                |                                                           
tttCTTGTTCGGCTGATCTTCGCCTTTCGGGTTACCCAGCTGACCCTGCAGATAAGCTTCGTTCGCtt                                                          >  1:1144551/1‑68 (MQ=255)
tttCTTGTTCGGCTGATCTTCGCCTTTCGGGTTACCCAGCTGACCCTGCAGATAAGCTTCGTTCGCtt                                                          >  1:2225171/1‑68 (MQ=255)
tttCTTGTTCGGCTGATCTTCGCCTTTCGGGTTACCCAGCTGACCCTGCAGATAAGCTTCGTTCGCtt                                                          >  1:2209122/1‑68 (MQ=255)
tttCTTGTTCGGCTGATCTTCGCCTTTCGGGTTACCCAGCTGACCCTGCAGATAAGCTTCGTTCGCtt                                                          >  1:2589353/1‑68 (MQ=255)
tttCTTGTTCGGCTGATCTTCGCCTTTCGGGTTACCCAGCTGACCCTGCAGATAAGCTTCGTTCGCtt                                                          >  1:3364083/1‑68 (MQ=255)
tttCTTGTTCGGCTGATCTTCGCCTTTCGGGTTACCCAGCTGACCCTGCAGATAAGCTTCGTTCGCtt                                                          >  1:1772682/1‑68 (MQ=255)
tttCTTGTTCGGCTGATCTTCGCCTTTCGGGTTACCCAGCTGACCCTGCAGATAAGCTTCGTTCGCtt                                                          >  1:1622108/1‑68 (MQ=255)
tttCTTGTTCGGCTGATCTTCGCCTTTCGGGTTACACAGCTGACCCTGCAGATAAGCTTCGTTCGCtt                                                          >  1:1439971/1‑68 (MQ=255)
  tCTTGTTCGGCTGATCTTCGCCTTTCGGGTTACCCAGCTGACCCTGCAGATAAGCTTCGTTCGGTTTg                                                        <  1:2959711/68‑1 (MQ=255)
  tCTTGTTCGGCTGATCTTCGCCTTTCGGGTTACCCAGCTGACCCTGCAGATAAGCTTCGTTCGCTTTg                                                        <  1:1816890/68‑1 (MQ=255)
  tCTTGTTCGGCTGATCTTCGCCTTTCGGGTTACCCAGCTGACCCTGCAGATAAGCTTCGTTCGCTTTg                                                        <  1:2396293/68‑1 (MQ=255)
  tCTTGTTCGGCTGATCTTCGCCTTTCGGGTTACCCAGCTGACCCTGCAGATAAGCTTCGTTCGCTTTg                                                        <  1:1846286/68‑1 (MQ=255)
  tCTTGTTCGGCTGATCTTCGCCTTTCGGGTTACCCAGCTGACCCTGCAGATAAGCTTCGTTCGCTTTg                                                        <  1:1687158/68‑1 (MQ=255)
                 cTTCGCCTTTCGGGTTTCCCAGCTGACCCTGCAGATAAGCTTCGTTCGCTTTGTAGTAGTTCAGAACg                                         <  1:276956/68‑1 (MQ=255)
                 cTTCGCCTTTCGGGTTACCCAGCTGACCCTGCAGATAAGCTTCGTTCGCtt                                                          <  1:2374934/51‑1 (MQ=255)
                 cTTCGCCTTTCGGGTTACCCAGCTGACCCTGCAGATAAGCTTCGTTCGCtt                                                          <  1:1169303/51‑1 (MQ=255)
                 cTTCGCCTTTCGGGTTACCCAGCTGACCCTGCAGATAAGCTTCGTTCGCtt                                                          <  1:489486/51‑1 (MQ=255)
                 cTTCGCCTTTCGGGTTACCCAGCTGACCCTGCAGATAAGCTTCGTTCGCtt                                                          <  1:453229/51‑1 (MQ=255)
                 cTTCGCCTTTCGGGTTACCCAGCTGACCCTGCAGATAAGCTTCGTTCGCtt                                                          <  1:255068/51‑1 (MQ=255)
                 cTTCGCCTTTCGGGTTACCCAGCTGACCCTGCAGATAAGCTTCGTTCGCtt                                                          <  1:780283/51‑1 (MQ=255)
                 cTTCGCCTTTCGGGTTACCCAGCTGACCCTGCAGATAAGCTTCGTTCGCtt                                                          <  1:1898485/51‑1 (MQ=255)
                 cTTCGCCTTTCGGGTTACCCAGCTGACCCTGCAGATAAGCTTCGTTCGCtt                                                          <  1:2967412/51‑1 (MQ=255)
                 cTTCGCCTTTCGGGTTACCCAGCTGACCCTGCAGATAAGCTTCGTTCGCtt                                                          <  1:2957927/51‑1 (MQ=255)
                 cTTCGCCTTTCGGGTTACCCAGCTGACCCTGCAGATAAGCTTCGTTCGCTTTgt                                                       >  1:1239236/1‑54 (MQ=255)
                 cTTCGCCTTTCGGGTTACCCAGCTGACCCTGCAGATAAGCTTCGTTCGCTTTgt                                                       >  1:2734838/1‑54 (MQ=255)
                 cTTCGCCTTTCGGGTTACCCAGCTGACCCTGCAGATAAGCTTCGTTCGCTTTGTAGTAGTTCAGAACg                                         <  1:726915/68‑1 (MQ=255)
                 cTTCGCCTTTCGGGTTACCCAGCTGACCCTGCAGATAAGCTTCGTTCGCTTTGTAGTAGTTCAGAACg                                         <  1:380759/68‑1 (MQ=255)
                 cTTCGCCTTTCGGGTTACCCAGCTGACCCTGCAGATAAGCTTCGTTCGCTTTGTAGTAGTTCAGAACg                                         <  1:1626444/68‑1 (MQ=255)
                                        tGACCCTGCAGATAAGCTTCGTTCGCTTTGTAGTAGTTCAGAACg                                         <  1:210232/45‑1 (MQ=255)
                                               gCAGGTGAGCTTCGTTTTCTTTGTAGTAGTTCAGAACACCTTCCCAGGTTGCCCa                        >  1:7250/1‑55 (MQ=255)
                                               gCAGGTAAGCTTCGTTTTCTTTTTAGTAGTTCAGAACACCTTCCCAGGTTGCCCa                        >  1:2980771/1‑55 (MQ=255)
                                               gCAGGTAAGCTTCGTTTTCTTTGTAGTAGTTCAGAACACCTTCCCAGGTTGCCCa                        >  1:973566/1‑55 (MQ=255)
                                               gCAGGTAAGCTTCGTTTTCTTTGTAGTAGTTCAGAACACCTTCCCAGGTTGCCCa                        >  1:89837/1‑55 (MQ=255)
                                               gCAGGTAAGCTTCGTTTTCTTTGTAGTAGTTCAGAACACCTTCCCAGGTTGCCCa                        >  1:2929410/1‑55 (MQ=255)
                                               gCAGGTAAGCTTCGTTTTCTTTGTAGTAGTTCAGAACACCTTCCCAGGTTGCCCa                        >  1:2787165/1‑55 (MQ=255)
                                               gCAGGTAAGCTTCGTTTTCTTTGTAGTAGTTCAGAACACCTTCCCAGGTTGCCCa                        >  1:999224/1‑55 (MQ=255)
                                               gCAGGTAAGCTTCGTTTTCTTTGTAGTAGTTCAGAACACCTTCCCAGGTTGCCCa                        >  1:3434956/1‑55 (MQ=255)
                                               gCAGGTAAGCTTCGTTTTCTTTGTAGTAGTTCAGAACACCTTCCCAGGTTGCCCa                        >  1:733216/1‑55 (MQ=255)
                                               gCAGGTAAGCTTCGTTTTCTTTGTAGTAGTTCAGAACACCTTCCCAGGTTGCCCa                        >  1:2529852/1‑55 (MQ=255)
                                               gCAGGTAAGCTTCGTTTTCTTTGTAGTAGTTCAGAACACCTTCCCAGGTTGCCCa                        >  1:982606/1‑55 (MQ=255)
                                               gCAGGTAAGCTTCGTTTTCTTTGTAGTAGTTCAGAACACCTTCCCAGGTTGCCCa                        >  1:1042252/1‑55 (MQ=255)
                                               gCAGGTAAGCTTCGTTTTCTTTGTAGTAGTTCAGAACACCTTCCCAGGTTGCCCa                        >  1:2278199/1‑55 (MQ=255)
                                               gCAGGTAAGCTTCGTTTTCTTTGTAGTAGTTCAGAACACCTTCCCAGGTTGCCCa                        >  1:2233008/1‑55 (MQ=255)
                                               gCAGGTAAGCTTCGTTTTCTTTGTAGTAGTTCAGAACACCTTCCCAGGTTGCCCa                        >  1:1313681/1‑55 (MQ=255)
                                               gCAGGTAAGCTTCGTTTTCTTTGTAGTAGTTCAGAACACCTTCCCAGGTTGCCCa                        >  1:1735115/1‑55 (MQ=255)
                                               gCAGGTAAGCTTCGTTTTCTTTGTAGTAGTTCAGAACACCTTCCCAGGTTGCCCa                        >  1:1771902/1‑55 (MQ=255)
                                               gCAGGTAAGCTTCGTTTTCTTTGTAGTAGTTCAGAACACCTTCCCAGGTTGCCCa                        >  1:1818750/1‑55 (MQ=255)
                                               gCAGGTAAGCTTCGTTTTCTTTGTAGTAGTTCAGAACACCTTCCCAGGTTGCCCa                        >  1:1838481/1‑55 (MQ=255)
                                               gCAGGTAAGCTTCGTTTTCTTTGTAGTAGTTCAGAACACCTTCCCAGGTTGCCCa                        >  1:1855275/1‑55 (MQ=255)
                                               gCAGGTAAGCTTCGTTTTCTTTGTAGTAGTTCAGAACACCTTCCCAGGTTGCCCa                        >  1:1905179/1‑55 (MQ=255)
                                               gCAGGTAAGCTTCGTTTTCTTTGTAGTAGTTCAGAACACCTTCCCAGGTTGCCCa                        >  1:200322/1‑55 (MQ=255)
                                               gCAGGTAAGCTTCGTTTTCTTTGTAGTAGTTCAGAACACCTTCCCAGGTTGCCCa                        >  1:2029302/1‑55 (MQ=255)
                                               gCAGGTAAGCTTCGTTTTCTTTGTAGTAGTTCAGAACACCTTCCCAGGTTGCCCa                        >  1:2261950/1‑55 (MQ=255)
                                                cAGGTAAGCTTCGTTTTCTTTGTAGTAGTTCAGAACACCTTCCCAGGTTGCCCAtt                      <  1:218172/56‑1 (MQ=255)
                                                  gATAAGCTTCGTTCGCTTTGTAGTAGTTCAGAACGCCTTCCCAGGTTGCCCATTGggtatcggtatc         <  1:293455/67‑1 (MQ=255)
                                                      aGCTTCGTTCGCTTTGTAGTAGTTCAGAACGCCTTCCCAGGTTGCCCATTGGGTATCGGTATCGATGtt   >  1:289341/1‑69 (MQ=255)
                                                      aGCTTCGTTCGCTTTGTAGTAGTTCAGAACGCCTTCCCAGGTTGCCCATTGGGTATCGGTATCGATGtt   >  1:1259209/1‑69 (MQ=255)
                                                      aGCTTCGTTCGCTTTGTAGTAGTTCAGAACGCCTTCCCAGGTTGCCCATTGGGTATCGGTATCGATGtt   >  1:1822309/1‑69 (MQ=255)
                                                      aGCTTCGTTCGCTTTGTAGTAGTTCAGAACGCCTTCCCAGGTTGCCCATTGGGTATCGGTATCGATGtt   >  1:1791243/1‑69 (MQ=255)
                                                        cTTCGTTTTCTTTGTAGTAGTTCAGAACACCTTCCCAGGTTGCCCATTGGGTATCGGTATCGATGTTc  <  1:1230458/68‑1 (MQ=255)
                                                        cTTCGTTTTCTTTGTAGTAGTTCAGAACACCTTCCCAGGTTGCCCATTGGGTATCGGTATCGATGTTc  <  1:2600841/68‑1 (MQ=255)
                                                        cTTCGTTTTCTTTGTAGTAGTTCAGAACACCTTCCCAGGTTGCCCATTGGGTATCGGTATCGATGTTc  <  1:836738/68‑1 (MQ=255)
                                                        cTTCGTTTTCTTTGTAGTAGTTCAGAACACCTTCCCAGGTTGCCCATTGGGTATCGGTATCGATGTTc  <  1:1397035/68‑1 (MQ=255)
                                                        cTTCGTTTTCTTTGTAGTAGTTCAGAACACCTTCCCAGGTTGCCCATTGGGTATCGGTATCGATGTTc  <  1:1526647/68‑1 (MQ=255)
                                                        cTTCGTTTTCTTTGTAGTAGTTCAGAACACCTTCCCAGGTTGCCCATTGGGTATCGGTATCGATGTTc  <  1:1594137/68‑1 (MQ=255)
                                                        cTTCGTTTTCTTTGTAGTAGTTCAGAACACCTTCCCAGGTTGCCCATTGGGTATCGGTATCGATGTTc  <  1:1714090/68‑1 (MQ=255)
                                                        cTTCGTTTTCTTTGTAGTAGTTCAGAACACCTTCCCAGGTTGCCCATTGGGTATCGGTATCGATGTTc  <  1:618753/68‑1 (MQ=255)
                                                        cTTCGTTTTCTTTGTAGTAGTTCAGAACACCTTCCCAGGTTGCCCATTGGGTATCGGTATCGATGTTc  <  1:1717966/68‑1 (MQ=255)
                                                        cTTCGTTTTCTTTGTAGTAGTTCAGAACACCTTCCCAGGTTGCCCATTGGGTATCGGTATCGATGTTc  <  1:1179747/68‑1 (MQ=255)
                                                        cTTCGTTTTCTTTGTAGTAGTTCAGAACACCTTCCCAGGTTGCCCATTGGGTATCGGTATCGATGTTc  <  1:2445671/68‑1 (MQ=255)
                                                        cTTCGTTTTCTTTGTAGTAGTTCAGAACACCTTCCCAGGTTGCCCATTGGGTATCGGTATCGATGTTc  <  1:3422099/68‑1 (MQ=255)
                                                        cTTCGTTTTCTTTGTAGTAGTTCAGAACACCTTCCCAGGTTGCCCATTGGGTATCGGTATCGATGTTc  <  1:2090839/68‑1 (MQ=255)
                                                        cTTCGTTTTCTTTGTAGTAGTTCAGAACACCTTCCCAGGTTGCCCATTGGGTATCGGTATCGATGTTc  <  1:2083953/68‑1 (MQ=255)
                                                        cTTCGTTTTCTTTGTAGTAGTTCAGAACACCTTCCCAGGTTGCCCATTGGGTATCGGTATCGATGTTc  <  1:1845109/68‑1 (MQ=255)
                                                        cTTCGTTTTCTTTGTAGTAGTTCAGAACACCTTCCCAGGTTGCCCATTGGGTATCGGTATCGATGTTc  <  1:2584257/68‑1 (MQ=255)
                                                        cTTCGTTTTCTTTGTAGTAGTTCAGAACACCTTCCCAGGTTGCCCATTGGGTATCGGTATCGATGTTc  <  1:18601/68‑1 (MQ=255)
                                                        cTTCGTTTTCTTTGTAGTAGTTCAGAACACCTTCCCAGGTTGCCCATTGGGTATCGGTATCGATGTTc  <  1:29207/68‑1 (MQ=255)
                                                        cTTCGTTTTCTTTGTAGTAGTTCAGAACACCTTCCCAGGTTGCCCATTGGGTATCGGTATCGATGTTc  <  1:2792009/68‑1 (MQ=255)
                                                        cTTCGTTTTCTTTGTAGTAGTTCAGAACACCTTCCCAGGTTGCCCATTGGGTATCGGTATCGATGTTc  <  1:1893435/68‑1 (MQ=255)
                                                        cTTCGTTTTCTTTGTAGTAGTTCAGAACACCTTACCAGGTTGCCCATTGGGTATCGGTATCGATGTTc  <  1:3290422/68‑1 (MQ=255)
                                                                |                                                           
TTTCTTGTTCGGCTGATCTTCGCCTTTCGGGTTACCCAGCTGACCCTGCAGATAAGCTTCGTTCGCTTTGTAGTAGTTCAGAACGCCTTCCCAGGTTGCCCATTGGGTATCGGTATCGATGTTC  >  minE/1959090‑1959213

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: