Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,252,614 +A coding (311/366 nt) yhaH → putative inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,252,6131.A100.0% 214.7 / NA 60Q104K (CAG→AAG) yhaHputative inner membrane protein
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base A (31/29);  total (31/29)

CTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGCCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGTTTATCGGCTGGCTTATTATCATCGTCTTTAACTGCCAGGCGGGTACGCCAGGCGAAAACCGCTTTGGACCAGATCCAAAACTGGAACCTTAAAACAAATCTGCTCGCCGGATGTCTGACTCATCCGGCTTTTATCATTATGTTGATGAATGAATAATTTACTAAT  >  NC_000913/3252477‑3252742
                                                                                                                                         |                                                                                                                                 
cTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGCCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGTTTATCGGCTGGCTTATTATCATCGTCTTTAACTGCCaag                                                                                                                                <  2:428961/140‑3 (MQ=255)
         tCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGCCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGTTTATCGGCTGGCTTATTATCATCGTCTTTAACTGCCAAGGCGGGTACg                                                                                                                        >  1:4194/1‑140 (MQ=255)
         tCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGCCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGTTTATCGGCTGGCTTATTATCATCGTCTTTAACTGCCaag                                                                                                                                <  1:1033215/131‑3 (MQ=255)
         tCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGCCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGTTTATCGGCTGGCTTATTATCATCGTCTTTAACTGCCaag                                                                                                                                >  2:1033215/1‑129 (MQ=255)
             gtgtTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGCCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGTTTATCGGCTGGCTTATTATCATCGTCTTTAACTGCCAAGGCGGGTACGCCAg                                                                                                                    >  1:133771/1‑140 (MQ=255)
                tttttACCGTGGTGGGCGGTTCAGTTCCGCCGCCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGTTTATCGGCTGGCTTATTATCATCGTCTTTAACTGCCAAGGCGGGTACGCCAGGCg                                                                                                                 >  2:298400/1‑140 (MQ=255)
                  tttACCGTGGTGGGCGGTTCAGTTCCGCCGCCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGTTTATCGGCTGGCTTATTATCATCGTCTTTAACTGCCAAGGCGGGTACGCCAGGCGaa                                                                                                               >  1:595657/1‑140 (MQ=255)
                    tACCGTGGTGGGCGGTTCAGTTCCGCCGCCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGTTTATCGGCTGGCTTATTATCATCGTCTTTAACTGCCAAGGCGGGTACGCCAGGCGaaaa                                                                                                             <  1:298400/140‑1 (MQ=255)
                     aCCGTGGTGGGCGGTTCAGTTCCGCCGCCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGTTTATCGGCTGGCTTATTATCATCGTCTTTAACTGCCAAGGCGGGTACGCCAGGCGAAAAc                                                                                                            <  2:4194/140‑1 (MQ=255)
                       cGTGGTGGGCGGTTCAGTTCCGCCGCCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGTTTATCGGCTGGCTTATTATCATCGTCTTTAACTGCCAAGGCGGGTACGc                                                                                                                       <  1:103113/127‑1 (MQ=255)
                       cGTGGTGGGCGGTTCAGTTCCGCCGCCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGTTTATCGGCTGGCTTATTATCATCGTCTTTAACTGCCAAGGCGGGTACGc                                                                                                                       >  2:103113/1‑127 (MQ=255)
                          ggtggGCGGTTCAGTTCCGCCGCCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGTTTATCGGCTGGCTTATTATCATCGTCTTTAACTGCCAAGGCGGGTACGCCAgg                                                                                                                   >  1:936702/1‑128 (MQ=255)
                          ggtggGCGGTTCAGTTCCGCCGCCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGTTTATCGGCTGGCTTATTATCATCGTCTTTAACTGCCAAGGCGGGTACGCCAgg                                                                                                                   <  2:936702/128‑1 (MQ=255)
                          ggtggGCGGTTCAGTTCCGCCGCCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGTTTATCGGCTGGCTTATTATCATCGTCTTTAACTGCCAAGGCGGGTACGCCAGGCGAAAAc                                                                                                            >  1:385976/1‑135 (MQ=255)
                          ggtggGCGGTTCAGTTCCGCCGCCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGTTTATCGGCTGGCTTATTATCATCGTCTTTAACTGCCAAGGCGGGTACGCCAGGCGAAAAc                                                                                                            <  2:385976/135‑1 (MQ=255)
                           gtggGCGGTTCAGTTCCGCCGCCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGTTTATCGGCTGGCTTATTATCATCGTCTTTAACTGCCAAGGCGGGTACGCCAGGCGAAAACCGCttt                                                                                                      >  1:598273/1‑140 (MQ=255)
                                cGGTTCAGTTCCGCCGCCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGTTTATCGGCTGGCTTATTATCATCGTCTTTAACTGCCAAGGCGGGTACGCCAGGCGAAAACCGCTTTGGAcc                                                                                                 <  1:774613/140‑1 (MQ=255)
                                    tCAGTTCCGCCGCCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGTTTATCGGCTGGCTTATTATCATCGTCTTTAACTGCCAAGGCGGGGACGCCAGGCGAAAAACGCCTTTGACCAGaa                                                                                             >  2:1054175/1‑139 (MQ=255)
                                             ccgccTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGTTTATCGGCTGGCTTATTATCATCGTCTTTAACTGCCAAGGCGGGTACGCCAGGCGAAAACCGCTTTGGACCAGATCCAAAACTg                                                                                    >  1:893417/1‑140 (MQ=255)
                                              cgccTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGTTTATCGGCTGGCTTATTATCATCGTCTTTAACTGCCAAGGCGGGTACGCCAGGCGAAAACCGCTTTGGACCAGATCCAAAACTgg                                                                                   <  2:598273/140‑1 (MQ=255)
                                              cgccTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGTTTATCGGCTGGCTTATTATCATCGTCTTTAACTGCCAAGGCGGGTACGCCAGGCGAAAACCGCTTTGGACCAGATCCAAAACTgg                                                                                   <  2:547419/140‑1 (MQ=255)
                                                              accgCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGTTTATCGGCTGGCTTATTATCATCGTCTTTAACTGCCAAGGCGGGTCCGCCCAGCGAAAACCGCTTTGGACCAGATCCAAAAATGGAAACTTaaaaaaaaac                                                                   >  2:430244/1‑138 (MQ=255)
                                                              accgCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGTTTATCGGCTGGCTTATTATCATCGTCTTTAACTGCCAAGGCGGGTACGCCAGGCGAAAACCGCTTTGGACCAGATCCAAAACTGGAACCTTAAAACAAATc                                                                   >  1:576374/1‑140 (MQ=255)
                                                              accgCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGTTTATCGGCTGGCTTATTATCATCGTCTTTAACTGCCAAGGCGGGTACGCCAGGCGAAAACCGCTTTGGACCAGATCCAAAACTGGAACCTTAAAACAAATc                                                                   >  2:1013458/1‑140 (MQ=255)
                                                              accgCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGTTTATCGGCTGGCTTATTATCATCGTCTTTAACTGCCAAGGCGGGGACGCCAGGCGAAAACCGCTTTGGACCAGATCCAAAACTGGAACCTTAAAACAAATc                                                                   >  1:572620/1‑140 (MQ=255)
                                                                cgCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGTTTATCGGCTGGCTTATTATCATCGTCTTTAACTGCCAAGGCGGGGACGCCAGGCGAAAACCGCTTTTGACCAGATCCAAAACTGGAAACTTAAAACAAAACTg                                                                 >  1:538253/1‑140 (MQ=255)
                                                                  cTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGTTTATCGGCTGGCTTATTATCATCGCCTTTAACTGCCAAGGCGGGTACGCCAGGCGAAAACCGCTTTGGACCAGATCCAAAACTGGAACCTTAAAACAAATCTGCt                                                               <  1:126755/140‑1 (MQ=255)
                                                                   tCGGCGTGGTGGGCACTGCTATTCTTAATCCCGTTTATCGGCTGGCTTATTATCATCGTCTTTAACTGCCAAGGCGGGTACGCCAGGCGAAAACCGCTTTGGACCAGATCCAAAACTGGAACCTTAAAACAAATCTGCTc                                                              >  1:1045963/1‑140 (MQ=255)
                                                                   tCGGCGTGGTGGGCACTGCTATTCTTAATCCCGTTTATCGGCTGGCTTATTATCATCGTCTTTAACTGCCAAGGCGGGTACGCCAGGCGAAAACCGCTTTGGACCAGATCCAAAACTGGAACCTTAAAACAAATCTGCTc                                                              >  2:481649/1‑140 (MQ=255)
                                                                       cGTGGTGGGCACTGCTATTCTTAATCCCGTTTATCGGCTGGCTTATTATCATCGTCTTTAACTGCCAAGGCGGGTACGCCAGGCGAAAACCGCTTTGGACCAGATCCAAAACTGGAACCTTAAAACAAATCTGCTCGCCg                                                          <  2:133771/140‑1 (MQ=255)
                                                                        gtggtgGGCACTGCTATTCTTAATCCCGTTTATCGGCTGGCTTATTATCATCGTCTTTAACTGCCAAGGCGGGTACGCCAGGCGAAAACCGCTTTGGACCAGATCCAAAACTGGAACCTTAAAACAAATCTGCTCGCCgg                                                         <  2:579273/140‑1 (MQ=255)
                                                                            tggGCACTGCTATTCTTAATCCCGTTTATCGGCTGGCTTATTATCATCGTCTTTAACTGCCAAGGCGGGTACGCCAGGCGAAAACCGCTTTGGACCAGATCCAAAACTGGAACCTTAAAACAAATCTGCTCGCCGGATGt                                                     <  1:430244/140‑1 (MQ=255)
                                                                               gCACTGCTATTCTTAATCCCGTTTATCGGCTGGCTTATTATCATCGTCTTTAACTGCCAAGGCGGGTACGCCAGGCGAAAACCGCTTTGCACCAGATCCAAAACTGGAACCTTAAAACAAATCTGCTCGCCGGATGTCTg                                                  >  1:278140/1‑140 (MQ=255)
                                                                                cACTGCTATTCTTAATCCCGTTTATCGGCTGGCTTATTATCATCGTCTTTAACTGCCAAGGCGGGTACGCCAGGCGAAAACCGCTTTGGACCAGATCCAAAACTGGAACCTTAAAACAAATCTGCTCGCCGGATGTCTGa                                                 <  1:1013458/140‑1 (MQ=255)
                                                                                cACTGCTATTCTTAATCCCGTTTATCGGCTGGCTTATTATCATCGTCTTTAACTGCCAAGGCGGGTACGCCAGGCGAAAACCGCTTTGGACCAGATCCAAAACTGGAACCTTAAAACAAATCTGCTCGCCGGATGTCTGa                                                 <  1:1054175/140‑1 (MQ=255)
                                                                                      tATTCTTAATCCCGTTTATCGGCTGGCTTATTATCATCGTCTTTAACTGCCAAGGCGGGTACGCCAGGCGAAAACCGCTTTGGACCAGATCCAAAACTGGAACCTTAAAACAAATCTGCTCGCCGGATGTCTGACTCATc                                           <  1:481649/140‑1 (MQ=255)
                                                                                      tATTCTTAATCCCGTTTATCGGCTGGCTTATTATCATCGTCTTTAACTGCCAAGGCGGGTACGCCAGGCGAAAACCGCTTTGGACCAGATCCAAAACTGGAACCTTAAAACAAATCTGCTCGCCGGATGTCTGACTCATc                                           <  2:73428/140‑1 (MQ=255)
                                                                                      tATTCTTAATCCCGTTTATCGGCTGGCTTATTATCATCGTCTTTAACTGCCAAGGCGGGTACGCCAGGCGAAAACCGCTTTGCACCAGATCCAAAACTGGAACCTTAAAACAAATCTGCTCGCCGGATGTCTGACTCATc                                           <  2:278140/140‑1 (MQ=255)
                                                                                        ttCTTAATCCCGTTTATCGGCTGGCTTATTATCATCGTCTTTAACTGCCAAGGCGGGTACGCCAGGCGAAAACCGCTTTGGACCAGATCCAAAACTGGAACCTTAAAACAAATCTGCTCGCCGGATGTCTGACTCATCCg                                         <  1:245846/140‑1 (MQ=255)
                                                                                             aaTCCCGTTTATCGGCTGGCTTATTATCATTGTCTTTAACTGCCAAGGCGGGTACGCCAGGCGAAAACCGCTTTGGACCAGATCCAAAACTGGAACCTTAAAACAAATCTGCTCGCCGGATGTCTGa                                                 >  1:919995/1‑127 (MQ=255)
                                                                                             aaTCCCGTTTATCGGCTGGCTTATTATCATTGTCTTTAACTGCCAAGGCGGGTACGCCAGGCGAAAACCGCTTTGGACCAGATCCAAAACTGGAACCTTAAAACAAATCTGCTCGCCGGATGTCTGa                                                 <  2:919995/127‑1 (MQ=255)
                                                                                             aaTCCCGTTTATCGGCTGGCTTATTATCATCGTCTTTAACTGCCAAGGCGGGTACGCCAGGCGAAAACCGCTTTGGACCAGATCCAAAACTGGAACCTTAAAACAAATCTGCTCGCCGGATGTCTGa                                                 <  2:486846/127‑1 (MQ=255)
                                                                                             aaTCCCGTTTATCGGCTGGCTTATTATCATCGTCTTTAACTGCCAAGGCGGGTACGCCAGGCGAAAACCGCTTTGGACCAGATCCAAAACTGGAACCTTAAAACAAATCTGCTCGCCGGATGTCTGa                                                 >  1:486846/1‑127 (MQ=255)
                                                                                             aaTCCCGTTTATCGGCTGGCTTATTATCATCGTCTTTAACTGCCAAGGCGGGTACGCCAGGCGAAAACCGCTTTGGACCAGATCCAAAACTGGAACCTTAAAACAAATCTGCTCGCCGGATGTCTGACTCATCCGGCttt                                    <  2:884512/140‑1 (MQ=255)
                                                                                             aaTCCCGTTTATCGGCTGGCTTATTATCATCGTCTTTAACTGCCAAGGCGGGTACGCCAGGCGAAAACCGCTTTGGACCAGATCCAAAACTGGAACCTTAAAACAAATCTGCTCGCCGGATGTCTGACTCATCCGGCttt                                    >  2:28817/1‑140 (MQ=255)
                                                                                             aaTCCCGTTTATCGGCTGGCTTATTATCATCGTCTTTAACTGCCAAGGCGGGGACGCCAGGCGAAAACCGCTTTGGACCAGATCCAAAACTGGAACCTTAAAACAAATCTGCTCGCCGGATGTCTGACTCATCCGGCttt                                    >  2:434018/1‑140 (MQ=255)
                                                                                                cccGTTTATCGGCTGGCTTATTATCATCGTCTTTAACTGCCAAGGCGGGTACGCCAGGCGAAAACCGCTTTGGACCAGATCCAAAACTGGAACCTTAAAACAAATCTGCTCGCCGGATGTCTGACTCATCCGGCTTTTAt                                 >  1:932763/1‑140 (MQ=255)
                                                                                                  cGTTTATCGGCTGGCTTATTATCATCGTCTTTAACTGCCAAGGCGGGTACGCCAGGCGAAAACCGCTTTGGACCAGATCCAAAACTGGAACCTTAAAACAAATCTGCTCGCCGGATGTCTGACTCATCCGGCTTTTATCa                               >  1:165015/1‑140 (MQ=255)
                                                                                                  cGTTTATCGGCTGGCTTATTATCATCGTCTTTAACTGCCAAGGCGGGTACGCCAGGCGAAAACCGCTTTGGACCAGATCCAAAACTGGAACCTTAAAACAAATCTGCTCGCCGGATGTCTGACTCATCCGGCTTTTATCa                               <  2:1045963/140‑1 (MQ=255)
                                                                                                  cGTTTATCGGCTGGCTTATTATCATCGTCTTTAACTGCCAAGGCGGGTACGCCAGGCGAAAACCGCTTTGGACCAGATCCAAAACTGGAACCTTAAAACAAATCTGCTCGCCGGATGTCTGACTCATCCGGCTTTTATCa                               <  1:667318/140‑1 (MQ=255)
                                                                                                    tttATCGGCTGGCTTATTATCATCGTCTTTAACTGCCAAGGCGGGTACGCCAGGCGAAAACCGCTTTGGACCAGATCCAAAACTGGAACCTTAAAACAAATCTGCTCGCCGGATGTCTGACTCATCCGGCTTTTATCAtt                             >  1:728666/1‑140 (MQ=255)
                                                                                                     ttATCGGCTGGCTTATTATCATCGTCTTTAACTGCCAAGGCGGGTACGCCAGGCGAAAACCGCTTTGGACCAGATCCAAAACTGGAACCTTAAAACAAATCTGCTCGCCGGATGTCTGACTCATCCGGCTTTTATCAtt                             <  2:405999/139‑1 (MQ=255)
                                                                                                     ttATCGGCTGGCTTATTATCATCGTCTTTAACTGCCAAGGCGGGTACGCCAGGCGAAAACCGCTTTGGACCAGATCCAAAACTGGAACCTTAAAACAAATCTGCTCGCCGGATGTCTGACTCATCCGGCTTTTATCAtt                             >  1:405999/1‑139 (MQ=255)
                                                                                                      tATCGGCTGGCTTATTATCATCGTCTTTAACTGCCAAGGCGGGTACGCCAGGCGAAAACCGCTTTGGACCAGATCCAAAACTGGAACCTTAAAACAAATCTGCTCGCCGGATGTCTGACTCATCCGGCTTTTATCATTAt                           <  2:595657/140‑1 (MQ=255)
                                                                                                       aTCGGCTGGCTTATTATCATCGTCTTTAACTGCCAAGGCGGGGACGCCAGGCGAAAACCGCTTTGGACCAGATCCAAAACTGGAACGTTAAAACAAATCTGCTCGCCGGATGTCTGACTCATCCGGCTTTTATCATTATg                          >  2:953411/1‑140 (MQ=255)
                                                                                                            ctggctTATTATCATCGTCTTTAACTGCCAAGGCGGGTACGCCAGGCGAAAACCGCTTTGGACCAGATCCAAAACTGGAACCTTAAAACAAATCTGCTCGCCGGATGTCTGACTCATCCGGCTTTTATCAtt                             <  2:592818/132‑1 (MQ=255)
                                                                                                            ctggctTATTATCATCGTCTTTAACTGCCAAGGCGGGTACGCCAGGCGAAAACCGCTTTGGACCAGATCCAAAACTGGAACCTTAAAACAAATCTGCTCGCCGGATGTCTGACTCATCCGGCTTTTATCAtt                             >  1:592818/1‑132 (MQ=255)
                                                                                                              ggctTATTATCATCGTCTTTAACTGCCAAGGCGGGTACGCCAGGCGAAAACCGCTTTGGACCAGATCCAAAACTGGAACCTTAAAACAAATCTGCTCGCCGGATGTCTGACTCATCCGGCTTTTATCATTATGTtgatga                   <  1:59717/140‑1 (MQ=255)
                                                                                                                         atcGTCTTTAACTGCCAAGGCGGGTACGCCAGGCGAAAACCGCTTTGGACCAGATCCAAAACTGGAACCTTAAAACAAATCTGCTCGCCGGATGTCTGACTCATCCGGCTTTTATCATTATGTTGatgaatg                >  2:340997/1‑132 (MQ=255)
                                                                                                                         atcGTCTTTAACTGCCAAGGCGGGTACGCCAGGCGAAAACCGCTTTGGACCAGATCCAAAACTGGAACCTTAAAACAAATCTGCTCGCCGGATGTCTGACTCATCCGGCTTTTATCATTATGTTGatgaatg                <  1:340997/132‑1 (MQ=255)
                                                                                                                         atcGTCTTTAACTGCCAAGGCGGGTACGCCAGGCGAAAACCGCTTTGGACCAGATCCAAAACTGGAACCTTAAAACAAATCTGCTCGCCGGATGTCTGACTCATCCGGCTTTTATCATTATGTTGATGAATGAATAAttt        >  2:896135/1‑140 (MQ=255)
                                                                                                                            gTCTTTAACTGCCAAGGCGGGTACGCCAGGCGAAAACCGCTTTGGACCAGATCCAAAACTGGAACCTTAAAACAAATCTGCTCGCCGGATGTCTGACTCATCCGGCTTTTATCATTATGTTGATGAATGAATAATTTACt     <  2:932763/140‑1 (MQ=255)
                                                                                                                               tttAACTGCCAAGGCGGGTACGCCAGGCGAAAACCGCTTTGGACCAGATCCAAAACTGGAACCTTAAAACAAATCTGCTCGCCGGATGTCTGACTCATCCGGCTTTTATCATTATGTTGATGAATGAATAATTTACTAAt  >  2:140624/1‑140 (MQ=255)
                                                                                                                                         |                                                                                                                                 
CTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGCCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGTTTATCGGCTGGCTTATTATCATCGTCTTTAACTGCCAGGCGGGTACGCCAGGCGAAAACCGCTTTGGACCAGATCCAAAACTGGAACCTTAAAACAAATCTGCTCGCCGGATGTCTGACTCATCCGGCTTTTATCATTATGTTGATGAATGAATAATTTACTAAT  >  NC_000913/3252477‑3252742

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: