Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 130,742 G→A 52.9% V308I (GTA→ATA)  yadE → predicted polysaccharide deacetylase lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE130,7420GA52.9% 3.5 / 72.0 84V308I (GTA→ATA) yadEpredicted polysaccharide deacetylase lipoprotein
Reads supporting (aligned to +/- strand):  ref base G (22/17);  new base A (14/30);  total (36/48)
Fisher's exact test for biased strand distribution p-value = 2.85e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGTTTCTGAGCTTAACGAAATTCGCGATGTATTTGATTTCCAGTCACATACCCATTTTTTGCATCGGGTAGATGGTTATCGCCGACCCATATTACTGAGCCGTAGTGAGCACAATATTCTGTTTG  >  minE/130675‑130799
                                                                   |                                                         
cGTTTCTGAGCTTAACGAAATTCGCGATGTATTTGATTTCCAGTCACATACCCATTTTTTGCATCggg                                                           <  1:1569650/68‑1 (MQ=255)
cGTTTCTGAGCTTAACGAAATTCGCGATGTATTTGATTTCCAGTCACATACCCATTTTTTGCATCggg                                                           <  1:595455/68‑1 (MQ=255)
cGTTTCTGAGCTTAACGAAATTCGCGATGTATTTGATTTCCAGTCACATACCCATTTTTTGCATCggg                                                           <  1:2679868/68‑1 (MQ=255)
cGTTTCTGAGCTTAACGAAATTCGCGATGTATTTGATTTCCAGTCACATACCCATTTTTTGCATCggg                                                           <  1:3241396/68‑1 (MQ=255)
 gTTTCTGAGCTTAACGAAATTCGCGATGTATTTGATTTCCAGTCACATACCCATTTTTTGCTTCGGGTa                                                         <  1:2985925/69‑1 (MQ=255)
 gTTTCTGAGCTTAACGAAATTCGCGATGTATTTGATTTCCAGTCACATACCCATTTTTTGCATCGGGTa                                                         <  1:2350739/69‑1 (MQ=255)
 gTTTCTGAGCTTAACGAAATTCGCGATGTATTTGATTTCCAGTCACATACCCATTTTTTGCATCGGGTa                                                         <  1:2098699/69‑1 (MQ=255)
 gTTTCTGAGCTTAACGAAATTCGCGATGTATTTGATTTCCAGTCACATACCCATTTTTTGCATCGGGTa                                                         <  1:600879/69‑1 (MQ=255)
 gTTTCTGAGCTTAACGAAATTCGCGATGTATTTGATTTCCAGTCACATACCCATTTTTTGCATCGGGTa                                                         <  1:2015829/69‑1 (MQ=255)
 gTTTCTGAGCTTAACGAAATTCGCGATGTATTTGATTTCCAGTCACATACCCATTTTTTGCATCGGGTa                                                         <  1:580024/69‑1 (MQ=255)
 gTTTCTGAGCTTAACGAAATTCGCGATGTATTTGATTTCCAGTCACATACCCATTTTTTGCATCGGGTa                                                         <  1:1645799/69‑1 (MQ=255)
 gTTTCTGAGCTTAACGAAATTCGCGATGTATTTGATTTCCAGTCACATACCCATTTTTTGCATCGGGTa                                                         <  1:1394758/69‑1 (MQ=255)
 gTTTCTGAGCTTAACGAAATTCGCGATGTATTTGATTTCCAGTCACATACCCATTTTTTGCATCGGGTa                                                         <  1:3002125/69‑1 (MQ=255)
  tttCTGATCTTAACGAAATTCGCGATGTATTTGATTTCCAGTCACATACCCATTTTTTGCATCGGTTa                                                         <  1:2599708/68‑1 (MQ=255)
    tCTGAGCTTAACGAAATTCGCGATGTATTTGATTTCCAGTCACATACCCATTTTTTGCATCGGGTAGAt                                                      >  1:1591929/1‑69 (MQ=255)
    tCTGAGCTTAACGAAATTCGCGATGTATTTGATTTCCAGTCACATACCCATTTTTTGCATCGGGTAGAt                                                      >  1:1630050/1‑69 (MQ=255)
    tCTGAGCTTAACGAAATTCGCGATGTATTTGATTTCCAGTCACATACCCATTTTTTGCATCGGGTAGAt                                                      >  1:3151831/1‑69 (MQ=255)
    tCTGAGCTTAACGAAATTCGCGATGTATTTGATTTCCAGTCACATACCCATTTTTTGCATCGGGTAGAt                                                      >  1:421563/1‑69 (MQ=255)
    tCTGAGCTTAACGAAATTCGCGATGTATTTGATTTCCAGTCACATACCCATTTTTTGCATCGGGTAGAt                                                      >  1:483996/1‑69 (MQ=255)
    tCTGAGCTTAACGAAATTCGCGATGTATTTGATTTCCAGTCACATACCCATTTTTTGCATCGGGTAGAt                                                      >  1:766881/1‑69 (MQ=255)
          cTTAACGAAATTCGCGATGTATTTGATTTCCAGTCACATACCCATTTTTTGCATCGGGTAGATGGTTAt                                                >  1:786711/1‑69 (MQ=255)
          cTTAACGAAATTCGCGATGTATTTGATTTCCAGTCACATACCCATTTTTTGCATCGGGTAGATGGTTAt                                                >  1:1562810/1‑69 (MQ=255)
          cTTAACGAAATTCGCGATGTATTTGATTTCCAGTCACATACCCATTTTTTGCATCGGGTAGATGGTTAt                                                >  1:1520027/1‑69 (MQ=255)
           ttAACGAAATTCGCGATGTATTTGATTTCCAGTCACATACTCATTTTTTGCATCGGATAGATGGTTATc                                               <  1:3016733/69‑1 (MQ=255)
           ttAACGAAATTCGCGATGTATTTGATTTCCAGTCACATACTCATTTTTTGCATCGGATAGATGGTTATc                                               <  1:3122126/69‑1 (MQ=255)
           ttAACGAAATTCGCGATGTATTTGATTTCCAGTCACATACTCATTTTTTGCATCGGATAGATGGTTATc                                               <  1:3425511/69‑1 (MQ=255)
           ttAACGAAATTCGCGATGTATTTGATTTCCAGTCACATACTCATTTTTTGCATCGGATAGATGGTTATc                                               <  1:354193/69‑1 (MQ=255)
           ttAACGAAATTCGCGATGTATTTGATTTCCAGTCACATACTCATTTTTTGCATCGGATAGATGGTTATc                                               <  1:2473046/69‑1 (MQ=255)
           ttAACGAAATTCGCGATGTATTTGATTTCCAGTCACATACTCATTTTTTGCATCGGATAGATGGTTATc                                               <  1:2427474/69‑1 (MQ=255)
           ttAACGAAATTCGCGATGTATTTGATTTCCAGTCACATACTCATTTTTTGCATCGGATAGATGGTTATc                                               <  1:1172731/69‑1 (MQ=255)
           ttAACGAAATTCGCGATGTATTTGATTTCCAGTCACATACTCATTTTTTGCATCGGATAGATGGTTATc                                               <  1:2086689/69‑1 (MQ=255)
           ttAACGAAATTCGCGATGTATTTGATTTCCAGTCACATACTCATTTTTTGCATCGGATAGATGGTTATc                                               <  1:2294374/69‑1 (MQ=255)
               cGAAATTCGCGATGTATTTGATTTCCAGTCACATACCCATTTTTTGCATCGGGTa                                                         <  1:532787/55‑1 (MQ=255)
               cGAAATTCGCGATGTATTTGATTTCCAGTCACATACCCATTTTTTGCATCGGGTAGATGGTTATCGCCg                                           >  1:2869332/1‑69 (MQ=255)
               cGAAATTCGCGATGTATTTGATTTCCAGTCACATACCCATTTTTTGCATCGGGTAGATGGTTATCGCCg                                           >  1:635098/1‑69 (MQ=255)
               cGAAATTCGCGATGTATTTGATTTCCAGTCACATACCCATTTTTTGCATCGGGTAGATGGTTATCGCCg                                           >  1:2405231/1‑69 (MQ=255)
               cGAAATTCGCGATGTATTTGATTTCCAGTCACATACCCATTTTTTGCATCGGGTAGATGGTTATCGCCg                                           >  1:995058/1‑69 (MQ=255)
               cGAAATTCGCGATGTATTTGATTTCCAGTCACATACCCATTTTTTGCATCGGGTAGATGGTTATCGCCg                                           >  1:2148885/1‑69 (MQ=255)
               cGAAATTCGCGATGTATTTGATTTCCAGTCACATACCCATTTTTTGCATCGGGTAGATGGTTATCGCCg                                           >  1:1257080/1‑69 (MQ=255)
               cGAAATTCGCGATGTATTTGATTTCAGTCACATACCCATTTTTTGCATCGGGTAGATGGTTATCGCCGa                                          >  1:3050963/1‑69 (MQ=255)
                   aTTCGCGATGTATTTGATTTCCAGTCACATACCCATTTTTTGCATGGGGTAGATGGTTATCGCCGAccc                                       <  1:2979932/69‑1 (MQ=255)
                    ttCGCGATGTATTTGATTTCCAGTCACATACCCATTTTTTGCATCGGGTAGATGGTTATCGCCGACCCa                                      <  1:1259238/69‑1 (MQ=255)
                    ttCGCGATGTATTTGATTTCCAGTCACATACCCATTTTTTGCATCGGGTAGATGGTTATCGCCGACCCa                                      <  1:1348653/69‑1 (MQ=255)
                    ttCGCGATGTATTTGATTTCCAGTCACATACCCATTTTTTGCATCGGGTAGATGGTTATCGCCGACCCa                                      <  1:1433423/69‑1 (MQ=255)
                    ttCGCGATGTATTTGATTTCCAGTCACATACCCATTTTTTGCATCGGGTAGATGGTTATCGCCGACCCa                                      <  1:1573818/69‑1 (MQ=255)
                         gATGTATTTGATTTCCAGTCACATACCCATTTTTTGCATCGGGTAGATGGTTATCGCCGACCCa                                      <  1:1915540/64‑1 (MQ=255)
                               tttGATTTCCAGTCACATACCCATTTTTTGCATCGGGTAGATGGTTATCGCCGACCCATATTACTGAGc                           >  1:2620451/1‑69 (MQ=255)
                                   aTTTCCAGTCACATACCCATTTTTTGCATCGGGTAGATGGTTATCGCCGACCCATATTACTGAGCCGt                        <  1:3066171/68‑1 (MQ=255)
                                                     aTTTTTTGCATCGGATAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTAGTGAGCACAATATTCTg      >  1:2452447/1‑68 (MQ=255)
                                                     aTTTTTTGCATCGGATAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTAGTGAGCACAATATTCTGt     >  1:1466739/1‑69 (MQ=255)
                                                     aTTTTTTGCATCGGATAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTAGTGAGCACAATATTCTGt     >  1:309454/1‑69 (MQ=255)
                                                     aTTTTTTGCATCGGATAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTAGTGAGCACAATATTCTGt     >  1:2388706/1‑69 (MQ=255)
                                                     aTTTTTTGCATCGGATAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTAGTGAGCACAATATTCTGt     >  1:197627/1‑69 (MQ=255)
                                                     aTTTTTTGCATCGGATAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTAGTGAGCACAATATTCTGt     >  1:1049846/1‑69 (MQ=255)
                                                     aTTTTTTGCATCGGATAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTAGTGAGCACAATATTCTGt     >  1:1300204/1‑69 (MQ=255)
                                                     aTTTTTTGCATCGGATAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTAGTGAGCACAATATTCTGt     >  1:1940560/1‑69 (MQ=255)
                                                     aTTTTTTGCATCGGATAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTAGTGAGCACAATATTCTGt     >  1:1761935/1‑69 (MQ=255)
                                                     aTTTTTTGCATCGGATAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTAGTGAGCACAATATTCTGt     >  1:3256514/1‑69 (MQ=255)
                                                     aTTTTTTGCATCGGATAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTAGTGAGCACAATATTCTGt     >  1:3279268/1‑69 (MQ=255)
                                                     aTTTTTTGCATCGGATAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTAGTGAGCACAATATTCTGt     >  1:1342739/1‑69 (MQ=255)
                                                     aTTTTTTGCATCGGATAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTAGTGAGCACAATATTCTGt     >  1:1079912/1‑69 (MQ=255)
                                                     aTTTTTTGCATAGGATAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTAGTGAGCACAATATTCTGt     >  1:3074273/1‑69 (MQ=255)
                                                        ttttGCATCGGGTAGATGGTTATCGCCGACCCATATTACTGa                             >  1:1123344/1‑42 (MQ=255)
                                                        ttttGCATCGGGTAGATGGTTATCGCCGACCCATATTACTGAGCCGTAGTGAGCACAATATTCTGTTTg  >  1:1659862/1‑69 (MQ=255)
                                                        ttttGCATCGGGTAGATGGTTATCGCCGACCCATATTACTGAGCCGTAGTGAGCACAATATTCTGTTTg  >  1:473200/1‑69 (MQ=255)
                                                        ttttGCATCGGGTAGATGGTTATCGCCGACCCATATTACTGAGCCGTAGTGAGCACAATATTCTGTTTg  >  1:550340/1‑69 (MQ=255)
                                                        ttttGCATCGGGTAGATGGTTATCGCCGACCCATATTACTGAGCCGTAGTGAGCACAATATTCTGTTTg  >  1:2822523/1‑69 (MQ=255)
                                                        ttttGCATCGGATAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTa                       <  1:1328291/48‑1 (MQ=255)
                                                        ttttGCATCGGATAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTa                       <  1:3055156/48‑1 (MQ=255)
                                                        ttttGCATCGGATAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTa                       <  1:93002/48‑1 (MQ=255)
                                                        ttttGCATCGGATAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTa                       <  1:108989/48‑1 (MQ=255)
                                                        ttttGCATCGGATAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTa                       <  1:714824/48‑1 (MQ=255)
                                                        ttttGCATCGGATAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTa                       <  1:1049094/48‑1 (MQ=255)
                                                        ttttGCATCGGATAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTa                       <  1:1385321/48‑1 (MQ=255)
                                                        ttttGCATCGGATAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTa                       <  1:3413015/48‑1 (MQ=255)
                                                        ttttGCATCGGATAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTa                       <  1:1513064/48‑1 (MQ=255)
                                                        ttttGCATCGGATAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTa                       <  1:229415/48‑1 (MQ=255)
                                                        ttttGCATCGGATAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTa                       <  1:3001774/48‑1 (MQ=255)
                                                        ttttGCATCGGATAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTa                       <  1:1739439/48‑1 (MQ=255)
                                                        ttttGCATCGGATAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTa                       <  1:2930983/48‑1 (MQ=255)
                                                        ttttGCATCGGATAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTa                       <  1:2920247/48‑1 (MQ=255)
                                                        ttttGCATCGGATAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTa                       <  1:2783278/48‑1 (MQ=255)
                                                        ttttGCATCGGATAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTa                       <  1:1876398/48‑1 (MQ=255)
                                                        ttttGCATCGGATAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTa                       <  1:2677249/48‑1 (MQ=255)
                                                        ttttGCATCGGATAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTa                       <  1:2676678/48‑1 (MQ=255)
                                                        ttttGCATCGGATAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTa                       <  1:2506743/48‑1 (MQ=255)
                                                        ttttGCATCGGATAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTa                       <  1:216183/48‑1 (MQ=255)
                                                        ttttGCATCGGATAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTa                       <  1:2244063/48‑1 (MQ=255)
                                                                   |                                                         
CGTTTCTGAGCTTAACGAAATTCGCGATGTATTTGATTTCCAGTCACATACCCATTTTTTGCATCGGGTAGATGGTTATCGCCGACCCATATTACTGAGCCGTAGTGAGCACAATATTCTGTTTG  >  minE/130675‑130799

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: