Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,178,850 G→A 47.8% A100A (GCG→GCA ansA → cytoplasmic L‑asparaginase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,178,8500GA47.8% 7.2 / 82.1 94A100A (GCG→GCAansAcytoplasmic L‑asparaginase I
Reads supporting (aligned to +/- strand):  ref base G (12/37);  new base A (19/26);  total (31/63)
Fisher's exact test for biased strand distribution p-value = 8.15e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTG  >  minE/1178786‑1178917
                                                                |                                                                   
aTGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCt                                                                   >  1:215505/1‑67 (MQ=255)
aTGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGt                                                                 >  1:2486516/1‑69 (MQ=255)
            atgGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGt                                                              <  1:1639198/60‑1 (MQ=255)
            atgGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGt                                                              <  1:2257998/60‑1 (MQ=255)
            atgGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGt                                                              <  1:2036724/60‑1 (MQ=255)
            atgGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGt                                                              <  1:2634533/60‑1 (MQ=255)
            atgGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGt                                                              <  1:2973523/60‑1 (MQ=255)
            atgGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGt                                                              <  1:557825/60‑1 (MQ=255)
            atgGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGt                                                              <  1:1435850/60‑1 (MQ=255)
                 tttGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGaa                                                 <  1:2758563/68‑1 (MQ=255)
                 tttGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGaa                                                 <  1:993890/68‑1 (MQ=255)
                 tttGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGaa                                                 <  1:850055/68‑1 (MQ=255)
                 tttGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGaa                                                 <  1:2028356/68‑1 (MQ=255)
                 tttGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGaa                                                 <  1:230257/68‑1 (MQ=255)
                 tttGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGaa                                                 <  1:1176718/68‑1 (MQ=255)
                 tttGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGaa                                                 <  1:637998/68‑1 (MQ=255)
                 tttGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGaa                                                 <  1:3395173/68‑1 (MQ=255)
                  ttGTCATTCTGCACGGTACAGACACGATGGCGTATACCGCCTCTGCACTGTCGTTCATGCTCGAGAAtc                                               >  1:495347/1‑69 (MQ=255)
                  ttGTCATTCTGCACGGTACAGACACGATGGCGTATACCGCCTCTGCACTGTCGTTCATGCTCGAGAAtc                                               >  1:1498042/1‑69 (MQ=255)
                  ttGTCATTCTGCACGGTACAGACACGATGGCGTATACCGCCTCTGCACTGTCGTTCATGCTCGAGAAtc                                               >  1:804464/1‑69 (MQ=255)
                  ttGTCATTCTGCACGGTACAGACACGATGGCGTATACCGCCTCTGCACTGTCGTTCATGCTCGAGAAtc                                               >  1:1963260/1‑69 (MQ=255)
                  ttGTCATTCTGCACGGTACAGACACGATGGCGTATACCGCCTCTGCACTGTCGTTCATGCTCGAGAAtc                                               >  1:1948663/1‑69 (MQ=255)
                  ttGTCATTCTGCACGGTACAGACACGATGGCGTATACCGCCTCTGCACTGTCGTTCATGCTCGAGAAtc                                               >  1:819074/1‑69 (MQ=255)
                  ttGTCATTCTGCACGGTACAGACACGATGGCGTATACCGCCTCTGCACTGTCGTTCATGCTCGAGAAtc                                               >  1:243729/1‑69 (MQ=255)
                  ttGTCATTCTGCACGGTACAGACACGATGGCGTATACCGCCTCTGCACTGTCGTTCATGCTCGAGAAtc                                               >  1:2483896/1‑69 (MQ=255)
                  ttGTCATTCTGCACGGTACAGACACGATGGCGTATACCGCCTCTGCACTGTCGTTCATGCTCGAGAAtc                                               >  1:2559906/1‑69 (MQ=255)
                  ttGTCATTCTGCACGGTACAGACACGATGGCGTATACCGCCTCTGCACTGTCGTTCATGCTCGAGAAtc                                               >  1:120273/1‑69 (MQ=255)
                  ttGTCATTCTGCACGGTACAGACACGATGGCGTATACCGCCTCTGCACTGTCGTTCATGCTCGAGAAtc                                               >  1:1619155/1‑69 (MQ=255)
                  ttGTCATTCTGCACGGTACAGACACGATGGCGTATACCGCCTCTGCACTGTCGTTCATGCTCGAGAAtc                                               >  1:838764/1‑69 (MQ=255)
                  ttGTCATTCTGCACGGTACAGACACGATGGCGTATACCGCCTCTGCACTGTCGTTCATGCTCGAGAAtc                                               >  1:3141321/1‑69 (MQ=255)
                  ttGTCATTCTGCACGGTACAGACACGATGGCGTATACCGCCTCTGCACTGTCGTTCATGCTCGAGAAtc                                               >  1:1353197/1‑69 (MQ=255)
                  ttGTCATTCTGCACGGTACAGACACGATGGCGTATACCGCCTCTGCACTGTCGTTCATGCTCGAGAAt                                                >  1:2197769/1‑68 (MQ=255)
                  ttGTCATTCTGCACGGTACAGACACGATGGCGTATACCGCCTCTGAACTGTCGTTCATGCTCGAGAAt                                                >  1:3285111/1‑68 (MQ=255)
                  ttGTCATTCTGCACGGTACAGACACGATGGCGTATACCGCATCTGCACTGTCGTTCATGCTCGAGAAtc                                               >  1:2070818/1‑69 (MQ=255)
                  ttGTCATTATGCACGGTACAGACACGATGGCGTATACCGCCTCTGCACTGTCGTTCATGCTCGAGaa                                                 >  1:2565139/1‑67 (MQ=255)
                  atGTCATTCTGCACGGTACAGACACGATGGCGTATACCGCCTCTGCACTGTCGTTCATGc                                                        >  1:2305838/2‑60 (MQ=255)
                      cATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGt                                                              <  1:173979/50‑1 (MQ=255)
                               cGGTACAGACACGATGGCGTATACCGCCTCTGCACTGTCGTTCATGCTCGAGAATCTCGGTAAACCg                                    <  1:514743/67‑1 (MQ=255)
                               cGGTACAGACACGATGGCGTATACCGCCTCTGCACTGTCGTTCATGCTCGAGAATCTCGGTAAACCg                                    <  1:3050909/67‑1 (MQ=255)
                               cGGTACAGACACGATGGCGTATACCGCCTCTGCACTGTCGTTCATGCTCGAGAATCTCGGTAAACCg                                    <  1:1727794/67‑1 (MQ=255)
                               cGGTACAGACACGATGGCGTATACCGCCTCTGCACTGTCGTTCATGCTCGAGAATCTCGGTAAACCg                                    <  1:3362331/67‑1 (MQ=255)
                               cGGTACAGACACGATGGCGTATACCGCCTCTGCACTGTCGTTCATGCTCGAGAATCTCGGTAAACCg                                    <  1:3379424/67‑1 (MQ=255)
                               cGGTACAGACACGATGGCGTATACCGCCTCTGCACTGTCGTTCATGCTCGAGAATCTCGGTAAACCg                                    <  1:351666/67‑1 (MQ=255)
                               cGGTACAGACACGATGGCGTATACCGCCTCTGCACTGTCGTTCATGCTCGAGAATCTCGGTAAACCg                                    <  1:353679/67‑1 (MQ=255)
                               cGGTACAGACACGATGGCGTATACCGCCTCTGCACTGTCGTTCATGCTCGAGAATCTCGGTAAACCg                                    <  1:395790/67‑1 (MQ=255)
                               cGGTACAGACACGATGGCGTATACCGCCTCTGCACTGTCGTTCATGCTCGAGAATCTCGGTAAACCg                                    <  1:502669/67‑1 (MQ=255)
                               cGGTACAGACACGATGGCGTATACCGCCTCTGCACTGTCGTTCATGCTCGAGAATCTCGGTAAACCg                                    <  1:2169870/67‑1 (MQ=255)
                               cGGTACAGACACGATGGCGTATACCGCCTCTGCACTGTCGTTCATGCTCGAGAATCTCGGTAAACCg                                    <  1:531040/67‑1 (MQ=255)
                               cGGTACAGACACGATGGCGTATACCGCCTCTGCACTGTCGTTCATGCTCGAGAATCTCGGTAAACCg                                    <  1:541920/67‑1 (MQ=255)
                               cGGTACAGACACGATGGCGTATACCGCCTCTGCACTGTCGTTCATGCTCGAGAATCTCGGTAAACCg                                    <  1:157458/67‑1 (MQ=255)
                               cGGTACAGACACGATGGCGTATACCGCCTCTGCACTGTCGTTCATGCTCGAGAATCTCGGTAAACCg                                    <  1:605519/67‑1 (MQ=255)
                               cGGTACAGACACGATGGCGTATACCGCCTCTGCACTGTCGTTCATGCTCGAGAATCTCGGTAAACCg                                    <  1:1072775/67‑1 (MQ=255)
                               cGGTACAGACACGATGGCGTATACCGCCTCTGCACTGTCGTTCATGCTCGAGAATCTCGGTAAACCg                                    <  1:899626/67‑1 (MQ=255)
                               cGGTACAGACACGATGGCGTATACCGCCTCTGCACTGTCGTTCATGCTCGAGAATCTCGGTAAACCg                                    <  1:105127/67‑1 (MQ=255)
                               cGGTACAGACACGATGGCGTATACCGCCTCTGCACTGTCGTTCATGCTCGAGAATCTCGGTAAACCg                                    <  1:2190353/67‑1 (MQ=255)
                               cGGTACAGACACGATGGCGTATACCGCCTCTGCACTGTCGTTCATGCTCGAGAATCTCGGTAAACCg                                    <  1:236701/67‑1 (MQ=255)
                               cGGTACAGACACGATGGCGTATACCGCCTCTGCACTGTCGTTCATGCTCGAGAATCTCGGTAAACCg                                    <  1:1919565/67‑1 (MQ=255)
                               cGGTACAGACACGATGGCGTATACCGCCTCTGCACTGTCGTTCATGCTCGAGAATCTCGGTAAACCg                                    <  1:248783/67‑1 (MQ=255)
                               cGGTACAGACACGATGGCGTATACCGCCTCTGCACTGTCGTTCATGCTCGAGAATCTCGGTAAACCg                                    <  1:1883401/67‑1 (MQ=255)
                               cGGTACAGACACGATGGCGTATACCGCCTCTGCACTGTCGTTCATGCTCGAGAATCTCGGTAAACCg                                    <  1:2629084/67‑1 (MQ=255)
                               cGGTACAGACACGATGGCGTATACCGCCTCTGCACTGTCGTTCATGCTCGAGAATCTCGGTAAACCg                                    <  1:2645369/67‑1 (MQ=255)
                               cGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGt                                  >  1:1936082/1‑69 (MQ=255)
                                  cACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCAt                               <  1:369061/69‑1 (MQ=255)
                                  cACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCAt                               <  1:1267678/69‑1 (MQ=255)
                                  cACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCAt                               <  1:1505281/69‑1 (MQ=255)
                                  cACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCAt                               <  1:1845956/69‑1 (MQ=255)
                                  cACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCAt                               <  1:3310876/69‑1 (MQ=255)
                                  cACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCAt                               <  1:2949516/69‑1 (MQ=255)
                                  cACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCAt                               <  1:2767495/69‑1 (MQ=255)
                                   aCCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCAt                               <  1:269187/68‑1 (MQ=255)
                                                  tataCCGCCTCTGCACTGTCGTTCATGCTCGAGAATCTCGGTAAACCg                                    <  1:1813262/48‑1 (MQ=255)
                                                      ccGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATcc           >  1:2661291/1‑69 (MQ=255)
                                                      ccGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATcc           >  1:1498439/1‑69 (MQ=255)
                                                      ccGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATcc           >  1:43005/1‑69 (MQ=255)
                                                      ccGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATcc           >  1:1596406/1‑69 (MQ=255)
                                                      ccGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATc            >  1:2323726/1‑68 (MQ=255)
                                                          ctctTCACTGTCGTTCATGCTCGAGAATCTCGGTAAACCg                                    <  1:3385929/40‑1 (MQ=38)
                                                             tGCGCTGTCTTTCATGCTCGAGAATCTCGGTACACCGGt                                  <  1:441874/39‑1 (MQ=37)
                                                             tGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGt                                  <  1:1387520/39‑1 (MQ=255)
                                                             tGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGt                                  <  1:1890368/39‑1 (MQ=255)
                                                             tGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGt                                  <  1:1886441/39‑1 (MQ=255)
                                                             tGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGt                                  <  1:3026054/39‑1 (MQ=255)
                                                             tGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGt                                  <  1:2549456/39‑1 (MQ=255)
                                                             tGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGt                                  <  1:391623/39‑1 (MQ=255)
                                                             tGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGt                                  <  1:2767888/39‑1 (MQ=255)
                                                             tGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGt                                  <  1:1622862/39‑1 (MQ=255)
                                                             tGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGt                                  <  1:3118564/39‑1 (MQ=255)
                                                             tGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGt                                  <  1:2037762/39‑1 (MQ=255)
                                                             tGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTcaca                >  1:3152890/1‑57 (MQ=255)
                                                             tGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTcaca                >  1:2732548/1‑57 (MQ=255)
                                                             tGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTcaca                >  1:1025974/1‑57 (MQ=255)
                                                             tGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTcaca                >  1:3209998/1‑57 (MQ=255)
                                                             tGCGCTGTAGTTCATGCTCGAGAATCTCGGTAAACCGGt                                  <  1:181879/39‑1 (MQ=255)
                                                             tGCGCTGGCGTTCATGCTCGAGAATCTCGGTAAACCGGt                                  <  1:398099/39‑1 (MQ=255)
                                                               cgcTGTCGTTCATGCTCGATAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCgctggctg  >  1:298175/1‑69 (MQ=255)
                                                               cgcTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCgctggctg  >  1:60775/1‑69 (MQ=255)
                                                               cgcTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCgctggctg  >  1:651637/1‑69 (MQ=255)
                                                                |                                                                   
ATGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTG  >  minE/1178786‑1178917

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