Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 910,937 C→T 25.1% G59G (GGC→GGT ycjV → predicted sugar transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE910,9370CT25.1% 118.2 / 32.8 95G59G (GGC→GGTycjVpredicted sugar transporter subunit
Reads supporting (aligned to +/- strand):  ref base C (37/34);  new base T (18/6);  total (55/40)
Fisher's exact test for biased strand distribution p-value = 5.86e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

TCGACCACCCTGCGCATGATTGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGACGGCAAACGAATGAATGACGTTCCAGCCAAAGCACGCAA  >  minE/910890‑910990
                                               |                                                     
tCGACCACCCTGCGCATGATTGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGACGGCaa                                   >  1:815931/1‑68 (MQ=255)
tCGACCACCCTGCGCATGATTGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGACGGCaa                                   >  1:2076118/1‑68 (MQ=255)
tCGACCACCCTGCGCATGATTGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGACGGCaa                                   >  1:2588544/1‑68 (MQ=255)
tCGACCACCCTGCGCATGATTGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGACGGCaa                                   >  1:3096705/1‑68 (MQ=255)
  gACCACCCTGCGCATGATTGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGACGGCAAACg                                <  1:1244326/69‑1 (MQ=255)
  gACCACCCTGCGCATGATTGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGACGGCAAACg                                <  1:1286227/69‑1 (MQ=255)
  gACCACCCTGCGCATGATTGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGACGGCAAACg                                <  1:3428678/69‑1 (MQ=255)
  gACCACCCTGCGCATGATTGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGACGGCAAACg                                <  1:424757/69‑1 (MQ=255)
  gACCACCCTGCGCATGATTGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGACGGCAAACg                                <  1:347193/69‑1 (MQ=255)
  gACCACCCTGCGCATGATTGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGACGGCAAACg                                <  1:424207/69‑1 (MQ=255)
  gACCACCCTGCGCATGATTGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGACGGCAAACg                                <  1:422827/69‑1 (MQ=255)
     caccCTGCGCATGATTGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGACGGCAAACgaa                              >  1:2126451/1‑68 (MQ=255)
     caccCTGCGCATGATTGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGACGGCAAACgaa                              >  1:216658/1‑68 (MQ=255)
     caccCTGCGCATGATTGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGACGGCAAACgaa                              >  1:3317872/1‑68 (MQ=255)
     caccCTGCGCATGATTGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGACGGCAAACgaa                              >  1:3415156/1‑68 (MQ=255)
     caccCTGCGCATGATTGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGACGGCAAACgaa                              >  1:1559758/1‑68 (MQ=255)
     caccCTGCGCATGATTGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGACGGCAAACgaa                              >  1:139455/1‑68 (MQ=255)
     caccCTGCGCATGATTGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGACGGCAAACgaa                              >  1:1121471/1‑68 (MQ=255)
     caccATGCGCATGATTGACGGGATTGAGGAGATCAGAGGCGGCGATCTGTTGATCGACGGCAAACgaa                              >  1:3037827/1‑68 (MQ=255)
      accCTGCGCATGATTGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGACGGCAAACgaatg                            <  1:3013393/69‑1 (MQ=255)
       cccTGCGCATGATTGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGACGGCAAACgaatg                            <  1:1430811/68‑1 (MQ=255)
       cccTGCGCATGATTGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGACGGCAAACgaatg                            <  1:2397327/68‑1 (MQ=255)
            cgtaTGATTGCCGGGCTTGAGGAGATCAGCGGCGGTGATCTGTTGATCGACGGCAAACgaatgaat                         >  1:2941780/4‑66 (MQ=255)
            cgtaTGATTGCCGGGCTTGAGGAGATCAGCGGCGGTGATCTGTTGATCGACGGCAAACgaatgaat                         >  1:2347645/4‑66 (MQ=255)
            cgtaTGATTGCCGGGCTTGAGGAGATCAGCGGCGGTGATCTGTTGATCGACGGCAAACGAATGAATGAc                      >  1:325250/4‑69 (MQ=255)
            cgtaTGATTGCCGGGCTTGAGGAGATCAGCGGCGGTGATCTGTTGATCGACGGCAAACGAATGAATGAc                      >  1:611504/4‑69 (MQ=255)
            cgtaTGATTGCCGGGCTTGAGGAGATCAGCGGCGGTGATCTGTTGATCGACGGCAAACGAATGAATGAc                      >  1:3421893/4‑69 (MQ=255)
            cgtaTGATTGCCGGGCTTGAGGAGATCAGCGGCGGTGATCTGTTGATCGACGGCAAACGAATGAATGAc                      >  1:2975418/4‑69 (MQ=255)
            cgtaTGATTGCCGGGCTTGAGGAGATCAGCGGCGGTGATCTGTTGATCGACGGCAAACGAATGAATGAc                      >  1:2784923/4‑69 (MQ=255)
            cgtaTGATTGCCGGGCTTGAGGAGATCAGCGGCGGTGATCTGTTGATCGACGGCAAACGAATGAATGAc                      >  1:71076/4‑69 (MQ=255)
            cgtaTGATTGCCGGGCTTGAGGAGATCAGCGGCGGTGATCTGTTGATCGACGGCAAACGAATGAATGAc                      >  1:596707/4‑69 (MQ=255)
            cgtaTGATTGCCGGGCTTGAGGAGATCAGCGGCGGTGATCTGTTGATCGACGGCAAACGAATGAATGAc                      >  1:724218/4‑69 (MQ=255)
            cgtaTGATTGCCGGGCTTGAGGAGATCAGCGGCGGTGATCTGTTGATCGACGGCAAACGAATGAATGAc                      >  1:976347/4‑69 (MQ=255)
            cgtaTGATTGCCGGGCTTGAGGAGATCAGCGGCGGTGATCTGTTGATCGACGGCAAACGAATGAATGAc                      >  1:2662114/4‑69 (MQ=255)
            cgtaTGATTGCCGGGCTTGAGGAGATCAGCGGCGGTGATCTGTTGATCGACGGCAAACGAATGAATGAc                      >  1:1031854/4‑69 (MQ=255)
            cgtaTGATTGCCGGGCTTGAGGAGATCAGCGGCGGTGATCTGTTGATCGACGGCAAACGAATGAATGAc                      >  1:1071633/4‑69 (MQ=255)
            cgtaTGATTGCCGGGCTTGAGGAGATCAGCGGCGGTGATCTGTTGATCGACGGCAAACGAATGAATGAc                      >  1:1193557/4‑69 (MQ=255)
            cgtaTGATTGCCGGGCTTGAGGAGATCAGCGGCGGTGATCTGTTGATCGACGGCAAACGAATGAATGAc                      >  1:1720163/4‑69 (MQ=255)
            cgtaTGATTGCCGGGCTTGAGGAGATCAGCGGCGGTGATCTGTTGATCGACGGCAAACGAATGAATGAc                      >  1:1947525/4‑69 (MQ=255)
            cgtaTGATTGCCGGGCTTGAGGAGATCAGCGGCGGTGATCTGTTGATCGACGGCAAACGAATGAATGAc                      >  1:195305/4‑69 (MQ=255)
                  aTTGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGACGGCAAACGAatgaatga                       >  1:3325030/1‑62 (MQ=255)
                  aTTGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGACGGCAAACGAatgaatga                       >  1:581033/1‑62 (MQ=255)
                  aTTGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGACGGCAAACGAatgaatga                       >  1:484350/1‑62 (MQ=255)
                  aTTGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGACGGCAAACGAatgaatga                       >  1:412698/1‑62 (MQ=255)
                  aTTGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGACGGCAAACGAatgaatga                       >  1:2119095/1‑62 (MQ=255)
                  aTTGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGACGGCAAACGAatgaatga                       >  1:3113340/1‑62 (MQ=255)
                  aTTGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGACGGCAAACGAatgaatga                       >  1:3089588/1‑62 (MQ=255)
                  aTTGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGACGGCAAACGAatgaatga                       >  1:3049715/1‑62 (MQ=255)
                  aTTGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGACGGCAAACGAatgaatga                       >  1:1976214/1‑62 (MQ=255)
                  aTTGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGACGGCAAACGAatgaatga                       >  1:2087138/1‑62 (MQ=255)
                  aTTGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGACGGCAAACGAatgaatga                       >  1:2274300/1‑62 (MQ=255)
                  aTTGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGACGGCAAACGAatgaatga                       >  1:2554427/1‑62 (MQ=255)
                  aTTGACGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGACGGCAAACGAatgaatga                       >  1:2676724/1‑62 (MQ=255)
                    tGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGACGGCAAACTAATGAATGACGTTCCAg               >  1:2729791/1‑68 (MQ=255)
                    tGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGACGGCAAACGAATGAATGACGTTCCAg               >  1:1942916/1‑68 (MQ=255)
                    tGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGACGGCAAACGAATGAATGACGTTCCAg               >  1:429281/1‑68 (MQ=255)
                    tGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGACGGCAAACGAATGAATGACGTTCCAg               >  1:591024/1‑68 (MQ=255)
                    tGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGACGGCAAACGAATGAATGACGTTCCAg               >  1:2933323/1‑68 (MQ=255)
                    tGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGACGGCAAACGAATGAATGACGTTCCAg               >  1:707937/1‑68 (MQ=255)
                    tGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGACGGCAAACGAATGAATGACGTTCCAg               >  1:724063/1‑68 (MQ=255)
                    tGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGACGGCAAACGAATGAATGACGTTCCAg               >  1:734456/1‑68 (MQ=255)
                    tGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGACGGCAAACGAATGAATGACGTTCCAg               >  1:741906/1‑68 (MQ=255)
                    tGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGACGGCAAACGAATGAATGACGTTCCAg               >  1:2008045/1‑68 (MQ=255)
                    tGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGACGGCAAACGAATGAATGACGTTCCAg               >  1:832555/1‑68 (MQ=255)
                    tGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGACGGCAAACGAATGAATGACGTTCCAg               >  1:891415/1‑68 (MQ=255)
                         ggCTTGAGGAGATCAGCGGCGGTGATCTGTTGATCGACGGCAAACGAATGAATGACGTTc                  <  1:1042202/60‑1 (MQ=255)
                         ggCTTGAGGAGATCAGCGGCGGTGATCTGTTGATCGACGGCAAACGAATGAATGACGTTc                  <  1:181498/60‑1 (MQ=255)
                         ggCTTGAGGAGATCAGCGGCGGTGATCTGTTGATCGACGGCAAACGAATGAATGACGTTc                  <  1:1660416/60‑1 (MQ=255)
                         ggCTTGAGGAGATCAGCGGCGGTGATCTGTTGATCGACGGCAAACGAATGAATGACGTTc                  <  1:2637312/60‑1 (MQ=255)
                         ggCTTGAGGAGATCAGCGGCGGTGATCTGTTGATCGACGGCAAACGAATGAATGACGTTc                  <  1:1905120/60‑1 (MQ=255)
                         ggCTTGAGGAGATCAGCGGCGGTGATCTGTTGATCGACGGCAAACGAATGAATGACGTTc                  <  1:2041495/60‑1 (MQ=255)
                                ggagATCAGCGGCGGCGATCTGTTGATCGACGGCAAACGAATGAATGACGTTCCAGCCAAAGCACGCa   <  1:1092183/68‑1 (MQ=255)
                                ggagATCAGCGGCGGCGATCTGTTGATCGACGGCAAACGAATGAATGACGTTCCAGCCAAAGCACGCa   <  1:1396786/68‑1 (MQ=255)
                                ggagATCAGCGGCGGCGATCTGTTGATCGACGGCAAACGAATGAATGACGTTCCAGCCAAAGCACGCa   <  1:470618/68‑1 (MQ=255)
                                ggagATCAGCGGCGGCGATCTGTTGATCGACGGCAAACGAATGAATGACGTTCCAGCCAAAGCACGCa   <  1:545725/68‑1 (MQ=255)
                                ggagATCAGCGGCGGCGATCTGTTGATCGACGGCAAACGAATGAATGACGTTCCAGCCAAAGCACGCa   <  1:1416224/68‑1 (MQ=255)
                                ggagATCAGCGGCGGCGATCTGTTGATCGACGGCAAACGAATGAATGACGTTCCAGCCAAAGCACGCa   <  1:158094/68‑1 (MQ=255)
                                ggagATCAGCGGCGGCGATCTGTTGATCGACGGCAAACGAATGAATGACGTTCCAGCCAAAGCACGCa   <  1:2009396/68‑1 (MQ=255)
                                ggagATCAGCGGCGGCGATCTGTTGATCGACGGCAAACGAATGAATGACGTTCCAGCCAAAGCACGCa   <  1:1778652/68‑1 (MQ=255)
                                ggagATCAGCGGCGGCGATCTGTTGATCGACGGCAAACGAATGAATGACGTTCCAGCCAAAGCACGCa   <  1:2609838/68‑1 (MQ=255)
                                ggagATCAGCGGCGGCGATCTGTTGATCGACGGCAAACGAATGAATGACGTTCCAGCCAAAGCACGCa   <  1:3288268/68‑1 (MQ=255)
                                 gagaTCAGCGGCGGCGATCTGTTGATCGACGGCAAACgaatg                            <  1:1719764/42‑1 (MQ=255)
                                 gagaTCAGCGGCGGCGATCTGTTGATCGACGGCAAACGAATGAATGACGTTCCAGCCAAAgcaagcaa  <  1:1850527/68‑1 (MQ=255)
                                 gagaTCAGCGGCGGCGATCTGTTGATCGACGGCAAACGAATGAATGACGTTCCAGCCAAAGCACGCaa  <  1:1000973/68‑1 (MQ=255)
                                 gagaTCAGCGGCGGCGATCTGTTGATCGACGGCAAACGAATGAATGACGTTCCAGCCAAAGCACGCaa  <  1:2722484/68‑1 (MQ=255)
                                 gagaTCAGCGGCGGCGATCTGTTGATCGACGGCAAACGAATGAATGACGTTCCAGCCAAAGCACGCaa  <  1:956672/68‑1 (MQ=255)
                                 gagaTCAGCGGCGGCGATCTGTTGATCGACGGCAAACGAATGAATGACGTTCCAGCCAAAGCACGCaa  <  1:949488/68‑1 (MQ=255)
                                 gagaTCAGCGGCGGCGATCTGTTGATCGACGGCAAACGAATGAATGACGTTCCAGCCAAAGCACGCaa  <  1:923968/68‑1 (MQ=255)
                                 gagaTCAGCGGCGGCGATCTGTTGATCGACGGCAAACGAATGAATGACGTTCCAGCCAAAGCACGCaa  <  1:880012/68‑1 (MQ=255)
                                 gagaTCAGCGGCGGCGATCTGTTGATCGACGGCAAACGAATGAATGACGTTCCAGCCAAAGCACGCaa  <  1:1833798/68‑1 (MQ=255)
                                 gagaTCAGCGGCGGCGATCTGTTGATCGACGGCAAACGAATGAATGACGTTCCAGCCAAAGCACGCaa  <  1:237873/68‑1 (MQ=255)
                                 gagaTCAGCGGCGGCGATCTGTTGATCGACGGCAAACGAATGAATGACGTTCCAGCCAAAGCACGCaa  <  1:218094/68‑1 (MQ=255)
                                 gagaTCAGCGGCGGCGATCTGTTGATCGACGGCAAACGAATGAATGACGTTCCAGCCAAAGCACGCaa  <  1:2479101/68‑1 (MQ=255)
                                 gagaTCAGCGGCGGCGATCTGTTGATCGACGGCAAACGAATGAATGACGTTCCAGCCAAAGCACGCaa  <  1:1458191/68‑1 (MQ=255)
                                 gagaTCAGCGGCGGCGATCTGTTGATCGACGGCAAACGAATGAATGACGTTCCAGCCAAAGCACGCaa  <  1:2637975/68‑1 (MQ=255)
                                               |                                                     
TCGACCACCCTGCGCATGATTGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGACGGCAAACGAATGAATGACGTTCCAGCCAAAGCACGCAA  >  minE/910890‑910990

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: