Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 2,203,709 T→C 12.9% K442K (AAA→AAG)  murF ← UDP‑N‑acetylmuramoyl‑tripeptide‑‑D‑alanyl‑D‑ alanine ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007302,203,7090TC12.9% 70.0 / 4.4 31K442K (AAA→AAG) murFUDP‑N‑acetylmuramoyl‑tripeptide‑‑D‑alanyl‑D‑ alanine ligase
Reads supporting (aligned to +/- strand):  ref base T (12/15);  new base C (2/2);  total (14/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.18e-01

AAAAGAAAAAGGCAATTTTGTATTTACACACAAAACCACCTCTTCTTTGTAATCTTTTAGTATTAATTCAAATCATTTTTATAAGTAAAAGGTCCCTCGACTAATCTCTATGAAATTAAAGCATTTACCACTTCTTCTAATTTCATACCACGTGATCCTTTAACTAATACACGGTCATGCGCTTTTAAATCGTTTATTAAAACTTCTATCATATCGTCTTTAGAATTGAAGTGTTGTGCTTTTTCGACATGTTGCTGGCCCGAATCATAAATATACT  >  CP000730/2203569‑2203845
                                                                                                                                            |                                                                                                                                        
aaaagaaaaagGCAATTTTGTATTTACACACAAAACCACCTCTTCTTTGTAATCTTTTAGTATTAATTCAAATCATTTTTATAAGTAAAAGGTCCCTCGACTAATCTCTATGAAATTAAAGCATTTACCACTTCTTCTAAt                                                                                                                                          >  2:366074/1‑141 (MQ=255)
          ggCAATTTTGTATTTACACACAAAACCACCTCTTCTTTGTAATCTTTTAGTATTAATTCAAATCATTTTTATAAGTAAAAGGTCCCTCGACTAATCTCTATGAAATTAAAGCATTTACCACTTCTTCTAATTTCATACCAc                                                                                                                                >  2:169001/1‑141 (MQ=255)
                     aTTTACACACAAAACCACCTCTTCTTTGTAATCTTTTAGTATTAATTCAAATCATTTTTATAAGTAAAAGGTCCCTCGACTAATCTCTATGAAATTAAAGCATTTACCACTTCTTCTAATTTCATACCACGTGATCCTTTa                                                                                                                     <  2:153817/141‑1 (MQ=255)
                           acacaAAACCACCTCTTCTTTGTAATCTTTTAGTATTAATTCAAATCATTTTTATAAGTAAAAGGTCGCTCGACTAATCTCTATGAAATTAAAGCATTTACAACTTCTTCTAACTTCATACCACGTGATCCTTTAACTAAt                                                                                                               >  1:357092/1‑141 (MQ=255)
                            cacaAACCCACCTCTTGTTTGTAATCTGTTAGTATTAATTCAAATCATTTTTAGAAGTAAAAGGTCCCTCTACTAATCTCTATGAAATTAAAGGATTTACCTCTTCTACTAATTTCATACCACGTGGTCCTTTAACTAAtg                                                                                                              <  2:48261/141‑2 (MQ=255)
                                   ccaccGCTTCTTTGTAATCTTTTAGTATTAATTCAAATCATTTTTATAAGTAAAAGGTCCCTCGACTAATCTCTATGAAATTAAAGCATTTACCACTTCTTCTAATTTCATACCACGTGATCCTTTAACTAATACACGGTc                                                                                                       <  1:332497/141‑1 (MQ=255)
                                      ccTCTTCTTTGTAATCTTTTAGTATTAATTCAAATCATTTTTATAAGTAAAAGGTCCCTCGACTAATCTCTATGAAATTAAAGCATTTACCACTTCTTCTAATTTCATACCACGTGATCCTTTAACTAATACACGGTCATg                                                                                                    <  2:10780/141‑1 (MQ=255)
                                              ttGTAATCTTTTAGTATTAATTCAAATCATTTTTATAAGTAAAAGGTCCCTCGACTAATCTCTATGAAATTAAAGCATTTACCACTTCTTCTAATTTCATACCAAGTGATCCTTTAACTAATACACGGTCATGCGCTTTTa                                                                                            <  1:409856/141‑1 (MQ=255)
                                                gTAATCTTTTAGTATTAATTCAAATCATTTTTATAAGTAAAAGGTCCCTCGACTAATCTCTATGAAATTAAAGCATTTACCACTTCTTCTAATTTCATACCACGTGATCCTTTAACTAATACACGGTCATGCGCTTTTaaa                                                                                          >  1:101850/1‑141 (MQ=255)
                                                     cTTTTAGTATTAATTCAAATCATTTTTATAAGTAAAAGGTCCCTCGACTAATCTCTATGAAATTAAAGCATTTACCACTTCTTCTAATTTCATACCACGTGATCATTTAACTAATACACGGTCATGCGCTTTTAAATCGtt                                                                                     <  2:397386/141‑1 (MQ=255)
                                                                aaTTCAAATCATTTTTATAAGTAAAAGGTCCCTCGACTAATCTCTATGAAATTAAAGCATTTACCACTTCTTCTAATTTCATACCACGTGATCCTTTAACTAATACACGGTCATGCGCTTTTAAATCGTTTATTAAAACtt                                                                          >  2:63915/1‑141 (MQ=255)
                                                                 aTTCAAATCATTTTTATAAGTAAAAGGTGNCTCAACNANTCTCTAGGAAATTAAAGCATTTACCACTTCTTCTAATTTCATACTACGTGATCATGTAACTAATACACGTTCATGCGCGTTTGAATCGTTTATTAAAACTTc                                                                         >  1:194679/1‑141 (MQ=255)
                                                                   tCAAATCATTTTTATAAGTAAAAGGTCCCTCGACTAATCTCTATGAAATTAAAGCATTTACCACTTCTTCTAATTTCATACCACGTGATCCTTTAACTAATACACGGTCATGCGCTTTTAAATCGTTTATTAAAACTTCTa                                                                       <  1:102946/141‑1 (MQ=255)
                                                                     aaaTCATTTTTATAAGTAAAAGGTCCCTCGACTAATCTCTATGAAATTAAAGCATTTACCACTTCTTCTAATTTCATACCACGTGATCCTTTAACTAATACACGGTCATGCGCTTTTAAATCGTTTATTAAAACTTCTATc                                                                     <  1:176803/141‑1 (MQ=255)
                                                                       aTCATTTTTATAAGTAAAAGGTCCCTCGACTAATCTCTATGAAATTAAAGCATTTACCACTTCTTCTAATTTCATACCACGTGATCCTTTAACTAATACACGGTCATGCGCTTTTAAATCGTTTATTAAAACTTCTACCat                                                                   >  1:352170/1‑141 (MQ=255)
                                                                                  aaGTAAAAGGTCCCTCGACTAATCTCTATGAAATTAAAGCATTTACCACTTCTTCTAATTTCATACCACGTGATCCTTTAACTAATACACGGTCATGCGCTTTTAAATCGTTTATTAAAACTTCTATCATATCGTCTTTAg                                                        >  2:152342/1‑141 (MQ=255)
                                                                                       aaaGGTCCCTCGACTAATCTCTATGAAATTAAAGCATTTACCACTTCTTCTAATTTCATACCACGTGATCCTTTAACTAATACACGGTCATGCGCTTTTAAATCGTTTATTAAAACTTCTATCATATCGTCTTTAGAATTg                                                   >  2:351313/1‑141 (MQ=255)
                                                                                         aGGACCCTCGACTAATCTCTATGAAATTAAAGCAATTACCACTTCTTCTAATGTCAGACCACTTGAGCCTTGAACTAATACACGGTCATGCGCTTTTAAATCGCTGAATAAAACTTCTAACAAATGGTCTTTAGAATTGag                                                 <  2:153376/141‑2 (MQ=255)
                                                                                               cTCGACTAATCTCTATGAAATTAAAGCATTTACCACTTCTTCTAATTTCATACCACGTGATCCTTTAACTAATACACGGTCATGCGCTTTTAAATCGTATATTAAAACTTCTATCATATCGTCTTTAGAATTGAAGtgttg                                           <  2:41277/141‑1 (MQ=255)
                                                                                                     tAATCTCTATGAAATTAAAGCATTTACCACTTCTTCTAATTTCATACCACGTGATCCTTTAACTAATACACGGTCATGCGCTTTTAAATCGTTTATTAAAACTTCTATCATATCGTCTTTAGAATTGAAGTGTTGTGCttt                                     <  2:422164/141‑1 (MQ=255)
                                                                                                       aTCTCTATGAAATTAAAGCATTTACCACTTCTTCTAATTTCATACCACGTGATCCTTTAACTAATACACGGTCATGCGCTTTTAAATCGTTTATTAAAACTTCTATCATATCGTCTTTAGAATTGAAGTGTTGTGCttttt                                   <  1:366074/141‑1 (MQ=255)
                                                                                                        tctcTATGAAATTAAAGCATTTACCACTTCTTCTAATTTCATACCACGTGATCCTTTAACTAATACACGGTCATGCGCTTTTAAATCGTTTATTAAAGCTTCTATCATATCGTCTTTAGAATTGAAGTTTTTTGCTTTTTc                                  >  2:424285/1‑141 (MQ=255)
                                                                                                            tATGAAATTAAAGCATTTACCACTTCTTCTAATTTCATACCACGTGATCCTTTAACTAATACACGGTCATGCGCTTTTAAATCGTTTATTAAAACTTCTATCATATCGTCTTTAGAATTGAAGTGTTGTGCTTTTTCGACa                              <  1:187140/141‑1 (MQ=255)
                                                                                                                 aaTTAAAGCATTTACCACTTCTTCTAATTTCATACCACGTGATCCTTTAACTAATACACGGTCATGCGCTTTTAAATCGTTTATTAAAACTTCTATCATATCGTCTTTAGAATTGAAGTGTTGTGCTTTTTCGACATGTTg                         >  2:238986/1‑141 (MQ=255)
                                                                                                                     aaaGCATTTACCACTTCTTCTAATTTCATACCACGTGATCCTTTAACTAATACACGGTCATGCGCTTTTAAATCGTTTATTAAAACTTCTATCATATCGTCTTTAGAATTGAAGTGTTGTGCTTTTTCGACATGTTGCTgg                     <  1:161264/141‑1 (MQ=255)
                                                                                                                        gCATTTACCACTTCTTCTAATTTCATACCACGTGATCCTTTAACTAATACACGGTCATGCGCTTTTAAATCGTTTATTAAAACTTCTATCATATCGTCTTTAGAATTGAAGTGTTGTGCTTTTTCGACATGTTGCTGGccc                  <  2:278620/141‑1 (MQ=255)
                                                                                                                           tttACCACTTCTTCTAATTTCATACCACGTGATCCTTTAACTAATACACGGTCATGCGCTTTTAAATCGTTTATTAAAACTTCTATCATATCGTCTTTAGAATTGAAGTGTTGTGCTTTTTCGACATGTTGCTGGCCCGaa               >  1:344157/1‑141 (MQ=255)
                                                                                                                           tttACCACTTCTTCTAATTTCATACCACGTGATCCTTTAACTAATACACGGTCATGCGCTTTTAAATCGTTTATTAAAACTTCTATCATATCGTCTTTAGAATTGAAGTGTTGTGCTTTTTCGACATGTTGCTGGCCCGaa               >  1:343386/1‑141 (MQ=255)
                                                                                                                               caaCTTCTTCTAACTTCATACCACGTGATCCTTTAACTAATACACGGTCATGCGCTTTTAAATCGCTTGTTAAAACTTCTATCATATCGTCTTTAGAATTGAAGTGTTGTGCTTTTTCGACATGTTGCTGGCCCGAATCAt           >  2:27733/3‑141 (MQ=255)
                                                                                                                                   ttcttcTAACTTCATACCACGTGATCCTTTAACTAATACACGGTCATGCGCTTTTAAATCGCTTGTTAAAACTTCTATCATATCGTCTTTAGAATTGAAGTGTTGTGCTTTTTCGACATGTTGCTGGCCCGAATCATAAat       <  2:357092/141‑1 (MQ=255)
                                                                                                                                        ctAATTTCATACCACGTGATCCTTTAACTAATACACGGTCATGCGCTTTTAAATCGTTTATTAAAACTTCTATCATATCGTCTTTAGAATTGAAGTGTTGTGCTTTTTCGACATGTTGCTGGCCCGAATCATAAATATACt  >  1:167188/1‑141 (MQ=255)
                                                                                                                                        ctAACTTCATACCACGTGATCCTTTAACTAATACACGGTCATGCGCTTTTAAATCGCTTGTTAAAACTTCTATCATATCGTCTTTAGAATTGAAGTGTTGTGCTTTTTCGACATGTTGCTGGCCCGAATCATAAATATACt  <  1:76513/141‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                        
AAAAGAAAAAGGCAATTTTGTATTTACACACAAAACCACCTCTTCTTTGTAATCTTTTAGTATTAATTCAAATCATTTTTATAAGTAAAAGGTCCCTCGACTAATCTCTATGAAATTAAAGCATTTACCACTTCTTCTAATTTCATACCACGTGATCCTTTAACTAATACACGGTCATGCGCTTTTAAATCGTTTATTAAAACTTCTATCATATCGTCTTTAGAATTGAAGTGTTGTGCTTTTTCGACATGTTGCTGGCCCGAATCATAAATATACT  >  CP000730/2203569‑2203845

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