Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_007793 1,862,563 T→C E287E (GAA→GAG)  SAUSA300_RS09220 ← cell division protein FtsK

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0077931,862,5630TC100.0% 230.7 / NA 84E287E (GAA→GAG) SAUSA300_RS09220cell division protein FtsK
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (6/78);  total (6/78)

CTTGCTCTTCGTTTTTATGTTCTTCACGCTTTTGTTGCAACGCTTTTTTAAATCGACGCTTTTGAAGCACTTTACGTTCTCGTTCACGTCTAATTTCTTCAACAATTTGTGAAGCATAAATATTTTCAATTTTGATAGTATTATCAACTTTTGAATAATTAGGCATTGATTTTGTTGTTGATGATAC  >  NC_007793/1862466‑1862652
                                                                                                 |                                                                                         
tccgatctTCGTTTTTATGTTCTTCACGCTTTTGTTGCAACGCTTTTTTAAATCGACGCTTTTGAAGCACTTTACGTTCTCGTTCACGTCTAATttcctc                                                                                         <  1:2432890/95‑1 (MQ=255)
  cgatctTCGTTTTTATGTTCTTCACGCTTTTGTTGCAACGCTTTTTTAAATCGACGCTTTTGAAGCACTTTACGTTCTCGTTCACGTCTAATTTCCTcaa                                                                                       <  1:6399825/97‑1 (MQ=255)
    cttccgatctttGTGTTCTTCACGCTTTTGTTGCAACGCTTTTTTAAATCGACGCTTTTGAAGCACTTTACGTTCTCGTTCACGTCTAATTTCCTcaaca                                                                                     <  1:4696836/91‑1 (MQ=255)
      ctTCGTTTTTATGTTCTTCACGCTTTTGTTGCAACGCTTTTTTAAATCGACGCTTTTGAAGCACTTTACGTTCTCGTTCACGTCTAATTTCCTCAACAAt                                                                                   >  1:4499327/1‑100 (MQ=255)
       tccgatcttATGTTCTTCGCGCTTTTGTTGCAACGCTTTTTTAAATCGACGCTTTTGAAGCACTTTACGTTCTCGTTCACGTCTAATTTCCTcaacaagc                                                                                  <  1:5013176/93‑3 (MQ=255)
        tCGCTTTTATGTTCTTCACGCTTTTGTTGCAACGCTTTTTTAAATCGACGCTTTTGAAGCACTTTACGTTCTCGTTCACGTCTAATTTCCTcaacaccac                                                                                 <  1:3952734/100‑5 (MQ=255)
         cGTTCTTATGTTCTTCACGCTTTTGTTGCAACGCTTTTTTAAATCGACGCTTTTGAAGCACTTTACGTTCTCGTTCACGTCTAATTTCCTCAACAAttta                                                                                <  1:6454025/100‑2 (MQ=255)
           tttttATGTTCTTCACGCTTTTGTTGCAACGCTTTTTTAAATCGACGCTTTTGAAGCACTTTACGTTCTCGTTCACGTCTAATTTCCTCAACAATTtgtg                                                                              <  1:4656074/100‑1 (MQ=255)
           tttttATGTTCTTCACGCTTTTGTTGCAACGCTTTTTTAAATCGACGCTTTTGAAGCACTTTACGTTCTCGTTCACGTCTAATTTCCTCAACAATTtgtg                                                                              <  1:4358679/100‑1 (MQ=255)
           tctttATGCTCTTCACGCTTTTGTTGCAACGCTTTTTTAAATCGACGCTTTTGAAGCACTTTACGTTCTCGTTCACGTCTAATTTCCTCAACAATTtgtg                                                                              <  1:764605/98‑1 (MQ=255)
           ctcttATGTTCTTCACGCTTTTGTTGCAACGCTTTTTTAAATCGACGCTTTTGAAGCACTTTACGTTCTCGTTCACGTCTAATTTCCTCAACAATTtgtg                                                                              <  1:1772869/97‑1 (MQ=255)
            tcttATGTTCTTCACGCTTTTGTTGCAACGCTTTTTTAAATCGACGCTTTTGAAGCACTTTACGTTCTCGTTCACGTCTAATTTCCTCAACAATTTGTGa                                                                             <  1:188496/98‑1 (MQ=255)
               tccgatctttACGCTTTTGTTGCAACGCTTTTTTAAATCGACGCTTTTGAAGCACTTTACGTTCTTGTTGGCGTCTAATTTCCTCAACAATTtgtggcta                                                                          <  1:3365566/95‑5 (MQ=255)
               tccgatcttcACGCTTTTGTTGCAACGCTTTTTTAAATCGACGCTTTTGAAGCACTTTACGTTCTCGTTCACGTCTAATTTCCTCAACAATTTGTGAAgg                                                                          <  1:810636/95‑2 (MQ=255)
                 gatctctcACGCTTCTGTTGCAACGCTTTTTTAAATCGACGCTTTTGAAGCACTTTACGTTCTCGTTCACGTCTAATTTCCTCAACAATTTGTGAAGCAt                                                                        <  1:4886162/94‑1 (MQ=255)
                  gatctccACGCTTTTGTTGCAACGCTCTTTTAAATCGACGCTTTTGAAGCACTTTACGTTCTCGTTCACGTCTAATTTCCTCAACAATTTGTGAAGCATa                                                                       <  1:57715/98‑1 (MQ=255)
                   ttcttcACGCTTTTGTTGCAACGCTTTTTTAAATCGACGCTTTTGAAGCACTTTACGTTCTCGTTCACGTCTAATTTCCTCAACAATTTGTGAAGCATat                                                                      <  1:2362815/100‑2 (MQ=255)
                    tcttcACGCTTTTGTTGCAACGCTTTTTTAAATCGACGCTTTTGAAGCACTTTACGTTCTCGTTCACGTCTAATTTCCTCAACAATTTGTGAAGCATagt                                                                     <  1:2219110/100‑3 (MQ=255)
                       tcACGCTTTTGTTGCAACGCTTTTTTAAATCGACGCTTTTGAAGCACTTTACGTTCTCGTTCACGTCTAATTTCCTCAACAATTTGTGAAGCATAAatag                                                                  <  1:3231934/100‑2 (MQ=255)
                         aCGCTTTTGTTGCAACGCTTTTTTAAATCGACGCTTTTGAAGCACTTTACGTTCTCGTTCACGTCTAATTTCCTCAACAATTTGTGAAGCATAAATAttt                                                                >  1:4260201/1‑100 (MQ=255)
                                  ttgCAACGCTTTTTTAAATCGACGCTTTTGAAGCACTTTACGTTCTCGTTCACGTCTAATTTCCTCAACAATTTGTGAAGCATAAATATTCTCAATTTTg                                                       <  1:3727626/100‑1 (MQ=255)
                                   tgCAACGCTTTTTTAAATCGACGCTTTTGAAGCACTTTACGTTCTCGTTCACGTCTAATTTCCTCAACAATTTGTGAAGCATAAATATTTTCAATTTTGa                                                      >  1:2357378/1‑100 (MQ=255)
                                    gCAACGCTTTTTTAAATCGACGCTTTTGAAGCACTTTACGTTCTCGTTCACGTCTGATTTCCTCAACAATGTGTGAAGCATAAATATTTTCAATGTTGAt                                                     <  1:5261423/100‑1 (MQ=255)
                                    gCAACGCTTTTTTAAATCGACGCTTTTGAAGCACTTTACGTTCTCGTTCACGTCTAATTTCCTCAACAATTTGTGAAGCATAAATATTTTCAATTTTGAt                                                     <  1:1143541/100‑1 (MQ=255)
                                    gCAACGCTTTTTTAAATCGACGCTTTTGAAGCACTTTACGTTCTCGTTCACGTCTAATTTCCTCAACAATTTGTGAAGCATAAATATTTTCAATTTTGAt                                                     <  1:1512748/100‑1 (MQ=255)
                                    gCAACGCTTTTTTAAATCGACGCTTTTGAAGCACTTTACGTTCTCGTTCACGTCTAATTTCCTCAACAATTTGTGAAGCATAAATATTTTCAATTTTGAt                                                     <  1:3949039/100‑1 (MQ=255)
                                    gCAACGCTTTTTTAAATCGACGCTTTTGAAGCACTTTACGTTCTCGTTCACGTCTAATTTCCTCAACAATTTGTGAAGCATAAATATTTTCAATTTTGAt                                                     <  1:498670/100‑1 (MQ=255)
                                    gCAACGCTTTTTTAAATCGACGCTTTTGAAGCACTTTACGTTCTCGTTCACGTCTAATTTCCTCAACAATTTGTGAAGCATAAATATTTTCAATTTTGAt                                                     <  1:5137250/100‑1 (MQ=255)
                                    gCAACGCTTTTTTAAATCGACGCTTTTGAAGCACTTTACGTTCTCGTTCACGTCTAATTTCCTCAACAATTTGTGAAGCATAAATATTTTCAATTTTGAt                                                     <  1:2425817/100‑1 (MQ=255)
                                    gCAACGCTTTTTTAAATCGACGCTTTTGAAGCACTTTACGTTCTCGTTCACGTCTAATTTCCTCAACAATTTGTGAAGCATAAATATTTTCAATTTTGAt                                                     <  1:5946704/100‑1 (MQ=255)
                                    gCAACGCTTTTTTAAATCGACGCTTTTGAAGCACTTTACGTTCTCGTTCACGTCTAATTTCCTCAACAATTTGTGAAGCATAAATATTTTCAATTTTGAt                                                     <  1:606902/100‑1 (MQ=255)
                                    gCAACGCTTTTTTAAATCGACGCTTTTGAAGCACTTTACGTTCTCGTTCACGTCTAATTTCCTCAACAATTTGTGAAGCATAAATATTTTCAATTTTGAt                                                     <  1:2441385/100‑1 (MQ=255)
                                      aaCGCTTTTTTAAATCGACGCTTTTGAAGCACTTTACGTTCTCGTTCACGTCTAATTTCCTCAACAATTTGTGAAGCATAAATATTTTCAATTTTGATAg                                                   <  1:2163420/100‑1 (MQ=255)
                                      aaCGCTTTTTTAAATCGACGCTTTTGAAGCACTTTACGTTCTCGTTCACGTCTAATTTCCTCAACAATTGGTGAAGCATAAATATTTTCAATTTTGATAg                                                   <  1:6641433/100‑1 (MQ=255)
                                       aCGCTTTTTTAAATCGACGCTTTTGAAGCACTTTACGTTCTCGTTCACGTCTAATTTCCTCAACAATTTGTGAAGCATAAATATTTTCAATTTTGATAGt                                                  <  1:1013549/100‑1 (MQ=255)
                                       aCGCTTTTTTAAATCGACGCTTTTGAAGCACTTTACGTTCTCGTTCACGTCTAATTTCCTCAACAATTTGTGAAGCATAAATATTTTCAATTTTGATAGt                                                  <  1:2239807/100‑1 (MQ=255)
                                         gCTTTTTTAAATCGACGCTTTTGAAGCACTTTACGTTCTCGTTCACGTCTAATTTCCTCAACAATTTGTGAAGCATAAATATTTTCAATTTTGATAGtag                                                <  1:5995063/100‑2 (MQ=255)
                                          cTTTTTTAAATCGACGCTTTTGAAGCACTTTACGTTCTCGTTCACGTCTAATTTCCTCAACAATTTGTGAAGCATAAATATTTTCAATTTTGATAGtatt                                               <  1:5462825/100‑1 (MQ=255)
                                          cTTTTTTAAATCGACGCTTTTGAAGCACTTTACGTTCTCGTTCACGTCTAATTTCCTCAACAATTTGTGAAGCATAAATATTTTCAATTTTGATAGtatt                                               <  1:5983220/100‑1 (MQ=255)
                                           ttttttAAATCGACGCTTTTGAAGCACTTTACGTTCTCGTTCACGTCTAATTTCCTCAACAATTTGTGAAGCATAAATATTTTCAATTTTGATAGtattt                                              <  1:3070653/100‑2 (MQ=255)
                                                            tttGAAGCACTTTATGTTCTCGTTCACGTCTAATTTCCTCAACAATTTGTGAAGCATAAATATTTTCAATTTTGATAGTATTATCAACTTTTGAATAAtt                             <  1:4871672/100‑1 (MQ=255)
                                                                aaGCACTTTACGTTCTCGTTCCCGTCTAATTTCCTCAACAATTTGTGAAGCATAAATATTTTCAATTTTGATAGTATTATCAACTTTTGAATAATTAGGc                         <  1:1571133/100‑1 (MQ=255)
                                                                aaGCACTTTACGTTCTCGTTCACGTCTAATTTCCTCAACAATTTGTGAAGCATAAATATTTTCAATTTTGATAGTATTATCAACTTTTGAATAATTAGGc                         <  1:1039332/100‑1 (MQ=255)
                                                                aaGCACTTTACGTTCTCGTTCACGTCTAATTTCCTCAACAATTTGTGAAGCATAAATATTTTCAATTTTGATAGTATTATCAACTTTTGAATAATTAGGc                         <  1:1127248/100‑1 (MQ=255)
                                                                aaGCACTTTAAGTTCTCGTTCACGTCTAATTTCCTTAACAATTTGTGAAGCATAAATATTTTCAATTTTGATAGTATTATCAACTTTTGAATAATTAGGc                         <  1:4243154/100‑1 (MQ=255)
                                                                 aGCACTTTACGTTCTCGTTCACGTCTAATTTCCTCAACAATTTGTGAAGCATAAATATTTTCAATTTTGATAGTATTATCAACTTTTGAATAATTAgggg                        <  1:4035485/100‑3 (MQ=255)
                                                                 aGCACTTTACGTTCTCGTTCACGTCTAATTTCCTCAACAATTTGTGAAGCATAAATATTTTCAATTTTGATAGTATTATCAACTTTTGAATAATTAGGCa                        <  1:2253110/100‑1 (MQ=255)
                                                                  gCACTTTACGTTCTCGTTCACGTCTAATTTCCTCAACAATTTGTGAAGCATAAATATTTTCAATTTTGATAGTATTATCAACTTTTGAATAATTAGGCAt                       <  1:4526556/100‑1 (MQ=255)
                                                                   cACTTTACGTTCTCGTTCACGTCTAATTTCCTCAACAATTTGTGAAGCATAAATATTTTCAATTTTGATAGTATTATCAACTTTTGAATAATTAGGCAtt                      <  1:3900036/100‑1 (MQ=255)
                                                                   cACTTTACGTTCTCGTTCACGTCTAATTTCCTCAACAATTTGTGAAGCATAAATATTTTCAATTTTGATAGTATTATCAACTTTTGAATAATTAGGCAtt                      <  1:5415017/100‑1 (MQ=255)
                                                                   cACTTTACGTTCTCGTTCACGTCTAATTTCCTCAACAATTTGTGAAGCATAAATATTTTCAATTTTGATAGTATTATCAACTTTTGAATAATTAGGCAtt                      <  1:2586602/100‑1 (MQ=255)
                                                                   cACTTTACGTTCTCGTTCACGTCTAATTTCCTCAACAATTTGTGAAGCATAAATATTTTCAATTTTGATAGTATTATCAACTTTTGAATAATTAGGCAtt                      <  1:5077127/100‑1 (MQ=255)
                                                                   cACTTTACGTTCTCGTTCACGTCTAATTTCCTCAACAATTTGTGAAGCATAAATATTTTCAATTTTGATAGTATTATCAACTTTTGAATAATTAGGCAtt                      <  1:2234257/100‑1 (MQ=255)
                                                                   cACTTTACGTTCTCGTTCACGTCTAATTTCCTCAACAATTTGTGAAGCATAAATATTTTCAATTTTGATAGTATTATCAACTTTTGAATAATTAGGCAtt                      <  1:2450560/100‑1 (MQ=255)
                                                                   cACTTTACGTTCTCGTTCACGTCTAATTTCCTCAACAATTTGTGAAGCATAAATATTTTCAATTTTGATAGTATTATCAACTTTTGAATAATTAGGCAtt                      <  1:4810882/100‑1 (MQ=255)
                                                                   cACTTTACGTTCTCGTTCACGTCTAATTTCCTCAACAATTTGTGAAGCATAAATATTTTCAATTTTGATAGTATTATCAACTTTTGAATAATTAGGCAtt                      <  1:6482782/100‑1 (MQ=255)
                                                                   cACTTTACGTTCTCGTTCACGTCTAATTTCCTCAACAATTTGTGAAGCATAAATATTTTCAATTTTGATAGTATTATCAACTTTTGAATAATTAGGCAtg                      <  1:6234141/100‑2 (MQ=255)
                                                                   cACTTTACGTTCTCGTTCACGTCTAATTTCCTCAACAATTTGTGAAGCATAAATATTTTCAATTTTGATAGTATTATCAACTTTTGAATAATTAGGCAtg                      <  1:4710141/100‑2 (MQ=255)
                                                                   cACTTTACGTTCTCGTTCACGTCTAATTTCCTCAACAATTTGTGAAGCATAAATATTTTCAATTTTGATAGTATTATCAACTTTTGAATAAGTAGGCAtt                      <  1:611567/100‑1 (MQ=255)
                                                                   cACTTTACGTTCTCGTTCACGACTAATTTCCTCAACAATTTGTGAAGCATAAATATTTTCAATTTTGATAGTATTATCAACTTTTGAATAATTAGGCAtt                      <  1:1627208/100‑1 (MQ=255)
                                                                   cACTTTACATTCTCGTTCACGTCTAATTTCCTCAACAATTTGTGAAGAATAAATATTTTCAATTTTGATAGTATTATCAACTTTTGAATAATTAGGCAtg                      <  1:4875291/100‑2 (MQ=255)
                                                                    gcTTTACGTTCTCGTTCACGTCTAATTTCCTCAGCAAATTGTGAAGCATAAATATTTTCAATTTTGATAGTATTATCAACATTTGAATAATTAGGCATTg                     <  1:5171825/99‑1 (MQ=255)
                                                                    aCTTTACGTTCTCGTTCACGTCTAATTTCCTCAACAATTTGTGAAGCATAAATATTTTCAATTTTGATAGTATTATCAACTTTTGAATAATTAGGCAtta                     <  1:3098238/100‑2 (MQ=255)
                                                                    aCTTTACGTTCTCGTTCACGTCTAATTTCCTCAACAATTTGTGAAGCATAAATATTTTCAATTTTGATAGTATTATCAACTTTTGAATAATTAGGCAtgc                     <  1:5369745/100‑3 (MQ=255)
                                                                    aCTTTACGTTCTCGTTCACGTCTAATTTCCTCAACAATTTGTGAAGCATAAATATTTTCAATTTTGATAGTATTATCAACTTTTGAATAATTAGGCAtgc                     <  1:5543086/100‑3 (MQ=255)
                                                                    aCTTTACGTTCTCGTTCACGTCTAATTTCCTCAACAATTTGTGAAGCATAAATATTTTCAATTTTGATAGTATTATCAACTTTTGAATAATTAGGCATTg                     <  1:3897156/100‑1 (MQ=255)
                                                                    aCTTTACGTTCTCGTTCACGTCTAATTTCCTCAACAATTTGTGAAGCATAAATATTTTCAATTTTGATAGTATTATCAACTTTTGAATAATTAGGCATTg                     <  1:436144/100‑1 (MQ=255)
                                                                    aCTTTACGTTCTCGTTCACGTCTAATTTCCTCAACAATTTGTGAAGCATAAATATTTTCAATTTTGATAGTATTATCAACTTTTGAATAATTAGGCATTg                     <  1:458969/100‑1 (MQ=255)
                                                                    aCTTTACGTTCTCGTTCACGTCTAATTTCCTCAACAATTTGTGAAGCATAAATATTTTCAATTTTGATAGTATTATCAACTTTTGAATAATTAGGCATTg                     <  1:4642944/100‑1 (MQ=255)
                                                                    aCTTTACGTTCTCGTTCACGTCTAATTTCCTCAACAATTTGTGAAGCATAAATATTTTCAATTTTGATAGTATTATCAACTTTTGAATAATTAGGCATTg                     <  1:4088359/100‑1 (MQ=255)
                                                                    aCTTTACGTTCTCGTTCACGTCTAATTTCCTCAACAATTTGTGAAGCATAAATATTTTCAATTTTGATAGTATTATCAACTTTTGAATAATTAGGCATTg                     <  1:3904638/100‑1 (MQ=255)
                                                                    aCTTTACGTTCTCGTTCACGTCTAATTTCCTCAACAATTTGTGAAGCATAAATATTTTCAATTTTGATAGTATTATCAACTTTTGAATAATTAGGCATTg                     <  1:6380844/100‑1 (MQ=255)
                                                                    aCTTTACGTTCTCGTTCACGTCTAATTTCCTCAACAATTTGTGAAGCATAAATATTTTCAATTTTGATAGTATTATCAACTTTTGAATAATTAGGCATTg                     <  1:14248/100‑1 (MQ=255)
                                                                    aCTTTACGTTCTCGTTCACGTCTAATTTCCTCAACAATTTGTGAAGCATAAATATTTTCAATTTTGATAGTATTATCAACTTTTGAATAATTAGGCATTg                     >  1:269629/1‑100 (MQ=255)
                                                                    aCTTTACGTTCTCGTTCACGTCTAATTTCCTCAACAATTTGTGAAGCATAAATATTTTCAATTTTGATAGTATTATCAACTTTTGAATAATTAGGCATTg                     <  1:2825771/100‑1 (MQ=255)
                                                                    aCTTTACGTTCTCGTTCACGTCTAATTTCCTCAACAATTTGTGAAGCATAAATATTTTCAATTTTGATAGTATTATCAACTTTTGAATAATTAGGCATTg                     >  1:2925449/1‑100 (MQ=255)
                                                                    aCTTTACGTTCTCGTTCACGTCTAATTTCCTCAACAATTTGTGAAGCATAAATATTTTCAATTTTGATAGTATTATCAACTTTTGAATAATTAGGCATTg                     <  1:3652282/100‑1 (MQ=255)
                                                                    aCTTTACGTTCTCGTTCACGTCTAATTTCCTCAACAATTTGTGAAGCATAAATATTTTCAATTTTGATAGTATTATCAACTTTTGAATAATTAGGCATTg                     <  1:555787/100‑1 (MQ=255)
                                                                    aCTTTACGTTCTCGTTCACGTCTAACTTCCTCAACAATTTGTGAAGCATAAATATTTTCAATTTTGATAGTATTATCAACTTTTGAATAATTAGGCATTg                     <  1:3573708/100‑1 (MQ=255)
                                                                     cTTTACGTTCTCGTTCACGTCTAATTTCCTCAACAATTTGTGAAGCATAAATATTTTCAATTTTGATAGTATTATCAACTTTTGAATAATTAGGCATTgt                    <  1:5503970/100‑2 (MQ=255)
                                                                              ctcGTTCACGTCTAATTTCCTCAACAATTTGTGAAGCATAAATATTTTCAATTTTGATAGTATTATCAACTTTTGAATAATTAGGCATTGATTttgttgt           <  1:518610/100‑1 (MQ=255)
                                                                                cGTTCACGTCTAATTTCCTCAACTATTTGTGAAGCATAAATATTTTCAATTTTGATAGTATTATCAACTTTTGAATAATTAGGCATTGATTttgttgttg         >  1:6035069/1‑100 (MQ=255)
                                                                                cGTTCACGTCTAATTTCCTCAACAATTTGTGAAGCATAAATATTTTCAATTTTGATAGTATTATCAACTTTTGAATAATTAGGCATTGATTttgttgttg         <  1:636034/100‑1 (MQ=255)
                                                                                       gTCTAATTTCCTCAACAATTTGTGAAGCATAAATATTTTCAATTTTGATAGTATTATCAACTTTTGAATAATTAGGCATTGATTTTGTTGTTGATGATAc  <  1:1987937/100‑1 (MQ=255)
                                                                                       gTCTAATTTCCTCAACAATTTGTGAAGCATAAATATTTTCAATTTTGATAGTATTATCAACTTTTGAATAATTAGGCATTGATTTTGTTGTTGATGATAc  <  1:4448431/100‑1 (MQ=255)
                                                                                                 |                                                                                         
CTTGCTCTTCGTTTTTATGTTCTTCACGCTTTTGTTGCAACGCTTTTTTAAATCGACGCTTTTGAAGCACTTTACGTTCTCGTTCACGTCTAATTTCTTCAACAATTTGTGAAGCATAAATATTTTCAATTTTGATAGTATTATCAACTTTTGAATAATTAGGCATTGATTTTGTTGTTGATGATAC  >  NC_007793/1862466‑1862652

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