Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 130,756 C→G 53.3% R312R (CGC→CGG yadE → predicted polysaccharide deacetylase lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE130,7560CG53.3% 11.5 / 86.9 75R312R (CGC→CGGyadEpredicted polysaccharide deacetylase lipoprotein
Reads supporting (aligned to +/- strand):  ref base C (22/13);  new base G (19/21);  total (41/34)
Fisher's exact test for biased strand distribution p-value = 2.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGAAATTCGCGATGTATTTGATTTCCAGTCACATACCCATTTTTTGCATCGGGTAGATGGTTATCGCCGACCCATATTACTGAGCCGTAGTGAGCACAATATTCTGTTTGATTTTGCACGTTCACGCCGCG  >  minE/130690‑130820
                                                                  |                                                                
cGAAATTCGCGATGTATTTGATTTCCAGTCACATACCCATTTTTTGCATCGGGTAGATGGTTATCGCCg                                                                >  1:995058/1‑69 (MQ=255)
cGAAATTCGCGATGTATTTGATTTCCAGTCACATACCCATTTTTTGCATCGGGTAGATGGTTATCGCCg                                                                >  1:2148885/1‑69 (MQ=255)
cGAAATTCGCGATGTATTTGATTTCCAGTCACATACCCATTTTTTGCATCGGGTAGATGGTTATCGCCg                                                                >  1:2405231/1‑69 (MQ=255)
cGAAATTCGCGATGTATTTGATTTCCAGTCACATACCCATTTTTTGCATCGGGTAGATGGTTATCGCCg                                                                >  1:2869332/1‑69 (MQ=255)
cGAAATTCGCGATGTATTTGATTTCCAGTCACATACCCATTTTTTGCATCGGGTAGATGGTTATCGCCg                                                                >  1:1257080/1‑69 (MQ=255)
cGAAATTCGCGATGTATTTGATTTCCAGTCACATACCCATTTTTTGCATCGGGTAGATGGTTATCGCCg                                                                >  1:635098/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:1761935/1‑69 (MQ=255)
                                      aTTTTTTGCATCGGATAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTAGTGAGCACAATATTCTGt                          >  1:309454/1‑69 (MQ=255)
                                      aTTTTTTGCATCGGATAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTAGTGAGCACAATATTCTGt                          >  1:1940560/1‑69 (MQ=255)
                                      aTTTTTTGCATCGGATAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTAGTGAGCACAATATTCTGt                          >  1:197627/1‑69 (MQ=255)
                                      aTTTTTTGCATCGGATAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTAGTGAGCACAATATTCTGt                          >  1:2388706/1‑69 (MQ=255)
                                      aTTTTTTGCATCGGATAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTAGTGAGCACAATATTCTGt                          >  1:3279268/1‑69 (MQ=255)
                                      aTTTTTTGCATCGGATAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTAGTGAGCACAATATTCTGt                          >  1:3256514/1‑69 (MQ=255)
                                      aTTTTTTGCATCGGATAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTAGTGAGCACAATATTCTGt                          >  1:1049846/1‑69 (MQ=255)
                                      aTTTTTTGCATCGGATAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTAGTGAGCACAATATTCTGt                          >  1:1342739/1‑69 (MQ=255)
                                      aTTTTTTGCATCGGATAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTAGTGAGCACAATATTCTGt                          >  1:1466739/1‑69 (MQ=255)
                                      aTTTTTTGCATCGGATAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTAGTGAGCACAATATTCTGt                          >  1:1300204/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:550340/1‑69 (MQ=255)
                                         ttttGCATCGGGTAGATGGTTATCGCCGACCCATATTACTGAGCCGTAGTGAGCACAATATTCTGTTTg                       >  1:473200/1‑69 (MQ=255)
                                         ttttGCATCGGGTAGATGGTTATCGCCGACCCATATTACTGAGCCGTAGTGAGCACAATATTCTGTTTg                       >  1:1659862/1‑69 (MQ=255)
                                         ttttGCATCGGGTAGATGGTTATCGCCGACCCATATTACTGAGCCGTAGTGAGCACAATATTCTGTTTg                       >  1:2822523/1‑69 (MQ=255)
                                         ttttGCATCGGATAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTa                                            <  1:3413015/48‑1 (MQ=255)
                                         ttttGCATCGGATAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTa                                            <  1:108989/48‑1 (MQ=255)
                                         ttttGCATCGGATAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTa                                            <  1:2783278/48‑1 (MQ=255)
                                         ttttGCATCGGATAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTa                                            <  1:714824/48‑1 (MQ=255)
                                         ttttGCATCGGATAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTa                                            <  1:3055156/48‑1 (MQ=255)
                                         ttttGCATCGGATAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTa                                            <  1:1328291/48‑1 (MQ=255)
                                         ttttGCATCGGATAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTa                                            <  1:3001774/48‑1 (MQ=255)
                                         ttttGCATCGGATAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTa                                            <  1:2930983/48‑1 (MQ=255)
                                         ttttGCATCGGATAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTa                                            <  1:2920247/48‑1 (MQ=255)
                                         ttttGCATCGGATAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTa                                            <  1:93002/48‑1 (MQ=255)
                                         ttttGCATCGGATAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTa                                            <  1:1739439/48‑1 (MQ=255)
                                         ttttGCATCGGATAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTa                                            <  1:1876398/48‑1 (MQ=255)
                                         ttttGCATCGGATAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTa                                            <  1:2244063/48‑1 (MQ=255)
                                         ttttGCATCGGATAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTa                                            <  1:2677249/48‑1 (MQ=255)
                                         ttttGCATCGGATAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTa                                            <  1:229415/48‑1 (MQ=255)
                                         ttttGCATCGGATAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTa                                            <  1:1049094/48‑1 (MQ=255)
                                         ttttGCATCGGATAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTa                                            <  1:216183/48‑1 (MQ=255)
                                         ttttGCATCGGATAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTa                                            <  1:1385321/48‑1 (MQ=255)
                                         ttttGCATCGGATAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTa                                            <  1:2506743/48‑1 (MQ=255)
                                         ttttGCATCGGATAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTa                                            <  1:2676678/48‑1 (MQ=255)
                                         ttttGCATCGGATAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTa                                            <  1:1513064/48‑1 (MQ=255)
                                                  ggatAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTAGTGAGCACAATATTCTGTTTGATTTTGCAc              >  1:46366/4‑69 (MQ=255)
                                                  ggatAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTAGTGAGCACAATATTCTGTTTGATTTTGCAc              >  1:1495543/4‑69 (MQ=255)
                                                  ggatAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTAGTGAGCACAATATTCTGTTTGATTTTGCAc              >  1:2067525/4‑69 (MQ=255)
                                                  ggatAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTAGTGAGCACAATATTCTGTTTGATTTTGCAc              >  1:3402962/4‑69 (MQ=255)
                                                  ggatAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTAGTGAGCACAATATTCTGTTTGATTTTGCAc              >  1:2162444/4‑69 (MQ=255)
                                                     tAGATGGTTATCGCCGACCCATATTACTGAGCCGTAGTGAGCACAATATTCTGTTTGATTTTGCACGtt           <  1:1335104/69‑1 (MQ=255)
                                                          ggTTATCGCCGACCCATATTACTGAGCCGTAGTGAGCACAATATTCTGTTTGATTTTg                 >  1:149634/1‑58 (MQ=255)
                                                          ggTTATCGCCGACCCATATTACTGAGCCGTAGTGAGCACAATATTCTGTTTGATTTTg                 >  1:3363384/1‑58 (MQ=255)
                                                          ggTTATCGCCGACCCATATTACTGAGCCGTAGTGAGCACAATATTCTGTTTGATTTTg                 >  1:848971/1‑58 (MQ=255)
                                                          ggTTATCGCCGACCCATATTACTGAGCCGTAGTGAGCACAATATTCTGTTTGATTTTg                 >  1:1215901/1‑58 (MQ=255)
                                                          ggTTATCGCCGACCCATATTACTGAGCCGTAGTGAGCACAATATTCTGTTTGATTTTg                 >  1:2276658/1‑58 (MQ=255)
                                                          ggTTATCGCCGACCCATATTACTGAGCCGTAGTGAGCACAATATTCTGTTTGATTTTg                 >  1:2426632/1‑58 (MQ=255)
                                                          ggTTATCGCCGACCCATATTACTGAGCCGTAGTGAGCACAATATTCTGTTTGATTTTg                 >  1:2575882/1‑58 (MQ=255)
                                                          ggTTATCGCCGACCCATATTACTGAGCCGTAGTGAGCACAATATTCTGTTTGATTTTg                 <  1:2861675/58‑1 (MQ=255)
                                                          ggTTATCGCCGACCCATATTACTGAGCCGTAGTGAGCACAATATTCTGTTTGATTTTg                 >  1:135888/1‑58 (MQ=255)
                                                          ggTTATCGCCGACCCATATTACTGAGCCGTAGTGAGCACAATATTCTGTTTGATTTTg                 >  1:2781839/1‑58 (MQ=255)
                                                          ggTTATCGCCGACCCATATTACTGAGCCGTAGTGAGCACAATATTCTGTTTGATTTTg                 <  1:1891845/58‑1 (MQ=255)
                                                              aTCGCCGACCCATATTACTGAGCCGTAGTGAGCACAATATTCTGTTTGATTTTGCACGTTCACGCcgcg  <  1:2245514/69‑1 (MQ=255)
                                                              aTCGCCGACCCATATTACTGAGCCGTAGTGAGCACAATATTCTGTTTGATTTTGCACGTTCACGCcgcg  <  1:273883/69‑1 (MQ=255)
                                                              aTCGCCGACCCATATTACTGAGCCGTAGTGAGCACAATATTCTGTTTGATTTTGCACGTTCACGCcgcg  <  1:2696447/69‑1 (MQ=255)
                                                                  |                                                                
CGAAATTCGCGATGTATTTGATTTCCAGTCACATACCCATTTTTTGCATCGGGTAGATGGTTATCGCCGACCCATATTACTGAGCCGTAGTGAGCACAATATTCTGTTTGATTTTGCACGTTCACGCCGCG  >  minE/130690‑130820

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: