Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,630,984 G→T T84T (ACG→ACT slt → soluble lytic murein transglycosylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,630,9840GT92.6% 64.3 / ‑3.2 27T84T (ACG→ACTsltsoluble lytic murein transglycosylase
Reads supporting (aligned to +/- strand):  ref base G (0/2);  new base T (13/12);  total (13/14)
Fisher's exact test for biased strand distribution p-value = 4.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.08e-01

CAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGTCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAGGATTATCCGCTTTATCCCTACCTGGAATACCGCCAGATCACCGATGATCTGATGAATCAACCGGCGGTGACGGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCGAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCT  >  NC_000913/4630744‑4631202
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                           
cAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGAAGCGTAGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGATACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGCGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACATGATGAATCAACCGGCGGTGACTGTCTCTAACt                                                                                                                                                                                                                   >  1:182575/1‑251 (MQ=255)
     aaGTTTCCTGGCGGCTGTTGTCGGCCGGTGTCTGGCTGCTGACGGTCAGCAGCGTGGCGCGAGCCTACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACGTTATCCCTACCTGGAATACCGTCAGATCACCGTTGGCCTGATGAATCAACCGGCGGTGATTGTCACGAACTTTGtt                                                                                                                                                                                                              <  2:192554/251‑1 (MQ=255)
                    tGTTGGCTGCCGGTGACTGTCTGCTGACGGTCAGAAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGGGACTGTCACTAACTTTGTTCGCGCTAAACCCACg                                                                                                                                                                                               >  1:109147/1‑251 (MQ=255)
                           tGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTtcctc                                                                                                                                                                                        >  1:291028/1‑251 (MQ=255)
                                   tctgtctgCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGc                                                                                                                                                                                >  1:33422/1‑251 (MQ=255)
                                             gACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCa                                                                                                                                                                        <  1:344879/249‑1 (MQ=255)
                                                  tCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCACTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGGCTGTCACTAACTTTGTTCGCGCAAACCCCACGCTTCCTCCCACTCGCACGCTGCAATCTCGt                                                                                                                                                                 >  1:114909/1‑251 (MQ=255)
                                                                 gAGCCGACGCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACtt                                                                                                                                                  <  1:235210/251‑2 (MQ=255)
                                                                              ggATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAAtctc                                                                                                                                                                   >  1:148927/1‑221 (MQ=255)
                                                                                aTGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGAc                                                                                                                                   <  1:146205/251‑1 (MQ=255)
                                                                                 tGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAAtctc                                                                                                                                                                   >  1:347843/1‑218 (MQ=255)
                                                                                 ggAGCAGCGTAGCCGTGACGCGCAAATCAAGCAGGCCTGGTATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCTCTTTATCCCTACCTGGAATACCGTCAGATCTCCGATGACCTGATGAAGCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAAtctc                                                                                                                                                                   <  2:347843/217‑1 (MQ=255)
                                                                                         gTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGTTAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCAACTTTATCCCTACCTGGAATACCATCAGATCACAGATGACCTGATGAATCAACCGGCGGTGACTGTCCCTAACTTTTTTCGCGCTAAACCCACGCTTCCCCACACCCGCCCTCCGAAATCCCCTTTCGTCACTGAACTGGCGCGGCGTGAAGaacagcgtggg                                                                                                                          >  2:160982/1‑241 (MQ=255)
                                                                                                          aTCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGccc                                                                                                         <  1:47889/251‑1 (MQ=255)
                                                                                                                        ggATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGgc                                                                                                                                                <  1:52515/198‑1 (MQ=255)
                                                                                                                                 aCAAATGGATGTGGTCGAACAAATGATGTCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCCTGTTAGCCTTTAGCCCGGGCAAAGCCGGAACTACCGa                                                                                   >  1:379634/1‑249 (MQ=255)
                                                                                                                                    aaTGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTTGCCTTTAGCCCGGGAAAGCCCGGAAATACCCAAgcgc                                                                              >  1:238817/1‑251 (MQ=255)
                                                                                                                                                  aacaaatatgtcttGGACTGAAGGATGTTTCGCTCTACCCTTACCTGGAATACCGCCAGTTCCCCGATGTTCTGATGAATCAACCGGTGGCGACGTTCACTAACTTTGTGCGCGCTAACCCCACGCTTCCTCCCGCTCGCAGGCTGCAATCTCGTTTCGTCAATGATCTGGCGCGGCGTGAAGACTGTCGTGTCTTGTTAGCCTTTATCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATtac                                                                   <  2:131987/235‑1 (MQ=255)
                                                                                                                                                    aaaaTGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATTCCGTCAGATCACCGAGGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTGTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACTCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTActacta                                                              <  1:160982/250‑1 (MQ=255)
                                                                                                                                                                ggACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACAGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCTATGTAATTACTACTATGCAAAATGGa                                                   <  1:339407/250‑1 (MQ=255)
                                                                                                                                                                   cTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAACCCGGAACTACCGAAGCGCCATGTAATTAATACTATGCAAAATGGaaaac                                               >  1:237233/1‑248 (MQ=255)
                                                                                                                                                                       aagATTTTCCACTTTATCCCTGCCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCGCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACCggg                                           <  2:76005/249‑1 (MQ=255)
                                                                                                                                                                               ccacTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACGTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTga                                   <  2:237233/248‑1 (MQ=255)
                                                                                                                                                                                caatgtagcccgaCCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCTCGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTgaa                                  <  2:152035/239‑1 (MQ=255)
                                                                                                                                                                                             ccTGGAATACCGCCAGATCACCGATGATCTGATGAATCAACCGGCGGTGACGGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCGAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAg                     <  1:166629/251‑1 (MQ=255)
                                                                                                                                                                                                           aGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAAAGGAAAACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGGGAAagagatg       >  1:147866/1‑248 (MQ=255)
                                                                                                                                                                                                                 ccGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGGTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCTTGGGTTGTTTGACTTTAGACCGGGAAAGACCGGGAATCCCGAAACGCAATGGTATTACTACTATGCAAAATTGAACACCGGGCAGCGGGAAGAAGCCTGGCCAGGGGGGGAAGAGCTGTGGCt  >  1:130496/1‑250 (MQ=255)
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                           
CAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGTCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAGGATTATCCGCTTTATCCCTACCTGGAATACCGCCAGATCACCGATGATCTGATGAATCAACCGGCGGTGACGGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCGAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCT  >  NC_000913/4630744‑4631202

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