Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,179,532 T→C 43.5% *328Q (TAA→CAA)  ansA → cytoplasmic L‑asparaginase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,179,5320TC43.5% 13.1 / 54.8 46*328Q (TAA→CAA) ansAcytoplasmic L‑asparaginase I
Reads supporting (aligned to +/- strand):  ref base T (18/8);  new base C (9/11);  total (27/19)
Fisher's exact test for biased strand distribution p-value = 1.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTAACCAAACTGCATTACCTGCTTAGCCAGGAACTGGATACTGAAACCATTCGCAAGGCCATGAGCTAAAACCTGCGCGGCGAACTGACGCTTTGCTCACATCTCACTTTAATCGTGCTCACAT  >  minE/1179465‑1179589
                                                                   |                                                         
gCTAACCAAACTGCATTACCTGCTTAGCTAGGAACTGGATACTGAAACCATTCGCAAGGCCATGAGCCaa                                                         >  1:58204/1‑70 (MQ=255)
gCTAACCAAACTGCATTACCTGCTTAGCCAGGAACTGGATACTGAAACCATTCGCAAGGCCATGAGCTaa                                                         >  1:1008528/1‑70 (MQ=255)
gCTAACCAAACTGCATTACCTGCTTAGCCAGGAACTGGATACTGAAACCATTCGCAAGGCCATGAGCTaa                                                         >  1:3641307/1‑70 (MQ=255)
    aCCAAACTGCATTACCTGCTTAGCCAGGAACTGGATACTGAAACCATTCGCAAGGCCATGAGCTAAAAcc                                                     <  1:605413/70‑1 (MQ=255)
    aCCAAACTGCATTACCTGCTTAGCCAGGAACTGGATACTGAAACCATTCGCAAGGCCATGAGCCaaa                                                        <  1:1412741/67‑1 (MQ=255)
    aCCAAACTGCATTACCTGCTTAGCCAGGAACTGGATACTGAAACCATTCGCAAGGCCATGAGCCAAAAc                                                      <  1:3546076/69‑1 (MQ=255)
      cAAACTGCATTACCTGCTTAGACAGGAACTGGATACTGAAACCATTCGCAAGGCCATGAGCCAAAACCTgc                                                  >  1:193328/1‑71 (MQ=255)
      cAAACTGCATTACATGCTTAGCCAGGAACTGGATACTGAAACCATTCGCAAGGCCATGAGCCaaaa                                                       >  1:3277624/1‑66 (MQ=255)
        aaCTGCATTACCTGCTTAGCCAGGAACTGGATACTGAAACCATTCGCAAGGCCATGAGCTaaaa                                                       >  1:1786279/1‑64 (MQ=255)
        aaCTGCATTACCTGCTTAGCCAGGAACTGGATACTGAAACCATTCGCAAGGCCATGAGCCAAAACCTgcgc                                                <  1:942166/71‑1 (MQ=255)
         aCTGCATTACCTGCTTAGCCAGGAACTGGATACTGAAACCATTCGCAAGGCCATGAGCCAAAACCTGcgcg                                               >  1:3244723/1‑71 (MQ=255)
             cATTACCTGCTTAGCCAGGAACTGGATACTGAAACCATTCGCAAGGCCATGAGCTAAAACCTGCgcggcg                                            <  1:368319/70‑1 (MQ=255)
               ttACCTGCTTAGCCAGGAACTGGATACTGAAACCATTCGCAAGGCCATGAGCTAAAACCTGCGCGGCGAAc                                         <  1:389046/71‑1 (MQ=255)
               ttACCTGCTTAGCCAGGAACTGGATACTGAAACCATTCGCAAGGCCATGAGCCAAAACCTGCGCGGCGAAc                                         >  1:3641139/1‑71 (MQ=255)
                 aCCTGCTTAGCCAGGAACTGGATACTGAAACCATTCGCAAGGCCATGAGCCAAAAcc                                                     >  1:2047378/1‑57 (MQ=255)
                  ccTGCTTAGCCAGGAACTGGATACTGAAACCATTCGCAAGGCCATGAGCCaaa                                                        <  1:3533053/53‑1 (MQ=255)
                  ccTGCTTAGCCAGGAACTGGATACTGAAACCATTCGCAAGGCCATGAGCCaaa                                                        <  1:1140260/53‑1 (MQ=255)
                  ccTGCTTAGCCAGGAACTGGATACTGAAACCATTCGCAAGGCCATGAGCCAAAACCTGCGCTGCGAACTg                                       <  1:2181000/70‑1 (MQ=255)
                  ccTGCTTAGCCAGGAACTGGATACTGAAACCATTCGCAAGGCCATGAGCCAAAACCTGCGCGGCGAACTGa                                      >  1:696434/1‑71 (MQ=255)
                  ccTGCTTAGCCAGGAACTGGATACTGAAACCATTCGCAAGGCCATGAGCCAAAACCTGCGCGGCGAACTGa                                      >  1:3119199/1‑71 (MQ=255)
                   cTGCTTAGCCAGGAACTGGATACTGAAACCATTCGCAAGGCCATGAGCCAAAACCTGCGCGGCGAACTGAc                                     <  1:1828396/71‑1 (MQ=255)
                      cTTAGCCAGGAACTGGATACTGAAACCATTCGCAAGGCCATGAGCCAAAACCTGCGCGGCGAACTGACg                                    <  1:2990964/69‑1 (MQ=255)
                       ttAGCCAGGAACTGGATACTGAAACCATTCGCAAGGCCATGAGCTAAAACCTGCGCGGCGAACTGACGCtt                                 >  1:2269378/1‑71 (MQ=255)
                             aGGAACTGGATACTGAAACCATTCGCAAGGCCATGAGCTAAAACCTGCGCGGCGAACTGACGCTTTGCTc                            >  1:3874412/1‑70 (MQ=255)
                              ggAACTGGATACTGAAACCATTCGCAAGGCCATGAGCCAAAACCTGCgcgg                                              >  1:179321/1‑51 (MQ=255)
                                    ggATACTGAAACCATTCGCAAGGCCATGAGCTAAAACCTGCGCGGCGAAc                                         >  1:3489145/1‑50 (MQ=255)
                                      aTACTGAAACCATTCGCAAGGCCATGAGCCAAAACCTgcg                                                 <  1:3080344/40‑1 (MQ=255)
                                      aTACTGAAACCATTCGCAAGGCCATGAGCCAAAACCTgcg                                                 <  1:3546221/40‑1 (MQ=255)
                                      aTACTGAAACCATTCGCAAGGCCATGAGCCAAAACCTGCGCGGAGAACTg                                       <  1:4812/50‑1 (MQ=255)
                                         cTGAAACCATTCGCAAGGCCATGAGCTAAAACCTGCGCGGCGAACTGACGCTTTGCTCACATCTCACttt                <  1:2603905/70‑1 (MQ=255)
                                         cTGAAACCATTCGCAAGGCCATGAGCTAAAACCTGCGCGGCGAACTGACGCTTTGCTCACATCTCACttt                <  1:3140001/70‑1 (MQ=255)
                                          tGAAACCATTCGCAAGGCCATGAGCTAAAACCTGCGCGGCGAACTGACGCTTTGCTCACATCTCACTTTaa              >  1:3764208/1‑71 (MQ=255)
                                          tGAAACCATTCGCAAGGCCATGAGCTAAAACCTGCGCGGCGAACTGACGCTTTGCTCACATCTCACTTTaa              >  1:155447/1‑71 (MQ=255)
                                           gAAACCATTCGCAAGGCCATGAGCTAAAACCTGCGCGGCGAACTGACGCTTTg                               >  1:1826578/1‑53 (MQ=255)
                                            aatccATTCGCAAGGCCATGAGCTAAAACCTGCGCGGCGAACTGACGCTTTGCTCACATCTCACTTTaa              <  1:2001524/66‑1 (MQ=255)
                                            aaaCCATTCGCAAGGCCATGAGCTAAAACCTGCGCGGCGAACTGACGCTTTGCTCACATCTCACTTTAATc            >  1:1262637/1‑71 (MQ=255)
                                            aaaCCATTCGCAAGGCCATGAGCTAAAACCTGCGCGGCGAACTGACGCTTTGCTCACATCTCACTTTAATc            >  1:2406704/1‑71 (MQ=255)
                                             aaCCATTCGCAAGGCCATGAGCTAAAACCTGCGCGGCGAACTGACGCTTTGCTCACATCTCACTTTAATCg           <  1:3302578/71‑1 (MQ=255)
                                             aaCCATTCGCAAGGCCATGAGCTAAAACCTGCGCGGCGAACTGACGCTTTGCTCACATCTCACTTTAATCg           >  1:2613304/1‑71 (MQ=255)
                                                   tCGCAAGGCCATGAGCTAAAACCTGCGCGGCGAACTGACGCTTTGCTCACATCTCACTTTAATCGTGCTca     >  1:684693/1‑71 (MQ=255)
                                                    cGCAAGGCCATGAGCTAAAACCTGCGCGGCGAACTGACGCTTTGCTCACATCTCACTTTAATCGTGCTca     >  1:153951/1‑70 (MQ=255)
                                                    cGCAAGGCCATGAGCTAAAACCTGCGCGGCGAACTGACGCTTTGCTCACATCTCACTTTAATCGTGCTca     >  1:1336465/1‑70 (MQ=255)
                                                    cGCAAGGCCATGAGCTAAAACCTGCGCGGCGAACTGACGCTTTGCTCACATCTCACTTTAATCGTGCTca     >  1:2291321/1‑70 (MQ=255)
                                                       aaGGCCATGAGCTAAAACCTGCGCGGCGAACTGACGCTTTGCTCACATCTCACTTTAATCGTGCTCACAt  >  1:671482/1‑70 (MQ=255)
                                                         ggCCATGAGCTAAAACCTGCGCGGCGAACTGACGCTTTGc                              >  1:1991086/1‑40 (MQ=255)
                                                           ccATGAGCTAAAACCTGCGCGGCGAACTGACGCTTTGCTCACATCTCACttt                <  1:336919/52‑1 (MQ=255)
                                                                   |                                                         
GCTAACCAAACTGCATTACCTGCTTAGCCAGGAACTGGATACTGAAACCATTCGCAAGGCCATGAGCTAAAACCTGCGCGGCGAACTGACGCTTTGCTCACATCTCACTTTAATCGTGCTCACAT  >  minE/1179465‑1179589

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: