Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 762,995 A→G 43.3% L235L (CTT→CTC potB ← polyamine transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE762,9950AG43.3% 20.5 / 64.2 83L235L (CTT→CTCpotBpolyamine transporter subunit
Reads supporting (aligned to +/- strand):  ref base A (22/25);  new base G (11/25);  total (33/50)
Fisher's exact test for biased strand distribution p-value = 1.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.67e-01

ATCAGGCCCATTACGATAGTCAGCGTAATGCTGGTAGCTGCACCAAACGGCCAGTCACGAATATTAAGGAACTGGACCTTGATGACGTTACCGATCAGCAGGTTTTTCGCACCGCCCATCAGGTCGGATA  >  minE/762930‑763059
                                                                 |                                                                
aTCAGGCCCATTACGATAGTCAGCGTAATGCTGGTAGCTGCACCAAACGGCCAGTCACGAATATTAAgg                                                               <  1:525168/69‑1 (MQ=255)
           tACGATAGTCAGCGTAATGCTGGTAGCTGCACCAAACGGCCAGTCACGAATATTAAGTAACTGGACCtt                                                    >  1:1873760/1‑69 (MQ=255)
           tACGATAGTCAGCGTAATGCTGGTAGCTGCACCAAACGGCCAGTCACGAATATTAAGGAACTGGACCtt                                                    >  1:1101715/1‑69 (MQ=255)
           tACGATAGTCAGCGTAATGCTGGTAGCTGCACCAAACGGCCAGTCACGAATATTAAGGAACTGGACCtt                                                    >  1:971326/1‑69 (MQ=255)
           tACGATAGTCAGCGTAATGCTGGTAGCTGCACCAAACGGCCAGTCACGAATATTAAGGAACTGGACCtt                                                    >  1:208059/1‑69 (MQ=255)
           tACGATAGTCAGCGTAATGCTGGTAGCTGCACCAAACGGCCAGTCACGAATATTAAGGAACTGGACCtt                                                    >  1:1321942/1‑69 (MQ=255)
           tACGATAGTCAGCGTAATGCTGGTAGCTGCACCAAACGGCCAGTCACGAATATTAAGGAACTGGACCtt                                                    >  1:2148092/1‑69 (MQ=255)
           tACGATAGTCAGCGTAATGCTGGTAGCTGCACCAAACGGCCAGTCACGAATATTAAGGAACTGGACCtt                                                    >  1:881198/1‑69 (MQ=255)
            aCGATAGTCAGCGTAATGCTGGTGGCTGCACCAAACGGCCAGTCACGAATATTGAGGAACTGTACCt                                                     >  1:1656989/1‑67 (MQ=255)
            aCGATAGTCAGCGTAATGCTGGTGGCTGCACCAAACGGCCAGTCACGAATATTGAGGAACTGTACCTtg                                                   >  1:2292157/1‑69 (MQ=255)
            aCGATAGTCAGCGTAATGCTGGTGGCTGCACCAAACGGCCAGTCACGAATATTGAGGAACTGTACCTtg                                                   >  1:427283/1‑69 (MQ=255)
            aCGATAGTCAGCGTAATGCTGGTGGCTGCACCAAACGGCCAGTCACGAATATTGAGGAACTGTACCTtg                                                   >  1:771409/1‑69 (MQ=255)
            aCGATAGTCAGCGTAATGCTGGTGGCTGCACCAAACGGCCAGTCACGAATATTGAGGAACTGTACCTtg                                                   >  1:1585307/1‑69 (MQ=255)
            aCGATAGTCAGCGTAATGCTGGTGGCTGCACCAAACGGCCAGTCACGAATATTGAGGAACTGTACCTtg                                                   >  1:2824321/1‑69 (MQ=255)
            aCGATAGTCAGCGTAATGCTGGTGGCTGCACCAAACGGCCAGTCACGAATATTGAGGAACTGTACCTtg                                                   >  1:2810771/1‑69 (MQ=255)
            aCGATAGTCAGCGTAATGCTGGTGGCTGCACCAAACGGCCAGTCACGAATATTGAGGAACTGTACCTtg                                                   >  1:2626528/1‑69 (MQ=255)
            aCGATAGTCAGCGTAATGCTGGTGGCTGCACCAAACGGCCAGTCACGAATATTGAGGAACTGTACCTtg                                                   >  1:1409288/1‑69 (MQ=255)
            aCGATAGTCAGCGTAATGCTGGTGGCTGCACCAAACGGCCAGTCACGAATATTGAGGAACTGTACCTtg                                                   >  1:2352868/1‑69 (MQ=255)
            aCGATAGTCAGCGTAATGCTGGTGGCTGCACCAAACGGCCAGTCACGAATATTGAGGAACTGTACCTtg                                                   >  1:1281796/1‑69 (MQ=255)
            aCGATAGTCAGCGTAATGCTGGTAGCTGCACCAAACGGCCAGTCACGAATATTAAGGAACTGGACCTtg                                                   >  1:1718994/1‑69 (MQ=255)
            aCGATAGTCAGCGTAATGCTGGTAGCTGCACCAAACGGCCAGTCACGAATATTAAGGAACTGGACCTtg                                                   >  1:1224415/1‑69 (MQ=255)
            aCGATAGTCAGCGTAATGCTGGTAGCTGCACCAAACGGCCAGTCACGAATATTAAGGAACTGGACCTtg                                                   >  1:1214308/1‑69 (MQ=255)
               aTAGTCAGCGTAATGCTGGTAGCTGCACCAAACGGCCAGTCACGAATATTAAGGAACTGGACCTtgat                                                 <  1:2443604/68‑1 (MQ=255)
                    cAGCGTAATGCTGGTAGCTGCACCAAACGGCCAGTCACGAATATTAAGGAACTGGACCTTGATGACGt                                            <  1:996980/68‑1 (MQ=255)
                           aTGCTGGTAGCTGCACCAAACGGCCAGTCACGAATATTAAGGAACTGGACCTtgatg                                                >  1:3329861/1‑57 (MQ=255)
                           aTGCTGGTAGCTGCACCAAACGGCCAGTCACGAATATTAAGGAACTGGACCTTGATGACGTTACCGATc                                    >  1:710423/1‑69 (MQ=255)
                             gCTGGTAGCTGCACCAAACGGCCAGTCACGAATATTAAGGAACTGGACCTTGATGACGTTACCGAt                                     >  1:1583644/1‑66 (MQ=255)
                             gCTGGTAGCTGCACCAAACGGCCAGTCACGAATATTAAGGAACTGGACCTTGATGACGTTACCGATc                                    >  1:1359435/1‑67 (MQ=255)
                             gCTGGTAGCTGCACCAAACGGCCAGTCACGAAAATTAAGGCACTGGACCTTGATGACGTTAGCGAt                                     >  1:148861/1‑66 (MQ=255)
                             gCTGGTAGCCGCACCAAACGGCCAGTCACGAATATTAAGGAACTGGACCTTGATGACGTTACCGATc                                    >  1:1540217/1‑67 (MQ=255)
                                      tGCACCAAACGGCCAGTCACGAATATTGAGGAACTTTACCTTGATGACGTTACCGATCAGCAGGtttt                          <  1:2581084/68‑1 (MQ=255)
                                      tGCACCAAACGGCCAGTCACGAATATTGAGGAACTGTACCTTGATGACGTTACCGATCAGCAGGtttt                          <  1:243389/68‑1 (MQ=255)
                                      tGCACCAAACGGCCAGTCACGAATATTGAGGAACTGTACCTTGATGACGTTACCGATCAGCAGGtttt                          <  1:46268/68‑1 (MQ=255)
                                      tGCACCAAACGGCCAGTCACGAATATTGAGGAACTGTACCTTGATGACGTTACCGATCAGCAGGtttt                          <  1:980477/68‑1 (MQ=255)
                                      tGCACCAAACGGCCAGTCACGAATATTGAGGAACTGTACCTTGATGACGTTACCGATCAGCAGGtttt                          <  1:246672/68‑1 (MQ=255)
                                      tGCACCAAACGGCCAGTCACGAATATTGAGGAACTGTACCTTGATGACGTTACCGATCAGCAGGtttt                          <  1:2584195/68‑1 (MQ=255)
                                      tGCACCAAACGGCCAGTCACGAATATTGAGGAACTGTACCTTGATGACGTTACCGATCAGCAGGtttt                          <  1:780524/68‑1 (MQ=255)
                                      tGCACCAAACGGCCAGTCACGAATATTGAGGAACTGTACCTTGATGACGTTACCGATCAGCAGGtttt                          <  1:274624/68‑1 (MQ=255)
                                      tGCACCAAACGGCCAGTCACGAATATTGAGGAACTGTACCTTGATGACGTTACCGATCAGCAGGtttt                          <  1:768763/68‑1 (MQ=255)
                                      tGCACCAAACGGCCAGTCACGAATATTGAGGAACTGTACCTTGATGACGTTACCGATCAGCAGGtttt                          <  1:2853360/68‑1 (MQ=255)
                                      tGCACCAAACGGCCAGTCACGAATATTGAGGAACTGTACCTTGATGACGTTACCGATCAGCAGGtttt                          <  1:3141984/68‑1 (MQ=255)
                                      tGCACCAAACGGCCAGTCACGAATATTGAGGAACTGTACCTTGATGACGTTACCGATCAGCAGGtttt                          <  1:349257/68‑1 (MQ=255)
                                      tGCACCAAACGGCCAGTCACGAATATTGAGGAACTGTACCTTGATGACGTTACCGATCAGCAGGtttt                          <  1:579471/68‑1 (MQ=255)
                                      tGCACCAAACGGCCAGTCACGAATATTGAGGAACTGTACCTTGATGACGTTACCGATCAGCAGGtttt                          <  1:202234/68‑1 (MQ=255)
                                      tGCACCAAACGGCCAGTCACGAATATTGAGGAACTGTACCTTGATGACGTTACCGATCAGCAGGtttt                          <  1:1388872/68‑1 (MQ=255)
                                      tGCACCAAACGGCCAGTCACGAATATTGAGGAACTGTACCTTGATGACGTTACCGATCAGCAGGtttt                          <  1:2128788/68‑1 (MQ=255)
                                      tGCACCAAACGGCCAGTCACGAATATTGAGGAACTGTACCTTGATGACGTTACCGATCAGCAGGtttt                          <  1:2152580/68‑1 (MQ=255)
                                      tGCACCAAACGGCCAGTCACGAATATTGAGGAACTGTACCTTGATGACGTTACCGATCAGCAGGtttt                          <  1:1459685/68‑1 (MQ=255)
                                      tGCACCAAACGGCCAGTCACGAATATTGAGGAACTGTACCTTGATGACGTTACCGATCAGCAGGtttt                          <  1:179824/68‑1 (MQ=255)
                                      tGCACCAAACGGCCAGTCACGAATATTGAGGAACTGTACCTTGATGACGTTACCGATCAGCAGGtttt                          <  1:1855843/68‑1 (MQ=255)
                                                 gCCAGTCACGAATATTGAGGAACTGTACCTTGATGACGTTACCGATCAGCAGGtttt                          <  1:1769707/57‑1 (MQ=255)
                                                  ccAGTCACGAATATTGAGGAACTGTACCTTGATGACGTTACCGATCAGCAGGTTTTTCGCACCGCCCAt             <  1:295305/69‑1 (MQ=255)
                                                  ccAGTCACGAATATTAAGGAACTGGACCTTGATGACGTTACCGATCAGCAGGTTTTTCGCACCGCCCAt             <  1:3336006/69‑1 (MQ=255)
                                                   cAGTCACGAATATTGAGGAACTGTACCTTGATGACGTTACCGATCAGCAGGTTTTTCGCACCGCCCAt             <  1:1154501/68‑1 (MQ=255)
                                                      tCACGAATATTGAGCAACTGTACCTTGATGACGTTACCGATCAGCAGGtttt                          <  1:2964990/52‑1 (MQ=255)
                                                       cACGAATATTAAGGAACTGGACCTTGATGACGTTACCGATCAGCAGGTTTTTCGCACCGCCCATCAGGt        >  1:2602473/1‑69 (MQ=255)
                                                       cACGAATATTAAGGAACTGGACCTTGATGACGTTACCGATCAGCAGGTTTTTCGCACCGCCCATCAGGt        >  1:2381597/1‑69 (MQ=255)
                                                       cACGAATATTAAGGAACTGGACCTTGATGACGTTACCGATCAGCAGGTTTTTCGCACCGCCCATCAGGt        >  1:3100910/1‑69 (MQ=255)
                                                       cACGAATATTAAGGAACTGGACCTTGATGACGTTACCGATCAGCAGGTTTTTCGCACCGCCCATCAGGt        >  1:249945/1‑69 (MQ=255)
                                                          gAATATTAAGGAACTGGACCTTGATGACGTTACCGATCAGCAGGTTTTTCGCACCGCCCATCAGGTCg      <  1:1838790/68‑1 (MQ=255)
                                                          gAATATTAAGGAACTGGACCTTGATGACGTTACCGATCAGCAGGTTTTTCGCACCGCCCATCAGGTCg      <  1:806390/68‑1 (MQ=255)
                                                          gAATATTAAGGAACTGGACCTTGATGACGTTACCGATCAGCAGGTTTTTCGCACCGCCCATCAGGTCg      <  1:1859984/68‑1 (MQ=255)
                                                          gAATATTAAGGAACTGGACCTTGATGACGTTACCGATCAGCAGGTTTTTCGCACCGCCCATCAGGTCg      <  1:1970421/68‑1 (MQ=255)
                                                          gAATATTAAGGAACTGGACCTTGATGACGTTACCGATCAGCAGGTTTTTCGCACCGCCCATCAGGTCg      <  1:3244024/68‑1 (MQ=255)
                                                          gAATATTAAGGAACTGGACCTTGATGACGTTACCGATCAGCAGGTTTTTCGCACCGCCCATCAGGTCg      <  1:2682553/68‑1 (MQ=255)
                                                          gAATATTAAGGAACTGGACCTTGATGACGTTACCGATCAGCAGGTTTTTCGCACCGCCCATCAGGTCg      <  1:2448113/68‑1 (MQ=255)
                                                          gAATATTAAGGAACTGGACCTTGATGACGTTACCGATCAGCAGGTTTTTCGCACCGCCCATCAGGTCg      <  1:2345247/68‑1 (MQ=255)
                                                          gAATATTAAGGAACTGGACCTTGATGACGTTACCGATCAGCAGGTTTTTCGCACCGCCCATCAGGTCg      <  1:2443325/68‑1 (MQ=255)
                                                          gAATATTAAGGAACTGGACCTTGATGACGTTACCGATCAGCAGGTTTTTCGCACCGCCCATCAGGTCg      <  1:2322956/68‑1 (MQ=255)
                                                          gAATATTAAGGAACTGGACCTTGATGACGTTACCGATCAGCAGGTTTTTCGCACCGCCCATCAGGTCg      <  1:2284584/68‑1 (MQ=255)
                                                          gAATATTAAGGAACTGGACCTTGATGACGTTACCGATCAGCAGGTTTTTCGCACCGCCCATCAGGTCg      <  1:2886690/68‑1 (MQ=255)
                                                          gAATATTAAGGAACTGGACCTTGATGACGTTACCGATCAGCAGGTTTTTCGCACCGCCCATCAGGTCg      <  1:3143075/68‑1 (MQ=255)
                                                          gAATATTAAGGAACTGGACCTTGATGACGTTACCGATCAGCAGGTTTTTCGCACCGCCCATCAGGACg      <  1:870835/68‑1 (MQ=255)
                                                           aaTATTAAGGAACTGGACCTTGATGACGTTACCGATCAGCAGGTTTTTCGc                      <  1:2377603/51‑1 (MQ=255)
                                                           aaTATTAAGGAACTGGACCTTGATGACGTTACCGATCAGCAGGTTTTTCGc                      <  1:2415294/51‑1 (MQ=255)
                                                           aaTATTAAGGAACTGGACCTTGATGACGTTACCGATCAGCAGGTTTTTCGc                      <  1:1266669/51‑1 (MQ=255)
                                                            atatTGAGGAACTGTACCTTGATGACGTTACCGATCAGCAGGtttt                          <  1:2028891/46‑1 (MQ=255)
                                                             tatTAAGGAACTGGACCTTGATGACGTTACCGATCAGCAGGTTTTTCGCACCGCCCATCAGGTCGGATa  <  1:2776768/69‑1 (MQ=255)
                                                             tatTAAGGAACTGGACCTTGATGACGTTACCGATCAGCAGGTTTTTCGCACCGCCCATCAGGTCGGATa  <  1:1691765/69‑1 (MQ=255)
                                                             tatTAAGGAACTGGACCTTGATGACGTTACCGATCAGCAGGTTTTTCGCACCGCCCATCAGGTCGGATa  <  1:1972669/69‑1 (MQ=255)
                                                             tatTAAGGAACTGGACCTTGATGACGTTACCGATCAGCAGGTTTTTCGCACCGCCCATCAGGTCGGATa  <  1:9963/69‑1 (MQ=255)
                                                               ttAAGGAACTGGACCTTGATGACGTTACCGATCAGCAGGTTTTTCGCACCGCCCATCAGGTCGGAt   >  1:1277580/1‑66 (MQ=255)
                                                               ttAAGGAACTGGACCTTGATGACGTTACCGATCAGCAGGTTTTTCGCACCGCCCATCAGGTCGGATa  >  1:446063/1‑67 (MQ=255)
                                                                 |                                                                
ATCAGGCCCATTACGATAGTCAGCGTAATGCTGGTAGCTGCACCAAACGGCCAGTCACGAATATTAAGGAACTGGACCTTGATGACGTTACCGATCAGCAGGTTTTTCGCACCGCCCATCAGGTCGGATA  >  minE/762930‑763059

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: