Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 130,768 A→G 50.0% L316L (TTA→TTG yadE → predicted polysaccharide deacetylase lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE130,7680AG50.0% ‑3.4 / 72.9 80L316L (TTA→TTGyadEpredicted polysaccharide deacetylase lipoprotein
Reads supporting (aligned to +/- strand):  ref base A (29/11);  new base G (19/21);  total (48/32)
Fisher's exact test for biased strand distribution p-value = 3.92e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTGATTTCCAGTCACATACCCATTTTTTGCATCGGGTAGATGGTTATCGCCGACCCATATTACTGAGCCGTAGTGAGCACAATATTCTGTTTGATTTTGCACGTTCACGCCGCGCTCTGGCGCAATTTAA  >  minE/130706‑130836
                                                              |                                                                    
tttGATTTCCAGTCACATACCCATTTTTTGCATCGGGTAGATGGTTATCGCCGACCCATATTACTGAGc                                                                >  1:2620451/1‑69 (MQ=255)
    aTTTCCAGTCACATACCCATTTTTTGCATCGGGTAGATGGTTATCGCCGACCCATATTACTGAGCCGt                                                             <  1:3066171/68‑1 (MQ=255)
                      aTTTTTTGCATCGGATAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTAGTGAGCACAATATTCTg                                           >  1:2452447/1‑68 (MQ=255)
                      aTTTTTTGCATCGGATAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTAGTGAGCACAATATTCTGt                                          >  1:1342739/1‑69 (MQ=255)
                      aTTTTTTGCATCGGATAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTAGTGAGCACAATATTCTGt                                          >  1:1761935/1‑69 (MQ=255)
                      aTTTTTTGCATCGGATAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTAGTGAGCACAATATTCTGt                                          >  1:1940560/1‑69 (MQ=255)
                      aTTTTTTGCATCGGATAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTAGTGAGCACAATATTCTGt                                          >  1:197627/1‑69 (MQ=255)
                      aTTTTTTGCATCGGATAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTAGTGAGCACAATATTCTGt                                          >  1:1466739/1‑69 (MQ=255)
                      aTTTTTTGCATCGGATAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTAGTGAGCACAATATTCTGt                                          >  1:309454/1‑69 (MQ=255)
                      aTTTTTTGCATCGGATAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTAGTGAGCACAATATTCTGt                                          >  1:3256514/1‑69 (MQ=255)
                      aTTTTTTGCATCGGATAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTAGTGAGCACAATATTCTGt                                          >  1:3279268/1‑69 (MQ=255)
                      aTTTTTTGCATCGGATAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTAGTGAGCACAATATTCTGt                                          >  1:1300204/1‑69 (MQ=255)
                      aTTTTTTGCATCGGATAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTAGTGAGCACAATATTCTGt                                          >  1:2388706/1‑69 (MQ=255)
                      aTTTTTTGCATCGGATAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTAGTGAGCACAATATTCTGt                                          >  1:1079912/1‑69 (MQ=255)
                      aTTTTTTGCATCGGATAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTAGTGAGCACAATATTCTGt                                          >  1:1049846/1‑69 (MQ=255)
                      aTTTTTTGCATAGGATAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTAGTGAGCACAATATTCTGt                                          >  1:3074273/1‑69 (MQ=255)
                         ttttGCATCGGGTAGATGGTTATCGCCGACCCATATTACTGa                                                                  >  1:1123344/1‑42 (MQ=255)
                         ttttGCATCGGGTAGATGGTTATCGCCGACCCATATTACTGAGCCGTAGTGAGCACAATATTCTGTTTg                                       >  1:473200/1‑69 (MQ=255)
                         ttttGCATCGGGTAGATGGTTATCGCCGACCCATATTACTGAGCCGTAGTGAGCACAATATTCTGTTTg                                       >  1:550340/1‑69 (MQ=255)
                         ttttGCATCGGGTAGATGGTTATCGCCGACCCATATTACTGAGCCGTAGTGAGCACAATATTCTGTTTg                                       >  1:2822523/1‑69 (MQ=255)
                         ttttGCATCGGGTAGATGGTTATCGCCGACCCATATTACTGAGCCGTAGTGAGCACAATATTCTGTTTg                                       >  1:1659862/1‑69 (MQ=255)
                         ttttGCATCGGATAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTa                                                            <  1:229415/48‑1 (MQ=255)
                         ttttGCATCGGATAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTa                                                            <  1:2506743/48‑1 (MQ=255)
                         ttttGCATCGGATAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTa                                                            <  1:2676678/48‑1 (MQ=255)
                         ttttGCATCGGATAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTa                                                            <  1:1049094/48‑1 (MQ=255)
                         ttttGCATCGGATAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTa                                                            <  1:2677249/48‑1 (MQ=255)
                         ttttGCATCGGATAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTa                                                            <  1:2783278/48‑1 (MQ=255)
                         ttttGCATCGGATAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTa                                                            <  1:2920247/48‑1 (MQ=255)
                         ttttGCATCGGATAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTa                                                            <  1:2930983/48‑1 (MQ=255)
                         ttttGCATCGGATAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTa                                                            <  1:3001774/48‑1 (MQ=255)
                         ttttGCATCGGATAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTa                                                            <  1:3055156/48‑1 (MQ=255)
                         ttttGCATCGGATAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTa                                                            <  1:3413015/48‑1 (MQ=255)
                         ttttGCATCGGATAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTa                                                            <  1:714824/48‑1 (MQ=255)
                         ttttGCATCGGATAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTa                                                            <  1:93002/48‑1 (MQ=255)
                         ttttGCATCGGATAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTa                                                            <  1:2244063/48‑1 (MQ=255)
                         ttttGCATCGGATAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTa                                                            <  1:1513064/48‑1 (MQ=255)
                         ttttGCATCGGATAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTa                                                            <  1:1385321/48‑1 (MQ=255)
                         ttttGCATCGGATAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTa                                                            <  1:216183/48‑1 (MQ=255)
                         ttttGCATCGGATAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTa                                                            <  1:108989/48‑1 (MQ=255)
                         ttttGCATCGGATAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTa                                                            <  1:1328291/48‑1 (MQ=255)
                         ttttGCATCGGATAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTa                                                            <  1:1739439/48‑1 (MQ=255)
                         ttttGCATCGGATAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTa                                                            <  1:1876398/48‑1 (MQ=255)
                                  ggatAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTAGTGAGCACAATATTCTGTTTGATTTTGCAc                              >  1:3402962/4‑69 (MQ=255)
                                  ggatAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTAGTGAGCACAATATTCTGTTTGATTTTGCAc                              >  1:1495543/4‑69 (MQ=255)
                                  ggatAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTAGTGAGCACAATATTCTGTTTGATTTTGCAc                              >  1:46366/4‑69 (MQ=255)
                                  ggatAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTAGTGAGCACAATATTCTGTTTGATTTTGCAc                              >  1:2067525/4‑69 (MQ=255)
                                  ggatAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTAGTGAGCACAATATTCTGTTTGATTTTGCAc                              >  1:2162444/4‑69 (MQ=255)
                                     tAGATGGTTATCGCCGACCCATATTACTGAGCCGTAGTGAGCACAATATTCTGTTTGATTTTGCACGtt                           <  1:1335104/69‑1 (MQ=255)
                                          ggTTATCGCCGACCCATATTACTGAGCCGTAGTGAGCACAATATTCTGTTTGATTTTg                                 >  1:2575882/1‑58 (MQ=255)
                                          ggTTATCGCCGACCCATATTACTGAGCCGTAGTGAGCACAATATTCTGTTTGATTTTg                                 <  1:1891845/58‑1 (MQ=255)
                                          ggTTATCGCCGACCCATATTACTGAGCCGTAGTGAGCACAATATTCTGTTTGATTTTg                                 >  1:1215901/1‑58 (MQ=255)
                                          ggTTATCGCCGACCCATATTACTGAGCCGTAGTGAGCACAATATTCTGTTTGATTTTg                                 >  1:3363384/1‑58 (MQ=255)
                                          ggTTATCGCCGACCCATATTACTGAGCCGTAGTGAGCACAATATTCTGTTTGATTTTg                                 >  1:2426632/1‑58 (MQ=255)
                                          ggTTATCGCCGACCCATATTACTGAGCCGTAGTGAGCACAATATTCTGTTTGATTTTg                                 >  1:2781839/1‑58 (MQ=255)
                                          ggTTATCGCCGACCCATATTACTGAGCCGTAGTGAGCACAATATTCTGTTTGATTTTg                                 >  1:135888/1‑58 (MQ=255)
                                          ggTTATCGCCGACCCATATTACTGAGCCGTAGTGAGCACAATATTCTGTTTGATTTTg                                 >  1:848971/1‑58 (MQ=255)
                                          ggTTATCGCCGACCCATATTACTGAGCCGTAGTGAGCACAATATTCTGTTTGATTTTg                                 <  1:2861675/58‑1 (MQ=255)
                                          ggTTATCGCCGACCCATATTACTGAGCCGTAGTGAGCACAATATTCTGTTTGATTTTg                                 >  1:2276658/1‑58 (MQ=255)
                                          ggTTATCGCCGACCCATATTACTGAGCCGTAGTGAGCACAATATTCTGTTTGATTTTg                                 >  1:149634/1‑58 (MQ=255)
                                              aTCGCCGACCCATATTACTGAGCCGTAGTGAGCACAATATTCTGTTTGATTTTGCACGTTCACGCcgcg                  <  1:2245514/69‑1 (MQ=255)
                                              aTCGCCGACCCATATTACTGAGCCGTAGTGAGCACAATATTCTGTTTGATTTTGCACGTTCACGCcgcg                  <  1:273883/69‑1 (MQ=255)
                                              aTCGCCGACCCATATTACTGAGCCGTAGTGAGCACAATATTCTGTTTGATTTTGCACGTTCACGCcgcg                  <  1:2696447/69‑1 (MQ=255)
                                                            ttACTGAGCCGTAGTGATCACAATATTCTGTTTGATTTTGCACGTTCACGCCgcgc                 <  1:1516020/56‑1 (MQ=255)
                                                            ttACTGAGCCGTAGTGAGCACAATATTCTGTTTGATTTTGCACGTTCACGCCgcgc                 <  1:3352177/56‑1 (MQ=255)
                                                            ttACTGAGCCGTAGTGAGCACAATATTCTGTTTGATTTTGCACGTTCACGCCgcgc                 <  1:2632034/56‑1 (MQ=255)
                                                            ttACTGAGCCGTAGTGAGCACAATATTCTGTTTGATTTTGCACGTTCACGCCgcgc                 <  1:462808/56‑1 (MQ=255)
                                                            ttACTGAGCCGTAGTGAGCACAATATTCTGTTTGATTTTGCACGTTCACGCCGCGCTCTGGCGCAAttt    >  1:1321249/1‑69 (MQ=255)
                                                            ttACTGAGCCGTAGTGAGCACAATATTCTGTTTGATTTTGCACGTTCACGCCGCGCTCTGGCGCAAttt    >  1:1026780/1‑69 (MQ=255)
                                                            ttACTGAGCCGTAGTGAGCACAATATTCTGTTTGATTTTGCACGTTCACGCCGCGCTCTGGCGCAAttt    >  1:8980/1‑69 (MQ=255)
                                                            ttACTGAGCCGTAGTGAGCACAATATTCTGTTTGATTTTGCACGTTCACGCCGCGCTCTGGCGCAAttt    >  1:2298937/1‑69 (MQ=255)
                                                            ttACTGAGCCGTAGTGAGCACAATATTCTGTTTGATTTTGCACGTTCACGCCGCGCTCTGGCGCAAttt    >  1:750912/1‑69 (MQ=255)
                                                            ttACTGAGCCGTAGTGAGCACAATATTCTGTTTGATTTTGCACGTTCACGCCGCGCTCTGGCGCAAttt    >  1:734313/1‑69 (MQ=255)
                                                            ttACTGAGCCGTAGTGAGCACAATATTCTGTTTGATTTTGCACGTTCACGCCGCGCTCTGGCGCAAttt    >  1:2390184/1‑69 (MQ=255)
                                                            ttACTGAGCCGTAGTGAGCACAATATTCTGTTTGATTTTGCACGTTCACGCCGCGCTCTGGCGCAAttt    >  1:2469907/1‑69 (MQ=255)
                                                            ttACTGAGCCGTAGTGAGCACAATATTCTGTTTGATTTTGCACGTTCACGCCGCGCTCTGGCGCAAttt    >  1:124193/1‑69 (MQ=255)
                                                            ttACTGAGCCGTAGTGAGCACAATATTCTGTTTGATTTTGCACGTTCACGCCGCGCTCTGGCGCAAttt    >  1:3349913/1‑69 (MQ=255)
                                                            ttACTGAGCCGTAGTGAGCACAATATTCTGTTTGATTTTGCACGTTCACGCCGCGCTCTGGCGCAAttt    >  1:2709307/1‑69 (MQ=255)
                                                            ttACTGAGCCGTAGTGAGCACAATATTCTGTTTGATTTTGCACGTTCACGCCGCGCTCTGGCGCAAttt    >  1:1662746/1‑69 (MQ=255)
                                                            ttACTGAGCCGTAGTGAGCACAATATTCTGTTTGATTTTGCACGTTCACGCCGCGCTCTGGCGCAAttt    >  1:1391328/1‑69 (MQ=255)
                                                            ttACTGAGCCGTAGTGAGCACAATATTCTGTTTGATTTTGCACGTTCACGCCGCGCTCTGGCGCAAttt    >  1:1419986/1‑69 (MQ=255)
                                                              aCTGAGCCGTAGTGAGCACAATATTCTGTTTGATTTTGCACGTTCACGCCGCGCTCTGGCGCAATTTaa  >  1:2908937/1‑69 (MQ=255)
                                                              aCTGAGCCGTAGTGAGCACAATATTCTGTTTGATTTTGCACGTTCACGCCGCGCTCTGGCGCAATTTa   >  1:2052744/1‑68 (MQ=255)
                                                              |                                                                    
TTTGATTTCCAGTCACATACCCATTTTTTGCATCGGGTAGATGGTTATCGCCGACCCATATTACTGAGCCGTAGTGAGCACAATATTCTGTTTGATTTTGCACGTTCACGCCGCGCTCTGGCGCAATTTAA  >  minE/130706‑130836

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: