Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 2,532,298 A→T 100% T15T (ACA→ACT USA300HOU_2392 → FemAB family peptidoglycan biosynthesis protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007302,532,2980AT100.0% 189.9 / NA 61T15T (ACA→ACTUSA300HOU_2392FemAB family peptidoglycan biosynthesis protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (23/38);  total (23/38)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GTGTATTTGTAATCAGTTTTTGATAAAATAAAAAAATAGCATTATTAATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACAACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTAAAGATGACAATGGTCAAGTGATTGCAGGATGCTTATTGACAGAAGCA  >  CP000730/2532162‑2532433
                                                                                                                                        |                                                                                                                                       
gtgtATTTGTAATTAGTTTTTGATAAAATAAAAAAATAGCATTATTAATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTActacaa                                                                                                                                     <  1:206707/141‑1 (MQ=255)
gtgtATTTGTAATCAGTTTTTGATAAAATAAAAAAATAGCATTATTAATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTAatacaa                                                                                                                                     >  1:665159/1‑141 (MQ=255)
gtgtATTTGTAATCAGTTTTTGATAAAATAAAAAAATAGCATTATTAATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTCTActacaa                                                                                                                                     >  2:551435/1‑141 (MQ=255)
   tATTTGTAATCAGTTTTTGATAAAATAAAAAAATAGCATTATTAATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTTAATGCGTTTACTACAAAGc                                                                                                                                  <  1:369919/141‑1 (MQ=255)
   tATTTGGAATCGGTTTTAGATAGAGTTAAAAAATAGCATTAGTAAGAAATTCAAAAATTATAGTACAATTCAGAGAGGGATGTGATAAAATGAACTTTGTAACGTTGACGTCGGAGGAGTCCAATGCGTTTACTACAAAGc                                                                                                                                  <  1:419324/141‑1 (MQ=255)
              aGTTTTTGATAAAATAAAAAAATAGCATTATTAATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACAt                                                                                                                       >  1:602558/1‑141 (MQ=255)
                ttttgGACAACATAATAAAATAGCATTATTAATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTa                                                                                                                     >  2:671324/1‑141 (MQ=255)
                   tttataaaataaaAAAATAGCATTATTAATAAATTCAAAAATTAGAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTGCAATGCGTTTACTACAAAGCATTTTTCACATTacac                                                                                                                  <  1:199186/138‑1 (MQ=255)
                   ttGATAAAATAAAAAAATAGCATTATTAATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTacac                                                                                                                  <  2:232363/141‑1 (MQ=255)
                    tGATAAAATAAAAAAATAGCATTATTAATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTacaca                                                                                                                 >  2:480651/1‑141 (MQ=255)
                    tGAGAAAATAAAAAAATAGCATTATTAATAAATTCAAAAATTATAGTTGAATTCAGTGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTacaca                                                                                                                 <  1:511927/141‑1 (MQ=255)
                        aaaataaaaAAATAGCATTATTAATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTCAAGCACTTTTTTCACATTACACACAAt                                                                                                             <  2:668572/141‑1 (MQ=255)
                         aaataaaaAAATAGCATTATTAATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATc                                                                                                            >  1:475634/1‑141 (MQ=255)
                         aaataaaaAAATAGCATTATTAATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACACTACACACAATc                                                                                                            >  2:270279/1‑141 (MQ=255)
                         aaataaaaAAATAGCATTATTAATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTAGCATAACACACaacc                                                                                                            >  1:87581/1‑139 (MQ=255)
                         aaataaaaAAATAGCATTATTAATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGATTACTACAAAGCATTTTTCACATTACACACAATc                                                                                                            >  1:317609/1‑141 (MQ=255)
                               aaaaaTAGCATTATTAATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATAGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGATCAATGCGTTTACTACAAAGCATTTTTCAAATTACACACAACCAGCTAt                                                                                                      >  1:267021/1‑141 (MQ=255)
                                   aTAGCATTATTAATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCAt                                                                                                  <  2:100441/141‑1 (MQ=255)
                                        attattAATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACaa                                                                                             >  2:207347/1‑141 (MQ=255)
                                              aataaATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATag                                                                                       <  2:385036/141‑1 (MQ=255)
                                                 aaaTTCAAAAGTTATAGTTCAATTCAGAGCGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTAGTACAAAGCATTTTTCACATTGCACAAAATCGGCGATTCATTAGAGTAATAGAGt                                                                                    >  1:187984/1‑141 (MQ=255)
                                                 aaaTTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGt                                                                                    >  2:57853/1‑141 (MQ=255)
                                                    ttCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGa                                                                                 <  1:251570/141‑1 (MQ=255)
                                                     tCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGAt                                                                                <  1:603856/141‑1 (MQ=255)
                                                          aaTTATCGTTCAACGCAGAGAGGGATGGGATAAAATGAACTTTCGAGCGTGGGCTGCGGATGAGTCCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTaaa                                                                           <  2:488823/141‑1 (MQ=255)
                                                          aaTTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTaaa                                                                           <  1:260849/141‑1 (MQ=255)
                                                              ataGTTCAATTCAGAGAGGGATGGGCTAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGc                                                                       <  1:158895/141‑1 (MQ=255)
                                                              ataGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGc                                                                       <  2:569885/141‑1 (MQ=255)
                                                               taGTGCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGGTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCg                                                                      <  1:45532/141‑1 (MQ=255)
                                                                aGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACCCAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGa                                                                     >  1:114011/1‑141 (MQ=255)
                                                                aGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGa                                                                     >  2:565308/1‑141 (MQ=255)
                                                                aGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGa                                                                     >  1:114017/1‑141 (MQ=255)
                                                                   tCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGAtgt                                                                  <  2:34100/141‑1 (MQ=255)
                                                                     aaTTCAGAGAGGGATGGGATAAAATGAACTTTGTGACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGc                                                                <  1:279995/141‑1 (MQ=255)
                                                                       ttCAGAGAGGGATGGGATAAAATGAACGTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCAt                                                              <  1:523311/141‑1 (MQ=255)
                                                                         cAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCt                                                            <  1:672678/141‑1 (MQ=255)
                                                                         cAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCt                                                            <  1:36370/141‑1 (MQ=255)
                                                                         cAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCt                                                            <  1:571369/141‑1 (MQ=255)
                                                                          agagagTGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCtt                                                           <  2:166625/141‑1 (MQ=255)
                                                                           gagagGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTg                                                          >  2:597111/1‑141 (MQ=255)
                                                                            agagGGATGGGAGAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATGACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGt                                                         >  1:490573/1‑141 (MQ=255)
                                                                                  atgggatAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGt                                                   <  1:287342/141‑1 (MQ=255)
                                                                                  atgggatAAAATGAACTTTGTAACGGTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGt                                                   <  2:426608/141‑1 (MQ=255)
                                                                                   tgggatAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGtt                                                  <  2:47798/141‑1 (MQ=255)
                                                                                   ggggatAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGtt                                                  <  2:269790/140‑1 (MQ=255)
                                                                                     ggatAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTTACTTTTTTTTCAATAAACTATTCATTACAATCATTGAGTTGATTTATAAAGCAATTTGCATTTTGTTTTGGTTaa                                                >  2:193169/1‑141 (MQ=255)
                                                                                          aaaTGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCTTTTACTACAAATCATTTTTCACATTACACACTATCAGCTATTCATTACAATCATAGAGTCGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTAAAGATg                                           <  1:22985/141‑1 (MQ=255)
                                                                                            aTGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTAAAGATGAc                                         >  1:198058/1‑141 (MQ=255)
                                                                                              gAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCCTTTTTCACATTACACACAATCAGCTATTCATTACAATGATAGAGTTGATTTAAAAGGCGATGTGCATATTGTAGGGCTTAAAGATGACaa                                       >  1:96861/1‑141 (MQ=255)
                                                                                               aaCTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTAAAGATGACAAt                                      <  1:231580/141‑1 (MQ=255)
                                                                                                     gTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTAAAGATGACAATGGTCaa                                >  1:106010/1‑141 (MQ=255)
                                                                                                         cGTTGACTTCGGATGAGTGCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTAAAGATGACAATGGTCAAGTGa                            <  2:47906/141‑1 (MQ=255)
                                                                                                            tGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTAAAGATGACAATGGTCAAGTGATTg                         <  2:569568/141‑1 (MQ=255)
                                                                                                                 tCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATATCAATATCTTTGTGCATCTTGTAGGGGTTAAAGATGACAATGGTCAAGTGATTGCGGGa                    >  1:193327/1‑141 (MQ=255)
                                                                                                                  cGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTAAAGATGACAATGGTCAAGTGATTGCAGGAt                   <  1:345212/141‑1 (MQ=255)
                                                                                                                   ggATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTAAAGATGACAATGGTCAAGTGATTGCAGGATg                  <  2:545936/141‑1 (MQ=255)
                                                                                                                    gATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTAAAGATGACAATGGTCAAGTGATTGCAGGATGc                 >  1:567366/1‑141 (MQ=255)
                                                                                                                    gATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTAAAGATGACAATGGTCAAGTGATTGCAGGATGc                 >  1:513757/1‑141 (MQ=255)
                                                                                                                       gAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTAAAGATGACAATGGTCAAGTGATTGCAGGATGCTTa              <  1:211620/141‑1 (MQ=255)
                                                                                                                        gggtctaTGTGTTTACTACAAAGCACTTTTCACAGTTCACACAATGAGTTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTAAAGATGACAATGGTCAAGTGATTGCAGGGTGCTTAt             <  2:586169/135‑1 (MQ=255)
                                                                                                                        aGTTCAATGCGTTGACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTAAAGATGACAATGGTCAAGTGATTGCAGGATGCTTAt             <  2:434427/141‑1 (MQ=255)
                                                                                                                                  ttttACTACAAAACATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTAAAGATGACAATGGTCAAGTGATTGCAGGATGCTTATTGACAGAAgc   <  2:326961/140‑1 (MQ=255)
                                                                                                                                   tttACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTAAAGATGACAATGGTCAAGTGATTGCAGGATGCTTATTGACAGAAgca  <  2:267021/141‑1 (MQ=255)
                                                                                                                                   tttACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTAAAGATGACAATGGTCAAGTGATTGCAGGATGCTTATTGACAGAAgca  <  2:75022/141‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                       
GTGTATTTGTAATCAGTTTTTGATAAAATAAAAAAATAGCATTATTAATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACAACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTAAAGATGACAATGGTCAAGTGATTGCAGGATGCTTATTGACAGAAGCA  >  CP000730/2532162‑2532433

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