Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 105,984 T→A 5.1% Y103F (TAC→TTC)  aroP ← aromatic amino acid transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE105,9840TA 5.1% 231.9 / 3.6 79Y103F (TAC→TTC) aroParomatic amino acid transporter
Reads supporting (aligned to +/- strand):  ref base T (42/33);  new base A (2/2);  total (44/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.08e-01

GGGATTTCCGGATACCAGAACTGAATGTATTTACCCACGGCAGTCAGCTCAGCCATGGCAACTAAAACGTACAGTACCCAGTAGTTCCAGCCAGAGGCGAAACCGGCAAAACTGCCCCAGTATTTATAAGCAAAGTGGC  >  minE/105915‑106053
                                                                     |                                                                     
gggATTTCCGGATACCAGAACTGAATGTATTTACCCACGGCAGTCAGCTCAGCCATGGCAACTAAAACGTa                                                                      <  1:3771609/71‑1 (MQ=255)
gggATTTCCGGATACCAGAACTGAATGTATTTACCCACGGCAGTCAGCTCAGCCATGGCAACTAAAACGTa                                                                      <  1:1794481/71‑1 (MQ=255)
 ggATTTCCGGATACCAGAACTGAATGTATTTACCCACGGCAGTCAGCTCAGCCATGGCAACTAAAACGTa                                                                      <  1:3707153/70‑1 (MQ=255)
   aTTTCCGGATACCAGAACTGAATGTATTTACCCACGGCAGTCAGCTCAGCCATGGCAACTAAAACGTACAg                                                                   >  1:3642768/1‑71 (MQ=255)
    tttCCGGATACCAGAACTGAATGTATTTACCCACGGCAGTCAGCTCAGCCATGGCAACTAAAACGTACAGt                                                                  <  1:762276/71‑1 (MQ=255)
    tttCCGGATACCAGAACTGAATGTATTTACCCACGGCAGTCAGCTCAGCCATGCCAACTAAAACGTACAg                                                                   >  1:2258200/1‑70 (MQ=255)
       ccGGATACCAGAACTGAATGTATTTACCCACGGCAGTCAGCTCAGCCATGGCAACTAAAACGTACAGTAc                                                                >  1:2371653/1‑70 (MQ=255)
       ccGGATACCAGAACTGAATGTATTTACCCACGGCAGTCAGCTCAGCCATGGCAACTAAAACGTACAGTAc                                                                <  1:1621756/70‑1 (MQ=255)
            tACAAGAACTGAATGTATTTACCCACGGCAGTCAGCTCAGCCATGGCAACTAAAACGTACAGTa                                                                 >  1:78304/1‑64 (MQ=255)
                aGAACTGAATGTATTTACCCACGGCAGTCAGCTCAGCCATGGCAACTAAAACGTACa                                                                    >  1:332819/1‑57 (MQ=255)
                aGAACTGAATGTATTTACCCACGGCAGTCAGCTCAGCCATGGCAACTAAAACGTACAGTACCCagtagta                                                       >  1:1036227/1‑69 (MQ=255)
                 gAACTGAATGTATTTACCCACGGCAGTCAGCTCAGCCATGGCAACTAAAACGTACAGTACCCag                                                            <  1:4089615/64‑1 (MQ=255)
                 gAACTGAATGTATTTACCCACGGCAGTCAGCTCAGCCATGGCAACTAAAACGTACAGTACCCAGTAGTTc                                                      >  1:1210586/1‑70 (MQ=255)
                     tGAATGTATTTACCCACGGCAGTCAGCTCAGCCATGGCAACTAAAACGTa                                                                      <  1:382850/50‑1 (MQ=255)
                       aaagTATTTACCCACGGCAGTCAGCTCAGCCATGGCAACTAAAACGTACAGTACCCAGTAGTTccagccag                                               <  1:3196029/68‑1 (MQ=255)
                        aTGTATTTACCCACGGCAGTCAGCTCAGCCATGGCAACTAAAACGt                                                                       <  1:1921153/46‑1 (MQ=255)
                         tGTATTTACCCACGGCAGTCAGCTCAGCCATGGCAACTAAAACGTa                                                                      <  1:1368323/46‑1 (MQ=255)
                         tGTATTTACCCACGGCAGTCAGCTCAGCCATGGCAACTAAAACGTACAGTACCCAGTAGTTCCAGCCagag                                             <  1:2176903/71‑1 (MQ=255)
                          gTATTTACCCACGGCAGTCAGCTCAGCCATGGCAACTAAAACGTAc                                                                     >  1:1701680/1‑46 (MQ=255)
                          gTATTTACCCACGGCAGTCAGCTCAGCCATGGCAACTAAAACGTACAGTAc                                                                <  1:434380/51‑1 (MQ=255)
                           tATTTACCCACGGCAGTCAGCTCAGCCATGGCAACTAAAACGTACa                                                                    >  1:2311322/1‑46 (MQ=255)
                           tATTTACCCACGGCAGTCAGCTCAGCCATGGCAACTAAAACGTACa                                                                    >  1:1990617/1‑46 (MQ=255)
                           tATTTACCCACGGCAGTCAGCTCAGCCATGGCAACTAAAACGTACa                                                                    >  1:135300/1‑46 (MQ=255)
                           tATTTACCCACGGCAGTCAGCTCAGCCATGGCAACTAAAACGTACAGTAc                                                                <  1:1434441/50‑1 (MQ=255)
                           tATTTACCCACGGCAGTCAGCTCAGCCATGGCAACTAAAACGTACAGTACCCAGTAGTTCCAGCCAGAGGc                                           >  1:1282508/1‑71 (MQ=255)
                           tATTTACCCACGGCAGTCAGCTCAGCCATGGCAACTAAAACGTACAGTACCCAGTAGTTCCAGCCAGAGGc                                           >  1:1307945/1‑71 (MQ=255)
                           tATTTACCCACGGCAGTCAGCTCAGCCATGGCAACTAAAACGTACAGTACCCAGTAGTTCAAGCCAGAGGc                                           >  1:3809937/1‑71 (MQ=255)
                            aTTTACCCACGGCAGTCAGCTCAGCCATGGCAACTAAAACGTa                                                                      >  1:2168005/1‑43 (MQ=255)
                             tttACCCACGGCAGTCAGCTCAGCCATGGCAACTAAAACGTACAGTACCCagta                                                          >  1:2671113/1‑54 (MQ=255)
                               tACCCACGGCAGTCAGCTCAGCCATGGCAACTAAAACGTACTGTACCCAGTAGTTCCAGCCAGAGGCGaaa                                       >  1:4048448/1‑71 (MQ=255)
                               tACCCACGGCAGTCAGCTCAGCCATGGCAACTAAAACGTACAGTAccc                                                              >  1:3524614/1‑48 (MQ=255)
                               tACCCACGGCAGTCAGCTCAGCCATGGCAACTAAAACGTACAGTAccc                                                              >  1:2281605/1‑48 (MQ=255)
                               tACCCACGGCAGTCAGCTCAGCCATGGCAACTAAAACGTACAGTAccc                                                              >  1:152817/1‑48 (MQ=255)
                               tACCCACGGCAGTCAGCTCAGCCATGGCAACTAAAACGTACAGTAccc                                                              >  1:157566/1‑48 (MQ=255)
                               tACCCACGGCAGTCAGCTCAGCCATGGCAACTAAAACGTACAGTACCCa                                                             >  1:1071948/1‑49 (MQ=255)
                               tACCCACGGCAGTCAGCTCAGCCATGGCAACTAAAACGTACAGTACCCAGTAGTTCCAGCCAGAGGCGaaa                                       >  1:1627064/1‑71 (MQ=255)
                               tACCCACGGCAGTCAGCTCAGCCATGGCAACTAAAACGTACAGTACCCAGTAGTTCCAGCCAGAGGCGaa                                        <  1:3544932/70‑1 (MQ=255)
                               tACCCACGGCAGTCAGCTCAGCCATGGCAACTAAAACGTACAGTACCCAGTAGTTCCAGCCAGAGGCGaa                                        <  1:1363538/70‑1 (MQ=255)
                                 cccACGGCAGTCAGCTCAGCCATGGCAACTAAAACGTACAGTACCCAGTAGTTccagc                                                  >  1:1188192/1‑58 (MQ=255)
                                  ccACGGCAGTCAGCTCAGCCATGGCAACTAAAACGTa                                                                      <  1:228981/37‑1 (MQ=255)
                                    aCGGCAGTCAGCTCAGCCATGGCAACTAAAACGTACAGTACCCa                                                             >  1:665807/1‑44 (MQ=255)
                                     cGGCAGTCAGCTCAGCCATGGCAAATAAAACGTACAg                                                                   >  1:2030201/1‑37 (MQ=255)
                                      ggCAGTCAGCTCAGCCATGGCAACTAAAACGTACAGTACCCAGTAGTTCCAGCCAGAGGCGAAATCGGCaa                                <  1:533082/71‑1 (MQ=255)
                                      ggCAGTCAGCTCAGCCATGGCAACTAAAACGTACAGTACCCAGTAGTTCCAGCCAGAGGCGAAACCGGCaa                                <  1:133881/71‑1 (MQ=255)
                                      ggCAGTCAGCTCAGCCATGGCAACTAAAACGTACAGTACCCAGTAGTTCCAGCCAGAGGCGAAACCGGCaa                                <  1:2715884/71‑1 (MQ=255)
                                       gCAGTCAGCTCAGCCATGGCAACTAAAACGTACAGTACCCAGTAGTTCCAGCCAGAGGCGAAACCGGCaa                                >  1:3114385/1‑70 (MQ=255)
                                         ggTCAGCTCTGCCATTCCCACCAGCACGAACATTACCCAGTAGTTccagcca                                                >  1:3868160/2‑52 (MQ=255)
                                         ggTCAGCTCTGCCATTCCCACCAGCACGAACATTACCCAGTAGTTccagcca                                                >  1:3834952/2‑52 (MQ=255)
                                         aGTCAGCTCAGCCATGGCAACTAAAACGTACAGTACCCAGTAGTTCCAGCCAGAGGCGAAACCGGCAAAAc                             <  1:3514581/71‑1 (MQ=255)
                                          gTCAGCTCTGCCATTCCCACCAGCACGAACATTACCCAGTAGTTCCAGCCaga                                              <  1:777500/53‑1 (MQ=255)
                                          gTCAGCTCAGCCATGGCAACTAAAACGTACAGTACCCAGTAGTTCCAGCCAGAGGCGAAACCGGCaaa                               <  1:2791479/68‑1 (MQ=255)
                                          gTCAGCTCAGCCATGGCAACTAAAACGTACAGTACCCAGTAGTTCCAGCCAGAGGCGAAACCGGCaaa                               <  1:385453/68‑1 (MQ=255)
                                           tcagctcagcCATGGCAACTAAAACGTACAGTACCCagt                                                           <  1:2434190/39‑1 (MQ=255)
                                           tcagctcagcCATGGCAACTAAAACGTACAGTACCCAGAAGTTCCAGCCAGAGGCGAAACCGGCaaa                               <  1:2172625/67‑1 (MQ=255)
                                             agctcagcCATGGCAACTAAAACGTACAGTACCCa                                                             >  1:140936/1‑35 (MQ=255)
                                             agctcagcCATGGCAACTAAAACGTACAGTACCCAGTAGTTc                                                      >  1:21301/1‑42 (MQ=255)
                                             agctcagcCATGGCAACTAAAACGTACAGTACCCAGTAGTTc                                                      >  1:55725/1‑42 (MQ=255)
                                             agctcagcCATGGCAACTAAAACGTACAGTACCCAGTAGTTCCAGCCAGAGGCGAAACCGGCAAAACTGcc                         >  1:1201327/1‑71 (MQ=255)
                                               ctcagcCATGGCAACTAAAACGTACAGTACCCAGTAGTTCCAGCCAGAGGCGAAACCGGCAAAACTGcccc                       <  1:74492/71‑1 (MQ=255)
                                                tcagcCATGGCAACTAAAACGTACAGTACCCAGTAGTTCCAGCCAGAGGCGAAACCGGCAAAACTGccc                        >  1:3515149/1‑69 (MQ=255)
                                                  agcCATGGCAACTAAAACGTACAGTACCCAGTAGTTCCAGCCAGAGGCGAAACCGGCAAAACTGCCCCAg                     >  1:1020786/1‑70 (MQ=255)
                                                    ccATGGCAACTAAAACGTACAGTACCCAGTAGTTCCAGCCAGAGGCGAAACCGGCaa                                <  1:3992337/57‑1 (MQ=255)
                                                      aTGGCAACTAAAACGTACAGTACCCAGTAGTTCCAGCCAGAGGCGAAACCGGCAAAACTGCCCCAGTAtt                 >  1:1199014/1‑70 (MQ=255)
                                                       tGGCAACTAAAACGTACAGTACCCAGTAGTTCCAGCCAGAGGCGAAACCGGCaaa                               <  1:4004223/55‑1 (MQ=255)
                                                       tGGCAACTAAAACGTACAGTACCCAGTAGTTCCAGCCAGAGGCGAAACCGGCAAAACTGCCCCATTAttt                >  1:203800/1‑70 (MQ=255)
                                                         gCAACTAAAACGTACAGTACCCAGTAGTTCCAGCCAGAGGCGAAACCGGCAAAACTGCCCCAg                     >  1:2303741/1‑63 (MQ=255)
                                                         gCAACTAAAACGTACAGTACCCAGTAGTTCCAGCCAGAGGCGAAACCGGCAAAACTGCCCCAGTATTtat              >  1:596156/1‑70 (MQ=255)
                                                         gCAACTAAAACGTACAGTACCCAGTAGTTCCAGCCAGAGGCGAAACCGGCAAAACTGCCCCAGTATTtat              >  1:2609106/1‑70 (MQ=255)
                                                           aaCTAAAACGTACAGTACCCAGTAGTTCCAGCCAGAGGCGAAACCGGCaa                                <  1:1822574/50‑1 (MQ=255)
                                                           aaCTAAAACGTACAGTACCCAGTAGTTCCAGCCAGAGGCGAAACCGGCaa                                <  1:3555752/50‑1 (MQ=255)
                                                             cTAAAACGTACAGTACCCAGTAGTTCCAGCCAGAGGCGAAACCGGCaa                                <  1:3610326/48‑1 (MQ=255)
                                                              tAAAACGTACAGTACCCAGTAGTTCCAGCCAGAGGCGAAACCGGCaa                                <  1:2029129/47‑1 (MQ=255)
                                                              tAAAACGTACAGTACCCAGTAGTTCCAGCCAGAGGCGAAACCGGCaa                                <  1:2499720/47‑1 (MQ=255)
                                                              tAAAACGTACAGTACCCAGTAGTTCCAGCCAGAGGCGAAACCGGCAAAACTGCCCCa                      >  1:2164309/1‑57 (MQ=255)
                                                                aaaCGTACAGTACCCAGTAGTTCCAGCCAGAGGCGAAACCGGCAAAACTGCCCCAGTATTTATAAGCAAAg      >  1:3323123/1‑71 (MQ=255)
                                                                 caCGAACATTACCCAGTAGTTCCAGCCAGAGAGGAAGCCCGCa                                 <  1:3176807/42‑1 (MQ=16)
                                                                 aaCGTACAGTACCCAGTAGTTCCAGCCAGAGGCGAAACCGGCaaaa                              >  1:1475965/1‑46 (MQ=255)
                                                                  aCGTACAGTACCCAGTAGTTCCAGCCAGAGGCGAAACCGGCAAAACTGCCCCAGTATTTATAAGCAAAGTg    <  1:106250/71‑1 (MQ=255)
                                                                  aCGTACAGTACCCAGTAGTTCCAGCCAGAGGCGAAACCGGCAAAACTGCCCCAGTATTTATAAGCAAAGTg    <  1:80708/71‑1 (MQ=255)
                                                                    gTACAGTACCCAGTAGTTCCAGCCAGAGGCGAAACCGGCAAAACTGCCCCAGTATTTATAAGCAAAGTGGc  >  1:4040760/1‑71 (MQ=255)
                                                                     |                                                                     
GGGATTTCCGGATACCAGAACTGAATGTATTTACCCACGGCAGTCAGCTCAGCCATGGCAACTAAAACGTACAGTACCCAGTAGTTCCAGCCAGAGGCGAAACCGGCAAAACTGCCCCAGTATTTATAAGCAAAGTGGC  >  minE/105915‑106053

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: