Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 904,059 C→A 48.0% R10R (CGC→CGA ycjP → predicted sugar transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE904,0590CA48.0% 16.1 / 80.0 75R10R (CGC→CGAycjPpredicted sugar transporter subunit
Reads supporting (aligned to +/- strand):  ref base C (24/15);  new base A (19/17);  total (43/32)
Fisher's exact test for biased strand distribution p-value = 4.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATTCTGCTGACCAGGAAACGGGTGAACCTCAATGGCAACAAATAAACGCACACTCAGTCGCATCGGTTTTTACTGCGGGCTGGCGCTGTTTCTCATCATCACGCTGTTTCCATTTTTTGTGATGCT  >  minE/903999‑904124
                                                            |                                                                 
aTTCTGCTGACCAGGAAACGGGTGAACCTCAATGGCAACAAATAAACGCACACTCAGTCGCATCGGtt                                                            >  1:626543/1‑68 (MQ=255)
aTTCTGCTGACCAGGAAACGGGTGAACCTCAATGGCAACAAATAAACGCACACTCAGTCGCATCGGtt                                                            >  1:1970961/1‑68 (MQ=255)
aTTCTGCTGACCAGGAAACGGGTGAACCTCAATGGCAACAAATAAACGCACACTCAGTCGCATCGGt                                                             <  1:1854982/67‑1 (MQ=255)
                    ggTGAACCTCAATGTCAACAAATAAACGCACACTCAGTCGCATCGGtttt                                                          >  1:647657/1‑50 (MQ=255)
                    ggTGAACCTCAATGGCAACAAATAAACGCACACTCAGTCGCATCGGtttt                                                          >  1:1526893/1‑50 (MQ=255)
                    ggTGAACCTCAATGGCAACAAATAAACGCACACTCAGTCGCATCGGtttt                                                          >  1:1834867/1‑50 (MQ=255)
                    ggTGAACCTCAATGGCAACAAATAAACGCACACTCAGTCGCATCGGtttt                                                          >  1:492016/1‑50 (MQ=255)
                    ggTGAACCTCAATGGCAACAAATAAACGCACACTCAGTCGCATCGGtttt                                                          >  1:496042/1‑50 (MQ=255)
                        aaCCTCAATGGCAACAAATAAAGGCACACTCAGTCGc                                                                   >  1:2410185/1‑37 (MQ=255)
                        aaCCTCAATGGCAACAAATAAACGCACACTCAGTCGc                                                                   >  1:989281/1‑37 (MQ=255)
                        aaCCTCAATGGCAACAAATAAACGCACACTCAGTCGc                                                                   >  1:600858/1‑37 (MQ=255)
                        aaCCTCAATGGCAACAAATAAACGCACACTCAGTCGc                                                                   >  1:3151112/1‑37 (MQ=255)
                        aaCCTCAATGGCAACAAATAAACGCACACTCAGTCGc                                                                   >  1:2216102/1‑37 (MQ=255)
                        aaCCTCAATGGCAACAAATAAACGCACACTCAGTCGc                                                                   >  1:2352726/1‑37 (MQ=255)
                        aaCCTCAATGGCAACAAATAAACGCACACTCAGTCGc                                                                   >  1:163486/1‑37 (MQ=255)
                             cAATGGCAACAAATAAACGCACACTCAGTCGCATCGGTTTTTACTGCGGGCTGGCGCt                                         >  1:172485/1‑58 (MQ=255)
                             cAATGGCAACAAATAAACGCACACTCAGTCGCATCGGTTTTTACTGCGGGCTGGCGCTGt                                       >  1:2447197/1‑60 (MQ=255)
                             cAATGGCAACAAATAAACGCACACTCAGTCGCATCGGTTTTTACTGCGGGCTGGCGCTGt                                       <  1:656376/60‑1 (MQ=255)
                             cAATGGCAACAAATAAACGCACACTCAGTCGCATCGGTTTTTACTGCGGGCTGGCGCTGt                                       >  1:2205986/1‑60 (MQ=255)
                             cAATGGCAACAAATAAACGCACACTCAGTCGCATCGGTTTTTACTGCGGGCTGGCGCTGt                                       >  1:2451087/1‑60 (MQ=255)
                             cAATGGCAACAAATAAACGCACACTCAGTCGCATCGGTTTTTACTGCGGGCTGGCGCTGt                                       >  1:3029606/1‑60 (MQ=255)
                             cAATGGCAACAAATAAACGCACACTCAGTCGCATCGGTTTTTACTGCGGGCTGGCGCTGt                                       <  1:2318117/60‑1 (MQ=255)
                             cAATGGCAACAAATAAACGCACACTCAGTCGCATCGGTTTTTACTGCGGGCTGGCGCTGt                                       <  1:3273051/60‑1 (MQ=255)
                             cAATGGCAACAAATAAACGCACACTCAGTCGCATCGGTTTTTACTGCGGGCTGGCGCTGt                                       <  1:334767/60‑1 (MQ=255)
                             cAATGGCAACAAATAAACGCACACTCAGTCGCATCGGTTTTTACTGCGGGCTGGCGCTGt                                       <  1:41532/60‑1 (MQ=255)
                             cAATGGCAACAAATAAACGCACACTCAGTCGCATCGGTTTTTACTGCGGGCTGGCGCTGt                                       <  1:454599/60‑1 (MQ=255)
                             cAATGGCAACAAATAAACGCACACTCAGTCGCATCGGTTTTTACTGCGGGCTGGCGCTGt                                       <  1:1815008/60‑1 (MQ=255)
                             cAATGGCAACAAATAAACGCACACTCAGTCGCATCGGTTTTTACTGCGGGCTGGCGCTGt                                       <  1:1582559/60‑1 (MQ=255)
                             cAATGGCAACAAATAAACGCACACTCAGTCGCATCGGTTTTTACTGCGGGCTGGCGCTGt                                       >  1:1525506/1‑60 (MQ=255)
                             cAATGGCAACAAATAAACGCACACTCAGTCGCATCGGTTTTTACTGCGGGCTGGCGCTGt                                       <  1:1321280/60‑1 (MQ=255)
                             cAATGGCAACAAATAAACGCACACTCAGTCGCATCGGTTTTTACTGCGGGCTGGCGCTGt                                       >  1:1282128/1‑60 (MQ=255)
                             cAATGGCAACAAATAAACGCACACTCAGTCGCATCGGTTTTACTGCGGGCTGGCGCTGt                                       <  1:2274608/59‑1 (MQ=255)
                              aaTGGCAACAAATAAACGCACACTCAGTCGCAGCGGTTTTTACTGCGGGCTGGCGCTGt                                       <  1:1113140/59‑1 (MQ=255)
                                  gCAACAAATAAACGCACACTCAGTCGCATCGGTTTTTACTGCGGGCTGGCGCTGTTTCTCATCATCACg                         >  1:3011531/1‑69 (MQ=255)
                                  gCAACAAATAAACGCACACTCAGTCGCATCGGTTTTTACTGCGGGCTGGCGCTGTTTCTCATCATCACg                         >  1:988357/1‑69 (MQ=255)
                                  gCAACAAATAAACGCACACTCAGTCGCATCGGTTTTTACTGCGGGCTGGCGCTGTTTCTCATCATCACg                         >  1:495495/1‑69 (MQ=255)
                                  gCAACAAATAAACGCACACTCAGTCGCATCGGTTTTTACTGCGGGCTGGCGCTGTTTCTCATCATCACg                         >  1:1319965/1‑69 (MQ=255)
                                  gCAACAAATAAACGCACACTCAGTCGCATCGGTTTTTACTGCGGGCTGGCGCTGTTTCTCATCATCACg                         >  1:1461278/1‑69 (MQ=255)
                                  gCAACAAATAAACGCACACTCAGTCGCATCGGTTTTTACTGCGGGCTGGCGCTGTTTCTCATCATCACg                         >  1:2358225/1‑69 (MQ=255)
                                  gCAACAAATAAACGCACACTCAGTCGCATCGGTTTTTACTGCGGGCTGGCGCTGTTTCTCATCATCACg                         >  1:3274634/1‑69 (MQ=255)
                                       aaaTAAACGCACACTCAGTCGAATCGGTTTTTACTGCGGGCTGGCGCTGTTTCTCATCATGACGCTGt                     <  1:2194469/68‑1 (MQ=255)
                                       aaaTAAACGCACACTCAGTCGAATCGGTTTTTACTGCGGGCTGGCGCTGTTTCTCATCATCACGCTGt                     <  1:2798268/68‑1 (MQ=255)
                                       aaaTAAACGCACACTCAGTCGAATCGGTTTTTACTGCGGGCTGGCGCTGTTTCTCATCATCACGCTGt                     <  1:2038905/68‑1 (MQ=255)
                                       aaaTAAACGCACACTCAGTCGAATCGGTTTTTACTGCGGGCTGGCGCTGTTTCTCATCATCACGCTGt                     <  1:2596675/68‑1 (MQ=255)
                                       aaaTAAACGCACACTCAGTCGAATCGGTTTTTACTGCGGGCTGGCGCTGTTTCTCATCATCACGCTGt                     <  1:1631911/68‑1 (MQ=255)
                                       aaaTAAACGCACACTCAGTCGAATCGGTTTTTACTGCGGGCTGGCGCTGTTTCTCATCATCACGCTGt                     <  1:1085684/68‑1 (MQ=255)
                                       aaaTAAACGCACACTCAGTCGAATCGGTTTTTACTGCGGGCTGGCGCTGTTTCTCATCATCACGCTGt                     <  1:24437/68‑1 (MQ=255)
                                       aaaTAAACGCACACTCAGTCGAATCGGTTTTTACTGCGGGCTGGCGCTGTTTCTCATCATCACGCTGt                     <  1:1677337/68‑1 (MQ=255)
                                       aaaTAAACGCACACTCAGTCGAATCGGTTTTTACTGCGGGCTGGCGCTGTTTCTCATCATCACGCTGt                     <  1:1077505/68‑1 (MQ=255)
                                       aaaTAAACGCACACTCAGTCGAATCGGTTTTTACTGCGGGCTGGCGCTGTTTCTCATCATCACGCTGt                     <  1:2254150/68‑1 (MQ=255)
                                       aaaTAAACGCACACTCAGTCGAATCGGTTTTTACTGCGGGCTGGCGCTGTTTCTCATCATCACGCTGt                     <  1:197231/68‑1 (MQ=255)
                                       aaaTAAACGCACACTCAGTCGAATCGGTTTTTACTGCGGGCTGGCGCTGTTTCTCATCATCACGCTGt                     <  1:1268084/68‑1 (MQ=255)
                                       aaaTAAACGCACACTCAGTCGAATCGGTTTTTACTGCGGGCTGGCGCTGTTTCTCATCATCACGCTGt                     <  1:607885/68‑1 (MQ=255)
                                       aaaTAAACGCACACTCAGTCGAATCGGTTTTTACTGCGGGCTGGCGCTGTTTCTCAACATCACGCTGt                     <  1:107798/68‑1 (MQ=255)
                                        aaTAAACGCACACTCAGTCGAATCGGTTTTTACTGCGGGCTGGCGCTGTTTCTCATCATCACGCTGt                     <  1:2391445/67‑1 (MQ=255)
                                          tAAACGCACACTCAGTCGAATCGGTTTTTACTGCGGGCTGGCGCTGTTTCTCATCATCACGCTGt                     <  1:2138786/65‑1 (MQ=255)
                                                   acTCAGTCGCATCGGTTTTTACTGCGGGCTGGCGCTGTTTCTCATCATCACGCTGTTTCCATTTTTtgt        <  1:2880855/69‑1 (MQ=255)
                                                   acTCAGTCGCATCGGTTTTTACTGCGGGCTGGCGCTGTTTCTCATCATCACGCTGTTTCCATTTTTtgt        <  1:3001546/69‑1 (MQ=255)
                                                   acTCAGTCGCATCGGTTTTTACTGCGGGCTGGCGCTGTTTCTCATCATCACGCTGTTTCCATTTTTtgt        <  1:2202703/69‑1 (MQ=255)
                                                   acTCAGTCGAATCGGTTTTTACTGCGGGCTGGCGCTGTTTCTCATCATCACGCTGt                     <  1:180883/56‑1 (MQ=255)
                                                   acTCAGTCGAATCGGTTTTTACTGCGGGCTGGCGCTGTTTCTCATCATCACGCTGTTTcc                 >  1:1699206/1‑60 (MQ=255)
                                                   acTCAGTCGAATCGGTTTTTACTGCGGGCTGGCGCTGTTTCTCATCATCACGCTGTTTCCATTTTTtgt        >  1:2330703/1‑69 (MQ=255)
                                                   acTCAGTCGAATCGGTTTTTACTGCGGGCTGGCGCTGTTTCTCATCATCACGCTGTTTCCATTTTTtgt        >  1:928682/1‑69 (MQ=255)
                                                   acTCAGTCGAATCGGTTTTTACTGCGGGCTGGCGCTGTTTCTCATCATCACGCTGTTTCCATTTTTtgt        >  1:76540/1‑69 (MQ=255)
                                                   acTCAGTCGAATCGGTTTTTACTGCGGGCTGGCGCTGTTTCTCATCATCACGCTGTTTCCATTTTTtgt        >  1:751379/1‑69 (MQ=255)
                                                   acTCAGTCGAATCGGTTTTTACTGCGGGCTGGCGCTGTTTCTCATCATCACGCTGTTTCCATTTTTtgt        >  1:1951516/1‑69 (MQ=255)
                                                   acTCAGTCGAATCGGTTTTTACTGCGGGCTGGCGCTGTTTCTCATCATCACGCTGTTTCCATTTTTtgt        >  1:1794232/1‑69 (MQ=255)
                                                   acTCAGTCGAATCGGTTTTTACTGCGGGCTGGCGCTGTTTCTCATCATCACGCTGTTTCCATTTTTtgt        >  1:1211635/1‑69 (MQ=255)
                                                   acTCAGTCGAATCGGTTTTTACTGCGGGCTGGCGCTGTTTCTCATCATCACGCTGTTTCCATTTTTtgt        >  1:2181551/1‑69 (MQ=255)
                                                   acTCAGTCGAATCGGTTTTTACTGCGGGCTGGCGCTGTTTCTCATCATCACGCTGTTTCCATTTTTtgt        >  1:1789873/1‑69 (MQ=255)
                                                   acTCAGTCGAATCGGTTTTTACTGCGGGCTGGCGCTGTTTCTCATCATCACGCTGTTTCCATTTTTtgt        >  1:1758546/1‑69 (MQ=255)
                                                   acTCAGTCGAATCGGTTTTTACTGCGGGCTGGCGCTGTTTCTCATCATCACGCTGTTTCCATTTTTtgt        >  1:2830527/1‑69 (MQ=255)
                                                   acTCAGTCGAATCGGTTTTTACTGCGGGCTGGCGCTGTTTCTCATCATCACGCTGTTTCCATTTTTtgt        >  1:1393771/1‑69 (MQ=255)
                                                   acTCAGTCGAATCGGTTTTTACTGCGGGCTGGCGCTGTTTCTCATCATCACGCTGTTTCCATTTTTtgt        >  1:391604/1‑69 (MQ=255)
                                                   acTCAGTCGAATCGGTTTTTACTGCGGGCTGGCGCTGTTTCTCATCATCACGCTGTTTCCATTTTTtgt        >  1:37842/1‑69 (MQ=255)
                                                   acTCAGTCGAATCGGTTTTTACTGCGGGCTGGCGCTGTTTCTCATCATCACGCTGTTTCCATTTTTtgt        >  1:163255/1‑69 (MQ=255)
                                                   acTCAGTCGAATCGGTTTTTACTGCGGGCTGGCGCTGTTTCTCATCATCACGCTGTTTCCATTTTTtgt        >  1:2714626/1‑69 (MQ=255)
                                                   acTCAGTCGAATCGGTTTTTACTGCGGGCTGGCGCTGTTTCTCATCATCACGCTGTTTCCATTTTTtgt        >  1:2971888/1‑69 (MQ=255)
                                                   acTCAGTCGAATCGGTTTTTACTGCGGGCTGGCGCTGTTTCTCATCATCACGCTGTTTCCATTTTTtgt        >  1:2913265/1‑69 (MQ=255)
                                                         tCGCATCGGTTTTTACTGCGGGCTGGCGCTGTTTCTCATCATCACGCTGTTTCCATTTTTTGTGATGCt  >  1:1507114/1‑69 (MQ=255)
                                                         tCGCATCGGTTTTTACTGCGGGCTGGCGCTGTTTCTCATCATCACGCTGTTTCCATTTTTTGTGATGCt  >  1:2768619/1‑69 (MQ=255)
                                                         tCGCATCGGTTTTTACTGCGGGCTGGCGCTGTTTCTCATCATCACGCTGTTTCCATTTTTTGTGATGCt  >  1:1821865/1‑69 (MQ=255)
                                                            |                                                                 
ATTCTGCTGACCAGGAAACGGGTGAACCTCAATGGCAACAAATAAACGCACACTCAGTCGCATCGGTTTTTACTGCGGGCTGGCGCTGTTTCTCATCATCACGCTGTTTCCATTTTTTGTGATGCT  >  minE/903999‑904124

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: