Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,179,174 C→G 56.6% T208T (ACC→ACG ansA → cytoplasmic L‑asparaginase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,179,1740CG56.6% 22.5 / 78.4 76T208T (ACC→ACGansAcytoplasmic L‑asparaginase I
Reads supporting (aligned to +/- strand):  ref base C (28/5);  new base G (25/18);  total (53/23)
Fisher's exact test for biased strand distribution p-value = 1.34e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTCGTCGTTTGAATACGCCACCCGCCCCGCACGGTGAAGGGGAATTAATCGTTCATCCAATCACCCCACAACCAATTGGCGTAGTGACGATTTATCCAGGGATTTCTGCTGACGTCGTGCGCAATTTT  >  minE/1179110‑1179237
                                                                |                                                               
ttCGTCGTTTGAATACGCCACCCGCCCCGCTCGGTGAAGGGGAATTAATCGTTCATCCAATCACCCcac                                                             >  1:141285/1‑69 (MQ=255)
ttCGTCGTTTGAATACGCCACCCGCCCCGCACGGTGAAGGGGAATTAATCGTTCATCCAATCACCCcac                                                             >  1:3331015/1‑69 (MQ=255)
ttCGTCGTTTGAATACGCCACCCGCCCCGCACGGTGAAGGGGAATTAATCGTTCATCCAATCACCCcac                                                             >  1:207845/1‑69 (MQ=255)
ttCGTCGTTTGAATACGCCACCCGCCCCGCACGGTGAAGGGGAATTAATCGTTCATCCAATCACCCcac                                                             >  1:2945645/1‑69 (MQ=255)
ttCGTCGTTTGAATACGCCACCCGCCCCGCACGGTGAAGGGGAATTAATCGTTCATCCAATCACCCcac                                                             >  1:3085129/1‑69 (MQ=255)
ttCGTCGTTTGAATACGCCACCCGCCCCGCACGGTGAAGGGGAATTAATCGTTCATCCAATCACCCcac                                                             >  1:2091867/1‑69 (MQ=255)
ttCGTCGTTTGAATACGCCACCCGCCCCGCACGGTGAAGGGGAATTAATCGTTCATCCAATCACCCcac                                                             >  1:3095967/1‑69 (MQ=255)
ttCGTCGTTTGAATACGCCACCCGCCCCGCACGGTGAAGGGGAATTAATCGTTCATCCAATCACCCcac                                                             >  1:2280430/1‑69 (MQ=255)
ttCGTCGTTTGAATACGCCACCCGCCCCGCACGGTGAAGGGGAATTAATCGTTCATCCAATCACCCcac                                                             >  1:370440/1‑69 (MQ=255)
ttCGTCGTTTGAATACGCCACCCGCCCCGCACGGTGAAGGGGAATTAATCGTTCATCCAATCACCCcac                                                             >  1:2750736/1‑69 (MQ=255)
ttCGTCGTTTGAATACGCCACCCGCCCCGCACGGTGAAGGGGAATTAATCGTTCATCCAATCACCCcac                                                             >  1:597804/1‑69 (MQ=255)
ttCGTCGTTTGAATACGCCACCCGCCCCGCACGGTGAAGGGGAATTAATCGTTCATCCAATCACCCcac                                                             >  1:610276/1‑69 (MQ=255)
ttCGTCGTTTGAATACGCCACCCGCCCCGCACGGTGAAGGGGAATTAATCGTTCATCCAATCACCCcac                                                             >  1:694236/1‑69 (MQ=255)
ttCGTCGTTTGAATACGCCACCCGCCCCGCACGGTGAAGGGGAATTAATCGTTCATCCAATCACCCcac                                                             >  1:1448115/1‑69 (MQ=255)
ttCGTCGTTTGAATACGCCACCCGCCCCGCACGGTGAAGGGGAATTAATCGTTCATCCAATCACCCcac                                                             >  1:781801/1‑69 (MQ=255)
ttCGTCGTTTGAATACGCCACCCGCCCCGCACGGTGAAGGGGAATTAATCGTTCATCCAATCACCCcac                                                             >  1:3422776/1‑69 (MQ=255)
ttCGTCGTTTGAATACGCCACCCGCCCCGCACGGTGAAGGGGAATTAATCGTTCATCCAATCACCCcac                                                             >  1:1217556/1‑69 (MQ=255)
     cgtTTGAATACGCCACCCGCCCCGCACGGTGAAGGGGAATTAATCGTTCATCCAATCACCCCACAAcc                                                         <  1:3219146/68‑1 (MQ=255)
         tGAATACGCCACCCGCCCCGCACGGTGAAGGGGAATTAATCGTTCATCCAATCACCCCACAACCAATTg                                                    >  1:2749567/1‑69 (MQ=255)
         tGAATACGCCACCCGCCCCGCACGGTGAAGGGGAATTAATCGTTCATCCAATCACCCCACAACCAATTg                                                    >  1:581399/1‑69 (MQ=255)
         tGAATACGCCACCCGCCCCGCACGGTGAAGGGGAATTAATCGTTCATCCAATCACCCCACAACCAATTg                                                    >  1:771047/1‑69 (MQ=255)
         tGAATACGCCACCCGCCCCGCACGGTGAAGGGGAATTAATCGTTCATCCAATCACCCCACAACCAATTg                                                    >  1:3110864/1‑69 (MQ=255)
         tGAATACGCCACCCGCCCCGCACGGTGAAGGGGAATTAATCGTTCATCCAATCACCCCACAACCAATTg                                                    >  1:2176797/1‑69 (MQ=255)
         tGAATACGCCACCCGCCCCGCACGGTGAAGGGGAATTAATCGTTCATCCAATCACCCCACAACCAATTg                                                    >  1:426595/1‑69 (MQ=255)
         tGAATACGCCACCCGCCCCGCACGGTGAAGGGGAATTAATCGTTCATCCAATCACCCCACAACCAATTg                                                    >  1:3398590/1‑69 (MQ=255)
             tACGCCACCCGCCCCGCACGGTGAAGGGGAATTAATCGTTCATCCAATCACCCCACAACCAATTGGCGt                                                <  1:2806031/69‑1 (MQ=255)
                   aCCCGTCCCGCCCGGTGAAGGGGCATTGATCGTTCATCCAATCACCCCGCAACCAATTg                                                    >  1:1928423/1‑59 (MQ=255)
                   aCCCGTCCCGCACGGTGAAGGGGCATTGATCGTTCATCCAATCACGCCGCAACCAATTGGCGTAGTGa                                           >  1:3196625/1‑68 (MQ=255)
                   aCCCGTCCCGCACGGTGAAGGGGCATTGATCGTTCATCCAATCACGCCGCAACCAATTGGCGTAGTGAc                                          >  1:2392907/1‑69 (MQ=255)
                   aCCCGTCCCGCACGGTGAAGGGGCATTGATCGTTCATCCAATCACGCCGCAACCAATTGGCGTAGTGAc                                          >  1:2708675/1‑69 (MQ=255)
                   aCCCGTCCCGCACGGTGAAGGGGCATTGATCGTTCATCCAATCACGCCGCAACCAATTGGCGTAGTGAc                                          >  1:2697765/1‑69 (MQ=255)
                   aCCCGTCCCGCACGGTGAAGGGGCATTGATCGTTCATCCAATCACGCCGCAACCAATTGGCGTAGTGAc                                          >  1:2691434/1‑69 (MQ=255)
                   aCCCGTCCCGCACGGTGAAGGGGCATTGATCGTTCATCCAATCACGCCGCAACCAATTGGCGTAGTGAc                                          >  1:1332834/1‑69 (MQ=255)
                   aCCCGTCCCGCACGGTGAAGGGGCATTGATCGTTCATCCAATCACGCCGCAACCAATTGGCGTAGTGAc                                          >  1:2563275/1‑69 (MQ=255)
                   aCCCGTCCCGCACGGTGAAGGGGCATTGATCGTTCATCCAATCACGCCGCAACCAATTGGCGTAGTGAc                                          >  1:1507378/1‑69 (MQ=255)
                   aCCCGTCCCGCACGGTGAAGGGGCATTGATCGTTCATCCAATCACGCCGCAACCAATTGGCGTAGTGAc                                          >  1:2069469/1‑69 (MQ=255)
                   aCCCGTCCCGCACGGTGAAGGGGCATTGATCGTTCATCCAATCACGCCGCAACCAATTGGCGTAGTGAc                                          >  1:2758716/1‑69 (MQ=255)
                   aCCCGTCCCGCACGGTGAAGGGGCATTGATCGTTCATCCAATCACGCCGCAACCAATTGGCGTAGTGAc                                          >  1:3146297/1‑69 (MQ=255)
                   aCCCGTCCCGCACGGTGAAGGGGCATTGATCGTTCATCCAATCACGCCGCAACCAATTGGCGTAGTGAc                                          >  1:1852521/1‑69 (MQ=255)
                   aCCCGTCCCGCACGGTGAAGGGGCATTGATCGTTCATCCAATCACGCCGCAACCAATTGGCGTAGTGAc                                          >  1:210540/1‑69 (MQ=255)
                   aCCCGTCCCGCACGGTGAAGGGGCATTGATCGTTCATCCAATCACGCCGCAACCAATTGGCGTAGTGAc                                          >  1:2910662/1‑69 (MQ=255)
                      cgccccgcACGGTGAAGGGGAATTAATCGTTCATCCAATCACCCCACAACCAATTGGCGTAGTGACGa                                        >  1:353751/1‑68 (MQ=255)
                      cgccccgcACGGTGAAGGGGAATTAATCGTTCATCCAATCACCCCACAACCAATTGGCGTAGTGACGa                                        >  1:3105302/1‑68 (MQ=255)
                                ggTGAAGGGGAATTAATCGTTCATCCAATCACCCCACAACCAATTGGCGTAGTGACGATTTATCCAggg                             <  1:3333376/69‑1 (MQ=255)
                                ggTGAAGGGGAATTAATCGTTCATCCAATCACCCCACAACCAATTGGCGTAGTGACGATTTATCCAggg                             <  1:1889677/69‑1 (MQ=255)
                                ggTGAAGGGGAATTAATCGTTCATCCAATCACCCCACAACCAATTGGCGTAGTGACGATTTATCCAggg                             <  1:1177715/69‑1 (MQ=255)
                                      ggggCATTGATCGTTCATCCAATCACGCCGCAACCaa                                                       >  1:2062139/1‑37 (MQ=21)
                                      ggggCATTGATCGTTCATCCAATCACGCCGCAACCAATTGGc                                                  >  1:1625688/1‑42 (MQ=25)
                                      ggggCATTGATCGTTCATCCAATCACGCCGCAACCAATTGGCg                                                 >  1:1006061/1‑43 (MQ=25)
                                      ggggCATTGATCGTTCATCCAATCACGCCGCAACCAATTGGCg                                                 >  1:2352389/1‑43 (MQ=25)
                                      ggggCATTGATCGTTCATCCAATCACGCCGCAACCAATTGGCg                                                 >  1:422942/1‑43 (MQ=25)
                                      ggggCATTGATCGTTCATCCAATCACGCCGCAACCAATTGGCg                                                 >  1:1996278/1‑43 (MQ=25)
                                      ggggCATTGATCGTTCATCCAATCACGCCGCAACCAATTGGCg                                                 >  1:1987074/1‑43 (MQ=25)
                                      ggggCATTGATCGTTCATCCAATCACGCCGCAACCAATTGGCg                                                 >  1:1704986/1‑43 (MQ=25)
                                      ggggCATTGATCGTTCATCCAATCACGCCGCAACCAATTGGCg                                                 >  1:52859/1‑43 (MQ=25)
                                      ggggCATTGATCGTTCATCCAATCACGCCGCAACCAATTGGCg                                                 >  1:1336958/1‑43 (MQ=25)
                                      ggggCATTGATCGTTCATCCAATCACGCCGCAACCAATTGGCg                                                 >  1:887807/1‑43 (MQ=25)
                                                  gTTCATCCAATCACCCCACAACCAATTGGCGTAGTGACGATTTATCCAGGGATTTCTGCTGACGTCGTg           >  1:2452893/1‑69 (MQ=255)
                                                           aTCACGCCGCAACCAATTGGCGTAGTGACGATTTATCCGGGGATTTCTGCTGACGTCGTGCGCAAtttt  <  1:495672/69‑1 (MQ=255)
                                                           aTCACGCCGCAACCAATTGGCGTAGTGACGATTTATCCGGGGATTTCTGCTGACGTCGTGCGCAAtttt  <  1:2538834/69‑1 (MQ=255)
                                                           aTCACGCCGCAACCAATTGGCGTAGTGACGATTTATCCGGGGATTTCTGCTGACGTCGTGCGCAAtttt  <  1:3362413/69‑1 (MQ=255)
                                                           aTCACGCCGCAACCAATTGGCGTAGTGACGATTTATCCGGGGATTTCTGCTGACGTCGTGCGCAAtttt  <  1:3286217/69‑1 (MQ=255)
                                                           aTCACGCCGCAACCAATTGGCGTAGTGACGATTTATCCGGGGATTTCTGCTGACGTCGTGCGCAAtttt  <  1:2577322/69‑1 (MQ=255)
                                                           aTCACGCCGCAACCAATTGGCGTAGTGACGATTTATCCGGGGATTTCTGCTGACGTCGTGCGCAAtttt  <  1:2517129/69‑1 (MQ=255)
                                                           aTCACGCCGCAACCAATTGGCGTAGTGACGATTTATCCGGGGATTTCTGCTGACGTCGTGCGCAAtttt  <  1:229445/69‑1 (MQ=255)
                                                           aTCACGCCGCAACCAATTGGCGTAGTGACGATTTATCCGGGGATTTCTGCTGACGTCGTGCGCAAtttt  <  1:1890136/69‑1 (MQ=255)
                                                           aTCACGCCGCAACCAATTGGCGTAGTGACGATTTATCCGGGGATTTCTGCTGACGTCGTGCGCAAtttt  <  1:1700683/69‑1 (MQ=255)
                                                           aTCACGCCGCAACCAATTGGCGTAGTGACGATTTATCCGGGGATTTCTGCTGACGTCGTGCGCAAtttt  <  1:1603909/69‑1 (MQ=255)
                                                           aTCACGCCGCAACCAATTGGCGTAGTGACGATTTATCCGGGGATTTCTGCTGACGTCGTGCGCAAtttt  <  1:1570382/69‑1 (MQ=255)
                                                           aTCACGCCGCAACCAATTGGCGTAGTGACGATTTATCCGGGGATTTCTGCTGACGTCGTGCGCAAtttt  <  1:1483324/69‑1 (MQ=255)
                                                           aTCACGCCGCAACCAATTGGCGTAGTGACGATTTATCCGGGGATTTCTGCTGACGTCGTGCGCAAtttt  <  1:1480041/69‑1 (MQ=255)
                                                           aTCACGCCGCAACCAATTGGCGTAGTGACGATTTATCCGGGGATTTCTGCTGACGTCGTGCGCAAtttt  <  1:1429643/69‑1 (MQ=255)
                                                           aTCACGCCGCAACCAATTGGCGTAGTGACGATTTATCCGGGGATTTCTGCTGACGTCGTGCGCAAtttt  <  1:1357580/69‑1 (MQ=255)
                                                           aTCACGCCGCAACCAATTGGCGTAGTGACGATTTATCCGGGGATTTCTGCTGACGTCGTGCGCAAtttt  <  1:1153427/69‑1 (MQ=255)
                                                            tCACGCCGCAACCAATTGGCGTAGTGACGATTTATCCGGGGATTTCTGCTGACGTCGTGCGCAAtttt  <  1:2001810/68‑1 (MQ=255)
                                                            tCACGCCGCAACCAATTGGCGTAGTGACGATTTATCCGGGGATTTCTGCTGACGTCGTGCGCAAtttt  <  1:1641983/68‑1 (MQ=255)
                                                                |                                                               
TTCGTCGTTTGAATACGCCACCCGCCCCGCACGGTGAAGGGGAATTAATCGTTCATCCAATCACCCCACAACCAATTGGCGTAGTGACGATTTATCCAGGGATTTCTGCTGACGTCGTGCGCAATTTT  >  minE/1179110‑1179237

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: