Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 105,988 G→T 7.4% L102M (CTG→ATG)  aroP ← aromatic amino acid transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE105,9880GT 7.4% 217.0 / 2.8 67L102M (CTG→ATG) aroParomatic amino acid transporter
Reads supporting (aligned to +/- strand):  ref base G (33/29);  new base T (2/3);  total (35/32)
Fisher's exact test for biased strand distribution p-value = 6.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.86e-01

ATTTCCGGATACCAGAACTGAATGTATTTACCCACGGCAGTCAGCTCAGCCATGGCAACTAAAACGTACAGTACCCAGTAGTTCCAGCCAGAGGCGAAACCGGCAAAACTGCCCCAGTATTTATAAGCAAAGTGGCTAA  >  minE/105918‑106056
                                                                      |                                                                    
aTTTCCGGATACCAGAACTGAATGTATTTACCCACGGCAGTCAGCTCAGCCATGGCAACTAAAACGTACAg                                                                      >  1:3642768/1‑71 (MQ=255)
 tttCCGGATACCAGAACTGAATGTATTTACCCACGGCAGTCAGCTCAGCCATGGCAACTAAAACGTACAGt                                                                     <  1:762276/71‑1 (MQ=255)
 tttCCGGATACCAGAACTGAATGTATTTACCCACGGCAGTCAGCTCAGCCATGCCAACTAAAACGTACAg                                                                      >  1:2258200/1‑70 (MQ=255)
    ccGGATACCAGAACTGAATGTATTTACCCACGGCAGTCAGCTCAGCCATGGCAACTAAAACGTACAGTAc                                                                   <  1:1621756/70‑1 (MQ=255)
    ccGGATACCAGAACTGAATGTATTTACCCACGGCAGTCAGCTCAGCCATGGCAACTAAAACGTACAGTAc                                                                   >  1:2371653/1‑70 (MQ=255)
         tACAAGAACTGAATGTATTTACCCACGGCAGTCAGCTCAGCCATGGCAACTAAAACGTACAGTa                                                                    >  1:78304/1‑64 (MQ=255)
             aGAACTGAATGTATTTACCCACGGCAGTCAGCTCAGCCATGGCAACTAAAACGTACAGTACCCagtagta                                                          >  1:1036227/1‑69 (MQ=255)
              gAACTGAATGTATTTACCCACGGCAGTCAGCTCAGCCATGGCAACTAAAACGTACAGTACCCag                                                               <  1:4089615/64‑1 (MQ=255)
              gAACTGAATGTATTTACCCACGGCAGTCAGCTCAGCCATGGCAACTAAAACGTACAGTACCCAGTAGTTc                                                         >  1:1210586/1‑70 (MQ=255)
                    aaagTATTTACCCACGGCAGTCAGCTCAGCCATGGCAACTAAAACGTACAGTACCCAGTAGTTccagccag                                                  <  1:3196029/68‑1 (MQ=255)
                      tGTATTTACCCACGGCAGTCAGCTCAGCCATGGCAACTAAAACGTACAGTACCCAGTAGTTCCAGCCagag                                                <  1:2176903/71‑1 (MQ=255)
                       gTATTTACCCACGGCAGTCAGCTCAGCCATGGCAACTAAAACGTACAGTAc                                                                   <  1:434380/51‑1 (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)
                          tttACCCACGGCAGTCAGCTCAGCCATGGCAACTAAAACGTACAGTACCCagta                                                             >  1:2671113/1‑54 (MQ=255)
                            tACCCACGGCAGTCAGCTCAGCCATGGCAACTAAAACGTACTGTACCCAGTAGTTCCAGCCAGAGGCGaaa                                          >  1:4048448/1‑71 (MQ=255)
                            tACCCACGGCAGTCAGCTCAGCCATGGCAACTAAAACGTACAGTAccc                                                                 >  1:157566/1‑48 (MQ=255)
                            tACCCACGGCAGTCAGCTCAGCCATGGCAACTAAAACGTACAGTAccc                                                                 >  1:3524614/1‑48 (MQ=255)
                            tACCCACGGCAGTCAGCTCAGCCATGGCAACTAAAACGTACAGTAccc                                                                 >  1:2281605/1‑48 (MQ=255)
                            tACCCACGGCAGTCAGCTCAGCCATGGCAACTAAAACGTACAGTAccc                                                                 >  1:152817/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)
                                 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:55725/1‑42 (MQ=255)
                                          agctcagcCATGGCAACTAAAACGTACAGTACCCAGTAGTTc                                                         >  1:21301/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:2609106/1‑70 (MQ=255)
                                                      gCAACTAAAACGTACAGTACCCAGTAGTTCCAGCCAGAGGCGAAACCGGCAAAACTGCCCCAGTATTtat                 >  1:596156/1‑70 (MQ=255)
                                                        aaCTAAAACGTACAGTACCCAGTAGTTCCAGCCAGAGGCGAAACCGGCaa                                   <  1:3555752/50‑1 (MQ=255)
                                                        aaCTAAAACGTACAGTACCCAGTAGTTCCAGCCAGAGGCGAAACCGGCaa                                   <  1:1822574/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)
                                                                   aCATTACCCAGTAGTTCCAGCCAGAGAGGAAGCCCGCAAACGGTCCCCAGTATTTATAGGCAAAGTgggca   <  1:962497/71‑4 (MQ=255)
                                                                    cAGTACCCAGTAGTTCCAGCCAGAGGCGAAACCGGCAAAACTGCCCCAGTATTTATAAGCAAAGTGGCTaa  <  1:2314086/71‑1 (MQ=255)
                                                                    cAGTACCCAGTAGTTCCAGCCAGAGGCGAAACCGGCAAAACTGCCCCAGTATTTATAAGCAAAGTGGCTaa  <  1:3070642/71‑1 (MQ=255)
                                                                      |                                                                    
ATTTCCGGATACCAGAACTGAATGTATTTACCCACGGCAGTCAGCTCAGCCATGGCAACTAAAACGTACAGTACCCAGTAGTTCCAGCCAGAGGCGAAACCGGCAAAACTGCCCCAGTATTTATAAGCAAAGTGGCTAA  >  minE/105918‑106056

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: