Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 105,158 G→A 40.7% S378S (AGC→AGT aroP ← aromatic amino acid transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE105,1580GA40.7% 31.7 / 52.1 76S378S (AGC→AGTaroParomatic amino acid transporter
Reads supporting (aligned to +/- strand):  ref base G (23/22);  new base A (14/17);  total (37/39)
Fisher's exact test for biased strand distribution p-value = 6.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.31e-01

CAGAGCAGGGAAGCGAGTTACCACGCCTTGTTCCTGCTTGGCGCGACGGAATTTCATATGCGCCAGGCTAATCATCGCCCAGTTGATTACCAGTGCAGATACCACCAGCGCCATT  >  minE/105092‑105206
                                                                  |                                                
cagggcagggAAGCGAGTTACCACGCCTTGTTCTTGTTTGGCACGACGGAATTTCATATGCGCCAGACt                                                >  1:985512/1‑69 (MQ=255)
cagggcagggAAGCGAGTTACCACGCCTTGTTCTTGTTTGGCACGACGGAATTTCATATGCGCCAGACt                                                >  1:2553230/1‑69 (MQ=255)
cagggcagggAAGCGAGTTACCACGCCTTGTTCTTGTTTGGCACGACGGAATTTCATATGCGCCAGACt                                                >  1:2557373/1‑69 (MQ=255)
cagggcagggAAGCGAGTTACCACGCCTTGTTCTTGTTTGGCACGACGGAATTTCATATGCGCCAGACt                                                >  1:2256391/1‑69 (MQ=255)
cagggcagggAAGCGAGTTACCACGCCTTGTTCTTGTTTGGCACGACGGAATTTCATATGCGCCAGACt                                                >  1:2239245/1‑69 (MQ=255)
cagggcagggAAGCGAGTTACCACGCCTTGTTCTTGTTTGGCACGACGGAATTTCATATGCGCCAGACt                                                >  1:2217730/1‑69 (MQ=255)
cagggcagggAAGCGAGTTACCACGCCTTGTTCTTGTTTGGCACGACGGAATTTCATATGCGCCAGACt                                                >  1:300396/1‑69 (MQ=255)
cagggcagggAAGCGAGTTACCACGCCTTGTTCTTGTTTGGCACGACGGAATTTCATATGCGCCAGACt                                                >  1:3133050/1‑69 (MQ=255)
cagggcagggAAGCGAGTTACCACGCCTTGTTCTTGTTTGGCACGACGGAATTTCATATGCGCCAGACt                                                >  1:3429309/1‑69 (MQ=255)
cagggcagggAAGCGAGTTACCACGCCTTGTTCTTGTTTGGCACGACGGAATTTCATATGCGCCAGACt                                                >  1:1670859/1‑69 (MQ=255)
cagggcagggAAGCGAGTTACCACGCCTTGTTCTTGTTTGGCACGACGGAATTTCATATGCGCCAGACt                                                >  1:395133/1‑69 (MQ=255)
cagggcagggAAGCGAGTTACCACGCCTTGTTCTTGTTTGGCACGACGGAATTTCATATGCGCCAGACt                                                >  1:505933/1‑69 (MQ=255)
cagggcagggAAGCGAGTTACCACGCCTTGTTCTTGTTTGGCACGACGGAATTTCATATGCGCCAGACt                                                >  1:663380/1‑69 (MQ=255)
cagggcagggAAGCGAGTTACCACGCCTTGTTCTTGTTTGGCACGACGGAATTTCATATGCGCCAGACt                                                >  1:1438798/1‑69 (MQ=255)
   agcagGGAAGCGAGTTACCACGCCTTGTTCCTGCTTGGCGCGACGGAATTTCATATGCGCCAGGCTAat                                             <  1:389966/69‑1 (MQ=255)
   agcagGGAAGCGAGTTACCACGCCTTGTTCCTGCTTGGCGCGACGGAATTTCATATGCGCCAGGCTAat                                             <  1:276570/69‑1 (MQ=255)
   agcagGGAAGCGAGTTACCACGCCTTGTTCCTGCTTGGCGCGACGGAATTTCATATGCGCCAGGCTAat                                             <  1:2867784/69‑1 (MQ=255)
   agcagGGAAGCGAGTTACCACGCCTTGTTCCTGCTTGGCGCGACGGAATTTCATATGCGCCAGGCTAat                                             <  1:2349422/69‑1 (MQ=255)
     cagGGAAGCGAGTTACCACGCCTTGTTCCTGCTTGGCGCGACGGAATTTCATATGCGCCAGGCTAatca                                           >  1:1770503/1‑69 (MQ=255)
     cagGGAAGCGAGTTACCACGCCTTGTTCCTGCTTGGCGCGACGGAATTTCATATGCGCCAGGCTAatca                                           >  1:3239351/1‑69 (MQ=255)
     cagGGAAGCGAGTTACCACGCCTTGTTCCTGCTTGGCGCGACGGAATTTCATATGCGCCAGGCTAatca                                           >  1:157331/1‑69 (MQ=255)
     cagGGAAGCGAGTTACCACGCCTTGTTCCTGCTTGGCGCGACGGAATTTCATATGCGACAGGCTAatca                                           >  1:322880/1‑69 (MQ=255)
        ggAAGCGAGTTACCACGCCTTGTTCCTGCTTGGCGCGACGGAATTTCATATGCGCCAGGCTAATCATCg                                        <  1:324001/69‑1 (MQ=255)
        ggAAGCGAGTTACCACGCCTTGTTCCTGCTTGGCGCGACGGAATTTCATATGCGCCAGGCTAATCATCg                                        <  1:2458481/69‑1 (MQ=255)
        ggAAGCGAGTTACCACGCCTTGTTCCTGCTTGGCGCGACGGAATTTCATATGCGCCAGGCTAATCATCg                                        <  1:2405033/69‑1 (MQ=255)
        ggAAGCGAGTTACCACGCCTTGTTCCTGCTTGGCGCGACGGAATTTCATATGCGCCAGGCTAATCATCg                                        <  1:2661706/69‑1 (MQ=255)
         gAAGCGTGTTACCACGCCTTGTTCCTGCTTGGCGCGACGGAATTTCATATGCGCCAGGCTAATCATCg                                        >  1:443248/1‑68 (MQ=255)
         gAAGCGAGTTACCACGCCTTGTTCCTGCTTGGCGCGACGGAATTTCATATGCGCCAGGCTAATCATCg                                        >  1:2708161/1‑68 (MQ=255)
         gAAGCGAGTTACCACGCCTTGTTCCTGCTTGGCGCGACGGAATTTCATATGCGCCAGGCTAATCATCg                                        >  1:2828677/1‑68 (MQ=255)
         gAAGCGAGTTACCACGCCTTGTTCCTGCTTGGCGCGACGGAATTTCATATGCGCCAGGCTAATCATCg                                        >  1:65772/1‑68 (MQ=255)
         gAAGCGAGTTACCACGCCTTGTTCCTGCTTGGCGCGACGGAATTTCATATGCGCCAGGCTAATCATCg                                        >  1:642834/1‑68 (MQ=255)
         gAAGCGAGTTACCACGCCTTGTTCCTGCTTGGCGCGACGGAATTTCATATGCGCCAGGCTAATCATCg                                        >  1:530296/1‑68 (MQ=255)
         gAAGCGAGTTACCACGCCTTGTTCCTGCTTGGCGCGACGGAATTTCATATGCGCCAGGCTAATCATCg                                        >  1:293280/1‑68 (MQ=255)
         gAAGCGAGTTACCACGCCTTGTTCCTGCTTGGCGCGACGGAATTTCATATGCGCCAGGCTAATCATCg                                        >  1:3024688/1‑68 (MQ=255)
         gAAGCGAGTTACCACGCCTTGTTCCTGCTTGGCGCGACGGAATTTCATATGCGCCAGGCTAATCATCg                                        >  1:3188176/1‑68 (MQ=255)
         gAAGCGAGTTACCACGCCTTGTTCCTGCTTGGCGCGACGGAATTTCATATGCGCCAGGCTAATCATCg                                        >  1:1953191/1‑68 (MQ=255)
         gAAGCGAGTTACCACGCCTTGTTCCTGCTTGGCGCGACGGAATTTCATATGCGCCAGGCTAATCATCg                                        >  1:1424792/1‑68 (MQ=255)
         gAAGCGAGTTACCACGCCTTGTTCCTGCTTGGCGCGACGGAATTTCATATGCGCCAGGCTAATCATCg                                        >  1:1688739/1‑68 (MQ=255)
         gAAGCGAGTTACCACGCCTTGTTCCTGCTTGGCGCGACGGAATTTCATATGCGCCAGGCTAATCATCg                                        >  1:1204387/1‑68 (MQ=255)
         gAAGCGAGTTACCACGCCTTGTTCCTGCTTGGCGCGACGGAATTTCATATGCGCCAGGCTAATCATCg                                        >  1:1936409/1‑68 (MQ=255)
         gAAGCGAGTTACCACGCCTTGTTCCTGCTTGGCGCGACGGAATTTCATATGCGCCAGGCTAATCATCg                                        >  1:2320139/1‑68 (MQ=255)
             cGAGTTACCACGCCTTGTTCCTGCTTGGCGCGACGGAATTTCATATGCGCCAGGCTAATCATCGCCCAg                                   <  1:3001362/69‑1 (MQ=255)
             cGAGTTACCACGCCTTGTTCCTGCTTGGCGCGACGGAATTTCATATGCGCCAGGCTAATCATCGCCCAg                                   <  1:2502304/69‑1 (MQ=255)
             cGAGTTACCACGCCTTGTTCCTGCTTGGCGCGACGGAATTTCATATGCGCCAGGCTAATCATCGCCCAg                                   <  1:3186164/69‑1 (MQ=255)
                gTTACCACGCCTTGTTCCTGCTTGGCGCGACGGAATTTCATATGCGCCAGGCTAATCATCGCCCAGTTg                                <  1:2965201/69‑1 (MQ=255)
                   aCCACGCCTTGTTCCTGCTTGGCGCGACGGAATTTCATATGCGCCAGGCTAATCATCGCCCAGTTGAtt                             >  1:203628/1‑69 (MQ=255)
                   aCCACGCCTTGTTCCTGCTTGGCGCGACGGAATTTCATATGCGCCAGGCTAATCATCGCCCAGTTGAtt                             >  1:1057318/1‑69 (MQ=255)
                   aCCACGCCTTGTTCCTGCTTGGCGCGACGGAATTTCATATGCGCCAGGCTAATCATCGCCCAGTTGAtt                             >  1:1878790/1‑69 (MQ=255)
                   aCCACGCCTTGTTCCTGCTTGGCGCGACGGAATTTCATATGCGCCAGGCTAATCATCGCCCAGTTGAtt                             >  1:1845618/1‑69 (MQ=255)
                            tGTTCCTGCTTGGCGCGACGGAATTTCATATGCGCCAGGCTAATCATCGCCCAGTTGATTACCAGTGCa                    <  1:2391600/69‑1 (MQ=255)
                            tGTTCCTGCTTGGCGCGACGGAATTTCATATGCGCCAGGCTAATCATCGCCCAGTTGATTACCAGTGCa                    <  1:57120/69‑1 (MQ=255)
                            tGTTCCTGCTTGGCGCGACGGAATTTCATATGCGCCAGGCTAATCATCGCCCAGTTGATTACCAGTGCa                    <  1:3065525/69‑1 (MQ=255)
                            tGTTCCTGCTTGGCGCGACGGAATTTCATATGCGCCAGGCTAATCATCGCCCAGTTGATTACCAGTGCa                    <  1:2343717/69‑1 (MQ=255)
                            tGTTCCTGCTTGGCGCGACGGAATTTCATATGCGCCAGGCTAATCATCGCCCAGTTGATTACCAGTGCa                    <  1:3285303/69‑1 (MQ=255)
                                  tGCTTGGCGCGACGGAATTTCATATGCGCCAGGCTAATCATCGCCCAGTTGATTACCAGTGCAGATac               <  1:657956/68‑1 (MQ=255)
                                   gCTTGGCGCGACGGAATTTCATATGCGCCAGGCTAATCATCg                                        <  1:88913/42‑1 (MQ=255)
                                   gCTTGGCGCGACGGAATTTCATATGCGCCAGGCTAATCATCGCCCAGTTGATTACCAGTGCAGATacc              <  1:741337/68‑1 (MQ=255)
                                          acgACGGAATTTCATATGCGCCAGACTGATCATCGCCCAGTTGATGACCAGTGCAGATACCACCAgcg       <  1:1286093/67‑1 (MQ=255)
                                          acgACGGAATTTCATATGCGCCAGACTGATCATCGCCCAGTTGATGACCAGTGCAGATACCACCAgcg       <  1:113931/67‑1 (MQ=255)
                                          acgACGGAATTTCATATGCGCCAGACTGATCATCGCCCAGTTGATGACCAGTGCAGATACCACCAgcg       <  1:1163476/67‑1 (MQ=255)
                                          acgACGGAATTTCATATGCGCCAGACTGATCATCGCCCAGTTGATGACCAGTGCAGATACCACCAgcg       <  1:567924/67‑1 (MQ=255)
                                          acgACGGAATTTCATATGCGCCAGACTGATCATCGCCCAGTTGATGACCAGTGCAGATACCACCAgcg       <  1:1357487/67‑1 (MQ=255)
                                          acgACGGAATTTCATATGCGCCAGACTGATCATCGCCCAGTTGATGACCAGTGCAGATACCACCAgcg       <  1:1513189/67‑1 (MQ=255)
                                          acgACGGAATTTCATATGCGCCAGACTGATCATCGCCCAGTTGATGACCAGTGCAGATACCACCAgcg       <  1:418488/67‑1 (MQ=255)
                                          acgACGGAATTTCATATGCGCCAGACTGATCATCGCCCAGTTGATGACCAGTGCAGATACCACCAgcg       <  1:1529851/67‑1 (MQ=255)
                                          acgACGGAATTTCATATGCGCCAGACTGATCATCGCCCAGTTGATGACCAGTGCAGATACCACCAgcg       <  1:1671095/67‑1 (MQ=255)
                                          acgACGGAATTTCATATGCGCCAGACTGATCATCGCCCAGTTGATGACCAGTGCAGATACCACCAgcg       <  1:367958/67‑1 (MQ=255)
                                          acgACGGAATTTCATATGCGCCAGACTGATCATCGCCCAGTTGATGACCAGTGCAGATACCACCAgcg       <  1:172190/67‑1 (MQ=255)
                                          acgACGGAATTTCATATGCGCCAGACTGATCATCGCCCAGTTGATGACCAGTGCAGATACCACCAgcg       <  1:3194644/67‑1 (MQ=255)
                                          acgACGGAATTTCATATGCGCCAGACTGATCATCGCCCAGTTGATGACCAGTGCAGATACCACCAgcg       <  1:3162080/67‑1 (MQ=255)
                                          acgACGGAATTTCATATGCGCCAGACTGATCATCGCCCAGTTGATGACCAGTGCAGATACCACCAgcg       <  1:2053284/67‑1 (MQ=255)
                                          acgACGGAATTTCATATGCGCCAGACTGATCATCGCCCAGTTGATGACCAGTGCAGATACCACCAgcg       <  1:2213128/67‑1 (MQ=255)
                                          acgACGGAATTTCATATGCGCCAGACTGATCATCGCCCAGTTGATGACCAGTGCAGATACCACCAgcg       <  1:2358946/67‑1 (MQ=255)
                                          acgACGGAATTTCATATGCGCCAGACTGATCATCGCCCAGTTGATGACCAGTGCAGATACCACCAgcg       <  1:2705610/67‑1 (MQ=255)
                                              cGGAATTTCATATGCGCCAGGCTAATCATCGCCCAGTTGATTACCAGTGCAGATACCACCAGCGCCAtt  <  1:63485/69‑1 (MQ=255)
                                              cGGAATTTCATATGCGCCAGGCTAATCATCGCCCAGTTGATTACCAGTGCAGATACCACCAGCGCCAtt  <  1:87605/69‑1 (MQ=255)
                                                                  |                                                
CAGAGCAGGGAAGCGAGTTACCACGCCTTGTTCCTGCTTGGCGCGACGGAATTTCATATGCGCCAGGCTAATCATCGCCCAGTTGATTACCAGTGCAGATACCACCAGCGCCATT  >  minE/105092‑105206

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: