Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 130,762 C→T 59.7% P314P (CCC→CCT yadE → predicted polysaccharide deacetylase lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE130,7620CT59.7% 23.0 / 47.6 67P314P (CCC→CCTyadEpredicted polysaccharide deacetylase lipoprotein
Reads supporting (aligned to +/- strand):  ref base C (15/12);  new base T (19/21);  total (34/33)
Fisher's exact test for biased strand distribution p-value = 6.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATTCGCGATGTATTTGATTTCCAGTCACATACCCATTTTTTGCATCGGGTAGATGGTTATCGCCGACCCATATTACTGAGCCGTAGTGAGCACAATATTCTGTTTGATTTTGCACGTTCACGCCGCG  >  minE/130694‑130820
                                                                    |                                                          
aTTCGCGATGTATTTGATTTCCAGTCACATACCCATTTTTTGCATGGGGTAGATGGTTATCGCCGAccc                                                            <  1:2979932/69‑1 (MQ=255)
 ttCGCGATGTATTTGATTTCCAGTCACATACCCATTTTTTGCATCGGGTAGATGGTTATCGCCGACCCa                                                           <  1:1573818/69‑1 (MQ=255)
 ttCGCGATGTATTTGATTTCCAGTCACATACCCATTTTTTGCATCGGGTAGATGGTTATCGCCGACCCa                                                           <  1:1259238/69‑1 (MQ=255)
 ttCGCGATGTATTTGATTTCCAGTCACATACCCATTTTTTGCATCGGGTAGATGGTTATCGCCGACCCa                                                           <  1:1433423/69‑1 (MQ=255)
 ttCGCGATGTATTTGATTTCCAGTCACATACCCATTTTTTGCATCGGGTAGATGGTTATCGCCGACCCa                                                           <  1:1348653/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:1342739/1‑69 (MQ=255)
                                  aTTTTTTGCATCGGATAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTAGTGAGCACAATATTCTGt                          >  1:1079912/1‑69 (MQ=255)
                                  aTTTTTTGCATCGGATAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTAGTGAGCACAATATTCTGt                          >  1:197627/1‑69 (MQ=255)
                                  aTTTTTTGCATCGGATAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTAGTGAGCACAATATTCTGt                          >  1:3279268/1‑69 (MQ=255)
                                  aTTTTTTGCATCGGATAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTAGTGAGCACAATATTCTGt                          >  1:1300204/1‑69 (MQ=255)
                                  aTTTTTTGCATCGGATAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTAGTGAGCACAATATTCTGt                          >  1:1466739/1‑69 (MQ=255)
                                  aTTTTTTGCATCGGATAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTAGTGAGCACAATATTCTGt                          >  1:3256514/1‑69 (MQ=255)
                                  aTTTTTTGCATCGGATAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTAGTGAGCACAATATTCTGt                          >  1:309454/1‑69 (MQ=255)
                                  aTTTTTTGCATCGGATAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTAGTGAGCACAATATTCTGt                          >  1:1049846/1‑69 (MQ=255)
                                  aTTTTTTGCATCGGATAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTAGTGAGCACAATATTCTGt                          >  1:1940560/1‑69 (MQ=255)
                                  aTTTTTTGCATCGGATAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTAGTGAGCACAATATTCTGt                          >  1:1761935/1‑69 (MQ=255)
                                  aTTTTTTGCATCGGATAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTAGTGAGCACAATATTCTGt                          >  1:2388706/1‑69 (MQ=255)
                                  aTTTTTTGCATAGGATAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTAGTGAGCACAATATTCTGt                          >  1:3074273/1‑69 (MQ=255)
                                     ttttGCATCGGGTAGATGGTTATCGCCGACCCATATTACTGa                                                  >  1:1123344/1‑42 (MQ=255)
                                     ttttGCATCGGGTAGATGGTTATCGCCGACCCATATTACTGAGCCGTAGTGAGCACAATATTCTGTTTg                       >  1:2822523/1‑69 (MQ=255)
                                     ttttGCATCGGGTAGATGGTTATCGCCGACCCATATTACTGAGCCGTAGTGAGCACAATATTCTGTTTg                       >  1:1659862/1‑69 (MQ=255)
                                     ttttGCATCGGGTAGATGGTTATCGCCGACCCATATTACTGAGCCGTAGTGAGCACAATATTCTGTTTg                       >  1:473200/1‑69 (MQ=255)
                                     ttttGCATCGGGTAGATGGTTATCGCCGACCCATATTACTGAGCCGTAGTGAGCACAATATTCTGTTTg                       >  1:550340/1‑69 (MQ=255)
                                     ttttGCATCGGATAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTa                                            <  1:2676678/48‑1 (MQ=255)
                                     ttttGCATCGGATAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTa                                            <  1:2677249/48‑1 (MQ=255)
                                     ttttGCATCGGATAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTa                                            <  1:2783278/48‑1 (MQ=255)
                                     ttttGCATCGGATAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTa                                            <  1:93002/48‑1 (MQ=255)
                                     ttttGCATCGGATAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTa                                            <  1:2920247/48‑1 (MQ=255)
                                     ttttGCATCGGATAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTa                                            <  1:2930983/48‑1 (MQ=255)
                                     ttttGCATCGGATAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTa                                            <  1:1049094/48‑1 (MQ=255)
                                     ttttGCATCGGATAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTa                                            <  1:3001774/48‑1 (MQ=255)
                                     ttttGCATCGGATAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTa                                            <  1:3055156/48‑1 (MQ=255)
                                     ttttGCATCGGATAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTa                                            <  1:714824/48‑1 (MQ=255)
                                     ttttGCATCGGATAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTa                                            <  1:3413015/48‑1 (MQ=255)
                                     ttttGCATCGGATAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTa                                            <  1:1876398/48‑1 (MQ=255)
                                     ttttGCATCGGATAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTa                                            <  1:108989/48‑1 (MQ=255)
                                     ttttGCATCGGATAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTa                                            <  1:1328291/48‑1 (MQ=255)
                                     ttttGCATCGGATAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTa                                            <  1:1385321/48‑1 (MQ=255)
                                     ttttGCATCGGATAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTa                                            <  1:1513064/48‑1 (MQ=255)
                                     ttttGCATCGGATAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTa                                            <  1:1739439/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)
                                     ttttGCATCGGATAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTa                                            <  1:229415/48‑1 (MQ=255)
                                              ggatAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTAGTGAGCACAATATTCTGTTTGATTTTGCAc              >  1:2162444/4‑69 (MQ=255)
                                              ggatAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTAGTGAGCACAATATTCTGTTTGATTTTGCAc              >  1:2067525/4‑69 (MQ=255)
                                              ggatAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTAGTGAGCACAATATTCTGTTTGATTTTGCAc              >  1:46366/4‑69 (MQ=255)
                                              ggatAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTAGTGAGCACAATATTCTGTTTGATTTTGCAc              >  1:3402962/4‑69 (MQ=255)
                                              ggatAGATGGTTATCGGCGGCCTATATTGCTGAGCCGTAGTGAGCACAATATTCTGTTTGATTTTGCAc              >  1:1495543/4‑69 (MQ=255)
                                                 tAGATGGTTATCGCCGACCCATATTACTGAGCCGTAGTGAGCACAATATTCTGTTTGATTTTGCACGtt           <  1:1335104/69‑1 (MQ=255)
                                                      ggTTATCGCCGACCCATATTACTGAGCCGTAGTGAGCACAATATTCTGTTTGATTTTg                 >  1:2575882/1‑58 (MQ=255)
                                                      ggTTATCGCCGACCCATATTACTGAGCCGTAGTGAGCACAATATTCTGTTTGATTTTg                 >  1:848971/1‑58 (MQ=255)
                                                      ggTTATCGCCGACCCATATTACTGAGCCGTAGTGAGCACAATATTCTGTTTGATTTTg                 >  1:1215901/1‑58 (MQ=255)
                                                      ggTTATCGCCGACCCATATTACTGAGCCGTAGTGAGCACAATATTCTGTTTGATTTTg                 >  1:2426632/1‑58 (MQ=255)
                                                      ggTTATCGCCGACCCATATTACTGAGCCGTAGTGAGCACAATATTCTGTTTGATTTTg                 >  1:135888/1‑58 (MQ=255)
                                                      ggTTATCGCCGACCCATATTACTGAGCCGTAGTGAGCACAATATTCTGTTTGATTTTg                 >  1:3363384/1‑58 (MQ=255)
                                                      ggTTATCGCCGACCCATATTACTGAGCCGTAGTGAGCACAATATTCTGTTTGATTTTg                 >  1:149634/1‑58 (MQ=255)
                                                      ggTTATCGCCGACCCATATTACTGAGCCGTAGTGAGCACAATATTCTGTTTGATTTTg                 >  1:2276658/1‑58 (MQ=255)
                                                      ggTTATCGCCGACCCATATTACTGAGCCGTAGTGAGCACAATATTCTGTTTGATTTTg                 <  1:1891845/58‑1 (MQ=255)
                                                      ggTTATCGCCGACCCATATTACTGAGCCGTAGTGAGCACAATATTCTGTTTGATTTTg                 >  1:2781839/1‑58 (MQ=255)
                                                      ggTTATCGCCGACCCATATTACTGAGCCGTAGTGAGCACAATATTCTGTTTGATTTTg                 <  1:2861675/58‑1 (MQ=255)
                                                          aTCGCCGACCCATATTACTGAGCCGTAGTGAGCACAATATTCTGTTTGATTTTGCACGTTCACGCcgcg  <  1:2696447/69‑1 (MQ=255)
                                                          aTCGCCGACCCATATTACTGAGCCGTAGTGAGCACAATATTCTGTTTGATTTTGCACGTTCACGCcgcg  <  1:273883/69‑1 (MQ=255)
                                                          aTCGCCGACCCATATTACTGAGCCGTAGTGAGCACAATATTCTGTTTGATTTTGCACGTTCACGCcgcg  <  1:2245514/69‑1 (MQ=255)
                                                                    |                                                          
ATTCGCGATGTATTTGATTTCCAGTCACATACCCATTTTTTGCATCGGGTAGATGGTTATCGCCGACCCATATTACTGAGCCGTAGTGAGCACAATATTCTGTTTGATTTTGCACGTTCACGCCGCG  >  minE/130694‑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: