Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 201,981 T→C 29.3% D108D (GAT→GAC gmhB → D,D‑heptose 1,7‑bisphosphate phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE201,9810TC29.3% 80.9 / 38.7 85D108D (GAT→GACgmhBD,D‑heptose 1,7‑bisphosphate phosphatase
Reads supporting (aligned to +/- strand):  ref base T (12/48);  new base C (7/18);  total (19/66)
Fisher's exact test for biased strand distribution p-value = 4.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.62e-01

ATCTATTATTGCCCGCATCATCCGCAGGGTAGTGTTGAAGAGTTTCGCCAGGTCTGCGATTGCCGCAAACCACATCCGGGGATGCTTTTGTCAGCACGCGATTATTTGCATATTGATA  >  minE/201922‑202039
                                                           |                                                          
aTCTATTATTGCCCGCATCATCCGCAGGGTAGTGTTGAAGAGTTTCGCCAGGTCTGCGATTGCCGCaa                                                    >  1:102928/1‑68 (MQ=255)
aTCTATTATTGCCCGCATCATCCGCAGGGTAGTGTTGAAGAGTTTCGCCAGGTCTGCGACTGCCGCaa                                                    <  1:935874/68‑1 (MQ=255)
aTCTATTATTGCCCGCATCATCCGCAGGGTAGTGTTGAAGAGTTTCGCCAGGTCTGCGACTGCCGCaa                                                    <  1:2332659/68‑1 (MQ=255)
aTCTATTATTGCCCGCATCATCCGCAGGGTAGTGTTGAAGAGTTTCGCCAGGTCTGCGACTGCCGCaa                                                    <  1:2049335/68‑1 (MQ=255)
aTCTATTATTGCCCGCATCATCCGCAGGGTAGTGTTGAAGAGTTTCGCCAGGTCTGCGACTGCCGCaa                                                    <  1:2007999/68‑1 (MQ=255)
aTCTATTATTGCCCGCATCATCCGCAGGGTAGTGTTGAAGAGTTTCGCCAGGTCTGCGACTGCCGCaa                                                    <  1:2927053/68‑1 (MQ=255)
aTCTATTATTGCCCGCATCATCCGCAGGGTAGTGTTGAAGAGTTTCGCCAGGTCTGCGACTGCCGCaa                                                    <  1:2608834/68‑1 (MQ=255)
aTCTATTATTGCCCGCATCATCCGCAGGGTAGTGTTGAAGAGTTTCGCCAGGTCTGCGACTGCCGCaa                                                    <  1:3266584/68‑1 (MQ=255)
aTCTATTATTGCCCGCATCATCCGCAGGGTAGTGTTGAAGAGTTTCGCCAGGTCTGCGACTGCCGCaa                                                    <  1:3369558/68‑1 (MQ=255)
aTCTATTATTGCCCGCATCATCCGCAGGGTAGTGTTGAAGAGTTTCGCCAGGTCTGCGACTGCCGCaa                                                    <  1:3387037/68‑1 (MQ=255)
aTCTATTATTGCCCGCATCATCCGCAGGGTAGTGTTGAAGAGTTTCGCCAGGTCTGCGACTGCCGCaa                                                    <  1:2890839/68‑1 (MQ=255)
aTCTATTATTGCCCGCATCATCCGCAGGGTAGTGTTGAAGAGTTTCGCCAGGTCTGCGACTGCCGCaa                                                    <  1:350485/68‑1 (MQ=255)
aTCTATTATTGCCCGCATCATCCGCAGGGTAGTGTTGAAGAGTTTCGCCAGGTCTGCGACTGCCGCaa                                                    <  1:1452668/68‑1 (MQ=255)
aTCTATTATTGCCCGCATCATCCGCAGGGTAGTGTTGAAGAGTTTCGCCAGGTCTGCGACTGCCGCaa                                                    <  1:2865838/68‑1 (MQ=255)
aTCTATTATTGCCCGCATCATCCGCAGGGTAGTGTTGAAGAGTTTCGCCAGGTCTGCGACTGCCGCaa                                                    <  1:1304447/68‑1 (MQ=255)
aTCTATTATTGCCCGCATCATCCGCAGGGTAGTGTTGAAGAGTTTCGCCAGGTCTGCGACTGCCGCaa                                                    <  1:115058/68‑1 (MQ=255)
aTCTATTATTGCCCGCATCATCCGCAGGGTAGTGTTGAAGAGTTTCGCCAGGTCTGCGACTGCCGCaa                                                    <  1:709033/68‑1 (MQ=255)
 tCTATTATTGCCCGCATCATCCGCAGGGTAGTGTTGAAGAGTTTCGCCAGGTCTGCGACTGCCGCaa                                                    <  1:532115/67‑1 (MQ=255)
  cTATTATTGCCCGCATCATCCGCAGGGTAGTGTTGAAGAGTTTCGCCAGGTCTGCGACTGCCGCaa                                                    <  1:1943900/66‑1 (MQ=255)
   tattatTGCCCGCATCATCCGCAGGGTAGTGTTGAAGAGTTTCGCCAGGTCTGCGATTGCCGCAAAc                                                  <  1:2139720/67‑1 (MQ=255)
   tattatTGCCCGCATCATCCGCAGGGTAGTGTTGAAGAGTTTCGCCAGGTCTGCGATTGCCGCAAAc                                                  <  1:1978667/67‑1 (MQ=255)
   tattatTGCCCGCATCATCCGCAGGGTAGTGTTGAAGAGTTTCGCCAGGTCTGCGATTGCCGCAAAc                                                  <  1:180358/67‑1 (MQ=255)
    attattGCCCGCATCATCCGCAGGGTAGTGTTGAAGAGTTTCGCCAGGTCTGCGATTGCCGCAAAc                                                  <  1:907752/66‑1 (MQ=255)
       attGCCCGCATCATCCGCAGGGTAGTGTTGAAGAGTTTCGCCAGGTCTGCGATTGCCGCAAACCACATc                                            >  1:2726243/1‑69 (MQ=255)
            ccGCATCATCCGCAGGGTAGTGTTGAAGAGTTTCGCCAGGTCTGCGATTGCCGCAAACCACATCCggg                                        <  1:3177008/68‑1 (MQ=255)
            ccGCATCATCCGCAGGGTAGTGTTGAAGAGTTTCGCCAGGTCTGCGATTGCCGCAAACCACATCCggg                                        <  1:2224482/68‑1 (MQ=255)
            ccGCATCATCCGCAGGGTAGTGTTGAAGAGTTTCGCCAGGTCTGCGATTGCCGCAAACCACATCCggg                                        <  1:2994168/68‑1 (MQ=255)
            ccGCATCATCCGCAGGGTAGTGTTGAAGAGTTTCGCCAGGTCTGCGATTGCCGCAAACCACATCCggg                                        <  1:3086696/68‑1 (MQ=255)
            ccGCATCATCCGCAGGGTAGTGTTGAAGAGTTTCGCCAGGTCTGCGATTGCCGCAAACCACATCCggg                                        <  1:1244050/68‑1 (MQ=255)
            ccGCATCATCCGCAGGGTAGTGTTGAAGAGTTTCGCCAGGTCTGCGATTGCCGCAAACCACATCCggg                                        <  1:1598821/68‑1 (MQ=255)
            ccGCATCATCCGCAGGGTAGTGTTGAAGAGTTTCGCCAGGTCTGCGATTGCCGCAAACCACATCCggg                                        <  1:819635/68‑1 (MQ=255)
            ccGCATCATCCGCAGGGTAGTGTTGAAGAGTTTCGCCAGGTCTGCGATTGCCGCAAACCACATCCggg                                        <  1:3368553/68‑1 (MQ=255)
            ccGCATCATCCGCAGGGTAGTGTTGAAGAGTTTCGCCAGGTCTGCGATTGCCGCAAACCACATCCggg                                        <  1:1289053/68‑1 (MQ=255)
              gCATCATCCGCAGGGTAGTGTTGAAGAGTTTCGCCAGGTCTGCGATTGCCGCAAACcaca                                              >  1:123003/1‑60 (MQ=255)
              gCATCATCCGCAGGGTAGTGTTGAAGAGTTTCGCCAGGTCTGCGATTGCCGCAAACCACATCCGGGGa                                      >  1:1156679/1‑68 (MQ=255)
              gCATCATCCGCAGGGTAGTGTTGAAGAGTTTCGCCAGGTCTGCGATTGCCGCAAACCACATCCGGGGAt                                     >  1:503894/1‑69 (MQ=255)
              gCATCATCCGCAGGGTAGTGTTGAAGAGTTTCGCCAGGTCTGCGATTGCCGCAAACCACATCCGGGGAt                                     >  1:2557200/1‑69 (MQ=255)
              gCATCATCCGCAGGGTAGTGTTGAAGAGTTTCGCCAGGTCTGCGATTGCCGCAAACCACATCCGGGGAt                                     >  1:2130822/1‑69 (MQ=255)
              gCATCATCCGCAGGGTAGTGTTGAAGAGTTTCGCCAGGTCTGCGATTGCCGCAAACCACATCCGGGGAt                                     >  1:2911331/1‑69 (MQ=255)
              gCATCATCCGCAGGGTAGTGTTGAAGAGTTTCGCCAGGTCTGCGATTGCCGCAAACCACATCCGGGGAt                                     >  1:3109255/1‑69 (MQ=255)
              gCATCATCCGCAGGGTAGTGTTGAAGAGTTTCGCCAGGTCTGCGATTGCCGCAAACCACATCCGGGGAt                                     >  1:1578985/1‑69 (MQ=255)
              gCATCATCCGCAGGGTAGTGTTGAAGAGTTTCGCCAGGTCTGCGATTGCCGCAAACCACATCCGGGGAt                                     >  1:503021/1‑69 (MQ=255)
                  catcCGCAGGGTAGTGTTGAAGAGTTTCGCCAGGTCTGCGATTGCCGCAAACCACATCCGGGGATGcc                                  <  1:2786430/68‑2 (MQ=255)
                       gCAGGGTAGTGTTGAAGAGTTTCGCCAGGGCTGCGATTGCCGCAAACCACATCCggg                                        <  1:586981/57‑1 (MQ=255)
                                    gAAGAGTTTCGCCAGGTCTGCGATTGCCGCAAACCACATCCGGGGATGCTTTTGTCAGCACGCGatta                <  1:345982/68‑1 (MQ=255)
                                    gAAGAGTTTCGCCAGGTCTGCGATTGCCGCAAACCACATCCGGGGATGCTTTTGTCAGCACGCGatta                <  1:456347/68‑1 (MQ=255)
                                       gagTTTCGCCAGGTCTGCGATTGCCGCAAACCACATCCGGGGATGCTTTTGTCAGCACGCGATTATTTg            <  1:3216456/69‑1 (MQ=255)
                                       gagTTTCGCCAGGTCTGCGATTGCCGCAAACCACATCCGGGGATGCTTTTGTCAGCACGCGATTATTTg            <  1:682470/69‑1 (MQ=255)
                                       gagTTTCGCCAGGTCTGCGATTGCCGCAAACCACATCCGGGGATGCTTTTGTCAGCACGCGATTATTTg            <  1:715156/69‑1 (MQ=255)
                                       gagTTTCGCCAGGTCTGCGATTGCCGCAAACCACATCCGGGGATGCTTTTGTCAGCACGCGATTATTTg            <  1:1464591/69‑1 (MQ=255)
                                       gagTTTCGCCAGGTCTGCGATTGCCGCAAACCACATCCGGGGATGCTTTTGTCAGCACGCGATTATTTg            <  1:1448423/69‑1 (MQ=255)
                                       gagTTTCGCCAGGTCTGCGATTGCCGCAAACCACATCCGGGGATGCTTTTGTCAGCACGCGATTATTTg            <  1:1134354/69‑1 (MQ=255)
                                       gagTTTCGCCAGGTCTGCGATTGCCGCAAACCACATCCGGGGATGCTTTTGTCAGCACGCGATTATTTg            <  1:1438074/69‑1 (MQ=255)
                                       gagTTTCGCCAGGTCTGCGATTGCCGCAAACCACATCCGGGGATGCTTTTGTCAGCACGCGATTATTTg            <  1:1389552/69‑1 (MQ=255)
                                       gagTTTCGCCAGGTCTGCGATTGCCGCAAACCACATCCGGGGATGCTTTTGTCAGCACGCGATTATTTg            <  1:1230375/69‑1 (MQ=255)
                                       gagTTTCGCCAGGTCTGCGATTGCCGCAAACCACATCCGGGGATGCTTTTGTCAGCACGCGATTATTTg            <  1:1305923/69‑1 (MQ=255)
                                       gagTTTCGCCAGGTCTGCGATTGCCGCAAACCACATCCGGGGATGCTTTTGTCAGCACGCGATTATTTg            <  1:1724501/69‑1 (MQ=255)
                                       gagTTTCGCCAGGTCTGCGATTGCCGCAAACCACATCCGGGGATGCTTTTGTCAGCACGCGATTATTTg            <  1:1601030/69‑1 (MQ=255)
                                       gagTTTCGCCAGGTCTGCGATTGCCGCAAACCACATCCGGGGATGCTTTTGTCAGCACGCGATTATTTg            <  1:1648546/69‑1 (MQ=255)
                                         gTTTCGCCAGGTCTGCGATTGCCGCAAACCACATCCGGGGATGCTTTTGTCAGCACGCGATTATTTGCa          >  1:2951052/1‑69 (MQ=255)
                                                cAGGTCTGCGATTGCCGCAAACCACATCCGGGGATGCTTTTGTCAGCACGCGATTATTTGCATATTGa    <  1:1679143/68‑1 (MQ=255)
                                                cAGGTCTGCGATTGCCGCAAACCACATCCGGGGATGCTTTTGTCAGCACGCGATTATTTGCATATTGa    <  1:1549667/68‑1 (MQ=255)
                                                cAGGTCTGCGATTGCCGCAAACCACATCCGGGGATGCTTTTGTCAGCACGCGATTATTTGCATATTGa    <  1:2834449/68‑1 (MQ=255)
                                                cAGGTCTGCGATTGCCGCAAACCACATCCGGGGATGCTTTTGTCAGCACGCGATTATTTGCATATTGa    <  1:3218961/68‑1 (MQ=255)
                                                cAGGTCTGCGATTGCCGCAAACCACATCCGGGGATGCTTTTGTCAGCACGCGATTATTTGCATATTGa    <  1:319362/68‑1 (MQ=255)
                                                cAGGTCTGCGATTGCCGCAAACCACATCCGGGGATGCTTTTGTCAGCACGCGATTATTTGCATATTGa    <  1:1831230/68‑1 (MQ=255)
                                                cAGGTCTGCGATTGCCGCAAACCACATCCGGGGATGCTTTTGTCAGCACGCGATTATTTGCATATTGa    <  1:2102311/68‑1 (MQ=255)
                                                cAGGTCTGCGATTGCCGCAAACCACATCCGGGGATGCTTTTGTCAGCACGCGATTATTTGCATATTGa    <  1:263202/68‑1 (MQ=255)
                                                cAGGTCTGCGATTGCCGCAAACCACATCCGGGGATGCTTTTGTCAGCACGCGATTATTTGCATATTGa    <  1:288143/68‑1 (MQ=255)
                                                cAGGTCTGCGATTGCCGCAAACCACATCCGGGGATGCTTTTGTCAGCACGCGATTATTTGCATATTGa    <  1:2694369/68‑1 (MQ=255)
                                                 aGGTCTGCGATTGCCGCAAACCACATCCGGGGATGCTTTTGTCAGCACGCGATTATTTGCATATTGa    <  1:3235941/67‑1 (MQ=255)
                                                 aGGTCTGCGACTGCCGCAAACCACATCCGGGGATGTTTTTGTCAGCACGCGATTATTTGCATATTGata  >  1:1473179/1‑69 (MQ=255)
                                                 aGGTCTGCGACTGCCGCAAACCACATCCGGGGATGTTTTTGTCAGCACGCGATTATTTGCATATTGata  >  1:1533733/1‑69 (MQ=255)
                                                 aGGTCTGCGACTGCCGCAAACCACATCCGGGGATGTTTTTGTCAGCACGCGATTATTTGCATATTGata  >  1:1730506/1‑69 (MQ=255)
                                                 aGGTCTGCGACTGCCGCAAACCACATCCGGGGATGTTTTTGTCAGCACGCGATTATTTGCATATTGata  >  1:1782392/1‑69 (MQ=255)
                                                 aGGTCTGCGACTGCCGCAAACCACATCCGGGGATGTTTTTGTCAGCACGCGATTATTTGCATATTGata  >  1:2419149/1‑69 (MQ=255)
                                                 aGGTCTGCGACTGCCGCAAACCACATCCGGGGATGTTTTTGTCAGCACGCGATTATTTGCATATTGat   >  1:1577668/1‑68 (MQ=255)
                                                 aGGTCTGCGACTGCCGCAAACCACATCCGGGGATGTTTTTGTCAGCACGCGATTATTTGCATATTGat   >  1:1341012/1‑68 (MQ=255)
                                                     cTGCGATTGCCGCAAACCACATCCGGGGATGCTTTTGTCAGCACGCGATTATTTGCATAtt      <  1:716454/61‑1 (MQ=255)
                                                     cTGCGATTGCCGCAAACCACATCCGGGGATGCTTTTGTCAGCACGCGATTATTTGCATAtt      <  1:112403/61‑1 (MQ=255)
                                                     cTGCGATTGCCGCAAACCACATCCGGGGATGCTTTTGTCAGCACGCGATTATTTGCATAtt      <  1:3300320/61‑1 (MQ=255)
                                                     cTGCGATTGCCGCAAACCACATCCGGGGATGCTTTTGTCAGCACGCGATTATTTGCATAtt      <  1:2907322/61‑1 (MQ=255)
                                                     cTGCGATTGCCGCAAACCACATCCGGGGATGCTTTTGTCAGCACGCGATTATTTGCATAtt      <  1:2706735/61‑1 (MQ=255)
                                                      tGCGATTGCCGCAAACCACATCCGGGGATGCTTTTGTCAGCACGCGATTATTTGCATAtt      <  1:2594017/60‑1 (MQ=255)
                                                          aTTGCCGCAAACCACATCCGGGGATGCTTTTGTCAGCACGCGATTATTTg            <  1:738222/50‑1 (MQ=255)
                                                           |                                                          
ATCTATTATTGCCCGCATCATCCGCAGGGTAGTGTTGAAGAGTTTCGCCAGGTCTGCGATTGCCGCAAACCACATCCGGGGATGCTTTTGTCAGCACGCGATTATTTGCATATTGATA  >  minE/201922‑202039

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: