Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,099,149 A→G 35.3% T130A (ACA→GCA)  ydiM → predicted transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,099,1490AG35.3% 30.1 / 31.6 34T130A (ACA→GCA) ydiMpredicted transporter
Reads supporting (aligned to +/- strand):  ref base A (8/14);  new base G (8/4);  total (16/18)
Fisher's exact test for biased strand distribution p-value = 1.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATGGCAAACAGTTTTCTCGATGCAGGCACTTATCCCAGTTTGATGGAAGCTTTTCCACGCTCACCTGGGACAGCCAATATTTTAATTAAAGCATTTGTTTCCAGCGGACAATTTTTATTACCGCTAATCATTAGCCTGTT  >  minE/1099080‑1099219
                                                                     |                                                                      
aTGGCAAACAGTTTTCTCGATGCAGGCACTTATCCCAGTTTGATGGAAGCTTTTCCACGCTCACCTGGGa                                                                        >  1:401463/1‑70 (MQ=255)
   gCAAACAGTTTTCTCGATGCAGGCACTTATCCCAGTTTGATGGAAGCTTTTCCACGCTCACCTGGGACAGc                                                                    >  1:1574037/1‑71 (MQ=255)
        cAGTTTTCTCGATGCAGGCACTTATCCCAGTTTGATGGAAGCTTTTCCACGCTCACCTGGGGCAGCCAata                                                               >  1:486019/1‑71 (MQ=255)
            tttCTCGATGCAGGCACTTATCCCAGTTTGATGGAAGCTTTTCCACGCTCACCTGGGa                                                                        >  1:1864908/1‑58 (MQ=255)
             ttCTCGATGCAGGCACTTATCCCAGTTTGATGGAAGCTTTTCCACGCTCACCTGGGGCAGCCAATAttt                                                            >  1:816728/1‑69 (MQ=255)
             ttCTCGATGCAGGCACTTATCCCAGTTTGATGGAAGCTTTTCCACGCTCACCTGGGACAGCCAATAtttt                                                           <  1:2015904/70‑1 (MQ=255)
             ttCTCGATGCAGGCACTTATCCCAGTTTGATGGAAGCTTTTCCACGCTCACCTGGGACAGCCAATAtttt                                                           >  1:788291/1‑70 (MQ=255)
                  gATGCAGGCACTTATCCCAGTTTGATGGAAGCTTTTCCACGCTCACCTGGGACAGCCAATATTttaatta                                                      <  1:1670328/70‑1 (MQ=255)
                  gATGCAGGCACTTATCCCAGTTTGATGGAAGCTTTTCCACGCTCACCTGGGACAGCCAATATTttaat                                                        <  1:2015627/68‑1 (MQ=255)
                      cAGGCACTTATCCCAGTTTGATGGAAGCTTTTCCACGCTCACCTGGGGCAGCCAATATTTTAATTAAAGCa                                                 >  1:1305080/1‑71 (MQ=255)
                             ttATCCCAGTTTGATGGAAGCTTTTCCACGCTCACCTGGGACAGCCAATAttt                                                            <  1:2113900/53‑1 (MQ=255)
                              tATCCCAGTTTGATGGAAGCTTTTCCACGCTCACCTGGGGCAGCCAATATTTTAATTAAAGCATTTGTTTc                                         >  1:1268993/1‑71 (MQ=255)
                               aTCCCAGTTTGATGGAAGCTTTTCCACGCTCACCTGGGa                                                                        >  1:1967590/1‑39 (MQ=255)
                               aTCCCAGTTTGATGGAAGCTTTTCCACGCTCACCTGGGACAGCCAATATTTTAATTAAAGCATTTGtt                                           >  1:1639268/1‑68 (MQ=255)
                                 cccAGTTTGATGGAAGCTTTTCCACGCTCACCTGGGACAGCCAATAtttt                                                           <  1:1225101/50‑1 (MQ=255)
                                   cAGTTTGATGGAAGCTTTTCCACGCTCACCTGGGGCAGCCAATATTttaat                                                        <  1:1223115/51‑1 (MQ=255)
                                   cAGTTTGATGGAAGCTTTTCCACGCTCACCTGGGACAGCCAATATTTTAATTAAAGCATTTGTTTCCAGCg                                    <  1:934314/71‑1 (MQ=255)
                                   cAGTTTGATGGAAGCTTTTCCACGCTCACCTGGGACAGCCAATATTTTAATTAAAGCATTTGTTTCCAGCg                                    <  1:959988/71‑1 (MQ=255)
                                        tGATGGAAGCTTTTCCACGCTCACCTGGGGCAGCCAATATTTTAATTAAAGCATTTGTTTCCAGCGGACaa                               >  1:972015/1‑71 (MQ=255)
                                         gATGGAAGCTTTTCCACGCTCACCTGGGACAGCCAATATTTTAATTAAAGCATTTGTTTCCAGCGGACaa                               <  1:1424589/70‑1 (MQ=255)
                                         gATGGAAGCTTTTCCACGCTCACCTGGGACAGCCAATATTTTAATTAAAGCATTTGTTTCCAGCGGACaa                               <  1:1013221/70‑1 (MQ=255)
                                             gAAGCTTTTCCACGCTCACCTGGGGCAGCCAATATTTTAATTAAAGCATTTGTTTCCAGCGGACAAttttt                          <  1:244531/71‑1 (MQ=255)
                                             gAAGCTTTTCCACGCTCACCTGGGACAGCCAATATTTTAATTAAAGCATTTGTTTCCAGCGGACAAtttt                           <  1:1682098/70‑1 (MQ=255)
                                              aaGCTTTTCCACGCTCACCTGGGACAGCCAATATTTTAATTAAAGCa                                                 <  1:640863/47‑1 (MQ=255)
                                                gCTTTTCCACGCTCACCTGGGACAGCCAATATTTTAATTAAAGCAtt                                               >  1:2023182/1‑47 (MQ=255)
                                                gCTTTTCCACGCTCACCTGGGACAGCCAATATTTTAATTAAAGCATTTGTTTCCAGCGGACAAttttt                          <  1:310566/68‑1 (MQ=255)
                                                gCTTTTCCACGCTCACCTGGGACAGCCAATATTTTAATTAAAGCATTTGTTTCCAGCGGACAATTTttatt                       <  1:1769502/71‑1 (MQ=255)
                                                 cTTTTCCACGCTCACCTGGGACAGCCAATATTTTAATTAAAGCATTTGTTTCCAGCGGACAATTTttatta                      >  1:2077529/1‑71 (MQ=255)
                                                  ttttCCACGCTCACCTGGGGCAGCCAATATTTTAATTAAAGCATTTGTTTc                                         >  1:212140/1‑51 (MQ=255)
                                                   tttCCACGCTCACCTGGGGCAGCCAATATTTTAATTAAAGCATTTGTTTCCAGCGGACAATTTTTATTAcc                    >  1:985851/1‑71 (MQ=255)
                                                   tttCCACGCTCACCTGGGACAGCCAATATTTTAATTAAAGCATTTGTTTCCAGCGGACAATTTttatta                      >  1:189786/1‑69 (MQ=255)
                                                     tCCACGCTCACCTGGGGCAGCCAATATTTTAATTAAAGCATTTGTTTCCAGCGGa                                  >  1:1578519/1‑55 (MQ=255)
                                                              aCCTGGGGCAGCCAATATTTTAATTAAAGCATTTGTTTCCAGCGGACAATTTTTATTACCGCTAATCAtt          <  1:1121448/70‑1 (MQ=255)
                                                                cTGGGGCAGCCAATATTTTAATTAAAGCATTTGTTTCCAGCGGACAATTTTTATTACCGCTAATCATTAGc       <  1:1552936/71‑1 (MQ=255)
                                                                cTGGGACAGCCAATATTTTAATTAAAGCATTTGTTTCCAGCGGACAAtttt                           <  1:245048/51‑1 (MQ=255)
                                                                 tGGGACAGCCAATATTTTAATTAAAGCATTTGTTTCCAGCGGACAAtt                             >  1:1796450/1‑48 (MQ=255)
                                                                 tGGGACAGCCAATATTTTAATTAAAGCATTTGTTTCCAGCGGACAAtt                             >  1:1203213/1‑48 (MQ=255)
                                                                     aCAGCCAATATTTTAATTAAAGCATTTGTTTCCAGCGGACAATTTTTATTACCGCTAATAATTAGCCTGtt  <  1:1840837/71‑1 (MQ=255)
                                                                     |                                                                      
ATGGCAAACAGTTTTCTCGATGCAGGCACTTATCCCAGTTTGATGGAAGCTTTTCCACGCTCACCTGGGACAGCCAATATTTTAATTAAAGCATTTGTTTCCAGCGGACAATTTTTATTACCGCTAATCATTAGCCTGTT  >  minE/1099080‑1099219

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: