Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,867,926 (C)6→7 coding (828/1140 nt) nlpD ← murein hydrolase activator NlpD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,867,9201.C96.7% 100.4 / ‑3.1 30coding (834/1140 nt)nlpDmurein hydrolase activator NlpD
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base C (17/12);  total (18/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.93e-01

ACTCAGGTAATCATCATTATGTTTGATGATAATCAGATTACCGTAGCCGCGCAGCGCGTTACCAGCATAAACAACGCGGCCATCTGCGGTCGCGATAATTGCCTGTCCTTTGCTGCCTGCGATATCAATCCCCTTGTTGCCCCCCTCAGAAGCGCCAAAGGTTTCGATCACTTTGCCCTCAGTCGGCCAGCGCCAGGTGGAGATAGGCGTACTGGTTGATGTACTGCTGACAGTCGGCTCGGTTGTGCTTGCTGTTGGTACCGTTACAGGCGC  >  NC_000913/2867782‑2868054
                                                                                                                                           |                                                                                                                                       
aCTCAGGTAATCATCATTATGTTTGATGATAATCAGATTACCGTAGCCGCGCAGCGCGTTACCAGAATAAACAACGCGGCCATCTGCGGTCGCGATAATTGCCTGTCCTTTGCTGCCTGCGATATCAATCCCCTTGTTGc                                                                                                                                        <  2:792511/140‑1 (MQ=255)
aCTCAGGTAATCATCATTATGTATGATGATAATCAGAATACCGTAGCCGCGCAGCGCGTTACCAGCATAAACAACGCGGCCATCTGCGGTCGCGATAATTGCCTGTCCTTTGCTGCCTGCGATATCAATCCCCTTGTTGc                                                                                                                                        <  1:324880/140‑1 (MQ=255)
           catcatTATGTTTGATGATAATCAGATTACCGTAGCCGCGCAGCGCGTTACCAGCATAAACAACGCGGCCATCTGCGGTCGCGATAATTGCCTGTCCTTTGCTGCCTGCGATATCAATCCCCTTGTTGCCCCCCCTCAGa                                                                                                                             <  2:554108/140‑1 (MQ=255)
                        gatgatAATCAGATTACCGTAGCCGCGCAGCGCGTTACCAGCATAAACAACGCGGCCATCTGCGGTCGCGATAATTGCCTGTCCTTTGCTGCCTGCGATATCAATCCCCTTGTTGCCCCCCCTCAGAAGCGCCAAAGGtt                                                                                                                >  1:833219/1‑140 (MQ=255)
                        gatgatAATCAGATTACCGTAGCCGCGCAGCGCGTTACCAGCATAAACAACGCGGCCATCTGCGGTCGCGATAATTGCCTGTCCTTTGCTGCCTGCGATATCAATCCCCTTGTTGCCCCCCCTCAGAAGCGCCAAAGGtt                                                                                                                >  1:601601/1‑140 (MQ=255)
                              aaTCAGATTACCGTAGCCGCGCAGCGCGTTACCAGCATAAACAACGCGGCCATCTGCGGTCGCGATAATTGCCTGTCCTTTGCTGCCTGCGATATCAATCCCCTTGTTGCCCCCCCTCAGAAGCGCCAAAGGTTTCGATc                                                                                                          <  1:594105/140‑1 (MQ=255)
                               aTCAGATTACCGTAGCCGCGCAGCGCGTTACCAGCATAAACAACGCGGCCATCTGCGGTCGCGATAATTGCCTGTCCTTTGCTGCCTGCGATATCAATCCCCTTGTTGCCCCCCCTCAGAAGCGCCAAAGGTTTCGATCa                                                                                                         >  1:653035/1‑140 (MQ=255)
                               aTCAGATTACCGTAGCCGCGCAGCGCGTTACCAGCATAAACAACGCGGCCATCTGCGGTCGCGATAATTGCCTGTCCTTTGCTGCCTGCGATATCAATCCCCTTGTTGCCCCCCGTCGGAAGCGCCAAAGGTTTCGATCa                                                                                                         >  1:339892/1‑140 (MQ=255)
                                        ccGTAGCCGCGCAGCGCGTTACCAGCATAAACAACGCGGCCATCTGCGGTCGCGATAATTGCCTGTCCTTTGCTGCCTGCGATATCAATCCCCTTGTTGCCCCCCCTCAGAAGCGCCAAAGGTTTCGATCACTTTGCCCt                                                                                                <  2:80835/140‑1 (MQ=255)
                                              ccGCGCAGCGCGTTACCAGCATAAACAACGCGGCCATCTGCGGTCGCGATAATTGCCTGTCCTTTGCTGCCTGCGATATCAATCCCCTTGTTGCCCCCCCTCAGAAGCGCCAAAGGTTTCGATCACTTTGCCCTCAGTCg                                                                                          >  1:735611/1‑140 (MQ=255)
                                                   cAGCGCGTTACCAGCATAAACAACGCGGCCATCTGCGGTCGCGATAATTGCCTGTCCTTTGCTGCCTGCGATATCAATCCCCTTGTTGCCCCCCCTCAGAAGCGCCAAAGGTTTCGATCACTTTGCCCTCAGTCGGCCAg                                                                                     >  1:700874/1‑140 (MQ=255)
                                                      cgcgTTACCAGCATAAACAACGCGGCCATCTGCGGTCGCGATAATTGCCTGTCCTTTGCTGCCTGCGATATCAATCCCCTTGTTGCCCCCCCTCAGAAGCGCCAAAGGTTTCGATCACTTTGCCCTCAGTCGGCCAgcgc                                                                                  >  2:776454/1‑140 (MQ=255)
                                                      cgcgTTACCAGCATAAACAACGCGGCCATCTGCGGTCGCGATAATTGCCTGTCCTTTGCTGCCTGCGATATCAATCCCCTTGTTGCCCCCCCTCAGAAGCGCCAAAGGTTTCGATCACTTTGCCCTCAGTCGGCCAgcgc                                                                                  <  2:601601/140‑1 (MQ=255)
                                                      cgcgTTACCAGCATAAACAACGCGGCCATCTGCGGTCGCGATAATTGCCTGTCCTTTGCTGCCTGCGATATCAATCCCCTTGTTGCCCCCCCTCAGAAGCGCCAAAGGTTTCGATCACTTTGCCCTCAGTCGGCCAgcgc                                                                                  >  2:465857/1‑140 (MQ=255)
                                                               aGCATAAACAACGCGGCCATCTGCGGTCGCGATAATTGCCTGTCCTTTGCTGCCTGCGATATCAATCCCCTTGTTGCCCCCCCTCAGAAGCGCCAAAGGTTTCGATCACTTTGCCCTCAGTCGGCCAGCGCCAGGTGgag                                                                         <  2:653035/140‑1 (MQ=255)
                                                                 cATAAACAACGCGGCCATCTGCGGTCGCGATAATTGCCTGTCCTTTGCTGCCTGCGATATCAATCCCCTTGTTGCCCCCCCTCAGAAGCGCCAAAGGTTTCGATCACTTTGCCCTCAGTCGGCCAGCGCCAGGTGGAGAt                                                                       <  2:735611/140‑1 (MQ=255)
                                                                          cgcgGCCATCTGCGGTCACGATCATTGCCTGTCCTTTGCAGCCTGCGACATCAATCCCCTTGTTGCCCCCCCTCAGAAGCGCCAAAGGTTTCGATCACTTTGCCCTCAGTCGGCCAGCGCCAGGTGGAGATAGGCGTACt                                                              <  2:324485/140‑1 (MQ=255)
                                                                            cgGCCATCTGCGGTCGCGATAACTGTCTGCCCTTTGCCGCCTGCGATATCAATCCCCTTGTTGCCCCCCCTCAGAAGCGCCAAAGGTTTCGATCACTTTGCCCTCAGTCGGCCAGCGCCAGGTGGAGATAGGCGTACTgg                                                            <  1:384986/140‑1 (MQ=255)
                                                                                  tCTGCGGTCGCGATAATTGCCTGTCCTTTGCTGCCTGCGATATCAATCCCCTTGTTGCCCCCCCTCAGAAGCGCCAAAGGTTTCGATCACTTTGCCCTCAGTCGGCCAGCGCCAGGTGGAGATAGGCGTACTGGTTGATg                                                      >  2:723535/1‑140 (MQ=255)
                                                                                                       tGTCCTTTGCTGCCTGCGATATCAATCCCCTTGTTGCCCCCCCTCAGAAGCGCCAAAGGTTTCGATCACTTTGCCCTCAGTCGGCCAGCGCCAGGTGGAGATAGGCGTACTGGTTGATGTACTGCTGACAGTCGGCTCgg                                 >  1:96482/1‑140 (MQ=255)
                                                                                                         tCCTTTGCTGCCTGCGATATCAATCCCCTTGTTGCCCCCCCTCAGAAGCGCCAAAGGTTTCGATCACTTTGCCCTCAGTCGGCCAGCGCCAGGTGGAGATAGTCGTACTGGTTGATGTACTGCTCACAGTCGGCTCGGtt                               >  1:380660/1‑140 (MQ=255)
                                                                                                         tCCTTTGCTGCCTGCGATATCAATCCCCTTGTTGCCCCCCCTCAGAAGCGCCAAAGGTTTCGATCACTTTGCCCTCAGTCGGCCAGCGCCAGGTGGAGATAGGCGTACTGGTTGATGTACTGCTGACAGTCGGCTCGGtt                               >  2:647380/1‑140 (MQ=255)
                                                                                                          ccTTTGCTGTCTGCGATATCAATCCACTTGGTGCCCCCCCTCAGAAGCGCCAAAGGTTTTGAGCACTTTGGCCTCAGTCGGCCTGAGCGGAGTGGAGATAGGCGGCCTGGTGGATGTACTGCTGATATACGGATCGGTtg                              >  2:537889/1‑140 (MQ=255)
                                                                                                               gctgcCTGCGATATCAATCCCCTTGTTGCCCCCCCTCAGAAGCGCCAAAGGTTTCGATCACTTTGCCCTCAGTCGGCCAGCGCCAGGTGGAGATAGGCGTACTGGTTGATGTACTGCTGACAGTCGGCTCGGTTGtgctt                         >  2:130079/1‑140 (MQ=255)
                                                                                                                   cctgcGATATCAATCCCCTTGTTGCCCCCCCTCAGAAGCGCCAAAGGTTTCGATCACTTTGCCCTCAGTCGGCCAGCGCCAGGTGGAGATAGGCGTACTGGTTGATGTACTGCTGACAGTCGGCTCGGTTGTGCTTGCTg                     >  1:65030/1‑140 (MQ=255)
                                                                                                                     tgcGATATCAATCCCCTTGTTGCCCCCCCTCAGAAGCGCCAAAGGTTTCGATCACTTTGCCCTCAGTCGGCCAGCGCCAGGTGGAGATAGGGGTACTGGTCGATCTACTGCTGACAGACGGCCCGGTTGTGCTTGCTGtt                   >  2:337442/1‑140 (MQ=255)
                                                                                                                     tgcGATATCAATCCCCTTGTTGCCCCCCCTCAGAAGCGCCAAAGGTTTCGATCACTTTGCCCTCAGTCGGCCAGCGCCAGGTGGAGATAGGCGTACTGGTTGATGTACTGCTGACAGTCGGCTCGGTTGTGCTTGCTGtt                   >  1:471325/1‑140 (MQ=255)
                                                                                                                      gcGATATCAATCCCCTTGTTGCCCCCCCTCAGAAGCGCCAAAGGTTTCGATCACTTTGCCCTCAGTCGGCCAGCGCCAGGTGGAGATAGGCGTACTGGTTGATGTACTGCTGACAGTCGGCTCGGTTGTGCTTGCTGTTg                  >  2:694390/1‑140 (MQ=255)
                                                                                                                         atatCAATCCCCTTGTTGCCCCCCCTCAGAAGCGCCAAAGGTTTCGATCACTTTGCCCTCAGTCGGCCAGCGCCAGGTGGAGATAGGCGTACTGGTTGATGTACTGCTGACAGTCGGCTCGGTTGTGCTTGCTGTTGGTa               <  1:465857/140‑1 (MQ=255)
                                                                                                                                 ccccTTGTTGCCCCCCCTCAGAAGCGCCAAAGGTTTCGATCACTTTGCCCTCAGTCGGCCAGCGCCAGGTGGAGATAGGCGTACTGGTTGATGTACTGCTGACAGTCGGCTCGGTTGTGCTTGCTGTTGGTACCGTTACa       <  2:96482/140‑1 (MQ=255)
                                                                                                                                   ccTTGTTGCCCCCCCTCAGAAGGGCAAAACGTTTCGATCACTTTGCCCTCAGTCGGCCAGCGCCAGGTGGAGATAGGCGTACTGGTTGATGTACTGCTGACAGTCGGCTCGGTTGTGCTTGCTGTTGGTACCGTTACAgg     <  1:337442/140‑1 (MQ=255)
                                                                                                                                   ccTTGTTGCCCCCCCTCAGAAGCGCCAAAGGTTTCGATCACTTTGCCCTCAGTCGGCCAGCGCCAGGTGGAGATAGGCGTACTGGTTGATGTACTGCTGACAGTCGGCTCGGTTGTGCTTGCTGTTGGTACCGTTACAgg     <  1:723535/140‑1 (MQ=255)
                                                                                                                                    cTTGTTGCCCCCCCTCAGAAGCGCCAAAGGTTTCGATCACTTTGCCCTCAGTCGGCCAGCGCCAGGTGGAGATAGGCGTACTGGTTGATGTGCTGCTGACGGTCGGCTGGGTTGTGCGTGCTGTTGGTACCGTTACAGgc    >  2:224977/1‑140 (MQ=255)
                                                                                                                                     ttgttgcCCCCCCTCAGAAGCGCCAAAGGTTTCGATCACTTTGCCCTCAGTCGGCCAGCGCCAGGTGGAGATAGGCGTACTGGTTGATGTACTGCTGACAGTCGGCTCGGTTGTGCTTGCTGTTGGTa               >  1:842615/1‑128 (MQ=255)
                                                                                                                                     ttgttgcCCCCCCTCAGAAGCGCCAAAGGTTTCGATCACTTTGCCCTCAGTCGGCCAGCGCCAGGTGGAGATAGGCGTACTGGTTGATGTACTGCTGACAGTCGGCTCGGTTGTGCTTGCTGTTGGTa               <  2:842615/128‑1 (MQ=255)
                                                                                                                                      tgttgcCCCCCCTCAGAAGCGCCAAAGGTTTCGATCACTTTGCCCTCAGTCGGCCAGCGCCAGGTGGAGATAGGCGTACTGGTTGATGTACTGCTGACAGTCGGGTCGTTTGTGCTTACTGTCGGTACCGTTACAGgcgc  >  2:324169/1‑140 (MQ=255)
                                                                                                                                      tgttgcCCCCCCTCAGAAGCGCCAAAGGTTTCGATCACTTTGCCCTCAGTCGGCCAGCGCCAGGTGGAGATAGGCGTACTGGTTGATGTACTGCTGACAGTCGGCTCGGTTGTGCTTGCTGTTGGTACCGTTACAGgcgc  >  1:103858/1‑140 (MQ=255)
                                                                                                                                           |                                                                                                                                       
ACTCAGGTAATCATCATTATGTTTGATGATAATCAGATTACCGTAGCCGCGCAGCGCGTTACCAGCATAAACAACGCGGCCATCTGCGGTCGCGATAATTGCCTGTCCTTTGCTGCCTGCGATATCAATCCCCTTGTTGCCCCCCTCAGAAGCGCCAAAGGTTTCGATCACTTTGCCCTCAGTCGGCCAGCGCCAGGTGGAGATAGGCGTACTGGTTGATGTACTGCTGACAGTCGGCTCGGTTGTGCTTGCTGTTGGTACCGTTACAGGCGC  >  NC_000913/2867782‑2868054

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