Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,839,562 A→G G144G (GGT→GGC nlpA ← lipoprotein‑28

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,839,5620AG100.0% 182.6 / NA 52G144G (GGT→GGCnlpAlipoprotein‑28
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (29/23);  total (29/23)

CGCGGGAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCGGTAGGTAATAACCCTTTTCCCTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGATTTTTTTGGAATAACCCGCCATAGGGAAGACAAACGTATTTCCCACCGCTACCAGTTTATAACCATGCG  >  NC_000913/3839425‑3839698
                                                                                                                                         |                                                                                                                                        
cgcgGGAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCGGTAGGTAATAACCCTTTTCCCTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGGcc                                                                                                                                        >  1:1013030/1‑140 (MQ=255)
cgcgGGAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCGGTAGGTAATAACCCTTTTCCCTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGGcc                                                                                                                                        <  2:1044981/140‑1 (MQ=255)
           cgcCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCGGTAGGTAATAACCCTTTTCCCTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGGCCAAGGTTggtcg                                                                                                                             >  2:297399/1‑140 (MQ=255)
           cgcCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCGGTAGGTAATAACCCTTTTCCCTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGGCCAAGGTTggtcg                                                                                                                             <  1:297399/140‑1 (MQ=255)
             cccccTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCGGTAGGTAATAACCCTTTTCCCTCTTTGAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGGCCAAGGTTggtcggg                                                                                                                           <  2:112882/140‑1 (MQ=255)
             cccccTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCGGTAGGTAATAACCCTTTTCCCTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGGCCAAGGTCggtcggg                                                                                                                           <  2:553050/140‑1 (MQ=255)
                ccTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCGGTAGGTAATAACCCTTTTCCCTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGGCCAAGGTTggtcgggtcg                                                                                                                        <  2:2191/140‑1 (MQ=255)
                                   cAAATGACGAGGATTATCGGTTATATCCAGCGCGGTAGGTAATAACCCTTTTCCCTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGGCCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTcg                                                                                                     >  2:1047423/1‑140 (MQ=255)
                                    aaaTGACGAGGATTATCGGTTATATCCAGCGCGGTAGGTAATAACCCTTTTCCCTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGGCCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTcgc                                                                                                    <  2:665123/140‑1 (MQ=255)
                                    aaaTGACGAGGATTATCGGTTATATCCAGCGCGGTAGGTAATAACCCTTTTCCCTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGGCCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTcgc                                                                                                    >  2:63185/1‑140 (MQ=255)
                                     aaTGACGAGGATTATCGGTTATATCCAGCGCGGTAGGTAATAACCCTTTTCCCTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGGCCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTcgcg                                                                                                   <  2:68844/140‑1 (MQ=255)
                                     aaTGACGAGGATTATCGGTTATATCCAGCGCGGTAGGTAATAACCCTTTTCCCTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGGCCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTcgcg                                                                                                   <  1:63185/140‑1 (MQ=255)
                                      aTGACGAGGATTATCGGTTATATCCAGCGCGGTAGGTAATAACCCTTTTCCCTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGGCCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCgcgc                                                                                                  <  2:855721/140‑1 (MQ=255)
                                        gACGAGGATTATCGGTTATATCCAGCGCGGTAGGTAATAACCCTTTTCCCTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGGCCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCct                                                                                                >  2:776353/1‑140 (MQ=255)
                                         aCGAGGATTATCGGTTATATCCAGCGCGGTAGGTAATAACCCTTTTCCCTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGGCCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCctt                                                                                               >  2:453944/1‑140 (MQ=255)
                                                        tataTCCAGCGCGGTAGGTAATAACCCTTTTCCCTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGGCCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCa                                                                                >  1:85839/1‑140 (MQ=255)
                                                                 cgcgGTAGGTAATAACCCTTTTCCCTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGGCCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGa                                                                       >  2:375056/1‑140 (MQ=255)
                                                                 cgcgGTAGGTAATAACCCTTTTCCCTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGGCCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGa                                                                       >  2:154149/1‑140 (MQ=255)
                                                                 cgcgGTAGGTAATAACCCTTTTCCCTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGGCCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGa                                                                       >  1:1036114/1‑140 (MQ=255)
                                                                 cgcgGTAGGTAATAACCCTTTTCCCTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGGCCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGa                                                                       >  2:61546/1‑140 (MQ=255)
                                                                 cgcgGTAGGTAATAACCCTTTTCCCTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGGCCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGa                                                                       >  2:730194/1‑140 (MQ=255)
                                                                 cgcgGTAGGTAATAACCCTTTTCCCTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGGCCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGa                                                                       >  1:207369/1‑140 (MQ=255)
                                                                      taggtaATAACCCTTTTCCCTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGGCCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGAttttt                                                                  >  2:274497/1‑140 (MQ=255)
                                                                      taggtaATAACCCTTTTCCCTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGGCCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGAttttt                                                                  <  2:1022331/140‑1 (MQ=255)
                                                                      taggtaATAACCCTTTTCCCTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGGCCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGAttttt                                                                  >  1:23759/1‑140 (MQ=255)
                                                                      taggtaATAACCCTTTTCCCTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGGCCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGAttttt                                                                  >  2:152337/1‑140 (MQ=255)
                                                                       aggtaATAACCCTTTTCCCTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGGCCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGAtttttt                                                                 <  1:776353/140‑1 (MQ=255)
                                                                          taataaCCCTTTTCCCTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGGCCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCCCGGTTTTGATTTTTTTgg                                                              <  2:85839/140‑1 (MQ=255)
                                                                          taataaCCCTTTTCCCTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGGCCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGATTTTTTTgg                                                              <  2:886979/140‑1 (MQ=255)
                                                                           aataaCCCTTTTCCCTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGGCCAAGGTTGGTCGGGTGGTGGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGATTTTTTTgaa                                                             <  2:803355/140‑3 (MQ=255)
                                                                            ataaCCCTTTTCCCTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGGCCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGATTTTTTTGGaa                                                            <  1:961175/140‑1 (MQ=255)
                                                                                    tttCCCTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGGCCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGATTTTTTTGGAATAACCCGc                                                    >  1:97762/1‑140 (MQ=255)
                                                                                        ccTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGGCCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGATTTTTTTGGAATACCCCGCCATa                                                <  2:376460/140‑1 (MQ=255)
                                                                                          tctTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGGCCAAGGTTGGTCGGGACGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGATTTTTTTGGAATAACCCGCCATAgg                                              >  1:121491/1‑140 (MQ=255)
                                                                                            tttCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGGCCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGATTTTTTTGGAATAACCCACCATAGGGa                                            <  1:1047423/140‑1 (MQ=255)
                                                                                             ttCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGGCCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGATTTTTTTGGAATAACCCGCCATAGGGaa                                           >  1:335036/1‑140 (MQ=255)
                                                                                                  agTAATTGTCTTCTTTTTTTGCACCAGTAAAAGTGCGCGGCCAGGGTTGGTGGGGTCGTTGGGGACCGCCTCGGTCTCGCCTTCTTTTATTTGCGCCGCGGTTTTGTTTTTTTGGGAATAGCCCGCCATAGGGAAGACaa                                      <  2:121491/140‑1 (MQ=255)
                                                                                                  agTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGGCCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGATTTTTTTGGAATAACCCGCCATAGGGAAGACaa                                      <  1:730194/140‑1 (MQ=255)
                                                                                                   gTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGGCCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGATTTTTTTGGAATAACCCGCCATAGGGAAGACaaa                                     <  2:97762/140‑1 (MQ=255)
                                                                                                      aTTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGGCCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGTTTTTGATTTTTTTGGAATAACCCGCCATAGGGAAGACAAACGt                                  <  2:207369/140‑1 (MQ=255)
                                                                                                      aTTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGGCCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGATTTTTTTGGAATAACCCGCCATAGGGAAGACAAACGt                                  >  1:112104/1‑140 (MQ=255)
                                                                                                              ctctTTTTGCAACAGTAAAAGTGCGCGGCCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGATTTTTTTGGAATAACCCGCCATAGGGAAGACAAACGTATTTCCCa                          >  2:831701/1‑140 (MQ=255)
                                                                                                               tctTTTTGCAACAGTAAAAGTGCGCGGCCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGATTTTTTTGGAATAACCCGCCATAGGGAAGACAAACGTATTTCCCAc                         <  1:274497/140‑1 (MQ=255)
                                                                                                                   tttGCAACAGTAAAAGTGCGCGGCCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGATTTTTTTGGAATAACCCGCCATAGGGAAGACAAACGTATTTCCCACCGCt                     >  2:311027/1‑140 (MQ=255)
                                                                                                                     tGCAACAGTAAAAGTGCGCGGCCAAGGTTGGTGGGGTCGTGGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGATTTTTTTGAAATAACCCGCCATAGGGAAGACAAACGTATTTCCCACCGCTAc                   <  1:314378/140‑1 (MQ=255)
                                                                                                                     tGCAACAGTAAAAGTGCGCGGCCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGATTTTTTTGGAATAACCCGCCATAGGGAAGACAAACGTATTTCCCACCGCTAc                   <  1:86342/140‑1 (MQ=255)
                                                                                                                     tGCAACAGTAAAAGTGCGCGGCCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGATTTTTTTGGAATAACCCGCCATAGGGAAGACAAACGTATTTCCCACCGCTAc                   <  1:857764/140‑1 (MQ=255)
                                                                                                                     tGCAACAGAAAAAGTGCGCGGCCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGATTTTTTTGGAATAACCCGCCATAGGGAAGACAAACGTATTTCCCACCGCTAc                   <  2:127053/140‑1 (MQ=255)
                                                                                                                       cAACAGTAAAAGTGCGCGGCCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGATTTTTTTGGAATAACCCGCCATAGGGAAGACAAACGTATTTCCCACCGCTACCa                 <  1:311027/140‑1 (MQ=255)
                                                                                                                           aGTAAAAGTGCGCGGCCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGATTTTTTTGGAATAACCCGCCATAGGGAAGACAAACGTATTTCCCACCGCTAc                   <  1:198033/134‑1 (MQ=255)
                                                                                                                           aGTAAAAGTGCGCGGCCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGATTTTTTTGGAATAACCCGCCATAGGGAAGACAAACGTATTTCCCACCGCTAc                   >  2:198033/1‑134 (MQ=255)
                                                                                                                             tAAAAGTGCGCGGCCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGATTTTTTTGGAATAACCCGCCATAGGGAAGACAAACGTATTTCCCACCGCTAc                   <  2:277971/132‑1 (MQ=255)
                                                                                                                             tAAAAGTGCGCGGCCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGATTTTTTTGGAATAACCCGCCATAGGGAAGACAAACGTATTTCCCACCGCTAc                   >  1:277971/1‑132 (MQ=255)
                                                                                                                               aaaGTGCGCGGCCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGATTTTTTTGGAATAACCCGCCATAGGGAAGACAAACGTATTTCCCACCGCTACCAGTTTATaa         >  1:216740/1‑140 (MQ=255)
                                                                                                                                    gcgcGGCCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGATTTTTTTGGAATAACCCGCCATAGGGAAGACAAACGTATTTCCCACCGCTACCAGTTTATAACCATg    >  2:70083/1‑140 (MQ=255)
                                                                                                                                     cgcgGCCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGATTTTTTTGGAATAACCCGCCATAGGGAAGACAAACGTATTTCCCACCGCTACCAGTTTATAACCATgc   >  1:1020661/1‑140 (MQ=255)
                                                                                                                                     cgcgGCCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGATTTTTTTGGAATAACCCGCCATAGGGAAGACAAACGTATTTCCCACCGCTACCAGTTTATAACCATgc   >  1:539306/1‑140 (MQ=255)
                                                                                                                                     cgcgGCCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGATTTTTTTGGAATAACCCGCCATAGGGAAGACAAACGTATTTCCCACCGCTACCAGTTTATAACCATgc   <  2:539306/140‑1 (MQ=255)
                                                                                                                                     cgcgGCCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGATTTTTTTGGAATAACCCGCCATAGGGAAGACAAACGTATTTCCCACCGCTACCAGTTTATAACCATgc   >  2:477886/1‑140 (MQ=255)
                                                                                                                                      gcgGCCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGATTTTTTTGGAATAACCCGCCATAGGGAAGACAAACGTATTTCCCACCGCTACCAGTTTATAACCATgcg  <  1:477886/140‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                        
CGCGGGAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCGGTAGGTAATAACCCTTTTCCCTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGATTTTTTTGGAATAACCCGCCATAGGGAAGACAAACGTATTTCCCACCGCTACCAGTTTATAACCATGCG  >  NC_000913/3839425‑3839698

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: