Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,888,190 G→A 62.0% R226R (CGC→CGT gudD ← (D)‑glucarate dehydratase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,888,1900GA62.0% 34.1 / 45.2 71R226R (CGC→CGTgudD(D)‑glucarate dehydratase 1
Reads supporting (aligned to +/- strand):  ref base G (13/14);  new base A (27/17);  total (40/31)
Fisher's exact test for biased strand distribution p-value = 3.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.33e-01

ATTTTAATCGCTTCGTTCAGCGACCAGGCACCGTTAGGATCGAGCGTAATACGCGCCTGCGGGAAGCGTTGCGCCAGTGCCACAATAGACTCGGCCTCTTCTTCCCCGGCCAGTACACCGCCCTTCA  >  minE/1888125‑1888251
                                                                 |                                                             
aTTTTAATCGCTTCGTTCAGCGACCAGGCACCGTTAGGATCGAGCGTAATACGCGCCTGCGGGAAACGt                                                            >  1:76382/1‑69 (MQ=255)
aTTTTAATCGCTTCGTTCAGCGACCAGGCACCGTTAGGATCGAGCGTAATACGCGCCTGCGGGAAACGt                                                            >  1:2595567/1‑69 (MQ=255)
aTTTTAATCGCTTCGTTCAGCGACCAGGCACCGTTAGGATCGAGCGTAATACGCGCCTGCGGGAAACGt                                                            >  1:2887166/1‑69 (MQ=255)
aTTTTAATCGCTTCGTTCAGCGACCAGGCACCGTTAGGATCGAGCGTAATACGCGCCTGCGGGAAACGt                                                            >  1:3038343/1‑69 (MQ=255)
aTTTTAATCGCTTCGTTCAGCGACCAGGCACCGTTAGGATCGAGCGTAATACGCGCCTGCGGGAAACGt                                                            >  1:3088916/1‑69 (MQ=255)
aTTTTAATCGCTTCGTTCAGCGACCAGGCACCGTTAGGATCGAGCGTAATACGCGCCTGCGGGAAACGt                                                            >  1:588251/1‑69 (MQ=255)
aTTTTAATCGCTTCGTTCAGCGACCAGGCACCGTTAGGATCGAGCGTAATACGCGCCTGCGGGAAACGt                                                            >  1:665755/1‑69 (MQ=255)
aTTTTAATCGCTTCGTTCAGCGACCAGGCACCGTTAGGATCGAGCGTAATACGCGCCTGCGGGAAACGt                                                            >  1:815349/1‑69 (MQ=255)
aTTTTAATCGCTTCGTTCAGCGACCAGGCACCGTTAGGATCGAGCGTAATACGCGCCTGCGGGAAACGt                                                            >  1:1713762/1‑69 (MQ=255)
aTTTTAATCGCTTCGTTCAGCGACCAGGCACCGTTAGGATCGAGCGTAATACGCGCCTGCGGGAAACGt                                                            >  1:1542418/1‑69 (MQ=255)
aTTTTAATCGCTTCGTTCAGCGACCAGGCACCGTTAGGATCGAGCGTAATACGCGCCTGCGGGAAACGt                                                            >  1:1223995/1‑69 (MQ=255)
aTTTTAATCGCTTCGTTCAGCGACCAGGCACCGTTAGGATCGAGCGTAATACGCGCCTGCGGGAAACGt                                                            >  1:1241129/1‑69 (MQ=255)
aTTTTAATCGCTTCGTTCAGCGACCAGGCACCGTTAGGATCGAGCGTAATACGCGCCTGCGGGAAACGt                                                            >  1:972570/1‑69 (MQ=255)
aTTTTAATCGCTTCGTTCAGCGACCAGGCACCGTTAGGATCGAGCGTAATACGCGCCTGCGGGAAACGt                                                            >  1:1311610/1‑69 (MQ=255)
aTTTTAATCGCTTCGTTCAGCGACCAGGCACCGTTAGGATCGAGCGTAATACGCGCCTGCGGGAAACGt                                                            >  1:1369268/1‑69 (MQ=255)
  tttAATCGCTTCGTTCAGCGACCAGGCACCGTTAGGATCGAGCGTAGTACGCGCGTGCGGGAAGCGTTg                                                          >  1:3145944/1‑69 (MQ=255)
  tttAATCGCTTCGTTCAGCGACCAGGCACCGTTAGGATCGAGCGTAATACGCGCCTGCGGGAAGCGTTg                                                          >  1:887547/1‑69 (MQ=255)
  tttAATCGCTTCGTTCAGCGACCAGGCACCGTTAGGATCGAGCGTAATACGCGCCTGCGGGAAGCGTTg                                                          >  1:1047675/1‑69 (MQ=255)
  tttAATCGCTTCGTTCAGCGACCAGGCACCGTTAGGATCGAGCGTAATACGCGCCTGCGGGAAGCGTTg                                                          >  1:313832/1‑69 (MQ=255)
  tttAATCGCTTCGTTCAGCGACCAGGCACCGTTAGGATCGAGCGTAATACGCGCCTGCGGGAAGCGTTg                                                          >  1:3243567/1‑69 (MQ=255)
  tttAATCGCTTCGTTCAGCGACCAGGCACCGTTAGGATCGAGCGTAATACGCGCCTGCGGGAAGCGTTg                                                          >  1:1999350/1‑69 (MQ=255)
  tttAATCGCTTCGTTCAGCGACCAGGCACCGTTAGGATCGAGCGTAATACGCGCCTGCGGGAAGCGTTg                                                          >  1:1699163/1‑69 (MQ=255)
  tttAATCGCTTCGTTCAGCGACCAGGCACCGTTAGGATCGAGCGTAATACGCGCCTGCGGGAAGCGTTg                                                          >  1:869379/1‑69 (MQ=255)
       tCGCTTCGTTCAGCGACCAGGCACCGTTAGGATCGAGCGTAATACGCGCCTGCGGGAAGCGTTGCGcc                                                      <  1:522727/68‑1 (MQ=255)
       tCGCTTCGTTCAGCGACCAGGCACCGTTAGGATCGAGCGTAATACGCGCCTGCGGGAAGCGTTGCGcc                                                      <  1:2519443/68‑1 (MQ=255)
         gCTTCGTTCAGCGACCAGGCACCGTTAGGATCGAGCGTAATACGCGCCTGCGGGAAGCGTTGCGCCAGt                                                   <  1:2335941/69‑1 (MQ=255)
                   gCGACCAGGCACCGTTAGGATTGAGCGTAATACGCGCCTGCGGGAAACGTTGCGCCAGTGCCACAATAg                                         >  1:3127169/1‑69 (MQ=255)
                   gCGACCAGGCACCGTTAGGATCGAGCGTAATACGCGCCTGCGGGAAACGTTGCGCCAGTGCCACAATAg                                         >  1:49252/1‑69 (MQ=255)
                   gCGACCAGGCACCGTTAGGATCGAGCGTAATACGCGCCTGCGGGAAACGTTGCGCCAGTGCCACAATAg                                         >  1:1029251/1‑69 (MQ=255)
                   gCGACCAGGCACCGTTAGGATCGAGCGTAATACGCGCCTGCGGGAAACGTTGCGCCAGTGCCACAATAg                                         >  1:97654/1‑69 (MQ=255)
                   gCGACCAGGCACCGTTAGGATCGAGCGTAATACGCGCCTGCGGGAAACGTTGCGCCAGTGCCACAATAg                                         >  1:893937/1‑69 (MQ=255)
                   gCGACCAGGCACCGTTAGGATCGAGCGTAATACGCGCCTGCGGGAAACGTTGCGCCAGTGCCACAATAg                                         >  1:2750291/1‑69 (MQ=255)
                   gCGACCAGGCACCGTTAGGATCGAGCGTAATACGCGCCTGCGGGAAACGTTGCGCCAGTGCCACAATAg                                         >  1:2143294/1‑69 (MQ=255)
                   gCGACCAGGCACCGTTAGGATCGAGCGTAATACGCGCCTGCGGGAAACGTTGCGCCAGTGCCACAATAg                                         >  1:2567635/1‑69 (MQ=255)
                   gCGACCAGGCACCGTTAGGATCGAGCGTAATACGCGCCTGCGGGAAACGTTGCGCCAGTGCCACAATAg                                         >  1:1734451/1‑69 (MQ=255)
                   gCGACCAGGCACCGTTAGGATCGAGCGTAATACGCGCCTGCGGGAAACGTTGCGCCAGTGCCACAATAg                                         >  1:225136/1‑69 (MQ=255)
                   gCGACCAGGCACCGTTAGGATCGAGCGTAATACGCGCCTGCGGGAAACGTTGCGCCAGTGCCACAATAg                                         >  1:1662693/1‑69 (MQ=255)
                   gCGACCAGGCACCGTTAGGATCGAGCGTAATACGCGCCTGCGGGAAACGTTGCGACAGTGCCACAATAg                                         >  1:2428051/1‑69 (MQ=255)
                             aCCGTTAGGATCGAGCGTAATACGCGCCTGCGGGAAGCGTTGCGCCAGTGCCACAATAGACTCGGCctc                               <  1:1645594/69‑1 (MQ=255)
                             aCCGTTAGGATCGAGCGTAATACGCGCCTGCGGGAAGCGTTGCGCCAGTGCCACAATAGACTCGGCctc                               <  1:2098724/69‑1 (MQ=255)
                             aCCGTTAGGATCGAGCGTAATACGCGCCTGCGGGAAGCGTTGCGCCAGTGCCACAATAGACTCGGCctc                               <  1:2012780/69‑1 (MQ=255)
                             aCCGTTAGGATCGAGCGTAATACGCGCCTGCGGGAAGCGTTGCGCCAGTGCCACAATAGACTCGGCctc                               <  1:2215382/69‑1 (MQ=255)
                             aCCGTTAGGATCGAGCGTAATACGCGCCTGCGGGAAGCGTTGCGCCAGTGCCACAATAGACTCGGCctc                               <  1:2843319/69‑1 (MQ=255)
                             aCCGTTAGGATCGAGCGTAATACGCGCCTGCGGGAAGCGTTGCGCCAGTGCCACAATAGACTCGGCctc                               <  1:2738846/69‑1 (MQ=255)
                             aCCGTTAGGATCGAGCGTAATACGCGCCTGCGGGAAGCGTTGCGCCAGTGCCACAATAGACTCGGCctc                               <  1:2778588/69‑1 (MQ=255)
                             aCCGTTAGGATCGAGCGTAATACGCGCCTGCGGGAAGCGTTGCGCCAGTGCCACAATAGACTCGGCctc                               <  1:2976730/69‑1 (MQ=255)
                             aCCGTTAGGATCGAGCGTAATACGCGCCTGCGGGAAGCGTTGCGCCAGTGCCACAATAGACTCGGCctc                               <  1:948655/69‑1 (MQ=255)
                                        cGAGCGTAATACGCGCCTGCGGGAAGCGTTGCGCCAGTGCCACAATAGACTCGGCCTCTTCTTCCCCg                     >  1:3081676/1‑68 (MQ=255)
                                        cGAGCGTAATACGCGCCTGCGGGAAGCGTTGCGCCAGTGCCACAATAGACTCGGCCTCTTCTTCCCCg                     >  1:7610/1‑68 (MQ=255)
                                        cGAGCGTAATACGCGCCTGCGGGAAGCGTTGCGCCAGTGCCACAATAGACTCGGCCTCTTCTTCCCCg                     >  1:3085709/1‑68 (MQ=255)
                                           gCGTAATACGCGCCTGCGGGAAGCGTTGCGCCAGTGCCACAATAGACTCGGCCTCTTCTTCCCCGGcc                  <  1:3223069/68‑1 (MQ=255)
                                           gCGTAATACGCGCCTGCGGGAAGCGTTGCGCCAGTGCCACAATAGACTCGGCCTCTTCTTCCCCGGcc                  <  1:2190605/68‑1 (MQ=255)
                                                   cgcgCCTGCGGGAAGCGTTGCGCCAGTGCCACAATAGACTCGGCCTCTTCTTCCCCGGCCAGTACACCg         >  1:2522184/1‑69 (MQ=255)
                                                   cgcgCCTGCGGGAAGCGTTGCGCCAGTGCCACAATAGACTCGGCCTCTTCTTCCCCGGCCAGTACACCg         >  1:58258/1‑69 (MQ=255)
                                                          gCGGGAAACGTTGCGCCAGTGCCACAATAGATTCGGCTTCTTCTTCCCCGGCCAGTACACCGCCTTTCa  <  1:1476181/69‑1 (MQ=255)
                                                          gCGGGAAACGTTGCGCCAGTGCCACAATAGATTCGGCTTCTTCTTCCCCGGCCAGTACACCGCCTTTCa  <  1:1510826/69‑1 (MQ=255)
                                                          gCGGGAAACGTTGCGCCAGTGCCACAATAGATTCGGCTTCTTCTTCCCCGGCCAGTACACCGCCTTTCa  <  1:80289/69‑1 (MQ=255)
                                                          gCGGGAAACGTTGCGCCAGTGCCACAATAGATTCGGCTTCTTCTTCCCCGGCCAGTACACCGCCTTTCa  <  1:1514735/69‑1 (MQ=255)
                                                          gCGGGAAACGTTGCGCCAGTGCCACAATAGATTCGGCTTCTTCTTCCCCGGCCAGTACACCGCCTTTCa  <  1:746257/69‑1 (MQ=255)
                                                          gCGGGAAACGTTGCGCCAGTGCCACAATAGATTCGGCTTCTTCTTCCCCGGCCAGTACACCGCCTTTCa  <  1:217060/69‑1 (MQ=255)
                                                          gCGGGAAACGTTGCGCCAGTGCCACAATAGATTCGGCTTCTTCTTCCCCGGCCAGTACACCGCCTTTCa  <  1:49082/69‑1 (MQ=255)
                                                          gCGGGAAACGTTGCGCCAGTGCCACAATAGATTCGGCTTCTTCTTCCCCGGCCAGTACACCGCCTTTCa  <  1:3385967/69‑1 (MQ=255)
                                                          gCGGGAAACGTTGCGCCAGTGCCACAATAGATTCGGCTTCTTCTTCCCCGGCCAGTACACCGCCTTTCa  <  1:3271018/69‑1 (MQ=255)
                                                          gCGGGAAACGTTGCGCCAGTGCCACAATAGATTCGGCTTCTTCTTCCCCGGCCAGTACACCGCCTTTCa  <  1:3252117/69‑1 (MQ=255)
                                                          gCGGGAAACGTTGCGCCAGTGCCACAATAGATTCGGCTTCTTCTTCCCCGGCCAGTACACCGCCTTTCa  <  1:3224267/69‑1 (MQ=255)
                                                          gCGGGAAACGTTGCGCCAGTGCCACAATAGATTCGGCTTCTTCTTCCCCGGCCAGTACACCGCCTTTCa  <  1:226813/69‑1 (MQ=255)
                                                          gCGGGAAACGTTGCGCCAGTGCCACAATAGATTCGGCTTCTTCTTCCCCGGCCAGTACACCGCCTTTCa  <  1:2341766/69‑1 (MQ=255)
                                                          gCGGGAAACGTTGCGCCAGTGCCACAATAGATTCGGCTTCTTCTTCCCCGGCCAGTACACCGCCTTTCa  <  1:2359783/69‑1 (MQ=255)
                                                          gCGGGAAACGTTGCGCCAGTGCCACAATAGATTCGGCTTCTTCTTCCCCGGCCAGTACACCGCCTTTCa  <  1:3074742/69‑1 (MQ=255)
                                                          gCGGGAAACGTTGCGCCAGTGCCACAATAGATTCGGCTTCTTCTTCCCCGGCCAGTACACCGCCTTTCa  <  1:2909454/69‑1 (MQ=255)
                                                          gCGGGAAACGTTGCGCCAGTGCCACAATAGATTCGGCTTCTTCTTCCCCGGCCAGTACACCGCCTTTCa  <  1:247733/69‑1 (MQ=255)
                                                                 |                                                             
ATTTTAATCGCTTCGTTCAGCGACCAGGCACCGTTAGGATCGAGCGTAATACGCGCCTGCGGGAAGCGTTGCGCCAGTGCCACAATAGACTCGGCCTCTTCTTCCCCGGCCAGTACACCGCCCTTCA  >  minE/1888125‑1888251

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: