Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 330,103 G→A 43.8% G535G (GGC→GGT ybaL ← predicted transporter with NAD(P)‑binding Rossmann‑fold domain

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE330,1030GA43.8% 22.0 / 62.2 82G535G (GGC→GGTybaLpredicted transporter with NAD(P)‑binding Rossmann‑fold domain
Reads supporting (aligned to +/- strand):  ref base G (32/14);  new base A (13/23);  total (45/37)
Fisher's exact test for biased strand distribution p-value = 3.59e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAAT  >  minE/330037‑330161
                                                                  |                                                          
cGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGcc                                                          <  1:2496977/69‑1 (MQ=255)
cGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGcc                                                          <  1:1665315/69‑1 (MQ=255)
cGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGcc                                                          <  1:1272331/69‑1 (MQ=255)
cGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGcc                                                          <  1:2568494/69‑1 (MQ=255)
cGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGcc                                                          <  1:138649/69‑1 (MQ=255)
cGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGcc                                                          <  1:143188/69‑1 (MQ=255)
 gTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGcc                                                          <  1:75724/68‑1 (MQ=255)
       accTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCAc                                                    >  1:1842343/1‑68 (MQ=255)
       accTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACt                                                   >  1:1975736/1‑69 (MQ=255)
       accTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACt                                                   >  1:286429/1‑69 (MQ=255)
       accTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACt                                                   >  1:2730085/1‑69 (MQ=255)
       accTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACt                                                   >  1:3034461/1‑69 (MQ=255)
       accTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACt                                                   >  1:1574012/1‑69 (MQ=255)
       accTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACt                                                   >  1:151846/1‑69 (MQ=255)
       accTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACt                                                   >  1:3122226/1‑69 (MQ=255)
       accTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACt                                                   >  1:1405140/1‑69 (MQ=255)
       accTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACt                                                   >  1:701007/1‑69 (MQ=255)
       accTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACt                                                   >  1:1142553/1‑69 (MQ=255)
       accTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACt                                                   >  1:1131793/1‑69 (MQ=255)
       accTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATAGTACGGGCAATTTCACGCTCACCCATCACt                                                   <  1:3372513/69‑1 (MQ=255)
       accTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATAGTACGGGCAATTTCACGCTCACCCATCACt                                                   <  1:2144532/69‑1 (MQ=255)
       accTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATAGTACGGGCAATTTCACGCTCACCCATCACt                                                   <  1:1533423/69‑1 (MQ=255)
       accTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATAGTACGGGCAATTTCACGCTCACCCATCACt                                                   <  1:158595/69‑1 (MQ=255)
       accTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATAGTACGGGCAATTTCACGCTCACCCATCACt                                                   <  1:73868/69‑1 (MQ=255)
       accTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATAGTACGGGCAATTTCACGCTCACCCATCACt                                                   <  1:94123/69‑1 (MQ=255)
        ccTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATAGTACGGGCAATTTCACGCTCACCCATCACt                                                   <  1:1526361/68‑1 (MQ=255)
         cTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCAt                                                       >  1:191140/1‑63 (MQ=255)
         cTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTAc                                                 >  1:2862289/1‑69 (MQ=255)
         cTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTAc                                                 >  1:2845257/1‑69 (MQ=255)
         cTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTAc                                                 >  1:2950079/1‑69 (MQ=255)
         cTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTAc                                                 >  1:1012258/1‑69 (MQ=255)
         cTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTAc                                                 >  1:630373/1‑69 (MQ=255)
         cTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTAc                                                 >  1:994613/1‑69 (MQ=255)
         cTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTAc                                                 >  1:2393704/1‑69 (MQ=255)
         cTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTAc                                                 >  1:1493923/1‑69 (MQ=255)
         cTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTAc                                                 >  1:15031/1‑69 (MQ=255)
            aCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATAGTACGGGCAATTTCACGCTCACCCATCACt                                                   <  1:679639/64‑1 (MQ=255)
                 ccgGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCg                                         <  1:3420763/69‑1 (MQ=255)
                 ccgGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCg                                         <  1:1228569/69‑1 (MQ=255)
                             ccAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCg                                         <  1:347468/57‑1 (MQ=255)
                                gcagTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGa                          <  1:1046596/69‑1 (MQ=255)
                                gcagTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGa                          <  1:2922782/69‑1 (MQ=255)
                                      ccAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATAc                     >  1:3386311/1‑68 (MQ=255)
                                      ccAGCATAGTACGGGCAATTTCACGCTCACCCATCACTACCTGAtt                                           >  1:1638598/1‑46 (MQ=38)
                                      ccAGCATAGTACGGGCAATTTCACGCTCACCCATCACTACCTGATTCGCACCACGTTCGGTGAtata                      >  1:2756057/1‑67 (MQ=255)
                                      ccAGCATAGTACGGGCAATTTCACGCTCACCCATCACTACCTGATTCGCACCACGTTCGGTGATATAc                     >  1:991835/1‑68 (MQ=255)
                                      ccAGCATAGTACGGGCAATTTCACGCTCACCCATCACTACCTGATTCGCACCACGTTCGGTGATATAc                     >  1:2543979/1‑68 (MQ=255)
                                      ccAGCATAGTACGGGCAATTTCACGCTCACCCATCACTACCTGATTCGCACCACGTTCGGTGATATAc                     >  1:130077/1‑68 (MQ=255)
                                      ccAGCATAGTACGGGCAATTTCACGCTCACCCATCACTACCTGATTCGCACCACGTTCGGTGATATAc                     >  1:285523/1‑68 (MQ=255)
                                      ccAGCATAGTACGGGCAATTTCACGCTCACCCATCACTACCTGATTCGCACCACGTTCGGTGATATAc                     >  1:201030/1‑68 (MQ=255)
                                      ccAGCATAGTACGGGCAATTTCACGCTCACCCATCACTACCTGATTCGCACCACGTTCGGTGATATAc                     >  1:1986958/1‑68 (MQ=255)
                                      ccAGCATAGTACGGGCAATTTCACGCTCACCCATCACTACCTGATTCGCACCACGTTCGGTGATATAc                     >  1:288622/1‑68 (MQ=255)
                                      ccAGCATAGTACGGGCAATTTCACGCTCACCCATCACTACCTGATTCGCACCACGTTCGGTGATATAc                     >  1:414879/1‑68 (MQ=255)
                                      ccAGCATAGTACGGGCAATTTCACGCTCACCCATCACTACCTGATTCGCACCACGTTCGGTGATATAc                     >  1:2616492/1‑68 (MQ=255)
                                      ccAGCATAGTACGGGCAATTTCACGCTCACCCATCACTACCTGATTCGCACCACGTTCGGTGATATAc                     >  1:3080263/1‑68 (MQ=255)
                                      ccAGCATAGTACGGGCAATTTCACGCTCACCCATCACTACCTGATTAGCACCACGTTCGGTGATATAc                     >  1:918687/1‑68 (MQ=255)
                                        aGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGc                   <  1:351258/68‑1 (MQ=255)
                                              gTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATAc                     >  1:2448074/1‑60 (MQ=255)
                                              gTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCGCTTCa            >  1:757166/1‑69 (MQ=255)
                                              gTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCGCTTCa            >  1:753197/1‑69 (MQ=255)
                                              gTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTc             >  1:736606/1‑68 (MQ=255)
                                              gTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCa            >  1:794266/1‑69 (MQ=255)
                                              gTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCa            >  1:1005477/1‑69 (MQ=255)
                                              gTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCa            >  1:2420887/1‑69 (MQ=255)
                                              gTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCa            >  1:2893188/1‑69 (MQ=255)
                                              gTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCa            >  1:2792814/1‑69 (MQ=255)
                                                     cGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCAt     <  1:953057/69‑1 (MQ=255)
                                                        tttCACGCTCACCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAAt  <  1:2746077/69‑1 (MQ=255)
                                                        tttCACGCTCACCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAAt  <  1:2493007/69‑1 (MQ=255)
                                                        tttCACGCTCACCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAAt  <  1:2299957/69‑1 (MQ=255)
                                                        tttCACGCTCACCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAAt  <  1:2227964/69‑1 (MQ=255)
                                                        tttCACGCTCACCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAAt  <  1:2676155/69‑1 (MQ=255)
                                                        tttCACGCTCACCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAAt  <  1:1191365/69‑1 (MQ=255)
                                                        tttCACGCTCACCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAAt  <  1:2688208/69‑1 (MQ=255)
                                                        tttCACGCTCACCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAAt  <  1:3284933/69‑1 (MQ=255)
                                                        tttCACGCTCACCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAAt  <  1:2024178/69‑1 (MQ=255)
                                                        tttCACGCTCACCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAAt  <  1:465577/69‑1 (MQ=255)
                                                        tttCACGCTCACCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAAt  <  1:1629143/69‑1 (MQ=255)
                                                        tttCACGCTCACCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAAt  <  1:3153606/69‑1 (MQ=255)
                                                        tttCACGCTCACCCATCACTACCTGATTCGCACCACGTTCGGGGATATACGCCACTTCATCGTCATAAt  <  1:1299079/69‑1 (MQ=255)
                                                         ttCACGCTCACCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAAt  <  1:20899/68‑1 (MQ=255)
                                                         ttCACGCTCACCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAAt  <  1:2799345/68‑1 (MQ=255)
                                                                  |                                                          
CGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAAT  >  minE/330037‑330161

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: