Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,586,707 G→A 63.6% intergenic (+269/‑115) cysK → / → ptsH cysteine synthase A, O‑acetylserine sulfhydrolase A subunit/phosphohistidinoprotein‑hexose phosphotransferase component of PTS system (Hpr)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,586,7070GA63.6% 42.4 / 44.9 74intergenic (+269/‑115)cysK/ptsHcysteine synthase A, O‑acetylserine sulfhydrolase A subunit/phosphohistidinoprotein‑hexose phosphotransferase component of PTS system (Hpr)
Reads supporting (aligned to +/- strand):  ref base G (12/15);  new base A (25/22);  total (37/37)
Fisher's exact test for biased strand distribution p-value = 6.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.40e-01

ATGTTTAATGACGTACGAAACGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTGCGCTCTTCGTGCGTCGCGTCTGTTAAAAACTGGCGCTAACAATACAGGCTAAAGTCGAACCGCC  >  minE/1586650‑1586771
                                                         |                                                                
aTGTTTAATGACGTACGAAACGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTGCGCTCTTCGTg                                                       >  1:416656/1‑69 (MQ=255)
  gTTTAATGACGTACGAAACGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTGCGCTCTTCGTgcg                                                     >  1:3108790/1‑69 (MQ=255)
  gTTTAATGACGTACGAAACGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTGCGCTCTTCGTgcg                                                     >  1:387580/1‑69 (MQ=255)
        tGACGTACGAAACGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTGCGCTCTTCGTgcgtcgcg                                                <  1:2083597/68‑1 (MQ=255)
         gACGTACGAAACGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTGCGCTCTTCGTgcgtcgcg                                                <  1:103124/67‑1 (MQ=255)
            gTACGAAACGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTGCGCTCTTCGTGCGTCGCGTCt                                             >  1:1459101/1‑67 (MQ=255)
            gTACGAAACGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTGCGCTCTTCGTGCGTCGCGTCTg                                            >  1:358674/1‑68 (MQ=255)
            gTACGAAACGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTGCGCTCTTCGTGCGTCGCGTCTg                                            >  1:3074572/1‑68 (MQ=255)
            gTACGAAACGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTGCGCTCTTCGTGCGTCGCGTCTg                                            >  1:1514797/1‑68 (MQ=255)
            gTACGAAACGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTGCGCTCTTCGTGCGTCGCGTCTg                                            >  1:3246262/1‑68 (MQ=255)
              aCGAAACGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTACGCTCTTCGTgcgtc                                                   <  1:328697/59‑1 (MQ=255)
              aCGAAACGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTACGCTCTTCGTgcgtc                                                   <  1:377120/59‑1 (MQ=255)
              aCGAAACGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTACGCTCTTCGTgcgtc                                                   <  1:1922685/59‑1 (MQ=255)
              aCGAAACGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTACGCTCTTCGTgcgtc                                                   <  1:3428910/59‑1 (MQ=255)
              aCGAAACGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTACGCTCTTCGTgcgtc                                                   <  1:3337771/59‑1 (MQ=255)
              aCGAAACGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTACGCTCTTCGTgcgtc                                                   <  1:2067389/59‑1 (MQ=255)
              aCGAAACGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTACGCTCTTCGTgcgtc                                                   <  1:2692853/59‑1 (MQ=255)
              aCGAAACGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTACGCTCTTCGTgcgtc                                                   <  1:3197262/59‑1 (MQ=255)
              aCGAAACGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTACGCTCTTCGTgcgtc                                                   <  1:3136994/59‑1 (MQ=255)
              aCGAAACGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTACGCTCTTCGTgcgtc                                                   <  1:2256393/59‑1 (MQ=255)
              aCGAAACGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTACGCTCTTCGTgcgtc                                                   <  1:2828849/59‑1 (MQ=255)
              aCGAAACGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTACGCTCTTCGTgcgtc                                                   <  1:2472048/59‑1 (MQ=255)
              aCGAAACGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTACGCTCTTCGTgcgtc                                                   <  1:389225/59‑1 (MQ=255)
              aCGAAACGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTACGCTCTTCGTgcgtc                                                   <  1:656873/59‑1 (MQ=255)
              aCGAAACGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTACGCTCTTCGTgcgtc                                                   <  1:577124/59‑1 (MQ=255)
              aCGAAACGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTACGCTCTTCGTgcgtc                                                   <  1:1284943/59‑1 (MQ=255)
              aCGAAACGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTACGCTCTTCGTgcgtc                                                   <  1:1342821/59‑1 (MQ=255)
              aCGAAACGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTACGCTCTTCGTgcgtc                                                   <  1:1108709/59‑1 (MQ=255)
              aCGAAACGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTACGCTCTTCGTgcgtc                                                   <  1:1101493/59‑1 (MQ=255)
               cGAAACGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTACGCTCTTCGTgcgtc                                                   <  1:969566/58‑1 (MQ=255)
                 aaaCGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTGCGCTCTTCGTGCGTCGCGTCTGTTaaa                                       <  1:424200/68‑1 (MQ=255)
                 aaaCGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTGCGCTCTTCGTGCGTCGCGTCTGTTaaa                                       <  1:1100082/68‑1 (MQ=255)
                 aaaCGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTGCGCTCTTCGTGCGTCGCGTCTGTTaaa                                       <  1:423056/68‑1 (MQ=255)
                 aaaCGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTGCGCTCTTCGTGCGTCGCGTCTGTTaaa                                       <  1:2212781/68‑1 (MQ=255)
                 aaaCGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTGCGCTCTTCGTGCGTCGCGTCTGTTaaa                                       <  1:864842/68‑1 (MQ=255)
                  aaCGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTACGCTCtt                                                           >  1:499536/1‑47 (MQ=255)
                  aaCGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTACGCTCTTCGTGCGTCGCGTCTGTTaaaaa                                     >  1:55511/1‑69 (MQ=255)
                  aaCGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTACGCTCTTCGTGCGTCGCGTCTGTTaaaaa                                     >  1:1055180/1‑69 (MQ=255)
                  aaCGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTACGCTCTTCGTGCGTCGCGTCTGTTaaaaa                                     >  1:3062674/1‑69 (MQ=255)
                  aaCGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTACGCTCTTCGTGCGTCGCGTCTGTTaaaaa                                     >  1:3033407/1‑69 (MQ=255)
                  aaCGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTACGCTCTTCGTGCGTCGCGTCTGTTaaaaa                                     >  1:871688/1‑69 (MQ=255)
                  aaCGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTACGCTCTTCGTGCGTCGCGTCTGTTaaaaa                                     >  1:881522/1‑69 (MQ=255)
                  aaCGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTACGCTCTTCGTGCGTCGCGTCTGTTaaaaa                                     >  1:1236706/1‑69 (MQ=255)
                  aaCGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTACGCTCTTCGTGCGTCGCGTCTGTTaaaaa                                     >  1:1679910/1‑69 (MQ=255)
                  aaCGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTACGCTCTTCGTGCGTCGCGTCTGTTaaaaa                                     >  1:1613176/1‑69 (MQ=255)
                  aaCGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTACGCTCTTCGTGCGTCGCGTCTGTTaaaaa                                     >  1:1525762/1‑69 (MQ=255)
                  aaCGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTACGCTCTTCGTGCGTCGCGTCTGTTaaaaa                                     >  1:1442605/1‑69 (MQ=255)
                  aaCGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTACGCTCTTCGTGCGTCGCGTCTGTTaaaaa                                     >  1:2009972/1‑69 (MQ=255)
                  aaCGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTACGCTCTTCGTGCGTCGCGTCTGTTaaaaa                                     >  1:2069607/1‑69 (MQ=255)
                  aaCGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTACGCTCTTCGTGCGTCGCGTCTGTTaaaaa                                     >  1:212808/1‑69 (MQ=255)
                  aaCGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTACGCTCTTCGTGCGTCGCGTCTGTTaaaaa                                     >  1:1290268/1‑69 (MQ=255)
                  aaCGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTACGCTCTTCGTGCGTCGCGTCTGTTaaaaa                                     >  1:1216071/1‑69 (MQ=255)
                  aaCGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTACGCTCTTCGTGCGTCGCGTCTGTTaaaaa                                     >  1:2272378/1‑69 (MQ=255)
                  aaCGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTACGCTCTTCGTGCGTCGCGTCTGTTaaaaa                                     >  1:2367004/1‑69 (MQ=255)
                  aaCGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTACGCTCTTCGTGCGTCGCGTCTGTTaaaaa                                     >  1:1199777/1‑69 (MQ=255)
                  aaCGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTACGCTCTTCGTGCGTCGCGTCTGTTaaaaa                                     >  1:1156537/1‑69 (MQ=255)
                  aaCGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTACGCTCTTCGTGCGTCGCGTCTGTTaaaaa                                     >  1:2410741/1‑69 (MQ=255)
                  aaCGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTACGCTCTTCGTGCGTCGCGTCTGTTaa                                        >  1:11218/1‑66 (MQ=255)
                  aaCGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATCACGCTCTTCGTGCGTCGCGTCTGTTaaaaa                                     >  1:982544/1‑69 (MQ=255)
                  aaCGTCAGAGGTCAACACCCGCCAGCAATGGACTGTATTACGCTCTTCGTGCGTCGCGTCTGTTaaaaa                                     >  1:1416272/1‑69 (MQ=255)
                        aGCGGTCAACACCCGCCAGCAATGGACTGTATTgcg                                                                >  1:1900253/1‑36 (MQ=255)
                        aGCGGTCAACACCCGCCAGCAATGGACTGTATTgcg                                                                >  1:849240/1‑36 (MQ=255)
                        aGCGGTCAACACCCGCCAGCAATGGACTGTATTgcg                                                                >  1:318570/1‑36 (MQ=255)
                        aGCGGTCAACACCCGCCAGCAATGGACTGTATTgcg                                                                >  1:2938310/1‑36 (MQ=255)
                          cGGTCAACACCCGCCAGCAATGGACTGTATTACGCTCTTCGTgcgtc                                                   <  1:209580/47‑1 (MQ=255)
                             tCAACACCCGCCAGCAATGGACTGTATTGCGCTCTTCGTGCGTCGCGTCTGTTAAAAACTGGCGCTAAc                          <  1:1967341/69‑1 (MQ=255)
                             tCAACACCCGCCAGCAATGGACTGTATTGCGCTCTTCGTGCGTCGCGTCTGTTAAAAACTGGCGCTAAc                          <  1:1465320/69‑1 (MQ=255)
                             tCAACACCCGCCAGCAATGGACTGTATTGCGCTCTTCGTGCGTCGCGTCTGTTAAAAACTGGCGCTAAc                          <  1:2188786/69‑1 (MQ=255)
                             tCAACACCCGCCAGCAATGGACTGTATTGCGCTCTTCGTGCGTCGCGTCTGTTAAAAACTGGCGCTAAc                          <  1:2258622/69‑1 (MQ=255)
                             tCAACACCCGCCAGCAATGGACTGTATTGCGCTCTTCGTGCGTCGCGTCTGTTAAAAACTGGCGCTAAc                          <  1:678065/69‑1 (MQ=255)
                             tCAACACCCGCCAGCAATGGACTGTATTGCGCTCTTCGTGCGTCGCGTCTGTTAAAAACTGGCGCTAAc                          <  1:2402726/69‑1 (MQ=255)
                             tCAACACCCGCCAGCAATGGACTGTATTGCGCTCTTCGTGCGTCGCGTCTGTTAAAAACTGGCGCTAAc                          <  1:2550003/69‑1 (MQ=255)
                               aaCACCCGCCAGCAATGGACTGTATTGCGCTCTTCGTGCGTCGCGTCTGTTAAAAACTGGCGCTAACaa                        <  1:3338921/69‑1 (MQ=255)
                                   cccGACAGCAATGGACTGTATTACGCTCTTCGTgcgtc                                                   <  1:1417323/38‑1 (MQ=38)
                                            aaTGGACTGTATTGCGCTCTTCGTGCGTCGCGTCTGTTAAAAACTGGCGCTAACAATACAGGCTAAAg            >  1:2995773/1‑68 (MQ=255)
                                                     tATTGCGCTCTTCGTGCGTCGCGTCTGTTAAAAACTGGCGCTAACAATACAGGCTAAAGTCGAACCGcc  >  1:2557817/1‑69 (MQ=255)
                                                     tATTGCGCTCTTCGTGCGTCGCGTCTGTTAAAAACTGGCGCTAACAATACAGGCTAAAGTCGAACCGcc  >  1:2831914/1‑69 (MQ=255)
                                                     tATTGCGCTCTTCGTGCGTCGCGTCTGTTAAAAACTGGCGCTAACAATACAGGCTAAAGTCGAACCGcc  >  1:1115758/1‑69 (MQ=255)
                                                         |                                                                
ATGTTTAATGACGTACGAAACGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTGCGCTCTTCGTGCGTCGCGTCTGTTAAAAACTGGCGCTAACAATACAGGCTAAAGTCGAACCGCC  >  minE/1586650‑1586771

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: