Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 668,204 T→C N5S (AAC→AGC)  mrdA ← peptidoglycan DD‑transpeptidase MrdA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913668,2040TC98.3% 209.7 / ‑4.6 60N5S (AAC→AGC) mrdApeptidoglycan DD‑transpeptidase MrdA
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (23/36);  total (23/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

CGGTAAAGCGAACAATTTGCAGATTATACAGGTTGGCGATAAGCACGCCGGTCAGCAGCAAAATCCCCAAAAAGGCGACCAGCGCCCGGCGCACAAACAGCGCGGACTCTGCCGTATAGTCGCGAAAAGAGTTCTGTAGTTTCATCCGCTGCGTTTTCTACTCAAAGCTCCCTTATCACTCACGGTGATAAGGATGGTTGGTGGTGATGCTCCACGCCCGGTACAGACTCTCTGCGACCAGCACGCGAACCAGCGGATGGG  >  NC_000913/668073‑668333
                                                                                                                                   |                                                                                                                                 
cGGTAAAGCGAACAATTTGCAGATTATACAGGTTGGCGATAAGCACGCCGGTCAGCAGCAAAATCCCCAAAAAGGCGACCAGCGCCCGGCGCACAAACAGCGCGGACTCTGCCGTATAGTCGCGAAAAGAGCTCTGTAGt                                                                                                                           >  1:7649/1‑140 (MQ=255)
     aaGCGAACAATTTGCAGATTATACAGGTTGGCGATAAGCACGCCGGTCAGCAGCAAAATCCCCAAAAAGGCGACCAGCGCCCGGCGCACAAACAGCGCGGACTCTGCCGTATAGTCGCGAAAAGAGCTCTGTAGTTTCAt                                                                                                                      >  2:750160/1‑140 (MQ=255)
          aaCAATTTGCAGATTATACAGGTTGGCGATAAGCACGCCGGTCAGCAGCAAAATCCCCAAAAAGGCGACCAGCGCCCGGCGCACAAACAGCGCGGACTCTGCCGTATAGTCGCGAAAAGAGCTCTGTAGTTTCATCCGCt                                                                                                                 <  2:537654/140‑1 (MQ=255)
            cAATTTGCAGATTATACAGGTTGGCGATAAGCACGCCGGTCAGCAGCAAAATCCCCAAAAAGGCGACCAGCGCCCGGCGCACAAACAGCGCGGACTCTGCCGTATAGTCGCGAAAAGAGCTCTGTAGTTTCATCCGCTGc                                                                                                               <  2:465554/140‑1 (MQ=255)
                 tGCAGATTATACAGGTTGGCGATAAGCACGCCGGTCAGCAGCAAAATCCCCAAAAAGGCGACCAGCGCCCGGCGCACAAACAGCGCGGACTCTGCCGTATAGTCGCGAAAAGAGCTCTGTAGTTTCATCCGCTGCGtttt                                                                                                          >  1:378632/1‑140 (MQ=255)
                     gATTATACAGGTTGGCGATAAGCACGCCGGTCAGCAGCAAAATCCCCAAAAAGGCGACCAGCGCCCGGCGCACAAACAGCGCGGACTCTGCCGTATAGTCGCGAAAAGAGCTCTGTAGTTTCATCCGCTGCGTTTTCTAc                                                                                                      <  2:319855/140‑1 (MQ=255)
                       ttATACAGGTTGGCGATAAGCACGCCGGTCAGCAGCAAAATCCCCAAAAAGGCGACCAGCGCCCGGCGCACAAACAGCGCGGACTCTGCCGTATAGTCGCGAAAAGAGCTCTGTAGTTTCATCCGCTGCGTTTTCTACTc                                                                                                    >  1:248870/1‑140 (MQ=255)
                        tataCAGGTTGGCGATAAGCACGCCGGTCAGCAGCAAAATCCCCAAAAAGGCGACCAGCGCCCGGCGCACAAACAGCGCGGACTCTGCCGTATAGTCGCGAAAAGAGCTCTGTAGTTTCATCCGCTGCGTTTTCTACTCa                                                                                                   <  1:939469/140‑1 (MQ=255)
                        tataCAGGTTGGCGATAAGCACGCCGGTCAGCAGCAAAATCCCCAAAAAGGCGACCAGCGCCCGGCGCACAAACAGCGCGGACTCTGCCGTATAGTCGCGAAAAGAGCTCTGTAGTTTCATCCGCTGCGTTTTCTACTCa                                                                                                   <  1:574047/140‑1 (MQ=255)
                            cAGGTTGGCGATAAGCACGCCGGTCAGCAGCAAAATCCCCAAAAAGGCGACCAGCGCCCGGCGCACAAACAGCGCGGACTCTGCCGTATAGTCGCGAAAAGAGCTCTGTAGTTTCATCCGCTGCGTTTTCTACTCAAAGc                                                                                               <  1:691555/140‑1 (MQ=255)
                             aGGTTGGCGATAAGCACGCCGGTCAGCAGCAAAATCCCCAAAAAGGCGACCAGCGCCCGGCGCACAAACAGCGCGGACTCTGCCGTATAGTCGCGAAAAGAGCTCTGTAGTTTCATCCGCTGCGTTTTCTACTCAAAGCt                                                                                              <  2:200560/140‑1 (MQ=255)
                               gTTGGCGATAAGCACGCCGGTCAGCAGCAAAATCCCCAAAAAGGCGACCAGCGCCCGGCGCACAAACAGCGCGGACTCTGCCGTATAGTCGCGAAAAGAGCTCTGTAGTTTCATCCGCTGCGTTTTCTACTCAAAGCTcc                                                                                            <  1:186832/140‑1 (MQ=255)
                               gTTGGCGATAAGCACGCCGGTCAGCAGCAAAATCCCCAAAAAGGCGACCAGCGCCCGGCGCACAAACAGCGCGGACTCTGCCGTATAGTCGCGAAAAGAGCTCTGTAGTTTCATCCGCTGCGTTTTCTACTCAAAGCTcc                                                                                            <  2:378632/140‑1 (MQ=255)
                                ttGGCGATAAGCACGCCGGTCAGCAGCAAAATCCCCAAAAAGGCGACCAGCGCCCGGCGCACAAACAGCGCGGACTCTGCCGTATAGTCGCGAAAAGAGCTCTGTAGTTTCATCCGCTGCGTTTTCTACTCAAAGCTccc                                                                                           >  2:547232/1‑140 (MQ=255)
                                 tGGCGATAAGCACGCCGGTCAGCAGCAAAATTCCCAAAAAGGCGACCAGCGCCCGGCGCACAAACAGCGCGGACTCTGCCGTATAGTCGCGAAAAGAGCTCTGTAGTTTCATCCGCTGCGTTTTCTACTCAAAGc                                                                                               >  2:147006/1‑135 (MQ=255)
                                 tGGCGATAAGCACGCCGGTCAGCAGCAAAATTCCCAAAAAGGCGACCAGCGCCCGGCGCACAAACAGCGCGGACTCTGCCGTATAGTCGCGAAAAGAGCTCTGTAGTTTCATCCGCTGCGTTTTCTACTCAAAGc                                                                                               <  1:147006/135‑1 (MQ=255)
                                  ggCGATAAGCACGCCGGTCAGCAGCAAAATCCCCAAAAAGGCGACCAGCGCCCGGCGCACAAACAGCGCGGACTCTGCCGTATAGTCGCGAAAAGAGCTCTGTAGTTTCATCCGCTGCGTTTTCTACTCAAAGCTCCCtt                                                                                         >  2:502782/1‑140 (MQ=255)
                                  ggCGATAAGCACGCCGGTCAGCAGCAAAATCCCCAAAAAGGCGACCAGCGCCCGGCGCACAAACAGCGCGGACTCTGCCGTATAGTCGCGAAAAGAGCTCTGTAGTTTCATCCGCTGCGTTTTCTACTCAAAGCTCCCtt                                                                                         >  2:467410/1‑140 (MQ=255)
                                   gCGATAAGCACGCCGGTCAGCAGCAAAATCCCCAAAAAGGCGACCAGCGCCCGGCGCACAAACAGCGCGGACTCTGCCGTATAGTCGCGAAAAGAGCTCTGTAGTTTCATCCGCTGCGTTTTCTACTCAAAGCTCCCTTa                                                                                        >  2:930621/1‑140 (MQ=255)
                                    cGATAAGCACGCCGGTCAGCAGCAAAATCCCCAAAAAGGCGACCAGCGCCCGGCGCACAAACAGCGCGGACTCTGCCGTATAGTCGCGAAAAGAGCTCTGTAGTTTCATCCGCTGCGTTTTCTACTCa                                                                                                   <  2:36499/128‑1 (MQ=255)
                                    cGATAAGCACGCCGGTCAGCAGCAAAATCCCCAAAAAGGCGACCAGCGCCCGGCGCACAAACAGCGCGGACTCTGCCGTATAGTCGCGAAAAGAGCTCTGTAGTTTCATCCGCTGCGTTTTCTACTCa                                                                                                   >  1:36499/1‑128 (MQ=255)
                                    cGATAAGCACGCCGGTCAGCAGCAAAATCCCCAAAAAGGCGACCAGCGCCCGGCGCACAAACAGCGCGGACTCTGCCGTATAGTCGCGAAAAGAGCTCTGTAGTTTCATCCGCTGCGTTTTCTACTCAAAGCTCCCTTAt                                                                                       >  1:96542/1‑140 (MQ=255)
                                    cGATAAGCACGCCGGTCAGCAGCAAAATCCCCAAAAAGGCGACCAGCGCCCGGCGCACAAACAGCGCGGACTCTGCCGTATAGTCGCGAAAAGAGCTCTGTAGTTTCATCCGCTGCGTTTTCTACTCAAAGCTCCCTTAt                                                                                       >  1:1031287/1‑140 (MQ=255)
                                    cGATAAGCACGCCGGTCAGCAGCAAAATCCCCAAAAAGGCGACCAGCGCCCGGCGCACAAACAGCGCGGACTCTGCCGTATAGTCGCGAAAAGAGCTCTGTAGTTTCATCCGCTGCGTTTTCTACTCAAAGCTCCCTTAt                                                                                       >  2:461585/1‑140 (MQ=255)
                                    cGATAAGCACGCCGGTCAGCAGCAAAATCCCCAAAAAGGCGACCAGCGCCCGGCGCACAAACAGCGCGGACTCTGCCGTATAGTCGCGAAAAGAGCTCTGTAGTTTCATCCGCTGCGTTTTCTACTCAAAGCTCCCTTAt                                                                                       >  1:47181/1‑140 (MQ=255)
                                        aaGCACGCCGGTCAGCAGCAAAATCCCCAAAAAGGCGACCAGCGCCCGGCGCACAAACAGCGCGGACTCTGCCGTATAGTCGCGAAAAGAGCTCTGTAGTTTCATCCGCTGCGTTTTCTACTCAAAGCTCCCTTAtcact                                                                                   <  1:502782/140‑1 (MQ=255)
                                          gCACGCCGGTCAGCAGCAAAATCCCCAAAAAGGCGACCAGCGCCCGGCGCACAAACAGCGCGGACTCTGCCGTATAGTCGCGAAAAGAGCTCTGTAGTTTCATCCGCTGCGTTTTCTACTCAAAGCTCCCTTAtcactca                                                                                 <  1:50212/140‑1 (MQ=255)
                                           cACGCCGGTCAGCAGCAAAATCCCCAAAAAGGCGACCAGCGCCCGGCGCACAAACAGCGCGGACTCTGCCGTATAGTCGCGAAAAGAGCTCTGTAGTTTCATCCGCTGCGTTTTCTACTCAAAGCTCCCTTAtcactcac                                                                                <  2:995052/140‑1 (MQ=255)
                                           cACGCCGGTCAGCAGCAAAATCCCCAAAAAGGCGACCAGCGCCCGGCGCACAAACAGCGCGGACTCTGCCGTATAGTCGCGAAAAGAGCTCTGTAGTTTCATCCGCTGCGTTTTCTACTCAAAGCTCCCTTAtcactcac                                                                                <  2:196767/140‑1 (MQ=255)
                                            aCGCCGGTCAGCAGCAAAATCCCCAAAAAGGCGACCAGCGCCCGGCGCACAAACAGCGCGGACTCTGCCGTATAGTCGCGAAAAGAGCTCTGTAGTTTCATCCGCTGCGTTTTCTACTCAAAGCTCCCTTATCACTCACg                                                                               <  1:785694/140‑1 (MQ=255)
                                            aCGCCGGTCAGCAGCAAAATCCCCAAAAAGGCGACCAGCGCCCGGCGCACAAACAGCGCGGACTCTGCCGTATAGTCGCGAAAAGAGCTCTGTAGTTTCATCCGCTGCGTTTTCTACTCAAAGCTCCCTTATCACTCACg                                                                               <  1:964411/140‑1 (MQ=255)
                                                  gTCAGCAGCAAAATCCCCAAAAAGGCGACCAGCGCCCGGCGCACAAACAGCGCGGACTCTGCCGTATAGTCGCGAAAAGAGCTCTGTAGTTTCATCCGCTGCGTTTTCTACTCAAAGCTCCCTTATCACTCACGGTGATa                                                                         >  1:1068525/1‑140 (MQ=255)
                                                  gTCAGCAGCAAAATCCCCAAAAAGGCGACCAGCGCCCGGCGCACAAACAGCGCGGACTCTGCCGTATAGTCGCGAAAAGAGCTCTGTAGTTTCATCCGCTGCGTTTTCTACTCAAAGCTCCCTTATCACTCACGGTGATa                                                                         <  2:1068525/140‑1 (MQ=255)
                                                   tCAGCAGCAAAATCCCCAAAAAGGCGACCAGCGCCCGGCGCACAAACAGCGCGGACTCTGCCGTATAGTCGCGAAAAGAGCTCTGTAGTTTCATCCGCTGCGTTTTCTACTCAAAGCTCCCTTAtcact                                                                                   >  1:376536/1‑129 (MQ=255)
                                                   tCAGCAGCAAAATCCCCAAAAAGGCGACCAGCGCCCGGCGCACAAACAGCGCGGACTCTGCCGTATAGTCGCGAAAAGAGCTCTGTAGTTTCATCCGCTGCGTTTTCTACTCAAAGCTCCCTTAtcact                                                                                   <  2:376536/129‑1 (MQ=255)
                                                        agcaAAATCCCCAAAAAGGCGACCAGCGCCCGGCGCACAAACAGCGCGGACTCTGCCGTATAGTCGCGAAAAGAGCTCTGTAGTTTCATCCGCTGCGTTTTCTACTCAAAGCTCCCTTATCACTCACGGTGATAAGGAtg                                                                   <  1:364432/140‑1 (MQ=255)
                                                         gcaAAATCCCCAAAAAGGCGACCAGCGCCCGGCGCACAAACAGCGCGGACTCTGCCGTATAGTCGCGAAAAGAGCTCTGTAGTTTCATCCGCTGCGTTTTCTACTCAAAGCTCCCTTATCACTCACGGTGATAAGGAtgg                                                                  <  1:750160/140‑1 (MQ=255)
                                                          caAAATCCCCAAAAAGGCGACCAGCGCCCGGCGCACAAACAGCGCGGACTCTGCCGTATAGTCGCGAAAAGAGCTCTGTAGTTTCATCCGCTGCGTTTTCTACTCAAAGCTCCCTTATCACTCACGGTGATAAGGAtggt                                                                 >  1:115588/1‑140 (MQ=255)
                                                          caAAATCCCCAAAAAGGCGACCAGCGCCCGGCGCACAAACAGCGCGGACTCTGCCGTATAGTCGCGAAAAGAGCTCTGTAGTTTCATCCGCTGCGTTTTCTACTCAAAGCTCCCTTATCACTCACGGTGATAAGGAtggt                                                                 <  2:1010103/140‑1 (MQ=255)
                                                             aaTCCCCAAAAAGGCGACCAGCGCCCGGCGCACAAACAGCGCGGACTCTGCCGTATAGTCGCGAAAAGAGCTCTGTAGTTTCATCCGCTGCGTTTTCTACTCAAAGCTCCCTTATCACTCACGGTGATAAGGAtggttgg                                                              <  2:734768/140‑1 (MQ=255)
                                                                   cAAAAAGGCGACCAGCGCCCGGCGCACAAACAGCGCGGACTCTGCCGTATAGTCGCGAAAAGAGCTCTGTAGTTTCATCCGCTGCGTTTTCTACTCAAAGCTCCCTTATCACTCACGGTGATAAGGATGGTTGGTGGTGa                                                        <  2:47181/140‑1 (MQ=255)
                                                                         ggCGACCAGCGCCCGGCGCACAAACAGCGCGGACTCTGCCGTATAGTCGCGAAAAGAGCTCTGTAGTTTCATCCGCTGCGTTTTCTACTCAAAGCTCCCTTATCACTCACGGTGATAAGGATGGTTGGTGGTGATGCTcc                                                  <  2:115588/140‑1 (MQ=255)
                                                                         ggCGACCAGCGCCCGGCGCACAAACAGCGCGGACTCTGCCGTATAGTCGCGAAAAGAGCTCTGTAGTTTCATCCGCTGCGTTTTCTACTCAAAGCTCCCTTATCACTCACGGTGATAAGGATGGTTGGTGGTGATGCTcc                                                  <  2:112270/140‑1 (MQ=255)
                                                                              ccAGCGCCCGGCGCACAAACAGCGCGGACTCTGCCGTATAGTCGCGAAAAGAGCTCTGTAGTTTCATCCGCTGCGTTTTCTACTCAAAGCTCCCTTATCACTCACGGTGATAAGGATGGTTGGTGGTGATGCTCCACGcc                                             <  2:764008/140‑1 (MQ=255)
                                                                              ccAGCGCCCGGCGCACAAACAGCGCGGACTCTGCCGTATAGTCGCGAAAAGAGCTCTGTAGTTTCATCCGCTGCGTTTTCTACTCAAAGCTCCCTTATCACTCACGGTGATAAGGATGGTTGGTGGTGATGCTCCACGcc                                             <  2:137830/140‑1 (MQ=255)
                                                                              ccAGCGCCCGGCGCACAAACAGCGCGGACTCTGCCGTATAGTCGCGAAAAGAGCTCTGTAGTTTCATCCGCTGCGTTTTCTACTCAAAGCTCCCTTATCACTCACGGTGATAAGGATGGTTGGTGGTGATGCTCCACGcc                                             <  2:1031287/140‑1 (MQ=255)
                                                                               cAGCGCCCGGCGCACAAACAGCGCGGACTCTGCCGTATAGTCGCGAAAAGAGTTCTGTAGTTTCATCCGCTGCGTTTTCTACTCAAAGCTCCCTTATCACTCACGGTAATAAGGATGGTTGGTGGTGATGCTCCACGccc                                            <  2:289130/140‑1 (MQ=255)
                                                                               cAGCGCCCGGCGCACAAACAGCGCGGACTCTGCCGTATAGTCGCGAAAAGAGCTCTGTAGTTTCATCCGCTGCGTTTTCTACTCAAAGCTCCCTTATCACTCACGGTGATAAGGATGGTTGGTGGTGATGCTCCACGccc                                            <  2:96542/140‑1 (MQ=255)
                                                                               cAGCGCCCGGCGCACAAACAGCGCGGACTCTGCCGTATAGTCGCGAAAAGAGCTCTGTAGTTTCATCCGCTGCGTTTTCTACTCAAAGCTCCCTTATCACTCACGGTGATAAGGATGGTTGGTGGTGATGCTCCACGccc                                            <  2:453747/140‑1 (MQ=255)
                                                                                   gcCCGGCGCACAAACAGCGCGGACTCTGCCGTATAGTCGCGAAAAGAGCTCTGTAGTTTCATCCGCTGCGTTTTCTACTCAAAGCTCCCTTATCACTCACGGTGATAAGGATGGTTGGTGGTGATGCTCCACGCCCGGTa                                        <  1:467410/140‑1 (MQ=255)
                                                                                        gcgcACAAACAGCGCGGACTCTGCCGTATAGTCGCGAAAAGAGCTCTGTAGTTTCATCCGCTGCGTTTTCTACTCAAAGCTCCCTTATCACTCACGGTGATAAGGATGGTTGGTGGTGATGCTCCACGCCCGGTACAGAc                                   <  1:912877/140‑1 (MQ=255)
                                                                                         cgcACAAACAGCGCGGACTCTGCCGTATAGTCGCGAAAAGAGCTCTGTAGTTTCATCCGCTGCGTTTTCTACTCAAAGCTCCCTTATCACTCACGGTGATAAGGATGGTTGGTGGTGATGCTCCACGCCCGGTACAGCct                                  >  2:179551/1‑140 (MQ=255)
                                                                                            acaAACAGCGCGGACTCTGCCGTATAGTCGCGAAAAGAGCTCTGTAGTTTCATCCGCTGCGTTTTCTACTCAAAGCTCCCTTATCACTCACGGTGATAAGGATGGTTGGTGGTGATGCTCCACGcc                                             <  1:402912/126‑1 (MQ=255)
                                                                                            acaAACAGCGCGGACTCTGCCGTATAGTCGCGAAAAGAGCTCTGTAGTTTCATCCGCTGCGTTTTCTACTCAAAGCTCCCTTATCACTCACGGTGATAAGGATGGTTGGTGGTGATGCTCCACGcc                                             >  2:402912/1‑126 (MQ=255)
                                                                                            acaAACAGCGCGGACTCTGCCGTATAGTCGCGAAAAGAGCTCTGTAGTTTCATCCGCTGCGTTTTCTACTCAAAGCTCCCTTATCACTCACGGTGATAAGGATGGTTGGTGGTGATGCTCCACGCCCGGTACAGActctc                               >  1:273061/1‑140 (MQ=255)
                                                                                               aaCAGCGCGGACTCTGCCGTATAGTCGCGAAAAGAGCTCTGTAGTTTCATCCGCTGCGTTTTCTACTCAAAGCTCCCTTATCACTCACGGTGATAAGGATGGTTGGTGGTGATGCTCCACGCCCGGTACAGACTCTCTGc                            >  1:593215/1‑140 (MQ=255)
                                                                                                  aGCGCGGACTCTGCCGTATAGTCGCGAAAAGAGCTCTGTAGTTTCATCCGCTGCGTTTTCTACTCAAAGCTCCCTTATCACTCACGGTGATAAGGATGGTTGGTGGTGATGCTCCACGCCCGGTACAGACTCTCTGCGAc                         <  2:229782/140‑1 (MQ=255)
                                                                                                    cgcgGACTCTGCCGTATAGTCGCGAAAAGAGCTCTGTAGTTTCATCCGCTGCGTTTTCTACTCAAAGCTCCCTTATCACTCACGGTGATAAGGATGGTTGGTGGTGATGCTCCACGCCCGGTACAGACTCTCTGCGACCa                       <  1:619569/140‑1 (MQ=255)
                                                                                                                         gcgAAAAGAGCTCTGTAGTTTCATCCGCTGCGTTTTCTACTCAAAGCTCCCTTATCACTCACGGTGATAAGGATGGTTGGTGGTGATGCTCCACGCCCGGTACAGACTCTCTGCGACCAGCACGCGAACCAGCGGATggg  >  1:310414/1‑140 (MQ=255)
                                                                                                                         gcgAAAAGAGCTCTGTAGTTTCATCCGCTGCGTTTTCTACTCAAAGCTCCCTTATCACTCACGGTGATAAGGATGGTTGGTGGTGATGCTCCACGCCCGGTACAGACTCTCTGCGACCAGCACGCGAACCAGCGGATggg  >  1:941110/1‑140 (MQ=255)
                                                                                                                                   |                                                                                                                                 
CGGTAAAGCGAACAATTTGCAGATTATACAGGTTGGCGATAAGCACGCCGGTCAGCAGCAAAATCCCCAAAAAGGCGACCAGCGCCCGGCGCACAAACAGCGCGGACTCTGCCGTATAGTCGCGAAAAGAGTTCTGTAGTTTCATCCGCTGCGTTTTCTACTCAAAGCTCCCTTATCACTCACGGTGATAAGGATGGTTGGTGGTGATGCTCCACGCCCGGTACAGACTCTCTGCGACCAGCACGCGAACCAGCGGATGGG  >  NC_000913/668073‑668333

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