Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,468,121:1 +A 8.6% coding (819/1047 nt) mltB ← membrane‑bound lytic murein transglycosylase B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,468,1211.A 8.6% 210.2 / 15.2 58coding (819/1047 nt)mltBmembrane‑bound lytic murein transglycosylase B
Reads supporting (aligned to +/- strand):  ref base . (26/27);  new base A (3/2);  total (29/29)
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.35e-01

AGTGATCGCGTAGAAGTTCTTCAGGCCCATCCAGTACTCGGGGCCGCTCTCGCCCTCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGCCGATGGCATCGTCAGGGTTGTTCCAGATATTGAT  >  NC_002947/5467974‑5468267
                                                                                                                                                    |                                                                                                                                                  
aGTGATCGCGTAGAAGTTCTTCAGGCCCATCCAGTACTCGGGGCCGCTCTCGCCCTCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTc                                                                                                                                                  >  1:1183803/1‑150 (MQ=255)
     tCGCGTAGAAGTTCTTCAGGCCCATCCAGTACTCGGGGCCGCTCTCGCCCTCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGcatcgagc                                                                                                                                              <  2:396393/150‑8 (MQ=255)
       gcgTAGAAGTTCTTCAGGCCCATCCAGTACTCGGGGCCGCTCTCGCCCTCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCg                                                                                                                                           >  1:1080833/1‑150 (MQ=255)
       gcgTAGAAGTTCTTCAGGCCCATCCAGTACTCGGGGCCGCTCTCGCCCTCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCg                                                                                                                                           <  2:1162024/150‑1 (MQ=255)
               gTTCTTCAGGCCCATCCAGTACTCGGGGCCGCTCTCGCCCTCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAg                                                                                                                                   <  2:1055684/150‑1 (MQ=255)
               gTTCTTCAGGCCCATCCAGTACTCGGGGCCGCTCTCGCCCTCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAg                                                                                                                                   <  1:47359/150‑1 (MQ=255)
                 tcttcAGGCCCATCCAGTACTCGGGGCCGCTCTCGCCCTCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGc                                                                                                                                 <  1:696267/150‑1 (MQ=255)
                        gCCCATCCAGTACTCGGGGCCGCTCTCGCCCTCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTCGGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAAGCCTtcgac                                                                                                                          >  2:303585/1‑148 (MQ=255)
                        gCCCATCCAGTACTCGGGGCCGCTCTCGCCCTCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTtcgtc                                                                                                                          >  1:888777/1‑150 (MQ=255)
                          ccATCCAGTACTCGGGGCCGCTCTCGCCCTCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCggtg                                                                                                                                        >  2:322382/1‑134 (MQ=255)
                          ccATCCAGTACTCGGGGCCGCTCTCGCCCTCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCggtg                                                                                                                                        <  1:322382/134‑1 (MQ=255)
                                aGTACTCGGGGCCGCTCTCGCCCTCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCt                                                                                                                  >  2:71569/1‑150 (MQ=255)
                                  tACTCGGGGCCGCTCTCGCCCTCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCa                                                                                                                >  2:67198/1‑150 (MQ=255)
                                    cTCGGGGCCGCTCTCGCCCTCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCAcc                                                                                                              >  2:1086189/1‑150 (MQ=255)
                                     tCGGGGCCGCTCTCGCCCTCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCt                                                                                                             >  1:1155303/1‑150 (MQ=255)
                                       ggggCCGCTCTCGCCCTCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTc                                                                                                           <  1:71569/150‑1 (MQ=255)
                                               tctcGCCCTCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGACAGGCCTTAGTCGGCGCGCTCACCTTCCACCCAgg                                                                                                   >  1:601485/1‑150 (MQ=255)
                                                 tcGCCCTCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAg                                                                                                                                   <  2:176063/116‑1 (MQ=255)
                                                 tcGCCCTCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAg                                                                                                                                   >  1:176063/1‑116 (MQ=255)
                                                    cccTCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGa                                                                                              <  2:1155303/150‑1 (MQ=255)
                                                      cTCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGAcg                                                                                            <  2:875482/150‑2 (MQ=255)
                                                        cGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGa                                                                                          <  1:316643/150‑1 (MQ=255)
                                                         gagaCGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGAc                                                                                         >  2:1053357/1‑150 (MQ=255)
                                                         gagaCGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTAACCTACCACCCAGGCGCAACTGAc                                                                                         <  2:838970/150‑1 (MQ=255)
                                                             cGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACcaca                                                                                     <  1:67198/150‑1 (MQ=255)
                                                                           cGGCAGATCACCGCGCAGCGAATTATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCAc                                                                       >  2:878139/1‑150 (MQ=255)
                                                                            ggCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCAcc                                                                      >  2:783916/1‑150 (MQ=255)
                                                                                  tcatcGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCg                                                                <  2:257130/150‑1 (MQ=255)
                                                                                     tcGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTgg                                                             <  2:1037489/150‑1 (MQ=255)
                                                                                     tcGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTgg                                                             <  2:66388/150‑1 (MQ=255)
                                                                                       gcgcAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCt                                                                                                             >  1:1409301/1‑100 (MQ=255)
                                                                                       gcgcAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCt                                                                                                             <  2:1409301/100‑1 (MQ=255)
                                                                                       gcgcAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGgcg                                                           >  2:1258258/1‑150 (MQ=255)
                                                                                                 tCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGCTGGTGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAg                                                 <  1:1086189/150‑1 (MQ=255)
                                                                                                 tCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAg                                                 <  1:1053357/150‑1 (MQ=255)
                                                                                                       aCCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTtcgtc                                                                                                                          <  2:1072543/71‑1 (MQ=255)
                                                                                                       aCCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTtcgtc                                                                                                                          >  1:1072543/1‑71 (MQ=255)
                                                                                                                  aaGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGCCGAt                                <  2:1240402/150‑1 (MQ=255)
                                                                                                                  aaGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGCCGAt                                >  2:50228/1‑150 (MQ=255)
                                                                                                                   aGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCTATGCCTGTGGTCAGGCCTTCGTCGGCTCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGCCGATg                               <  2:601485/150‑1 (MQ=255)
                                                                                                                        aCGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGCCGATGGCATc                          <  1:1343770/150‑1 (MQ=255)
                                                                                                                          gCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAg                                                 <  1:402582/125‑1 (MQ=255)
                                                                                                                          gCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAg                                                 >  2:402582/1‑125 (MQ=255)
                                                                                                                                 cccTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTg                                                       >  1:365682/1‑112 (MQ=255)
                                                                                                                                 cccTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTg                                                       <  2:365682/112‑1 (MQ=255)
                                                                                                                                  ccTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCCGCCTAGGCGCTTGAAGTAGCTGGCCACGCTGCCGATGGCATCGTCAGGgttg                >  1:87072/1‑150 (MQ=255)
                                                                                                                                   cTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGCCGATGGCATCGTCAGGgtttt               >  2:1363050/1‑148 (MQ=255)
                                                                                                                                       cGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGCCGATGGCATCGTCAGGGTTGTTCCa           >  2:668385/1‑150 (MQ=255)
                                                                                                                                        gTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGCCGATGGCATCGTCAGGGTTGTTCCAg          >  2:126985/1‑150 (MQ=255)
                                                                                                                                          gaccACCGGCATCGAGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGCCGATGGCATCGTCAGGGTTGTTCCAGa         >  1:706007/4‑149 (MQ=255)
                                                                                                                                          gaccACCGGCATCGAGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGCCGATGGCATCGTCAGGGTTGTTCCAGa         <  2:706007/146‑1 (MQ=255)
                                                                                                                                            ccACCGGCATCGAGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGCCGAt                                <  1:1376499/123‑1 (MQ=255)
                                                                                                                                            ccACCGGCATCGAGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGCCGAt                                >  2:1376499/2‑124 (MQ=255)
                                                                                                                                             cACCGGCATCGAGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGCCGATGGCATCGTCAGGGTTGTTCCAGATAtt     >  1:1209718/1‑150 (MQ=255)
                                                                                                                                               ccGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGCCGATGGCATCGTCAGGGTTGTTCCAGATATTGa   >  2:10078/1‑150 (MQ=255)
                                                                                                                                                cGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCCCCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGCCGATGGCATCGTCAGGGTTTTTCCAGATATTGat  <  1:1363050/150‑1 (MQ=255)
                                                                                                                                                cGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGCCGATGGCATCGTCAGGGTTGTTCCAGATATTGat  <  1:878139/150‑1 (MQ=255)
                                                                                                                                                  gCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGc                                    <  1:1419096/114‑1 (MQ=255)
                                                                                                                                                  gCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGc                                    >  2:1419096/1‑114 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
AGTGATCGCGTAGAAGTTCTTCAGGCCCATCCAGTACTCGGGGCCGCTCTCGCCCTCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGCCGATGGCATCGTCAGGGTTGTTCCAGATATTGAT  >  NC_002947/5467974‑5468267

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