Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,406,881 A→G 100% L59L (TTA→TTG yejE → predicted oligopeptide transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,406,8810AG75.2% 41.5 / 8.9 36L59L (TTA→TTGyejEpredicted oligopeptide transporter subunit
Reads supporting (aligned to +/- strand):  ref base A (0/9);  new base G (0/27);  total (0/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TTCTGAACTTATCGCCAACGATAAACCGTTGCTGGTGCGTTATGACGGCAGTTGGTATTTCCCGTTATTGAAAAACTACAGCGAAAGCGATTTTGGCGGCCCGCTGGCAAGTCAGGCTGATTATCAGGACCCGTG  >  minE/1406815‑1406949
                                                                  |                                                                    
ttCTGAACTTATCGCCAACGATAAACCGTTGCTGGTGCGTTATGACGGCAGTTGGTATTTCCCGTTAtt                                                                    <  1:1813507/69‑1 (MQ=255)
       cTTATCGCCAACGATAAACCGTTGCTGGTGCGTTATGACGGCAGTTGGTATTTCCCGTTATTGAAAAAc                                                             <  1:2736051/69‑1 (MQ=255)
        ttATCGCCAACGATAAACCGTTGCTGGTGCGTTATGACGGCAGTTGGTATTTCCCGTTATTGAAAAAc                                                             <  1:2176484/68‑1 (MQ=255)
            cGCCAACGATAAACCGTTGCTGGTGCGTTATGACGGCAGTTGGTATTTCCCGTTGTTGAAAAACTAc                                                          <  1:1241920/67‑1 (MQ=255)
            cGCCAACGATAAACCGTTGCTGGTGCGTTATGACGGCAGTTGGTATTTCCCGTTGTTGAAAAACTAc                                                          <  1:905139/67‑1 (MQ=255)
            cGCCAACGATAAACCGTTGCTGGTGCGTTATGACGGCAGTTGGTATTTCCCGTTGTTGAAAAACTAc                                                          <  1:853573/67‑1 (MQ=255)
            cGCCAACGATAAACCGTTGCTGGTGCGTTATGACGGCAGTTGGTATTTCCCGTTGTTGAAAAACTAc                                                          <  1:659918/67‑1 (MQ=255)
            cGCCAACGATAAACCGTTGCTGGTGCGTTATGACGGCAGTTGGTATTTCCCGTTGTTGAAAAACTAc                                                          <  1:605984/67‑1 (MQ=255)
            cGCCAACGATAAACCGTTGCTGGTGCGTTATGACGGCAGTTGGTATTTCCCGTTGTTGAAAAACTAc                                                          <  1:419519/67‑1 (MQ=255)
            cGCCAACGATAAACCGTTGCTGGTGCGTTATGACGGCAGTTGGTATTTCCCGTTGTTGAAAAACTAc                                                          <  1:408433/67‑1 (MQ=255)
            cGCCAACGATAAACCGTTGCTGGTGCGTTATGACGGCAGTTGGTATTTCCCGTTGTTGAAAAACTAc                                                          <  1:3345775/67‑1 (MQ=255)
            cGCCAACGATAAACCGTTGCTGGTGCGTTATGACGGCAGTTGGTATTTCCCGTTGTTGAAAAACTAc                                                          <  1:2913328/67‑1 (MQ=255)
            cGCCAACGATAAACCGTTGCTGGTGCGTTATGACGGCAGTTGGTATTTCCCGTTGTTGAAAAACTAc                                                          <  1:2744843/67‑1 (MQ=255)
            cGCCAACGATAAACCGTTGCTGGTGCGTTATGACGGCAGTTGGTATTTCCCGTTGTTGAAAAACTAc                                                          <  1:262744/67‑1 (MQ=255)
            cGCCAACGATAAACCGTTGCTGGTGCGTTATGACGGCAGTTGGTATTTCCCGTTGTTGAAAAACTAc                                                          <  1:2563216/67‑1 (MQ=255)
            cGCCAACGATAAACCGTTGCTGGTGCGTTATGACGGCAGTTGGTATTTCCCGTTGTTGAAAAACTAc                                                          <  1:222812/67‑1 (MQ=255)
            cGCCAACGATAAACCGTTGCTGGTGCGTTATGACGGCAGTTGGTATTTCCCGTTGTTGAAAAACTAc                                                          <  1:1998399/67‑1 (MQ=255)
            cGCCAACGATAAACCGTTGCTGGTGCGTTATGACGGCAGTTGGTATTTCCCGTTGTTGAAAAACTAc                                                          <  1:1868461/67‑1 (MQ=255)
            cGCCAACGATAAACCGTTGCTGGTGCGTTATGACGGCAGTTGGTATTTCCCGTTGTTGAAAAACTAc                                                          <  1:1355639/67‑1 (MQ=255)
            cGCCAACGATAAACCGTTGCTGGTGCGTTATGACGGCAGTTGGTATTTCCCGTTGTTGAAAAACTAc                                                          <  1:1509133/67‑1 (MQ=255)
            cGCCAACGATAAACCGTTGCTGGTGCGTTATGACGGCAGTTGGTATTTCCCGTTGTTGAAAAACTAc                                                          <  1:1836214/67‑1 (MQ=255)
            cGCCAACGATAAACCGTTGCTGGTGCGTTATGACGGCAGTTGGTATTTCCCGTTGTTGAAAAACTAc                                                          <  1:1741943/67‑1 (MQ=255)
            cGCCAACGATAAACCGTTGCTGGTGCGTTATGACGGCAGTTGGTATTTCCCGTTGTTGAAAAACTAc                                                          <  1:1514401/67‑1 (MQ=255)
            cGCCAACGATAAACCGTTGCTGGTGCGTTATGACGGCAGTTGGTATTTCCCGTTGTTGAAAAACTAc                                                          <  1:1615768/67‑1 (MQ=255)
            cGCCAACGATAAACCGTTGCTGGTGCGTTATGACGGCAGTTGGTATTTCCCGTTGTTGAAAAACTAc                                                          <  1:1611391/67‑1 (MQ=255)
             gCCAACGATAAACCGTTGCTGGTGCGTTATGACGGCAGTTGGTATTTCCCGTTGTTGAAAAACTAc                                                          <  1:2042907/66‑1 (MQ=255)
              ccAACGATAAACCGTTGCTGGTGCGTTATGACGGCAGTTGGTATTTCCCGTTGTTGAAAAACTAc                                                          <  1:674837/65‑1 (MQ=255)
                  cGATAAACCGTTGCTGGTGCGTTATGACGGCAGTTGGTATTTCCCGTTGTTGAAAAACTAc                                                          <  1:2349009/61‑1 (MQ=255)
                         ccGTTGCTGGTGCGTTATGACGGCAGTTGGTATTTCCCGTTGTTGAAAAACTAc                                                          <  1:831977/54‑1 (MQ=255)
                                   tGCGTTATGACGGCAGTTGGTATTTCCCGTTGTTGAAAAACTAc                                                          <  1:2798273/44‑1 (MQ=255)
                                                 aGTTGGTATTTCCCGTTATTGAAAAACTACAGCGAAAGCGATTTTGGCGGCCCGCTGGCAAGTCAGGCt                   <  1:1659718/69‑1 (MQ=255)
                                                 aGTTGGTATTTCCCGTTATTGAAAAACTACAGCGAAAGCGATTTTGGCGGCCCGCTGGCAAGTCAGGCt                   <  1:1877835/69‑1 (MQ=255)
                                                     ggTATTTCCCGTTATTGAAAAACTACAGCGAAAGCGATTTTGGCGGCCCGCTGGCAAGTCAg                      <  1:1267217/62‑1 (MQ=255)
                                                     ggTATTTCCCGTTATTGAAAAACTACAGCGAAAGCGATTTTGGCGGCCCGCTGGCAAGTCAg                      <  1:1182989/62‑1 (MQ=255)
                                                     ggTATTTCCCGTTATTGAAAAACTACAGCGAAAGCGATTTTGGCGGCCCGCTGGCAAGTCAGGCTGAtt               <  1:2491420/69‑1 (MQ=255)
                                                        aTTTCCCGTTATTGAAAAACTACAGCGAAAGCGATTTTGGCGGCCCGCTGGCAAGTCAGGCTGATTATc            <  1:264945/69‑1 (MQ=255)
                                                                  aTTGAAAAACTACAGCGAAAGCGATTTTGGCGGCCCGCTGGCAAGTCAGGCTGATTATCAGGACCCGt   >  1:2650434/1‑68 (MQ=255)
                                                                  aTTGAAAAACTACAGCGAAAGCGATTTTGGCGGCCCGCTGGCAAGTCAGGCTGATTATCAGGACCCGTg  >  1:3093302/1‑69 (MQ=255)
                                                                  |                                                                    
TTCTGAACTTATCGCCAACGATAAACCGTTGCTGGTGCGTTATGACGGCAGTTGGTATTTCCCGTTATTGAAAAACTACAGCGAAAGCGATTTTGGCGGCCCGCTGGCAAGTCAGGCTGATTATCAGGACCCGTG  >  minE/1406815‑1406949

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: