Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 1,191,016 C→G 12.5% T241T (ACC→ACG)  mraY → phospho‑N‑acetylmuramoyl‑pentapeptide‑ transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007301,191,0160CG12.5% 91.2 / 4.8 31T241T (ACC→ACG) mraYphospho‑N‑acetylmuramoyl‑pentapeptide‑ transferase
Reads supporting (aligned to +/- strand):  ref base C (19/8);  new base G (2/2);  total (21/10)
Fisher's exact test for biased strand distribution p-value = 5.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGAAACGGCAATTGGTATTTTCTGTATCATTATGTTGTTTGCACTTTTAGGATTTTTACCATATAACATTAACCCTGCTAAAGTGTTTATGGGAGATACAGGTAGCTTAGCTTTAGGTGGTATATTTGCTACCATTTCAATCATGCTTAATCAGGAATTATCATTAATTTTTATAGGTTTAGTATTCGTAATTGAAACATTATCTGTTATGTTACAAGTCGCTAGCTTTAAATTGACTGGAAAGCGTATATTTAAAATGAGTCCGATTCATCA  >  CP000730/1190884‑1191156
                                                                                                                                    |                                                                                                                                            
agaAACGGCAATTGGTATTTTCTGTATCATTATGTTGTTTGCACTTTTAGGATTTTTACCATATAACATTAACCCTGCTAAAGTGTTTATGGGAGATACAGGTAGCTTAGCTTTAGGTGGTATATTTGCTACGATTTCAAt                                                                                                                                      >  2:273986/1‑141 (MQ=255)
  agaCGGCAATTGGTATTTTCTGTATCATTATGTTGTTTGCACTTTTAGGATTTTTACCATATAACATTAACCCTGCTAAAGTGTTTATGGGAGATACAGGTAGCTTAGCTTTAGGTGGTATATTTGCTACGATTTCAATCa                                                                                                                                    <  1:280026/139‑1 (MQ=255)
  aaaCGGCAATTGGTATTTTCTGTATCATTATGTTGTTTGCACTTTTAGGATTTTTACCATATAACATTAACCCTGCTAAAGTGTTTATGGGAGATACAGGTAGCTTAGCTTTAGGTGGTATATTTGCTACCATTTCAATCa                                                                                                                                    <  2:136946/141‑1 (MQ=255)
    aCGGCAATTGGTATTTTCTGTATCATTATGTTGTTTGCACTTTTAGGATTTTTACCATATAACATTAACCCTGCTAAAGTGTTTATGGGAGATACAGGTAGCTTAGCTTTAGGTGGTATATTTGCTACCATTTCAATCATg                                                                                                                                  >  1:305825/1‑141 (MQ=255)
        cAATTGGTATTTTCTGTATCATTATGTTGTTTGCACTTTTAGGATTTTTACCATATAACATTAACCCTGCTAAAGTGTTTATGGGAGATACAGGTAGCTTAGCTTTAGGTGGTATATTTGCTACCATTTCAATCATGCTTa                                                                                                                              >  1:456545/1‑141 (MQ=255)
                      tGTATCATTATGTTGTTTGCACTTTTAGGATTTTTACCATATAACATTAACCCTGCTAAAGTGTTTATGGGAGATACAGGTAGCTTAGCTTTAGGTGGTATATTTGCTACCATTTCAATCATGCTTAATCAGGAATTATCa                                                                                                                <  1:193841/141‑1 (MQ=255)
                          tCATTATGTTGTTTGCACTTTTAGGATTTTTACCATATAACATTAACCCTGCTAAAGTGTTTATGGGAGATACAGGTAGCTTAGCTTTAGGTGGTATATTTGCTACCATTTCAATCATGCTTAATCAGGAATTATCATTaa                                                                                                            >  2:387610/1‑141 (MQ=255)
                                    gttTGCACTTTTAGGATTTTTACCATATAACATTAACCCTGCTAAAGTGTTTATGGGAGATACAGGTAGCTTAGCTTTAGGTGGTATATTTGCTACCGTTTCAATCATGCTTAATCAGGAATTATCATTAATTTTTATAgg                                                                                                  >  2:425497/1‑141 (MQ=255)
                                    gttTGCACTTTTAGGATTTTTACCATATAACATTAACCCTGCTAAAGTGTTTATGGGAGATACAGGTAGCTTAGCTTTAGGTGGTATATTTGCTACCATTTCAATCATGCTTAATCAGGAATTATCATTAATTTTTATAgg                                                                                                  >  1:95932/1‑141 (MQ=255)
                                      ttGCACTTTTAGGATTTTTACCATATAACATTAACCCTGCTAAAGTGTTTATGGGAGATACAGGTAGCTTAGCTTTAGGTGGTATATTTGCTACCATTTCAATCATGCTTAATCAGGAATTATCATTAATTTTTATAGGtt                                                                                                >  1:110202/1‑141 (MQ=255)
                                            ttttAGGATTTTTACCATATAACATTAACCCTGCTAAAGTGTTTATGGGAGATACAGGTAGCTTAGCTTTAGGTGGTATATTTGCTACCATTTCAATCATGCTTAATCAGGAATTATCATTAATTTTTATAGGTTTAGTAt                                                                                          <  2:88715/141‑1 (MQ=255)
                                                 ggATTTTTACCATATAACATTAACCCTGCTAAAGTGTTTATGGGAGATACAGGTAGCTTAGCTTTAGGTGGTATATTTGCTACCATTTCAATCATGCTTACTCAGGAATTATCATTAATTTTTATAGGTTTAGTATTCGTa                                                                                     >  2:338973/1‑141 (MQ=255)
                                                 ggATTTTTACCATATAACATTAACCCTGCTAAAGTGTTTATGGGAGATACAGGTAGCTTAGCTTTAGGTGGTATATTTGCTACCATTTCAATCATGCTTAATCAGGAATTATCATTAATTTTTATAGGTTTAGTATTCGTa                                                                                     >  2:298149/1‑141 (MQ=255)
                                                   aTTTTTACCATATAACATTAACCCTGCTAAAGTGTTTATGGGAGATACAGGTAGCTTAGCTTTAGGTGGTATATTTGCTACCATTTCAATCATGCTTAATCAGGAATTATCATTAATTTTTATAGGTTTAGTATTCGTAAt                                                                                   >  1:196403/1‑141 (MQ=255)
                                                                    ttAACCCTGCTAAAGTGTTTATGGGAGATACAGGTAGCTTAGCTTTAGGTGGTATATTTGCTACCATTTCAATCATGCTTAATCAGGAATTATCATTAATTTTTATAGGTTTAGTATTCGTAATTGAAACATTATCtgtta                                                                  >  1:107981/1‑141 (MQ=255)
                                                                     tAACCCTGCTAAAGTGTTTATGGGAGATACAGGTAGCTTAGCTTTAGGTGGTATATTTGCTACCATTTCAATCATGCTTAATCAGGAATTATCATTAATTTTTATAGGTTTAGTATTCGTAATTGAAACATTATCtgttat                                                                 >  1:261204/1‑141 (MQ=255)
                                                                     tAACCCTGCTAAAGTGTTTATGGGAGATACAGGTAGCTTAGCTTTAGGTGGTATATTTGCTACCATTTCAATCATGCTTAATCAGGAATTATCATTAATTTTTATAGGTTTAGTATTCGTAATTGAAACATTATCtgttat                                                                 >  1:34518/1‑141 (MQ=255)
                                                                         ccTGCTAAAGTGTTTATGGGAGATACAGGTAGCTTAGCTTTAGGTGGTATATTTGCTACCATTTCAATCATGCTTAATCAGGAATTATCATTAATTTTTATAGGTTTAGTATTCGTAATTGAAACATTATCtgttatgtta                                                             >  1:403585/1‑141 (MQ=255)
                                                                            gCTAAAGTGTTTATGGGAGATACAGGTAGCTTAGCTTTAGGTGGTATATTTGCTACCATTTCAATCATGCTTAATCAGGAATTATCATTAATTTTTATAGGTTTAGTATTCGTAATTGAAACATTATCTGTTATGTTACaa                                                          <  2:306082/141‑1 (MQ=255)
                                                                             cTAAAGTGTTTATGGGAGATACAGGTAGCTTAGCTTTAGGTGGTATATTTGCTACCATTTCAATCATGCTTAATCAGGAATTATCATTAATTTTTATAGGTTTAGTATTCGTAATTGAAACATTATCTGTTATGTTACAAg                                                         >  1:393166/1‑141 (MQ=255)
                                                                                    ttatATGGGAGATACAGGTAGCTTAGCTTTAGGTGGTATATTTGCTACGATTTCAATCATGCTTAATCAGGAATTATCATTAATTTTTATAGGTTTAGTATTCGTAATTGAAACATTATCTGTTATGTTACAAGTCGCTAg                                                  <  1:122638/138‑1 (MQ=255)
                                                                                    gtTTATGGGAGATACAGGTAGCTTAGCTTTAGGTGGTATATTTGCTACCATTTCAATCATGCTTAATCAGGAATTATCATTAATTTTTATAGGTTTAGTATTCGTAATTGAAACATTATCTGTTATGTTACAAGTCGCTAg                                                  >  2:277135/1‑141 (MQ=255)
                                                                                                  cAGGTGGCTAAGCATTAGGGGGTCTGTAGGCTACGATATAAATGGTACTTACTCCGGCCTTACCCTTCTGTTTTCTATTTTTCGTATTCGTACTTGAAACGTTATATGTTATGTTACAAGTCGCTAGCTTTGAAGTGACTg                                    >  2:443200/1‑141 (MQ=255)
                                                                                                   aGGTAGCTTAGCTTTAGGTGGTATATTTGCTACCATTTCAATCATGCTTAATCAGGAATTATCATTAATTTTTATAGGTTTAGTATTCGTAATTGAAACATTATCTGTTATGTTACAAGTCGCTAGCTTTAAATTGACTgg                                   <  2:289733/141‑1 (MQ=255)
                                                                                                   aGGTAGCTTAGCTTTAGGTGGTATATTTGCTACCATTTCAATCATGCTTAATCAGGAATTATCATTAATTTTTATAGGTTTAGTATTCGTAATTGAAACATTATCTGTTATGTTACAAGTCGCTAGCTTTAAATTGACTgg                                   <  2:290022/141‑1 (MQ=255)
                                                                                                             gcttTAGGTGGTATATTTGCTACCATTTCAATCATGCTTAATCAGGAATTATCATTAATTTTTCTAGGTTTAGTATTAGTACTTGCATACCTCTCCGTCATGTTACAAGTCGCCAGCTTTAAATTGACTGGAAAGCGtata                         >  1:385072/1‑141 (MQ=255)
                                                                                                              cttTAGGTGGTATATTTGCTACCATTTCAATCATGCTTAATCAGGAATTATCATTAATTTTTATAGGTTTAGTATTCGTAATTGAAACATTATCTGTTATGTTACAAGTCGCTAGCTTTAAATTGACTGGAAAGCGTatat                        >  1:439632/1‑141 (MQ=255)
                                                                                                                 tacgtggtATATTTGCTACCATTTCAATCATGCTTAATCAGGAATTATCATTAATTTTTATAGGTTTAGTATTCGTAATTGAAACATTATCTGTTATGTTACAAGTCGCTAGCTTTAAATTGACTGGAAAGCGTATATTTa                     <  2:393166/138‑1 (MQ=255)
                                                                                                                   ggtggtATATTTGCTACCATTTCAATCATGCTTAATCAGGAATTATCATTAATTTTTATAGGTTTAGTATTCGTAATTGAAACATTATCTGTTATGTTACAAGTCGCTAGCTTTAAATTGACTGGAAAGCGTATATTTaaa                   >  1:450298/1‑141 (MQ=255)
                                                                                                                    aggttcaatctcggaCCATTTCAATCATGCTTAATCAGGAATTATCATTAATTTTTATAGGTTTAGTATTCGTAATTGAAACATTATCTGTTATGTTACAAGTCGCTAGCTTTAAATTGACTGGAAAGCGTATATTTaaaa                  <  1:232072/127‑1 (MQ=255)
                                                                                                                              tGCTACCATTTCAATCATGCTTAATCAGGAATTATCATTAATTTTTATAGGTTTAGTATTCGTAATTGAAACATTATCTGTTATGTTACAAGTCGCTAGCTTTAAATTGACTGGAAAGCGTATATTTAAAATGAGTCCGAt        >  2:278832/1‑141 (MQ=255)
                                                                                                                                    cATTTCAATCATGCTTAATCAGGAATTATCATTAATTTTTATAGGTTTAGTATTCGTAATTGAAACATTATCTGTTATGTTACAAGTCGCTAGCTTTAAATTGACTGGAAAGCGTATATTTAAAATGAGTCCGATtcatca  >  2:273496/1‑141 (MQ=255)
                                                                                                                                    |                                                                                                                                            
AGAAACGGCAATTGGTATTTTCTGTATCATTATGTTGTTTGCACTTTTAGGATTTTTACCATATAACATTAACCCTGCTAAAGTGTTTATGGGAGATACAGGTAGCTTAGCTTTAGGTGGTATATTTGCTACCATTTCAATCATGCTTAATCAGGAATTATCATTAATTTTTATAGGTTTAGTATTCGTAATTGAAACATTATCTGTTATGTTACAAGTCGCTAGCTTTAAATTGACTGGAAAGCGTATATTTAAAATGAGTCCGATTCATCA  >  CP000730/1190884‑1191156

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