Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,719,058 G→A P182L (CCG→CTG)  anmK ← anhydro‑N‑acetylmuramic acid kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,719,0580GA93.8% 84.7 / ‑1.4 32P182L (CCG→CTG) anmKanhydro‑N‑acetylmuramic acid kinase
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (17/13);  total (18/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.04e-01

AGCAGTGGGAGAATAACTTTACCCGCCCGTGCCCACTCGGCATCTTTATCGTAAGGTTTACCGGCCTGACGCCAGATCCAGGCATCCATCAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCGGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATTCGTCGCTCGGTTGGGTGAGCCAGCAGCGCATGATGGAACGCAGGTACCAGCGGTGCGCCTTGTCCTCCCAAGGCA  >  NC_000913/1718932‑1719189
                                                                                                                              |                                                                                                                                   
agcagTGGGAGAATAACTTTACCCGCCCGTGCCCACTCGGCATCTTTATCGTAAGGTTTACCGGCCTGACGCCAGATCCAGGCATCCATCAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCGGCTGCCCAGGaa                                                                                                                         >  2:124691/1‑139 (MQ=255)
    gTGGGAGAATAACTTTACCCGCCCGTGCCCACTCGGCATCTTTATCGTAAGGTTTACCGGCCTGACGCCAGATCCAGGCATCCATCAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAg                                                                                                                     >  2:369466/1‑139 (MQ=255)
       ggAGAATAACTTTACCCGCCCGTGCCCACTCGGCATCTTTATCGTAAGGTTTACCGGCCTGACGCCAGATCCAGGCATCCATCAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAg                                                                                                                  >  1:112774/1‑139 (MQ=255)
                      ccgcccgTGCCCACTCGGCATCTTTATCGTAAGGTTTACCGGCCTGACGCCAGATCCAGGCATCCATCAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGAt                                                                                                   <  2:905286/139‑1 (MQ=255)
                        gcccgTGCCCACTCGGCATCTTTATCGTAAGGTTTACCGGCCTGACGCCAGATCCAGGCATCCATCAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATgc                                                                                                 <  2:1104858/139‑1 (MQ=255)
                             tGCCCACTCGGCATCTTTATCGTAAGGTTTACCGGCCTGACGCCAGATCCAGGCTTCCATCAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCCGTGACTGATTGGCGTTGCCGCCa                                                                                            >  2:1026949/1‑139 (MQ=255)
                                              tATCGTAAGGTTTACCGGCCTGACGCCAGATCCAGGCATCCATCAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATtcg                                                                           <  1:753593/139‑1 (MQ=255)
                                              tATCGTAAGGTTTACCGGCCTGACGCCAGATCCAGGCATCCATCAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATtcg                                                                           >  2:753593/1‑139 (MQ=255)
                                                tCGTAAGGTTTACCGGCCTGACGCCAGATCCAGGCATCCATCAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATtcgtc                                                                         <  1:369466/139‑1 (MQ=255)
                                                tCGTAAGGTTTACCGGCCTGACGCCAGATCCAGGCATCCATCAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATtcgtc                                                                         <  2:941209/139‑1 (MQ=255)
                                                  gTAAGGTTTACCGGCCTGACGCCAGATCCAGGCATCCATCAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATTCGTCGc                                                                       <  2:112774/139‑1 (MQ=255)
                                                          tACCGGCCTGACGCAAGATCCAGGCATCCATCAGCATGTTGCCAGGACCGGTATCGTAGCCCCCAACCAGCTACCCAGGAATGAGCAGTGACAGATTGGAGATGCCGCCAATATTGAGCACCATTCGTCGCTCGGTTgg                                                               >  1:676821/1‑139 (MQ=255)
                                                           aCCGGCCTGACGCCAGATCCAGGCATCCATCAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATTCGTTGCTCGGTTggg                                                              <  2:1000241/139‑1 (MQ=255)
                                                                 cTGACGCCAGATCCAGGCATCCATCAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATTCGTCGCTCGGTTGGGTGAGcc                                                        >  2:206858/1‑139 (MQ=255)
                                                                  tGACGCCAGATCCAGGCATCCATCAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGcc                                                                                                                               <  1:730373/67‑1 (MQ=255)
                                                                  tGACGCCAGATCCAGGCATCCATCAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGcc                                                                                                                               >  2:730373/1‑67 (MQ=255)
                                                                  tGACGCCAGATCCAGGCATCCATCAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATTCGTCGCTCGGTTGGGTGAGCca                                                       >  2:923154/1‑139 (MQ=255)
                                                                  tGACGCCAGATCCAGGCATCCATCAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATTCGTCGCTCGGTTGGGTGAGCca                                                       >  1:950664/1‑139 (MQ=255)
                                                                   gACGCCAGATCCAGGCATCCATCAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATTCGTCGCTCGGTTGGGTGAGCcag                                                      <  1:206858/139‑1 (MQ=255)
                                                                    aCGCCAGATCCAGGCATCCATCAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATTCGTCGCTCGGTTGGGTGAGCcagc                                                     >  1:1101344/1‑139 (MQ=255)
                                                                                    tccatcAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAg                                                                                                                     <  2:949509/59‑1 (MQ=255)
                                                                                    tccatcAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAg                                                                                                                     >  1:949509/1‑59 (MQ=255)
                                                                                    tccatcAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATTCGTCGCTCGGTTGGGTGAGCCAGCAGCGCATGATGGAACg                                     >  1:1102500/1‑139 (MQ=255)
                                                                                     ccatcAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCACTATTGa                                                                                     >  2:640864/1‑90 (MQ=255)
                                                                                     ccatcAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGa                                                                                     <  1:640864/90‑1 (MQ=255)
                                                                                      catcAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCGGCTGCCCAGGAATGAGCAGAGACAGATTGGCGATGCCGCCAATATTGAGCACCATTCGTCGCTCGGTTGGGTGAGCCAGCAGCGCATGATGGAACGCa                                   <  1:124691/139‑1 (MQ=255)
                                                                                       atcAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATtcgtcg                                                                        >  2:502109/1‑101 (MQ=255)
                                                                                       atcAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATtcgtcg                                                                        <  1:502109/101‑1 (MQ=255)
                                                                                                            gTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATTCGTCGCTCGGTTGGGTGAGCCAGCAGCGCatg                                             >  1:864027/1‑107 (MQ=255)
                                                                                                            gTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATTCGTCGCTCGGTTGGGTGAGCCAGCAGCGCatg                                             <  2:864027/107‑1 (MQ=255)
                                                                                                                       cccAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATTCGTCGCTCGGTTGGGTGAGCCAGCAGCGCATGATGGAACGCAGGTACCAGCGGTGCGCCTTGTCCTCCCAAGGCa  >  2:394337/1‑139 (MQ=255)
                                                                                                                       cccAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATTCGTCGCTCGGTTGGGTGAGCCAGCAGCGCATGATGGAACGCAGGTACCAGCGGTGCGCCTTGTCCTCCCAAGGCa  >  2:5865/1‑139 (MQ=255)
                                                                                                                              |                                                                                                                                   
AGCAGTGGGAGAATAACTTTACCCGCCCGTGCCCACTCGGCATCTTTATCGTAAGGTTTACCGGCCTGACGCCAGATCCAGGCATCCATCAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCGGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATTCGTCGCTCGGTTGGGTGAGCCAGCAGCGCATGATGGAACGCAGGTACCAGCGGTGCGCCTTGTCCTCCCAAGGCA  >  NC_000913/1718932‑1719189

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: