Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 331,222 A→G 24.9% V162V (GTT→GTC 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
*minE331,2220AG24.9% 102.8 / 31.3 88V162V (GTT→GTCybaLpredicted transporter with NAD(P)‑binding Rossmann‑fold domain
Reads supporting (aligned to +/- strand):  ref base A (18/48);  new base G (5/17);  total (23/65)
Fisher's exact test for biased strand distribution p-value = 7.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.86e-01

AGTGGCAAAGCCCACATCGCCCTGTTCCATCATTCCTGCCACTGCGGGCAGCAACACCAGCGTCAGAACCATTACCAGGTCTTCCACAATCAACCAACCGATGGCGATTTGCCCACGCTGACTGTCAATTA  >  minE/331156‑331286
                                                                  |                                                                
agTGGCAAAGCCCACATCGCCCTGTTCCATCATTCCTGCCACTGCGGGCAGCAACACCAGCGTCAGAAc                                                                <  1:1627926/69‑1 (MQ=255)
agTGGCAAAGCCCACATCGCCCTGTTCCATCATTCCTGCCACTGCGGGCAGCAACACCAGCGTCAGAAc                                                                <  1:1899591/69‑1 (MQ=255)
agTGGCAAAGCCCACATCGCCCTGTTCCATCATTCCTGCCACTGCGGGCAGCAACACCAGCGTCAGAAc                                                                <  1:2194958/69‑1 (MQ=255)
agTGGCAAAGCCCACATCGCCCTGTTCCATCATTCCTGCCACTGCGGGCAGCAACACCAGCGTCAGAAc                                                                <  1:1755538/69‑1 (MQ=255)
agTGGCAAAGCCCACATCGCCCTGTTCCATCATTCCTGCCACTGCGGGCAGCAACACCAGCGTCAGAAc                                                                <  1:2247438/69‑1 (MQ=255)
agTGGCAAAGCCCACATCGCCCTGTTCCATCATTCCTGCCACTGCGGGCAGCAACACCAGCGTCAGAAc                                                                <  1:1668880/69‑1 (MQ=255)
agTGGCAAAGCCCACATCGCCCTGTTCCATCATTCCTGCCACTGCGGGCAGCAACACCAGCGTCAGAAc                                                                <  1:1640162/69‑1 (MQ=255)
agTGGCAAAGCCCACATCGCCCTGTTCCATCATTCCTGCCACTGCGGGCAGCAACACCAGCGTCAGAAc                                                                <  1:2399584/69‑1 (MQ=255)
agTGGCAAAGCCCACATCGCCCTGTTCCATCATTCCTGCCACTGCGGGCAGCAACACCAGCGTCAGAAc                                                                <  1:384112/69‑1 (MQ=255)
agTGGCAAAGCCCACATCGCCCTGTTCCATCATTCCTGCCACTGCGGGCAGCAACACCAGCGTCAGAAc                                                                <  1:2782213/69‑1 (MQ=255)
agTGGCAAAGCCCACATCGCCCTGTTCCATCATTCCTGCCACTGCGGGCAGCAACACCAGCGTCAGAAc                                                                <  1:1417233/69‑1 (MQ=255)
agTGGCAAAGCCCACATCGCCCTGTTCCATCATTCCTGCCACTGCGGGCAGCAACACCAGCGTCAGAAc                                                                <  1:2950690/69‑1 (MQ=255)
agTGGCAAAGCCCACATCGCCCTGTTCCATCATTCCTGCCACTGCGGGCAGCAACACCAGCGTCAGAAc                                                                <  1:1260591/69‑1 (MQ=255)
agTGGCAAAGCCCACATCGCCCTGTTCCATCATTCCTGCCACTGCGGGCAGCAACACCAGCGTCAGAAc                                                                <  1:1128145/69‑1 (MQ=255)
agTGGCAAAGCCCACATCGCCCTGTTCCATCATTCCTGCCACTGCGGGCAGCAACACCAGCGTCAGAAc                                                                <  1:1106623/69‑1 (MQ=255)
 gTGGCAAAGCCCACATCGCCCTGTTCCATCATTCCTGCCACTGCGGGCAGCAACACCAGCGTCAGAAcc                                                               >  1:2421848/1‑69 (MQ=255)
 gTGGCAAAGCCCACATCGCCCTGTTCCATCATTCCTGCCACTGCGGGCAGCAACACCAGCGTCAGAAc                                                                <  1:561779/68‑1 (MQ=255)
   ggCAAAGCCCACATCGCCCTGTTCCATCATTCCTGCTACCGCGGGTAGCAATACCAGCGTCAGGACCAt                                                             >  1:2881164/1‑69 (MQ=255)
   ggCAAAGCCCACATCGCCCTGTTCCATCATTCCTGCTACCGCGGGTAGCAATACCAGCGTCAGGACCAt                                                             >  1:2107080/1‑69 (MQ=255)
   ggCAAAGCCCACATCGCCCTGTTCCATCATTCCTGCTACCGCGGGTAGCAATACCAGCGTCAGGACCAt                                                             >  1:3187285/1‑69 (MQ=255)
   ggCAAAGCCCACATCGCCCTGTTCCATCATTCCTGCTACCGCGGGTAGCAATACCAGCGTCAGGACCAt                                                             >  1:1674164/1‑69 (MQ=255)
   ggCAAAGCCCACATCGCCCTGTTCCATCATTCCTGCTACCGCGGGTAGCAATACCAGCGTCAGGACCAt                                                             >  1:22846/1‑69 (MQ=255)
      aaaGCCCACATCGCCCTGTTCCATCATTCCTGCCACTGCGGGCAGCAACACCAGCGTCAGAACCATTa                                                           <  1:3138667/68‑1 (MQ=255)
               aTCGCCCTGTTCCATCATTCCTGCCACTGCGGGCAGCAACACCAGCGTCAGAACCATTACCAGGTCTTc                                                 <  1:2912723/69‑1 (MQ=255)
               aTCGCCCTGTTCCATCATTCCTGCCACTGCGGGCAGCAACACCAGCGTCAGAACCATTACCAGGTCTTc                                                 <  1:511011/69‑1 (MQ=255)
               aTCGCCCTGTTCCATCATTCCTGCCACTGCGGGCAGCAACACCAGCGTCAGAACCATTACCAGGTCTTc                                                 <  1:107855/69‑1 (MQ=255)
                       gTTCCATCATTCCTGCCACTGCGGGCAGCAACACCAGCGTCAGAACCATTACCAGGTCTTCCACAATca                                         >  1:3159453/1‑69 (MQ=255)
                       gTTCCATCATTCCTGCCACTGCGGGCAGCAACACCAGCGTCAGAACCATTACCAGGTCTTCCACAATca                                         >  1:3172695/1‑69 (MQ=255)
                       gTTCCATCATTCCTGCCACTGCGGGCAGCAACACCAGCGTCAGAACCATTACCAGGTCTTCCACAATca                                         >  1:1442489/1‑69 (MQ=255)
                         tCCATCATTCCTGCCACTGCGGGCAGCAACACCAGCGTCAGAAc                                                                <  1:1725189/44‑1 (MQ=255)
                          ccATCATTCCTGCCACTGCGGGCAGCAACTCCAGCGTCAGAACCATTACCAGGTCTTCCACAATcaacc                                      <  1:2822889/69‑1 (MQ=255)
                          ccATCATTCCTGCCACTGCGGGCAGCAACACCAGCGTCAGAACCATTACCAGGTCTTCCTCAATcaacc                                      >  1:1100292/1‑69 (MQ=255)
                          ccATCATTCCTGCCACTGCGGGCAGCAACACCAGCGTCAGAACCATTACCAGGTCTTCCACAATcaagc                                      >  1:724166/1‑67 (MQ=255)
                          ccATCATTCCTGCCACTGCGGGCAGCAACACCAGCGTCAGAACCATTACCAGGTCTTCCACAATcaacc                                      <  1:1952123/69‑1 (MQ=255)
                          ccATCATTCCTGCCACTGCGGGCAGCAACACCAGCGTCAGAACCATTACCAGGTCTTCCACAATcaacc                                      >  1:275062/1‑69 (MQ=255)
                          ccATCATTCCTGCCACTGCGGGCAGCAACACCAGCGTCAGAACCATTACCAGGTCTTCCACAATcaacc                                      >  1:170691/1‑69 (MQ=255)
                          ccATCATTCCTGCCACTGCGGGCAGCAACACCAGCGTCAGAACCATTACCAGGTCTTCCACAATcaacc                                      <  1:2934221/69‑1 (MQ=255)
                          ccATCATTCCTGCCACTGCGGGCAGCAACACCAGCGTCAGAACCATTACCAGGTCTTCCACAATcaacc                                      <  1:51721/69‑1 (MQ=255)
                          ccATCATTCCTGCCACTGCGGGCAGCAACACCAGCGTCAGAACCATTACCAGGTCTTCCACAATcaacc                                      <  1:1296442/69‑1 (MQ=255)
                                     gCCACTGCGGTCAGCAACACCAGCGTCAGAACCATTACCAGGTCTTCCACAATCAACCAACCGATGGCg                           <  1:1970638/69‑1 (MQ=255)
                                     gCCACTGCGGGCAGCAACACCAGCGTCAGGACCATTACCAGGTCTTCCACAATCAACCAACCGATGGCg                           <  1:2313503/69‑1 (MQ=255)
                                     gCCACTGCGGGCAGCAACACCAGCGTCAGAACCATTACCAGGTCTTCCACAATCaaccaacc                                  <  1:3230003/62‑1 (MQ=255)
                                     gCCACTGCGGGCAGCAACACCAGCGTCAGAACCATTACCAGGTCTTCCACAATCaaccaacc                                  <  1:67764/62‑1 (MQ=255)
                                     gCCACTGCGGGCAGCAACACCAGCGTCAGAACCATTACCAGGTCTTCCACAATCaaccaacc                                  <  1:478540/62‑1 (MQ=255)
                                     gCCACTGCGGGCAGCAACACCAGCGTCAGAACCATTACCAGGTCTTCCACAATCaaccaacc                                  >  1:2078198/1‑62 (MQ=255)
                                     gCCACTGCGGGCAGCAACACCAGCGTCAGAACCATTACCAGGTCTTCCACAATCaaccaacc                                  <  1:2809508/62‑1 (MQ=255)
                                     gCCACTGCGGGCAGCAACACCAGCGTCAGAACCATTACCAGGTCTTCCACAATCaaccaacc                                  <  1:191972/62‑1 (MQ=255)
                                     gCCACTGCGGGCAGCAACACCAGCGTCAGAACCATTACCAGGTCTTCCACAATCaaccaacc                                  <  1:1089112/62‑1 (MQ=255)
                                     gCCACTGCGGGCAGCAACACCAGCGTCAGAACCATTACCAGGTCTTCCACAATCAACCAACCGATGGc                            <  1:3259317/68‑1 (MQ=255)
                                     gCCACTGCGGGCAGCAACACCAGCGTCAGAACCATTACCAGGTCTTCCACAATCAACCAACCGATGGc                            <  1:271006/68‑1 (MQ=255)
                                     gCCACTGCGGGCAGCAACACCAGCGTCAGAACCATTACCAGGTCTTCCACAATCAACCAACCGATGGCg                           <  1:1328558/69‑1 (MQ=255)
                                     gCCACTGCGGGCAGCAACACCAGCGTCAGAACCATTACCAGGTCTTCCACAATCAACCAACCGATGGCg                           <  1:1499726/69‑1 (MQ=255)
                                     gCCACTGCGGGCAGCAACACCAGCGTCAGAACCATTACCAGGTCTTCCACAATCAACCAACCGATGGCg                           <  1:2481332/69‑1 (MQ=255)
                                     gCCACTGCGGGCAGCAACACCAGCGTCAGAACCATTACCAGGTCTTCCACAATCAACCAACCGATGGCg                           <  1:2585690/69‑1 (MQ=255)
                                     gCCACTGCGGGCAGCAACACCAGCGTCAGAACCATTACCAGGTCTTCCACAATCAACCAACCGATGGCg                           <  1:256417/69‑1 (MQ=255)
                                     gCCACTGCGGGCAGCAACACCAGCGTCAGAACCATTACCAGGTCTTCCACAATCAACCAACCGATGGCg                           <  1:1794323/69‑1 (MQ=255)
                                     gCCACTGCGGGCAGCAACACCAGCGTCAGAACCATTACCAGGTCTTCCACAATCAACCAACCGATGGCg                           <  1:181793/69‑1 (MQ=255)
                                     gCCACTGCGGGCAGCAACACCAGCGTCAGAACCATTACCAGGTCTTCCACAATCAACCAACCGATGGCg                           <  1:1837688/69‑1 (MQ=255)
                                            cGGGCAGCAACACCAGCGTCAGAACCATTACCAGGTCTTCCACAATCAACCAACCGATgg                             >  1:2404113/1‑60 (MQ=255)
                                            cGGGCAGCAACACCAGCGTCAGAACCATTACCAGGTCTTCCACAATCAACCAACCGATGGCGATTTGcc                    >  1:1605932/1‑69 (MQ=255)
                                            cGGGCAGCAACACCAGCGTCAGAACCATTACCAGGTCTTCCACAATCAACCAACCGATGGCGATTTGcc                    >  1:730854/1‑69 (MQ=255)
                                            cGGGCAGCAACACCAGCGTCAGAACCATTACCAGGTCTTCCACAATCAACCAACCGATGGCGATTTGcc                    >  1:1631935/1‑69 (MQ=255)
                                            cGGGCAGCAACACCAGCGTCAGAACCATTACCAGGTCTTCCACAATCAACCAACCGATGGCGATTTGcc                    >  1:229411/1‑69 (MQ=255)
                                            cGGGCAGCAACACCAGCGTCAGAACCATTACCAGGTCTTCCACAATCAACCAACCGATGGCGATTTGcc                    >  1:225420/1‑69 (MQ=255)
                                            cGGGCAGCAACACCAGCGTCAGAACCATTACCAGGTCTTCCACAATCAACCAACCGATGGCGATTTGcc                    >  1:2967247/1‑69 (MQ=255)
                                                 agcaACACCAGCGTCAGAACCATTACCAGGTCTTCCACAATCAACCCACCGATGGCGATTTGCCCAc                 >  1:1745750/1‑67 (MQ=255)
                                                 agcaACACCAGCGTCAGAACCATTACCAGGTCTTCCACAATCAACCAACCGATGGCGATTTGCCCACg                >  1:1962801/1‑68 (MQ=255)
                                                  gcaATACCAGCGTCAGGACCATCACCATGTCTTCCACAATCAACCAACCGATGGCGATTTGCCCACGCt              <  1:2395847/69‑1 (MQ=255)
                                                  gcaATACCAGCGTCAGGACCATCACCAGGTCTTCCACAATCAACCAACCGATGGCGATTTGCCCACGCt              <  1:2222830/69‑1 (MQ=255)
                                                  gcaATACCAGCGTCAGGACCATCACCAGGTCTTCCACAATCAACCAACCGATGGCGATTTGCCCACGCt              <  1:926420/69‑1 (MQ=255)
                                                  gcaATACCAGCGTCAGGACCATCACCAGGTCTTCCACAATCAACCAACCGATGGCGATTTGCCCACGCt              <  1:865938/69‑1 (MQ=255)
                                                  gcaATACCAGCGTCAGGACCATCACCAGGTCTTCCACAATCAACCAACCGATGGCGATTTGCCCACGCt              <  1:1492920/69‑1 (MQ=255)
                                                  gcaATACCAGCGTCAGGACCATCACCAGGTCTTCCACAATCAACCAACCGATGGCGATTTGCCCACGCt              <  1:1588060/69‑1 (MQ=255)
                                                  gcaATACCAGCGTCAGGACCATCACCAGGTCTTCCACAATCAACCAACCGATGGCGATTTGCCCACGCt              <  1:1729718/69‑1 (MQ=255)
                                                  gcaATACCAGCGTCAGGACCATCACCAGGTCTTCCACAATCAACCAACCGATGGCGATTTGCCCACGCt              <  1:562040/69‑1 (MQ=255)
                                                  gcaATACCAGCGTCAGGACCATCACCAGGTCTTCCACAATCAACCAACCGATGGCGATTTGCCCACGCt              <  1:3143420/69‑1 (MQ=255)
                                                  gcaATACCAGCGTCAGGACCATCACCAGGTCTTCCACAATCAACCAACCGATGGCGATTTGCCCACGCt              <  1:2360230/69‑1 (MQ=255)
                                                  gcaATACCAGCGTCAGGACCATCACCAGGTCTTCCACAATCAACCAACCGATGGCGATTTGCCCACGCt              <  1:2361429/69‑1 (MQ=255)
                                                  gcaATACCAGCGTCAGGACCATCACCAGGTCTTCCACAATCAACCAACCGATGGCGATTTGCCCACGCt              <  1:2654218/69‑1 (MQ=255)
                                                  gcaATACCAGCGTCAGGACCATCACCAGGTCTTCCACAATCAACCAACCGATGGCGATTTGCCCACGCt              <  1:33472/69‑1 (MQ=255)
                                                  gcaATACCAGCGTCAGGACCATCACCAGGTCTTCCACAATCAACCAACCGATGGCGATTTGCCCACGCt              <  1:2909792/69‑1 (MQ=255)
                                                  gcaATACCAGCGTCAGGACCATCACCAGGTCTTCCACAATCAACCAACCGATGGCGATTTGCCCACGCt              <  1:3255524/69‑1 (MQ=255)
                                                  gcaATACCAGCGTCAGGACCATCACCAGGTCTTCCACAATCAACCAACCGATGGCGATTTGCCCACGCt              <  1:301900/69‑1 (MQ=255)
                                                   caACACCAGCGTCAGAACCATTACCAGGTCTTCCACAATCAACCAACCGATGGCGATTTGCCCACGCTg             <  1:660252/69‑1 (MQ=255)
                                                   caACACCAGCGTCAGAACCATTACCAGGTCTTCCACAATCAACCAACCGATGGCGATTTGCCCACGCTg             <  1:292556/69‑1 (MQ=255)
                                                              tCAGAACCATTACCAGGTCTTCCACAATCAACCAACCGATGGCGATTTGCCCACGCTGACTGTCaatta  <  1:2122848/69‑1 (MQ=255)
                                                              tCAGAACCATTACCAGGTCTTCCACAATCAACCAACCGATGGCGATTTGCCCACGCTGACTGTCaatta  <  1:580323/69‑1 (MQ=255)
                                                              tCAGAACCATTACCAGGTCTTCCACAATCAACCAACCGATGGCGATTTGCCCACGCTGACTGTCaatta  <  1:3310646/69‑1 (MQ=255)
                                                                  |                                                                
AGTGGCAAAGCCCACATCGCCCTGTTCCATCATTCCTGCCACTGCGGGCAGCAACACCAGCGTCAGAACCATTACCAGGTCTTCCACAATCAACCAACCGATGGCGATTTGCCCACGCTGACTGTCAATTA  >  minE/331156‑331286

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: