Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,062,3550TG24.0% 42.8 / 10.8 25intergenic (+20/+36)asnV/ldtAtRNA‑Asn/L,D‑transpeptidase linking Lpp to murein
Reads supporting (aligned to +/- strand):  ref base T (11/8);  new base G (2/4);  total (13/12)
Fisher's exact test for biased strand distribution p-value = 3.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TCACCATGCACTGCCATTAATATGCGCCCCGTTCACACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACGCAGCATTTGCGCGATTTACATCAACTTCTTGCCCGTTGATAAACGC  >  NC_000913/2062219‑2062498
                                                                                                                                        |                                                                                                                                               
tCACCATGCACTGCCATTAATATGCGCCCCGTTCACACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCttt                                                                                                                                     <  3:95373/149‑1 (MQ=255)
  aCCATGCACTGCCATTAATATGCGCCCCGTTCACACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCtttt                                                                                                                                    <  4:114898/148‑1 (MQ=255)
   ccATGCACTGCCATTAATATGCGCCCCGTTCACACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTAGAAAGCAGGCtttttt                                                                                                                                  <  1:73128/149‑2 (MQ=255)
   ccATGCACTGCAATTAATATGCGCCCCGTTCACACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCTGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGATCCAAATTCAAAAAAGCCTGCTAGAAAGCAGGCtttttt                                                                                                                                  <  1:70678/149‑2 (MQ=17)
    cATGCACTGCCATTACTATGCCCCCCTTTCACAAGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGGCTGTTAATCCGGATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGc                                                                                                                                 >  4:106902/1‑149 (MQ=18)
    cATGCACTGCCATTAATATGCGCCCCGTTCACACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGc                                                                                                                                 >  3:17974/1‑149 (MQ=255)
        cACTGCCATTAATATGCGCCCCGTTCACACGAGTCCTCTGTAGTTCAGTCGGTAGAAGGTCGGACTGGTAATCCGTATGTCACTGGGTCGAGTCCAGTCAGAGCAGCCAAATGCAAAAAAGCCTGCTATCTAGCAGGCTTTTTGCttt                                                                                                                              <  3:106902/148‑1 (MQ=18)
                                      gATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCa                                                                                                                                           >  4:72778/1‑105 (MQ=255)
                                                tAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGaa                                                                                     <  4:17974/149‑1 (MQ=255)
                                                  gTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAAc                                                                                    >  1:51490/1‑148 (MQ=255)
                                                     cAGTCGGTAGAACGGCTGACTGTTAATCCTTATGTCACTGGTTCGAGTACAGTCAGAGTAGCCAAATTCAAAAAAGCCTGATTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGATCTTAAAACATCTGTCTTTAACAAGAACTAAt                                                                                >  8:78220/1‑149 (MQ=255)
                                                              gAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAgg                                                                       >  6:14581/1‑149 (MQ=255)
                                                                                         aaggtttCGAGTCCAGTCAGATGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACgcagc                                            <  7:78220/144‑1 (MQ=255)
                                                                                         aCTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTACTTACCAACGCTCTTAACACATCTGTCTAGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACgcag                                             >  4:74508/1‑148 (MQ=255)
                                                                                         aCTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACgcag                                             >  6:157583/1‑148 (MQ=255)
                                                                                                       aGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACGCAGCATTTGCGCGATTTa                              >  3:125078/1‑149 (MQ=255)
                                                                                                       aGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACGCAGCATTTGCGCGATTTa                              >  4:96351/1‑149 (MQ=255)
                                                                                                           agagGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACGCAGCATTTGCGCGATTTACAt                           <  5:157583/148‑1 (MQ=255)
                                                                                                              ggagCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACGCAGCATTTGCGCGATTTACATCAAc                       <  5:14581/149‑1 (MQ=255)
                                                                                                                        agcaaaaaGCCTGCTAGAAAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACGCAGCATTTGCGCGATTTACa                            <  4:105406/131‑1 (MQ=255)
                                                                                                                        agcaaaaaGCCTGCTAGAAAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACGCAGCATTTGCGCGATTTACa                            >  3:105406/4‑134 (MQ=255)
                                                                                                                          aaaaaaGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACGCAGCATTTGCGCGATTTACATCAACTTCTTGCCCGtt           <  4:125078/149‑1 (MQ=255)
                                                                                                                            aaaaGCCTGCTAGAAAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACGCAGCATTTGCGCGATTTAc                             <  7:72533/129‑1 (MQ=255)
                                                                                                                            aaaaGCCTGCTAGAAAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACGCAGCATTTGCGCGATTTAc                             >  8:72533/1‑129 (MQ=255)
                                                                                                                                   tGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACGCAGCATTTGCGCGATTTACATCAACTTCTTGCCCGTTGATAAACGc  >  1:74037/1‑149 (MQ=255)
                                                                                                                                        |                                                                                                                                               
TCACCATGCACTGCCATTAATATGCGCCCCGTTCACACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACGCAGCATTTGCGCGATTTACATCAACTTCTTGCCCGTTGATAAACGC  >  NC_000913/2062219‑2062498

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