Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,178,823 C→A 45.7% T91T (ACC→ACA ansA → cytoplasmic L‑asparaginase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,178,8230CA45.7% 33.4 / 95.7 94T91T (ACC→ACAansAcytoplasmic L‑asparaginase I
Reads supporting (aligned to +/- strand):  ref base C (21/30);  new base A (19/24);  total (40/54)
Fisher's exact test for biased strand distribution p-value = 8.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.68e-01

GATATTAAAGCGCACTATGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCAT  >  minE/1178770‑1178888
                                                     |                                                                 
gATATTAAAGCGCACTATGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGtat                                                    >  1:2250596/1‑69 (MQ=255)
gATATTAAAGCGCACTATGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGtat                                                    >  1:2333695/1‑69 (MQ=255)
gATATTAAAGCGCACTATGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGtat                                                    >  1:2413705/1‑69 (MQ=255)
gATATTAAAGCGCACTATGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGtat                                                    >  1:2520734/1‑69 (MQ=255)
gATATTAAAGCGCACTATGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGtat                                                    >  1:3052050/1‑69 (MQ=255)
gATATTAAAGCGCACTATGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGtat                                                    >  1:440705/1‑69 (MQ=255)
gATATTAAAGCGCACTATGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGtat                                                    >  1:3366166/1‑69 (MQ=255)
      aaaGCGCACTATGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGc                                               <  1:2917763/68‑1 (MQ=255)
      aaaGCGCACTATGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGc                                               <  1:1299201/68‑1 (MQ=255)
      aaaGCGCACTATGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGc                                               <  1:134140/68‑1 (MQ=255)
      aaaGCGCACTATGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGc                                               <  1:2030140/68‑1 (MQ=255)
      aaaGCGCACTATGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGc                                               <  1:351305/68‑1 (MQ=255)
      aaaGCGCACTATGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGc                                               <  1:2100216/68‑1 (MQ=255)
                aTGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGacac                                                              >  1:1504711/1‑43 (MQ=255)
                aTGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGacac                                                              >  1:1508569/1‑43 (MQ=255)
                aTGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGacac                                                              >  1:621147/1‑43 (MQ=255)
                aTGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGacac                                                              >  1:10720/1‑43 (MQ=255)
                aTGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCt                                      >  1:215505/1‑67 (MQ=255)
                aTGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGt                                    >  1:2486516/1‑69 (MQ=255)
                        tATGATGGTTTTGTCATTCTGCACGGCACCGacac                                                              >  1:1752567/1‑35 (MQ=255)
                        tATGATGGTTTTGTCATTCTGCACGGCACCGacac                                                              >  1:2443668/1‑35 (MQ=255)
                        tATGATGGTTTTGTCATTCTGCACGGCACCGacac                                                              >  1:2663560/1‑35 (MQ=255)
                        tATGATGGTTTTGTCATTCTGCACGGCACCGacac                                                              >  1:2667456/1‑35 (MQ=255)
                        tATGATGGTTTTGTCATTCTGCACGGCACCGacac                                                              >  1:1832364/1‑35 (MQ=255)
                        tATGATGGTTTTGTCATTCTGCACGGCACCGacac                                                              >  1:3103171/1‑35 (MQ=255)
                        tATGATGGTTTTGTCATTCTGCACGGCACCGacac                                                              >  1:960295/1‑35 (MQ=255)
                            atgGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGt                                 <  1:2973523/60‑1 (MQ=255)
                            atgGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGt                                 <  1:2634533/60‑1 (MQ=255)
                            atgGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGt                                 <  1:2036724/60‑1 (MQ=255)
                            atgGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGt                                 <  1:1639198/60‑1 (MQ=255)
                            atgGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGt                                 <  1:1435850/60‑1 (MQ=255)
                            atgGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGt                                 <  1:557825/60‑1 (MQ=255)
                            atgGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGt                                 <  1:2257998/60‑1 (MQ=255)
                                 tttGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGaa                    <  1:2758563/68‑1 (MQ=255)
                                 tttGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGaa                    <  1:993890/68‑1 (MQ=255)
                                 tttGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGaa                    <  1:3395173/68‑1 (MQ=255)
                                 tttGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGaa                    <  1:637998/68‑1 (MQ=255)
                                 tttGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGaa                    <  1:850055/68‑1 (MQ=255)
                                 tttGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGaa                    <  1:230257/68‑1 (MQ=255)
                                 tttGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGaa                    <  1:1176718/68‑1 (MQ=255)
                                 tttGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGaa                    <  1:2028356/68‑1 (MQ=255)
                                  ttGTCATTCTGCACGGTACAGACACGATGGCGTATACCGCCTCTGCACTGTCGTTCATGCTCGAGAAtc                  >  1:3141321/1‑69 (MQ=255)
                                  ttGTCATTCTGCACGGTACAGACACGATGGCGTATACCGCCTCTGCACTGTCGTTCATGCTCGAGAAtc                  >  1:1619155/1‑69 (MQ=255)
                                  ttGTCATTCTGCACGGTACAGACACGATGGCGTATACCGCCTCTGCACTGTCGTTCATGCTCGAGAAtc                  >  1:1498042/1‑69 (MQ=255)
                                  ttGTCATTCTGCACGGTACAGACACGATGGCGTATACCGCCTCTGCACTGTCGTTCATGCTCGAGAAtc                  >  1:495347/1‑69 (MQ=255)
                                  ttGTCATTCTGCACGGTACAGACACGATGGCGTATACCGCCTCTGCACTGTCGTTCATGCTCGAGAAtc                  >  1:1948663/1‑69 (MQ=255)
                                  ttGTCATTCTGCACGGTACAGACACGATGGCGTATACCGCCTCTGCACTGTCGTTCATGCTCGAGAAtc                  >  1:1353197/1‑69 (MQ=255)
                                  ttGTCATTCTGCACGGTACAGACACGATGGCGTATACCGCCTCTGCACTGTCGTTCATGCTCGAGAAtc                  >  1:1963260/1‑69 (MQ=255)
                                  ttGTCATTCTGCACGGTACAGACACGATGGCGTATACCGCCTCTGCACTGTCGTTCATGCTCGAGAAtc                  >  1:120273/1‑69 (MQ=255)
                                  ttGTCATTCTGCACGGTACAGACACGATGGCGTATACCGCCTCTGCACTGTCGTTCATGCTCGAGAAtc                  >  1:2559906/1‑69 (MQ=255)
                                  ttGTCATTCTGCACGGTACAGACACGATGGCGTATACCGCCTCTGCACTGTCGTTCATGCTCGAGAAtc                  >  1:804464/1‑69 (MQ=255)
                                  ttGTCATTCTGCACGGTACAGACACGATGGCGTATACCGCCTCTGCACTGTCGTTCATGCTCGAGAAtc                  >  1:2483896/1‑69 (MQ=255)
                                  ttGTCATTCTGCACGGTACAGACACGATGGCGTATACCGCCTCTGCACTGTCGTTCATGCTCGAGAAtc                  >  1:819074/1‑69 (MQ=255)
                                  ttGTCATTCTGCACGGTACAGACACGATGGCGTATACCGCCTCTGCACTGTCGTTCATGCTCGAGAAtc                  >  1:243729/1‑69 (MQ=255)
                                  ttGTCATTCTGCACGGTACAGACACGATGGCGTATACCGCCTCTGCACTGTCGTTCATGCTCGAGAAtc                  >  1:838764/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:351666/67‑1 (MQ=255)
                                               cGGTACAGACACGATGGCGTATACCGCCTCTGCACTGTCGTTCATGCTCGAGAATCTCGGTAAACCg       <  1:2629084/67‑1 (MQ=255)
                                               cGGTACAGACACGATGGCGTATACCGCCTCTGCACTGTCGTTCATGCTCGAGAATCTCGGTAAACCg       <  1:248783/67‑1 (MQ=255)
                                               cGGTACAGACACGATGGCGTATACCGCCTCTGCACTGTCGTTCATGCTCGAGAATCTCGGTAAACCg       <  1:605519/67‑1 (MQ=255)
                                               cGGTACAGACACGATGGCGTATACCGCCTCTGCACTGTCGTTCATGCTCGAGAATCTCGGTAAACCg       <  1:2645369/67‑1 (MQ=255)
                                               cGGTACAGACACGATGGCGTATACCGCCTCTGCACTGTCGTTCATGCTCGAGAATCTCGGTAAACCg       <  1:541920/67‑1 (MQ=255)
                                               cGGTACAGACACGATGGCGTATACCGCCTCTGCACTGTCGTTCATGCTCGAGAATCTCGGTAAACCg       <  1:531040/67‑1 (MQ=255)
                                               cGGTACAGACACGATGGCGTATACCGCCTCTGCACTGTCGTTCATGCTCGAGAATCTCGGTAAACCg       <  1:514743/67‑1 (MQ=255)
                                               cGGTACAGACACGATGGCGTATACCGCCTCTGCACTGTCGTTCATGCTCGAGAATCTCGGTAAACCg       <  1:502669/67‑1 (MQ=255)
                                               cGGTACAGACACGATGGCGTATACCGCCTCTGCACTGTCGTTCATGCTCGAGAATCTCGGTAAACCg       <  1:2169870/67‑1 (MQ=255)
                                               cGGTACAGACACGATGGCGTATACCGCCTCTGCACTGTCGTTCATGCTCGAGAATCTCGGTAAACCg       <  1:395790/67‑1 (MQ=255)
                                               cGGTACAGACACGATGGCGTATACCGCCTCTGCACTGTCGTTCATGCTCGAGAATCTCGGTAAACCg       <  1:2190353/67‑1 (MQ=255)
                                               cGGTACAGACACGATGGCGTATACCGCCTCTGCACTGTCGTTCATGCTCGAGAATCTCGGTAAACCg       <  1:353679/67‑1 (MQ=255)
                                               cGGTACAGACACGATGGCGTATACCGCCTCTGCACTGTCGTTCATGCTCGAGAATCTCGGTAAACCg       <  1:236701/67‑1 (MQ=255)
                                               cGGTACAGACACGATGGCGTATACCGCCTCTGCACTGTCGTTCATGCTCGAGAATCTCGGTAAACCg       <  1:157458/67‑1 (MQ=255)
                                               cGGTACAGACACGATGGCGTATACCGCCTCTGCACTGTCGTTCATGCTCGAGAATCTCGGTAAACCg       <  1:3379424/67‑1 (MQ=255)
                                               cGGTACAGACACGATGGCGTATACCGCCTCTGCACTGTCGTTCATGCTCGAGAATCTCGGTAAACCg       <  1:1919565/67‑1 (MQ=255)
                                               cGGTACAGACACGATGGCGTATACCGCCTCTGCACTGTCGTTCATGCTCGAGAATCTCGGTAAACCg       <  1:3362331/67‑1 (MQ=255)
                                               cGGTACAGACACGATGGCGTATACCGCCTCTGCACTGTCGTTCATGCTCGAGAATCTCGGTAAACCg       <  1:899626/67‑1 (MQ=255)
                                               cGGTACAGACACGATGGCGTATACCGCCTCTGCACTGTCGTTCATGCTCGAGAATCTCGGTAAACCg       <  1:1072775/67‑1 (MQ=255)
                                               cGGTACAGACACGATGGCGTATACCGCCTCTGCACTGTCGTTCATGCTCGAGAATCTCGGTAAACCg       <  1:1883401/67‑1 (MQ=255)
                                               cGGTACAGACACGATGGCGTATACCGCCTCTGCACTGTCGTTCATGCTCGAGAATCTCGGTAAACCg       <  1:1727794/67‑1 (MQ=255)
                                               cGGTACAGACACGATGGCGTATACCGCCTCTGCACTGTCGTTCATGCTCGAGAATCTCGGTAAACCg       <  1:3050909/67‑1 (MQ=255)
                                               cGGTACAGACACGATGGCGTATACCGCCTCTGCACTGTCGTTCATGCTCGAGAATCTCGGTAAACCg       <  1:105127/67‑1 (MQ=255)
                                               cGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGt     >  1:1936082/1‑69 (MQ=255)
                                                  cACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCAt  <  1:1845956/69‑1 (MQ=255)
                                                  cACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCAt  <  1:1267678/69‑1 (MQ=255)
                                                  cACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCAt  <  1:369061/69‑1 (MQ=255)
                                                  cACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCAt  <  1:1505281/69‑1 (MQ=255)
                                                  cACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCAt  <  1:3310876/69‑1 (MQ=255)
                                                  cACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCAt  <  1:2767495/69‑1 (MQ=255)
                                                  cACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCAt  <  1:2949516/69‑1 (MQ=255)
                                                   aCCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCAt  <  1:269187/68‑1 (MQ=255)
                                                     |                                                                 
GATATTAAAGCGCACTATGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCAT  >  minE/1178770‑1178888

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: