Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,633,585 T→C 59.7% N201N (AAT→AAC livG → leucine/isoleucine/valine transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,633,5850TC59.7% 21.4 / 33.0 52N201N (AAT→AAClivGleucine/isoleucine/valine transporter subunit
Reads supporting (aligned to +/- strand):  ref base T (8/13);  new base C (22/9);  total (30/22)
Fisher's exact test for biased strand distribution p-value = 2.45e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCTTAACCCGAAAGAGACGAAAGAGCTGGATGAGCTGATTGCCGAACTGCGCAATCATCACAACACCACTATCTTGTTGATTGAACACGATATGAAGCTGGTGATGGGAATTTCGGACCGA  >  minE/2633531‑2633651
                                                      |                                                                  
tCTTAACCCGAAAGAGACGAAAGAGCTGGATGAGCTGATTGCCGAACTGCGCAATCATCACAAcaccac                                                      >  1:2710397/1‑69 (MQ=255)
 cTTAACCCGAAAGAGACGAAAGAGCTGGATGAGCTGATTGCCGAACTGCGCAatca                                                                  >  1:2462446/1‑56 (MQ=255)
 cTTAACCCGAAAGAGACGAAAGAGCTGGATGAGCTGATTGCCGAACTGCGCAATCATCACAACACCACt                                                     >  1:1641603/1‑69 (MQ=255)
 cTTAACCCGAAAGAGACGAAAGAGCTGGATGAGCTGATTGCCGAACTGCGCAATCATCACAACACCACt                                                     >  1:1167703/1‑69 (MQ=255)
      cccGAAAGAGACGAAAGAGCTGGATGAGCTGATTGCCGAACTGCGCAATCATCACAACACCACTATc                                                  <  1:3233050/67‑1 (MQ=255)
         gAAAGAGACGAAAGAGCTGGATGAGCTGATTGCCGAGCTGCGTAACCATCACAACACCAATATCttgtt                                             >  1:3113578/1‑69 (MQ=255)
         gAAAGAGACGAAAGAGCTGGATGAGCTGATTGCCGAACTGCGTAACCATCACAACACCACt                                                     >  1:2558345/1‑61 (MQ=255)
         gAAAGAGACGAAAGAGCTGGATGAGCTGATTGCCGAACTGCGTAACCATCACAACACCACTATCttgtt                                             >  1:3153607/1‑69 (MQ=255)
         gAAAGAGACGAAAGAGCTGGATGAGCTGATTGCCGAACTGCGTAACCATCACAACACCACTATCttgtt                                             >  1:921681/1‑69 (MQ=255)
         gAAAGAGACGAAAGAGCTGGATGAGCTGATTGCCGAACTGCGTAACCATCACAACACCACTATCttgtt                                             >  1:2567161/1‑69 (MQ=255)
         gAAAGAGACGAAAGAGCTGGATGAGCTGATTGCCGAACTGCGTAACCATCACAACACCACTATCttgtt                                             >  1:2580874/1‑69 (MQ=255)
         gAAAGAGACGAAAGAGCTGGATGAGCTGATTGCCGAACTGCGTAACCATCACAACACCACTATCttgtt                                             >  1:2272595/1‑69 (MQ=255)
         gAAAGAGACGAAAGAGCTGGATGAGCTGATTGCCGAACTGCGTAACCATCACAACACCACTATCttgtt                                             >  1:2684962/1‑69 (MQ=255)
         gAAAGAGACGAAAGAGCTGGATGAGCTGATTGCCGAACTGCGTAACCATCACAACACCACTATCttgtt                                             >  1:1148069/1‑69 (MQ=255)
         gAAAGAGACGAAAGAGCTGGATGAGCTGATTGCCGAACTGCGTAACCATCACAACACCACTATCttgtt                                             >  1:2188994/1‑69 (MQ=255)
         gAAAGAGACGAAAGAGCTGGATGAGCTGATTGCCGAACTGCGTAACCATCACAACACCACTATCttgtt                                             >  1:2804912/1‑69 (MQ=255)
         gAAAGAGACGAAAGAGCTGGATGAGCTGATTGCCGAACTGCGTAACCATCACAACACCACTATCttgtt                                             >  1:2171162/1‑69 (MQ=255)
         gAAAGAGACGAAAGAGCTGGATGAGCTGATTGCCGAACTGCGTAACCATCACAACACCACTATCttgtt                                             >  1:2153632/1‑69 (MQ=255)
         gAAAGAGACGAAAGAGCTGGATGAGCTGATTGCCGAACTGCGTAACCATCACAACACCACTATCttgtt                                             >  1:2822281/1‑69 (MQ=255)
         gAAAGAGACGAAAGAGCTGGATGAGCTGATTGCCGAACTGCGTAACCATCACAACACCACTATCttgtt                                             >  1:1905807/1‑69 (MQ=255)
         gAAAGAGACGAAAGAGCTGGATGAGCTGATTGCCGAACTGCGTAACCATCACAACACCACTATCttgtt                                             >  1:3231047/1‑69 (MQ=255)
         gAAAGAGACGAAAGAGCTGGATGAGCTGATTGCCGAACTGCGTAACCATCACAACACCACTATCttgtt                                             >  1:2909637/1‑69 (MQ=255)
         gAAAGAGACGAAAGAGCTGGATGAGCTGATTGCCGAACTGCGTAACCATCACAACACCACTATCttgtt                                             >  1:1507594/1‑69 (MQ=255)
         gAAAGAGACGAAAGAGCTGGATGAGCTGATTGCCGAACTGCGTAACCATCACAACACCACTATCttgtt                                             >  1:2972294/1‑69 (MQ=255)
         gAAAGAGACGAAAGAGCTGGATGAGCTGATTGCCGAACTGCGTAACCATCACAACACCACTATCttgtt                                             >  1:1449848/1‑69 (MQ=255)
         gAAAGAGACGAAAGAGCTGGATGAGCTGATTGCCGAACTGCGTAACCATCACAACACCACTATCttgtt                                             >  1:1398398/1‑69 (MQ=255)
         gAAAGAGACGAAAGAGCTGGATGAGCTGATTGCCGAACTGCGTAACCATCACAACACCACTATCttgat                                             >  1:2481566/1‑67 (MQ=255)
          aaaGAGACGAAAGAGCTGGATGAGCTGATTGCCGAACTGCGCAATCATCACAACACCACTATCttgttg                                            <  1:2321837/69‑1 (MQ=255)
          aaaGAGACGAAAGAGCTGGATGAGCTGATTGCCGAACTGCGCAATCATCACAACACCACTATCttgttg                                            <  1:1957926/69‑1 (MQ=255)
          aaaGAGACGAAAGAGCTGGATGAGCTGATTGCCGAACTGCGCAATCATCACAACACCACTATCttgttg                                            <  1:37593/69‑1 (MQ=255)
                     agagCTGGATGAGCTGATTGCCGAACTGCGCAATCATCACAACACCACTATCTTGTTGATTGAACACGa                                 <  1:792878/69‑1 (MQ=255)
                     agagCTGGATGAGCTGATTGCCGAACTGCGCAATCATCACAACACCACTATCTTGTTGATTGAACACGa                                 <  1:2996742/69‑1 (MQ=255)
                     agagCTGGATGAGCTGATTGCCGAACTGCGCAATCATCACAACACCACTATCTTGTTGATTGAACACGa                                 <  1:3069994/69‑1 (MQ=255)
                     agagCTGGATGAGCTGATTGCCGAACTGCGCAATCATCACAACACCACTATCTTGTTGATTGAACACGa                                 <  1:3181177/69‑1 (MQ=255)
                      gagCTGGATGAGCTGATTGCCGAACTGCGCAATCATCACAACACCACTATCTTGTTGATTGAACACGa                                 <  1:226850/68‑1 (MQ=255)
                         cTGGATGAGCTGATTGCCGAACTGCGCAATCATCACAACACCACTATCTTGTTGATTGAACACGa                                 <  1:1163530/65‑1 (MQ=255)
                                aGCTGATTGCCGAACTGCGCAATCATCACAACACCACTATCTTGTTGATTGAACACGa                                 >  1:1498751/1‑58 (MQ=255)
                                  cTGATTGCCGAACTGCGCAATCATCACAACACCACTATCTTGTTTATTGAACACGATATGAAGCTGGt                     <  1:2786019/68‑1 (MQ=255)
                                    gATTGCCGAACTGCGTAACCATCACAACACCACTATCttg                                               <  1:1071323/40‑1 (MQ=38)
                                    gATTGCCGAACTGCGTAACCATCACAACACCACTATCttg                                               <  1:1034024/40‑1 (MQ=38)
                                    gATTGCCGAACTGCGTAACCATCACAACACCACTATCttg                                               <  1:841884/40‑1 (MQ=38)
                                    gATTGCCGAACTGCGTAACCATCACAACACCACTATCttg                                               <  1:3305591/40‑1 (MQ=38)
                                    gATTGCCGAACTGCGTAACCATCACAACACCACTATCttg                                               <  1:3285861/40‑1 (MQ=38)
                                    gATTGCCGAACTGCGTAACCATCACAACACCACTATCttg                                               <  1:313241/40‑1 (MQ=38)
                                    gATTGCCGAACTGCGTAACCATCACAACACCACTATCttg                                               <  1:3025148/40‑1 (MQ=38)
                                    gATTGCCGAACTGCGTAACCATCACAACACCACTATCttg                                               <  1:1595995/40‑1 (MQ=38)
                                    gATTGCCGAACTGCGTAACCATCACAACACCACTATCttg                                               <  1:2182722/40‑1 (MQ=38)
                                       tGCCGAACTGCGCAATCATCACAACACCACTATCTTGTTGATTGaa                                      >  1:1199760/1‑46 (MQ=255)
                                       tGCCGAACTGCGCAATCATCACAACACCACTATCTTGTTGATTGAACACGATATGAAGCTGGTGATggg               >  1:449374/1‑69 (MQ=255)
                                       tGCCGAACTGCGCAATCATCACAACACCACTATCTTGTTGATTGAACACGATATGAAGCTGGTGATggg               >  1:2267485/1‑69 (MQ=255)
                                       tGCCGAACTGCGCAATCATCACAACACCACTATCTTGTTGATTGAACACGATATGAAGCTGGTGATggg               >  1:2284006/1‑69 (MQ=255)
                                            aaCTGCGCAATCATCACAACACCACTATCTTGTTGATTGAACACGATATGAAGCTGGTGATGGGAAttt          <  1:1413152/69‑1 (MQ=255)
                                                    aaTCATCACAACACCACTATCTTGTTGATTGAACACGATATGAAGCTGGTGATGGGAATTTCGGACCGa  <  1:1102617/69‑1 (MQ=255)
                                                      |                                                                  
TCTTAACCCGAAAGAGACGAAAGAGCTGGATGAGCTGATTGCCGAACTGCGCAATCATCACAACACCACTATCTTGTTGATTGAACACGATATGAAGCTGGTGATGGGAATTTCGGACCGA  >  minE/2633531‑2633651

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: