Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 614,741 G→A 51.6% A130A (GCG→GCA serC → 3‑phosphoserine/phosphohydroxythreonine aminotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE614,7410GA51.6% 1.2 / 141.2 152A130A (GCG→GCAserC3‑phosphoserine/phosphohydroxythreonine aminotransferase
Reads supporting (aligned to +/- strand):  ref base G (54/19);  new base A (19/59);  total (74/78)
Fisher's exact test for biased strand distribution p-value = 8.90e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 169 total aligned reads displayed.

AAAGAAGCGAAAAAATACTGCACGCCTAATGTCTTTGACGCCAAAGTGACTGTTGATGGTCTGCGCGCGGTTAAGCCAATGCGTGAATGGCAACTCTCTGATAATGCTGCTTATATGCATTATTGCCCGAA  >  minE/614673‑614803
                                                                    |                                                              
aaaGAAGCGAAAAAATACTGCACGCCTAATGTCTTTGACGCCAAAGTGACTGTTGATGGTCTGcgcgcg                                                                >  1:3137586/1‑69 (MQ=255)
                aCTGCACGCCTAATGTCTTTGACGCCAAAGTGACTGTTGATGGTCTGCGCGCGGTTAAGCCAATGCGTg                                                <  1:73318/69‑1 (MQ=255)
                            aTGTCTTTGACGCCAAAGTGACTGTTGATGGTCTGCGCGCGGTTAAGCCAATGCGTGAATGGCAActct                                    >  1:153209/1‑69 (MQ=255)
                            aTGTCTTTGAAGCCAAAGTGACTGTTGATGGTCTGCGCGCGGTTAAGCCAATGCGTGaa                                              >  1:3029024/1‑59 (MQ=255)
                               tCTTTGACGCCAAAGTGACTGTTGATGGTCTGCGCGCGGTTAAGCCAATGCGTg                                                >  1:544617/1‑54 (MQ=255)
                               tCTTTGACGCCAAAGTGACTGTTGATGGTCTGCGCGCGGTTAAGCCAATGCGTg                                                >  1:579829/1‑54 (MQ=255)
                               tCTTTGACGCCAAAGTGACTGTTGATGGTCTGCGCGCGGTTAAGCCAATGCGTg                                                >  1:833183/1‑54 (MQ=255)
                               tCTTTGACGCCAAAGTGACTGTTGATGGTCTGCGCGCGGTTAAGCCAATGCGTg                                                >  1:2986355/1‑54 (MQ=255)
                               tCTTTGACGCCAAAGTGACTGTTGATGGTCTGCGCGCGGTTAAGCCAATGCGTg                                                >  1:1210119/1‑54 (MQ=255)
                               tCTTTGACGCCAAAGTGACTGTTGATGGTCTGCGCGCGGTTAAGCCAATGCGTg                                                >  1:3426773/1‑54 (MQ=255)
                               tCTTTGACGCCAAAGTGACTGTTGATGGTCTGCGCGCGGTTAAGCCAATGCGTGAATGGCAActctct                                  >  1:3377772/1‑68 (MQ=255)
                               tCTTTGACGCCAAAGTGACTGTTGATGGTCTGCGCGCGGTTAAGCCAATGCGTGAATGGCAActctct                                  >  1:544214/1‑68 (MQ=255)
                               tCTTTGACGCCAAAGTGACTGTTGATGGTCTGCGCGCGGTTAAGCCAATGCGTGAATGGCAActctct                                  >  1:2133783/1‑68 (MQ=255)
                               tCTTTGACGCCAAAGTGACTGTTGATGGTCTGCGCGCGGTTAAGCCAATGCGTGAATGGCAActctct                                  >  1:529544/1‑68 (MQ=255)
                               tCTTTGACGCCAAAGTGACTGTTGATGGTCTGCGCGCGGTTAAGCCAATGCGTGAATGGCAActctct                                  >  1:2328773/1‑68 (MQ=255)
                               tCTTTGACGCCAAAGTGACTGTTGATGGTCTGCGCGCGGTTAAGCCAATGCGTGAATGGCAActctct                                  >  1:1432910/1‑68 (MQ=255)
                               tCTTTGACGCCAAAGTGACTGTTGATGGTCTGCGCGCGGTTAAGCCAATGCGTGAATGGCAActctct                                  >  1:727928/1‑68 (MQ=255)
                               tCTTTGACGCCAAAGTGACTGTTGATGGTCTGCGCGCGGTTAAGCCAATGCGTGAATGGCAActctct                                  >  1:1184896/1‑68 (MQ=255)
                               tCTTTGACGCCAAAGTGACTGTTGATGGTCTGCGCGCGGTTAAGCCAATGCGTGAATGGCAActctct                                  >  1:1075391/1‑68 (MQ=255)
                               tCTTTGACGCCAAAGTGACTGTTGATGGTCTGCGCGCGGTTAAGCCAATGCGTGAATGGCAActctct                                  >  1:921564/1‑68 (MQ=255)
                               tCTTTGACGCCAAAGTGACTGTTGATGGTCTGGCGCGGTTAAGCCAATGCGTGAATGGCAACTCTCTg                                 >  1:1875822/1‑68 (MQ=255)
                                              tGACCGTTGATGGTCTTCGCGCAGTTAAGCCAATGAGTGAATGGCAACTCTCTGATAATGCTGCTtata                  <  1:937678/69‑1 (MQ=255)
                                              tGACCGTTGATGGTCTTCGCGCAGTTAAGCCAATGAGTGAATGGCAACTCTCTGATAATGCTGCTtata                  <  1:2799155/69‑1 (MQ=255)
                                              tGACCGTTGATGGTCTTCGCGCAGTTAAGCCAATGAGTGAATGGCAACTCTCTGATAATGCTGCTtata                  <  1:2558120/69‑1 (MQ=255)
                                              tGACCGTTGATGGTCTTCGCGCAGTTAAGCCAATGAGTGAATGGCAACTCTCTGATAATGCTGCTtata                  <  1:2819795/69‑1 (MQ=255)
                                              tGACCGTTGATGGTCTTCGCGCAGTTAAGCCAATGAGTGAATGGCAACTCTCTGATAATGCTGCTtata                  <  1:461832/69‑1 (MQ=255)
                                              tGACCGTTGATGGTCTTCGCGCAGTTAAGCCAATGAGTGAATGGCAACTCTCTGATAATGCTGCTtata                  <  1:3121319/69‑1 (MQ=255)
                                              tGACCGTTGATGGTCTTCGCGCAGTTAAGCCAATGAGTGAATGGCAACTCTCTGATAATGCTGCTtata                  <  1:2271149/69‑1 (MQ=255)
                                              tGACCGTTGATGGTCTTCGCGCAGTTAAGCCAATGAGTGAATGGCAACTCTCTGATAATGCTGCTtata                  <  1:3050973/69‑1 (MQ=255)
                                              tGACCGTTGATGGTCTTCGCGCAGTTAAGCCAATGAGTGAATGGCAACTCTCTGATAATGCTGCTtata                  <  1:2229021/69‑1 (MQ=255)
                                              tGACCGTTGATGGTCTTCGCGCAGTTAAGCCAATGAGTGAATGGCAACTCTCTGATAATGCTGCTtata                  <  1:2183226/69‑1 (MQ=255)
                                              tGACCGTTGATGGTCTTCGCGCAGTTAAGCCAATGAGTGAATGGCAACTCTCTGATAATGCTGCTtata                  <  1:3121021/69‑1 (MQ=255)
                                              tGACCGTTGATGGTCTTCGCGCAGTTAAGCCAATGAGTGAATGGCAACTCTCTGATAATGCTGCTtata                  <  1:1132404/69‑1 (MQ=255)
                                              tGACCGTTGATGGTCTTCGCGCAGTTAAGCCAATGAGTGAATGGCAACTCTCTGATAATGCTGCTtata                  <  1:716903/69‑1 (MQ=255)
                                              tGACCGTTGATGGTCTTCGCGCAGTTAAGCCAATGAGTGAATGGCAACTCTCTGATAATGCTGCTtata                  <  1:1781588/69‑1 (MQ=255)
                                              tGACCGTTGATGGTCTTCGCGCAGTTAAGCCAATGAGTGAATGGCAACTCTCTGATAATGCTGCTtata                  <  1:654922/69‑1 (MQ=255)
                                              tGACCGTTGATGGTCTTCGCGCAGTTAAGCCAATGAGTGAATGGAACTCTCTGATAATGCTGCTTatat                 <  1:1745184/69‑1 (MQ=255)
                                               gACTGGTGATGGTCTGCGCGCGGTTAAGCCAATGCGTGAATGGCAACTCTCTGATAATGCTGCTtat                   >  1:3210558/1‑67 (MQ=255)
                                               gACCGTTGATGGTCTTCGCGCAGTTAAGCCAATGAGTGAATGGCAACTCTCTGATAATGCTGCTTatat                 <  1:1655569/69‑1 (MQ=255)
                                               gACCGTTGATGGTCTTCGCGCAGTTAAGCCAATGAGTGAATGGCAACTCTCTGATAATGCTGCTTatat                 <  1:3295648/69‑1 (MQ=255)
                                               gACCGTTGATGGTCTTCGCGCAGTTAAGCCAATGAGTGAATGGCAACTCTCTGATAATGCTGCTTatat                 <  1:3360740/69‑1 (MQ=255)
                                               gACCGTTGATGGTCTTCGCGCAGTTAAGCCAATGAGTGAATGGCAACTCTCTGATAATGCTGCTTatat                 <  1:2792464/69‑1 (MQ=255)
                                               gACCGTTGATGGTCTTCGCGCAGTTAAGCCAATGAGTGAATGGCAACTCTCTGATAATGCTGCTTatat                 <  1:1404641/69‑1 (MQ=255)
                                               gACCGTTGATGGTCTTCGCGCAGTTAAGCCAATGAGTGAATGGCAACTCTCTGATAATGCTGCTTatat                 <  1:2470835/69‑1 (MQ=255)
                                               gACCGTTGATGGTCTTCGCGCAGTTAAGCCAATGAGTGAATGGCAACTCTCTGATAATGCTGCTTatat                 <  1:639110/69‑1 (MQ=255)
                                               gACCGTTGATGGTCTTCGCGCAGTTAAGCCAATGAGTGAATGGCAACTCTCTGATAATGCTGCTTatat                 <  1:1366922/69‑1 (MQ=255)
                                               gACCGTTGATGGTCTTCGCGCAGTTAAGCCAATGAGTGAATGGCAACTCTCTGATAATGCTGCTTatat                 <  1:676973/69‑1 (MQ=255)
                                               gACCGTTGATGGTCTTCGCGCAGTTAAGCCAATGAGTGAATGGCAACTCTCTGATAATGCTGCTTatat                 <  1:843740/69‑1 (MQ=255)
                                               gACCGTTGATGGTCTTCGCGCAGTTAAGCCAATGAGTGAATGGCAACTCTCTGATAATGCTGCTTatat                 <  1:3068988/69‑1 (MQ=255)
                                                accgTTGATGGTCTTCGCGCAGTTAAGCCAATGAGTGAATGGCAACTCTCTGATAATGCTGCTTatat                 <  1:2421793/65‑1 (MQ=255)
                                                  tGTTGATGGTCTGCGCGCGGTTAAGCCAATGCGTGAATGGc                                          >  1:2066083/1‑41 (MQ=255)
                                                  tGTTGATGGTCTGCGCGCGGTTAAGCCAATGCGTGAATGGc                                          >  1:77382/1‑41 (MQ=255)
                                                  tGTTGATGGTCTGCGCGCGGTTAAGCCAATGCGTGAATGGc                                          >  1:114361/1‑41 (MQ=255)
                                                  tGTTGATGGTCTGCGCGCGGTTAAGCCAATGCGTGAATGGc                                          >  1:1046183/1‑41 (MQ=255)
                                                  tGTTGATGGTCTGCGCGCGGTTAAGCCAATGCGTGAATGGc                                          >  1:2621193/1‑41 (MQ=255)
                                                  tGTTGATGGTCTGCGCGCGGTTAAGCCAATGCGTGAATGGc                                          >  1:387446/1‑41 (MQ=255)
                                                  tGTTGATGGTCTGCGCGCGGTTAAGCCAATGCGTGAATGGc                                          >  1:2347028/1‑41 (MQ=255)
                                                  tGTTGATGGTCTGCGCGCGGTTAAGCCAATGCGTGAATGGc                                          >  1:568897/1‑41 (MQ=255)
                                                  tGTTGATGGTCTGCGCGCGGTTAAGCCAATGCGTGAATGGc                                          >  1:2209784/1‑41 (MQ=255)
                                                  tGTTGATGGTCTGCGCGCGGTTAAGCCAATGCGTGAATGGc                                          >  1:1921325/1‑41 (MQ=255)
                                                  tGTTGATGGTCTGCGCGCGGTTAAGCCAATGCGTGAATGGc                                          >  1:2041426/1‑41 (MQ=255)
                                                    ttGATGGTCTTCGCGCTGTTAAGCCAATGAGTGAATGGCAACTCTCTGATAATGCTGCTTATATGCatt            >  1:1889091/1‑69 (MQ=255)
                                                    ttGATGGTCTTCGCGCAGTTAAGCCAATGAGTGAATGGCAACTCTCTGATAATGCTGCTTATATGCatt            >  1:633776/1‑69 (MQ=255)
                                                    ttGATGGTCTTCGCGCAGTTAAGCCAATGAGTGAATGGCAACTCTCTGATAATGCTGCTTATATGCatt            >  1:2652157/1‑69 (MQ=255)
                                                    ttGATGGTCTTCGCGCAGTTAAGCCAATGAGTGAATGGCAACTCTCTGATAATGCTGCTTATATGCatt            >  1:2973641/1‑69 (MQ=255)
                                                    ttGATGGTCTTCGCGCAGTTAAGCCAATGAGTGAATGGCAACTCTCTGATAATGCTGCTTATATGCatt            >  1:1288326/1‑69 (MQ=255)
                                                    ttGATGGTCTTCGCGCAGTTAAGCCAATGAGTGAATGGCAACTCTCTGATAATGCTGCTTATATGCatt            >  1:2737569/1‑69 (MQ=255)
                                                    ttGATGGTCTTCGCGCAGTTAAGCCAATGAGTGAATGGCAACTCTCTGATAATGCTGCTTATATGCatt            >  1:1449468/1‑69 (MQ=255)
                                                    ttGATGGTCTTCGCGCAGTTAAGCCAATGAGTGAATGGCAACTCTCTGATAATGCTGCTTATATGCatt            >  1:2709024/1‑69 (MQ=255)
                                                    ttGATGGTCTTCGCGCAGTTAAGCCAATGAGTGAATGGCAACTCTCTGATAATGCTGCTTATATGCatt            >  1:270460/1‑69 (MQ=255)
                                                    ttGATGGTCTTCGCGCAGTTAAGCCAATGAGTGAATGGCAACTCTCTGATAATGCTGCTTATATGCatt            >  1:1828297/1‑69 (MQ=255)
                                                    ttGATGGTCTTCGCGCAGTTAAGCCAATGAGTGAATGGCAACTCTCTGATAATGCTGCTTATATGCatt            >  1:1970916/1‑69 (MQ=255)
                                                    ttGATGATCTTCGCGCAGTTAAGCCAATGAGTGAATGGCAACTCTCTGATAATGCTGCTTATATGCatt            >  1:2394/1‑69 (MQ=255)
                                                     tGATGGTCTGCGCGCGGTTAAGCCAATGCGTGAATGGCAACTCTCTGATaa                             <  1:104678/51‑1 (MQ=255)
                                                     tGATGGTCTGCGCGCGGTTAAGCCAATGCGTGAATGGCAACTCTCTGATaa                             <  1:1227169/51‑1 (MQ=255)
                                                     tGATGGTCTGCGCGCGGTTAAGCCAATGCGTGAATGGCAACTCTCTGATaa                             <  1:2458454/51‑1 (MQ=255)
                                                     tGATGGTCTGCGCGCGGTTAAGCCAATGCGTGAATGGCAACTCTCTGATaa                             <  1:2872521/51‑1 (MQ=255)
                                                     tGATGGTCTGCGCGCGGTTAAGCCAATGCGTGAATGGCAACTCTCTGATaa                             <  1:2597339/51‑1 (MQ=255)
                                                     tGATGGTCTGCGCGCGGTTAAGCCAATGCGTGAATGGCAACTCTCTGATaa                             <  1:423397/51‑1 (MQ=255)
                                                      gATGGTCTTCGCGCAGTTAAGCCAATGAGTGAATGGCAACTCTCTGATAATGCTGCtt                     <  1:3248767/58‑1 (MQ=255)
                                                      gATGGTCTTCGCGCAGTTAAGCCAATGAGTGAATGGCAACTCTCTGATAATGCTGCtt                     <  1:308188/58‑1 (MQ=255)
                                                      gATGGTCTTCGCGCAGTTAAGCCAATGAGTGAATGGCAACTCTCTGATAATGCTGCtt                     <  1:1055090/58‑1 (MQ=255)
                                                      gATGGTCTTCGCGCAGTTAAGCCAATGAGTGAATGGCAACTCTCTGATAATGCTGCtt                     <  1:883954/58‑1 (MQ=255)
                                                      gATGGTCTTCGCGCAGTTAAGCCAATGAGTGAATGGCAACTCTCTGATAATGCTGCtt                     <  1:788269/58‑1 (MQ=255)
                                                      gATGGTCTTCGCGCAGTTAAGCCAATGAGTGAATGGCAACTCTCTGATAATGCTGCtt                     <  1:1259279/58‑1 (MQ=255)
                                                      gATGGTCTTCGCGCAGTTAAGCCAATGAGTGAATGGCAACTCTCTGATAATGCTGCtt                     <  1:319447/58‑1 (MQ=255)
                                                      gATGGTCTTCGCGCAGTTAAGCCAATGAGTGAATGGCAACTCTCTGATAATGCTGCtt                     <  1:1525529/58‑1 (MQ=255)
                                                      gATGGTCTTCGCGCAGTTAAGCCAATGAGTGAATGGCAACTCTCTGATAATGCTGCtt                     <  1:504958/58‑1 (MQ=255)
                                                      gATGGTCTTCGCGCAGTTAAGCCAATGAGTGAATGGCAACTCTCTGATAATGCTGCtt                     <  1:489918/58‑1 (MQ=255)
                                                      gATGGTCTTCGCGCAGTTAAGCCAATGAGTGAATGGCAACTCTCTGATAATGCTGCtt                     <  1:2359473/58‑1 (MQ=255)
                                                       aTGGTCTTCGCGCAGTTAAGCCAATGAGTGAATGGCAACTCTCTGATAATGCTGCtt                     <  1:1938426/57‑1 (MQ=255)
                                                            cTGCGCGCGGTTAAGCCAATGCGTGAATGGCAACTCTCTGATAATGCTGCTTATATGCATTATTGCCCg    <  1:2273394/69‑1 (MQ=255)
                                                            cTGCGCGCGGTTAAGCCAATGCGTGAATGGCAACTCTCTGATAATGCTGCTTATATGCATTATTGCCCg    <  1:87773/69‑1 (MQ=255)
                                                             tGCGCGCGGTTAAGCCAATGCGTGAATGGCAACTCTCTGATAATGCTGCTTATATGCATTATTGCCCg    <  1:1340412/68‑1 (MQ=255)
                                                               cgcgcgGTTAAGCCAATGCGTTAATGGCAACTCTCTGATaa                             <  1:1261149/41‑1 (MQ=255)
                                                               cgcgcgGTTAAGCCAATGCGTGAATGGCAACTCTCTGATAATGCTGCTTATATGCATTATTGCCCGaa  >  1:589072/1‑68 (MQ=255)
                                                               cgcgcgGTTAAGCCAATGCGTGAATGGCAACTCTCTGATAATGCTGCTTATATGCATTATTGCCCGaa  >  1:1830128/1‑68 (MQ=255)
                                                               cgcgcgGTTAAGCCAATGCGTGAATGGCAACTCTCTGATAATGCTGCTTATATGCATTATTGCCCGaa  >  1:1138217/1‑68 (MQ=255)
                                                               cgcgcgGTTAAGCCAATGCGTGAATGGCAACTCTCTGATAATGCTGCTTATATGCATTATTGCCCGaa  >  1:1100968/1‑68 (MQ=255)
                                                                 cgcaGTTAAGCCAATGAGTGAATGGCAACTCTCTGATAATGCTGCTTatat                 <  1:1815409/51‑1 (MQ=255)
                                                                    |                                                              
AAAGAAGCGAAAAAATACTGCACGCCTAATGTCTTTGACGCCAAAGTGACTGTTGATGGTCTGCGCGCGGTTAAGCCAATGCGTGAATGGCAACTCTCTGATAATGCTGCTTATATGCATTATTGCCCGAA  >  minE/614673‑614803

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: