Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,329,743 G→T D261Y (GAT→TAT)  dacB → peptidoglycan DD‑endopeptidase DacB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,329,7430GT97.1% 107.0 / ‑3.6 35D261Y (GAT→TAT) dacBpeptidoglycan DD‑endopeptidase DacB
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base T (17/17);  total (17/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.43e-01

CTACAGTGCCCCAAAGCCTGGTGATATGGCTTTTATACGCGTGGCATCTTATTACCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAAGATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGATAGGCCATGCGCGCTTCAATGTGCCTGGAACATGGCGGGCCGGGTCGGACGCCGTGCGTCA  >  NC_000913/3329496‑3329975
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                             
ctACAGTGCCCCAAAGCCTGGTGATATGGCTTTTATACGCGTGGAATCTTATTACCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATg                                                                                                                                                                                                                                            >  1:53611/1‑251 (MQ=255)
  aCAGTGCCCCAAAGCCTGGTGATATGGCTTTTATACGCGTGGCATCTTATTACCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATg                                                                                                                                                                                                                                            >  2:14494/1‑249 (MQ=255)
             aaGCCTGGTGATATGGCTTTTATACGCGTGTCATCTTATTACCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGGTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGGGTTAAAACAGGc                                                                                                                                                                                                                               >  2:129523/1‑251 (MQ=255)
              aGCCGGGTGATATGCCTTTTATACGCGTGGCACCTTTGTCCCCCGTTACGTTGTCCAGCCAGTTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACGGGATGTGGTGCCAGGCGACTTGAACCGCTTTACGCTGACGGGATGCCTGCCAAAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCg                                                                                                                                                                                                                              <  2:69191/251‑1 (MQ=255)
               gcctggggaactGGCTTTTATACGCGGGGCATCTTAGTCCACCGTAACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACGGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCgg                                                                                                                                                                                                                             <  2:108911/240‑1 (MQ=255)
                            cgaactataCGCGTGGCATCTTATTACCCCGTTACGAGGTTCAGCCAGGTACGCACCCTCCCCCGGGGGTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGCCCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCGTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAg                                                                                                                                                                                                                <  2:302121/246‑1 (MQ=255)
                               tttATCCGCGTTGCATCTTATTACCCCGTTACGATGTTCAGCCAGGTACGCACCCCCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGc                                                                                                                                                                                                               <  2:217672/249‑1 (MQ=255)
                                  ataCGCGTGGCATCTTATTACCCCGTTACGATGTTCAGCCAGGTCCGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATTTGTTGCCAGGTGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGaa                                                                                                                                                                                                          <  1:459403/250‑1 (MQ=255)
                                  ataCGCGTGGCATCTTATTACCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGGATCACCTGGAg                                                                                                                                                                                                                >  1:302121/1‑245 (MQ=255)
                                                ttattcCCCCGTTACGATGTGCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACGGGATGTGGTGCCAGGCGACCTGACCCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAg                                                                                                                                                                                            <  2:286159/250‑1 (MQ=255)
                                                             ccgggtttCAGCCAGGTACGCCCCCTCCCCCGTGGTTCGGCCGAAGCGCAATACTGCGAACGGGAGGGGGGGCCAGGCGACCTGACCCTCTTGACGCTGACGGGAGGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCGAGCTATGCCGGTGCAAGTCTGAAATATGAGTTAAAACAGGCGGGGATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACg                                                                                                                                                                               <  2:471770/244‑1 (MQ=255)
                                                               gATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAACCAGGCGGGTAACCCCCGGGGCGGGACACAGCTGCCCCAGAATCAGGTTAACGAAc                                                                                                                                                                            >  1:191380/1‑251 (MQ=255)
                                                               gATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGCACACTGCTGCGCCAGACTCCGGTTAACGAAc                                                                                                                                                                            >  1:29548/1‑251 (MQ=255)
                                                                             tACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCCCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTCTCACCTGGAGCGGAACACTGCCGCGCCAGACGCAGGGTAACGAAACGGGAACgggaggt                                                                                                                                                              >  1:232943/1‑246 (MQ=255)
                                                                                         cccGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTCTCACCTGGAGCGGTAAACTCCTGCGCCAGAACCCGGTTACCGAACCTGGGACAGTCGTTGACAGGAAACa                                                                                                                                                   >  1:23849/1‑250 (MQ=255)
                                                                                                                          tGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGGACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGGTGCCAGTAAACAGGCGGCCCCCCTGCACGATCTGCTTAAGAttttg                                                                                                                 >  1:168432/1‑248 (MQ=255)
                                                                                                                                      cAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCg                                                                                                   <  1:225061/251‑1 (MQ=255)
                                                                                                                                             ccTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGGGCGGGACACTGCTGCGCCAGACTCAGGTTAACGAAACTGGAACGGTAGTTGCCAGTAAAAAGTCGGCCCCGCTGCACGGTCTGCTTAAGATTATGCTGAAAAAGTCggcaaac                                                                                             >  1:102390/1‑245 (MQ=255)
                                                                                                                                                gaccctcttcacgtggtcgtgATGCCTGCCACAACGTTTTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGGCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATg                                                                                          <  2:149860/230‑1 (MQ=255)
                                                                                                                                                gAACCGCTTTACGCTGACTGGAGGGCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCGTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGTAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGa                                                                                         <  1:189773/251‑1 (MQ=255)
                                                                                                                                                 aaCCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGGAATTTCCAGTAAACAGTTGGGCCCGCTGCACGATATGCTTAAGTTTAAGCTGACAAAGTCGGccaacctta                                                                                         >  1:133215/1‑246 (MQ=255)
                                                                                                                                                  aCCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATc                                                                                       <  2:168432/251‑1 (MQ=255)
                                                                                                                                                               gACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTTAAAAAGTCGGACAACATTATCGCCGATACGGttt                                                                          <  1:395536/251‑1 (MQ=255)
                                                                                                                                                                aCGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGGTTCCCCTGGAGCGGGACACTGCTGCGGCAGAACCAGGTTAACGAACCTGGGACGGTAGTTTCCAGTAAAAAGTCGGCCCCGCTGCCCGATCTGGTTAAGATTATGGTGAAAAAGTCGGACAAAATGATCGCCGATCAGGtttt                                                                         >  2:460063/1‑251 (MQ=255)
                                                                                                                                                                                                cccGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTGGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGATAGGCCATGCGCGCTTCAAtgtg                                        <  1:428996/250‑1 (MQ=255)
                                                                                                                                                                                                   gTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGGGCGGAACAATGCCGCGGCAGCATAAGGTTAACGAACCTGGAACGGTAGTTGGCCGTAAACCGTCCGCCCCGCTGCACGATCTGCTTAAGGTTAAGTTGGTAAAAGCGGACAAAAAGATCCGCGATACGGTTTTCCGCATGATAGGCCCTGGGCGCATCAAtggtcctg                                    >  2:458429/1‑245 (MQ=255)
                                                                                                                                                                                                     tGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGAt                                                                <  2:186335/223‑1 (MQ=255)
                                                                                                                                                                                                           ttGCTGTGCAGGAGGGAGCCAGCTATGACGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATTATCTGGAGCGTAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCTGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGGTGGACAACATGATCGCCGATACGGTTTTCCGCATGATAGGCCATGCGCGCTTCAATGTGCCTGGAACATg                             <  1:302961/250‑1 (MQ=255)
                                                                                                                                                                                                               cGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCCCCTGGAGGGGAACACAGCCGGGCCAGACTCAGGCTGAAGAAACTGGAACGGTAGTTGCCAGTAAAAAGACGGCCCCCCTCCACAAACTACTTAAGATTATGATGACAAAGCCGAACAAAATAATCGCCCTTAAGGGTTTTCGCATGATAGGCCCATGGCGCTTTAATATGCCTGGCACATGGcacg                        >  2:483718/1‑248 (MQ=255)
                                                                                                                                                                                                               cGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGGAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGGTTTCCGCATGATAGGCCATGCGCGCTTCAATGTGCCTGGAACATGGCggg                        >  1:142964/1‑251 (MQ=255)
                                                                                                                                                                                                               cGTGCAGGATGGAGCCAGCTAGGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGATAGGCCATGCGCGCTTCAATGTGCCTGGAACATGGCggg                        <  1:460063/251‑1 (MQ=255)
                                                                                                                                                                                                                   cAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTTGAGCGGGAACCTGCTGCGCCAGACTCAGGTTTACGAACCTGGGACGGTGGTTTCCAGTAAACAGTCGGCCCCGCTGCCCGATCTGCTTAAGATTTTGCTTGAAAAAGCCGAACAAATAATCGCCCTAACGGTTTTCCGCATGATTGGGCCTTCGCGCTTTCAAGTGCCTGGAACATGGCggggccg                    >  2:403576/1‑247 (MQ=255)
                                                                                                                                                                                                                    aggagggtgCCAGCTATGCCGGGGCAATTCTGAAATAGGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACGGCTGCGCCAGACTCAGGTTAGCGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCGGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCTGCATGATAGGCCATGCGCGCTGCAATGTGCCTGGAACATGGCGGGCCggg                   <  2:451994/243‑1 (MQ=255)
                                                                                                                                                                                                                           aGCCAGCTATGCCGGTGCAATTCTGAAAGATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGATAGGc                                                            <  1:181945/205‑1 (MQ=255)
                                                                                                                                                                                                                                       cGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAAAGAACCTGGAACGGGAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTTTGCTGGAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGATAAGGCATGCGCGCCTCAAAGTGCCTGGAACATGGAGGGCCGGGGCGGGCGcccggcggcc  >  1:142072/1‑241 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                             
CTACAGTGCCCCAAAGCCTGGTGATATGGCTTTTATACGCGTGGCATCTTATTACCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAAGATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGATAGGCCATGCGCGCTTCAATGTGCCTGGAACATGGCGGGCCGGGTCGGACGCCGTGCGTCA  >  NC_000913/3329496‑3329975

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