Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_007793 1,631,245 A→G F17F (TTT→TTC)  SAUSA300_RS08060 ← integrase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0077931,631,2450AG100.0% 235.6 / NA 89F17F (TTT→TTC) SAUSA300_RS08060integrase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (0/89);  total (0/89)

TTTGCATTTTCAATAGAACCACGTTGATATGGACAACCCAAATCACTAAAATAAATACTACAATTCAATTCATCTTCTATTTCTTCAAATAGAGCAAACTCCTTACCTCTATCCATAGTTATACTTTTGGGTAATAATGGTTTCATAATATTTATAATCTCATTACAAACTACATTGGATTTACGTGTATT  >  NC_007793/1631150‑1631340
                                                                                               |                                                                                               
tttACATTTTCAATGGAACCACGTTGATATGGACAACCCAAATCACTAAAATAAATACTACAATTCAATTCATCTTCTATTTGTTCAAATAGAGCGAAca                                                                                             <  1:85563/100‑2 (MQ=255)
    cATTTTCAATGGAATCACGTTGATATGGACAACCCAAATCACTAAAATAAATACTACAATTCAATTCATCTTCTATTTGTTCAAATAGAGCGAACTCCtt                                                                                         <  1:4798546/100‑1 (MQ=255)
    cATTTTCAATGGAACCACGTTGATATGGACAACCCAAATCACTAAAATAAATACTACAATTCAATTCATCTTCTATTTGTTCAAATAGAGCGAACTCCtt                                                                                         <  1:6146213/100‑1 (MQ=255)
     aTTTTCAATGGAACTACGTTGATATGGACAACCCAAATCACTAAAATAAATACTACAATTCAATTCATCTTCTATTTGTTCAAATAGAGCGAACTCCTTa                                                                                        <  1:4937813/100‑1 (MQ=255)
      ttttCAATGGAACCACGTTGATATGGACAACCCAAATCACTAAAATAAATACTACAATTCAATTCATCTTCTATTTGTTCAAATAGAGCGAACTCCTTAc                                                                                       <  1:3575429/100‑1 (MQ=255)
       tttCAATGGAACCACGTTGATATGGACAACCCAAATCACTAAAATAAATACTACAATTCAATTCATCTTCTATTTGTTCAAATAGAGCGAACTCCTTAcc                                                                                      <  1:3317466/100‑1 (MQ=255)
        ttCAATGGAACCACGTTGATATGGACAACCCAAATCACTAAAATAAATACTACAATTCAATTCATCTTCTATTTGTTCAAATAGAGCGAACTCCTTAccc                                                                                     <  1:623127/100‑2 (MQ=255)
        ttCAATGGAACCACGTTGATATGGACAACCCAAATCACTAAAATAAATACTACAATTCAATTCATCTTCTATTTGTTCAAATAGAGCGAACTCCTTACct                                                                                     <  1:5300910/100‑1 (MQ=255)
        ttCAATGGAACCACGTTGATATGGACAACCCAAATCACTAAAATAAATACTACAATTCAATTCATCTTCTATTTGTTCAAATAGAGCGAACTCCTTACct                                                                                     <  1:843346/100‑1 (MQ=255)
         tCAATGGAACCACGTTGATATGGACAACCCAAATCACTAAAATAAATACTACAATTCAATTCATCTTCTATTTGTTCAAATAGAGCGAACTCCTTAccag                                                                                    <  1:3918626/100‑3 (MQ=255)
         tCAATGGAACCACGTTGATATGGACAACCCAAATCACTAAAATAAATACTACAATTCAATTCATCTTCTATTTGTTCAAATAGAGCGAACTCCTTACctg                                                                                    <  1:999421/100‑2 (MQ=255)
         tCAATGGAACCACGTTGATATGGACAACCCAAATCACTAAAATAAATACTACAATTCAATTCATCTTCTATTTGTTCAAATAGAGCGAACTCCTTACctc                                                                                    <  1:1656270/100‑1 (MQ=255)
         tCAATGGAACCACGTTGATATGGACAACCCAAATCACTAAAATAAATACTACAATTCAATTCATCTTCTATTTGTTCAAATAGAGCGAACTCCTTACctc                                                                                    <  1:5725119/100‑1 (MQ=255)
         tCAATGGAACCACGTTGATATGGACAACCCAAATCACTAAAATAAATACTACAATTCAATTCATCTTCTATTTGTTCAAATAGAGCGAACTCCTTACcta                                                                                    <  1:4588255/100‑2 (MQ=255)
          cAATGGAACCACGTTGATATGGACAACCCAAATCACTAAAATAAATACTACAATTCAATTCATCTTCTATTTGTTCAAATAGAGCGAACTCCTTAccatg                                                                                   <  1:5361800/100‑4 (MQ=255)
          cAATGGAACCACGTTGATATGGACAACCCAAATCACTAAAATAAATACTACAATTCAATTCATCTTCTATTTGTTCAAATAGAGCGAACTCCTTACctct                                                                                   <  1:3880056/100‑1 (MQ=255)
          cAATGGAACCACGTTGATATGGACAACCCAAATCACTAAAATAAATACTACAATTCAATTCATCTTCTATTTGTTCAAATAGAGCGAACTCCTTACctct                                                                                   <  1:3658229/100‑1 (MQ=255)
          cAATGGAACCACGTTGATATGGACAACCCAAATCACTAAAATAAATACTACAATTCAATTCATCTTCTATTTGTTCAAATAGAGCGAACTCCTTACctct                                                                                   <  1:773519/100‑1 (MQ=255)
          cAATGGAACCACGTTGATATGGACAACCCAAATCACTAAAATAAATACTACAATTCAATTCATCTTCTATTTGTTCAAATAGAGCGAACTCCTTACctct                                                                                   <  1:1272068/100‑1 (MQ=255)
           aaTGGAACCACGTTGATATGGACAACCCAAATCACTAAAATAAATACTACAATTCAATTCATCTTCTATTTGTTCAAATAGAGCGAACTCCTTAcccaag                                                                                  <  1:5306849/100‑5 (MQ=255)
           aaTGGAACCACGTTGATATGGACAACCCAAATCACTAAAATAAATACTACAATTCAATTCATCTTCTATTTGTTCAAATAGAGCGAACTCCTTAccagct                                                                                  <  1:4436747/100‑5 (MQ=255)
            atggAACCACGTTGATATGGACAACCCAAATCTCTACAATAAATACTACAATTCAATTCATCTTCTATTTGTTCAAATAGAGCGAACTCCTTACCTCTag                                                                                 <  1:6062216/97‑2 (MQ=255)
            atggAACCACGTTGATATGGACAACCCAAATCACTAAAATAAATACTACAATTCAATTCATCTTCTATTTGTTCAAATAGAGCGAACTCCTTAcccgcag                                                                                 <  1:813573/97‑6 (MQ=255)
            atggAACCACGTTGATATGGACAACCCAAATCACTAAAATAAATACTACAATTCAATTCATCTTCTATTTGTTCAAATAGAGCGAACTCCTTACCTCTac                                                                                 <  1:1280461/97‑2 (MQ=255)
            atggAACCACGTTGATATGGACAACCCAAATCACTAAAATAAATACTACAATTCAATTCATCTTCTATTTGTTCAAATAGAGCGAACTCCTTACCTCTAt                                                                                 <  1:2663916/97‑1 (MQ=255)
             tggAACCACGTTGATATGGACAACCCAAATCACTAAAATAAATACTACAATTCAATTCATCTTCTATTTGTTCAAATAGAGCGAACTCCTTACCTCTaag                                                                                <  1:3687912/98‑3 (MQ=255)
             tggAACCACGTTGATATGGACAACCCAAATCACTAAAATAAATACTACAATTCAATTCATCTTCTATTTGTTCAAATAGAGCGAACTCCTTACCTCTAtg                                                                                <  1:1362690/98‑2 (MQ=255)
              ggAACCACGTTGATATGGACACCCCAAATCCCTAAAATAAATACTACAATTCAATTCATCTTCTATTTGTTCAAATAGAGCGAACTCCTTCCCTCTagaa                                                                               <  1:4696689/99‑4 (MQ=255)
              ggAACCACGTTGATATGGACAACCCAAATCACTAAAATAAATACTACAATTCAATTCATCTTCTATTTGTTCAAATAGAGCGAACTCCTTACCTCTATcc                                                                               <  1:4487191/99‑1 (MQ=255)
              ggAACCACGTTGATATGGACAACCCAAATCACTAAAATAAATACTACAATTCAATTCATCTTCTATTTGTTCAAATAGAGCGAACTCCTTACCTCTATcc                                                                               <  1:4313650/99‑1 (MQ=255)
              ggAACCACGTTGATATGGACAACCCAAATCACTAAAATAAATACTACAATTCAATTCATCTTCTATTTGTTCAAATAGAGCGAACTCCTTACCTCTATcc                                                                               <  1:4495080/99‑1 (MQ=255)
              ggAACCACGTTGATATGGACAACCCAAATCACTAAAATAAATACTACAATTCAATTCATCTTCTATTTGTTCAAATAGAGCGAACTCCTTACCTCTATcc                                                                               <  1:1192364/99‑1 (MQ=255)
              ggAACCACGTTGATATGGACAACCCAAATCACTAAAATAAATACTACAATTCAATTCATCTTCTATTTGTTCAAATAGAGCGAACTCCCTAccagcacac                                                                               <  1:3156012/99‑8 (MQ=255)
                aaCCACGTTGATATGGACAACCCAAATCGCTAAAATAAATACTACAATTCAATTCATCTTCTATGTGTTCAAATAGAGCGAACTCCTTACCTCTATCCAt                                                                             <  1:3668945/100‑1 (MQ=255)
                aaCCACGTTGATATGGACAACCCAAATCACTAAAATAAGTACTACAATTCAATTCATCTTCTATTTGTTCAAATAGAGCGAACTCCTTACCTCTATCCAt                                                                             <  1:1071928/100‑1 (MQ=255)
                aaCCACGTTGATATGGACAACCCAAATCACTAAAATAAATACTACAATTCAATTCATCTTCTATTTGTTCAAATAGAGCGAACTCCTTACCTCTATccct                                                                             <  1:1709010/100‑3 (MQ=255)
                aaCCACGTTGATATGGACAACCCAAATCACTAAAATAAATACTACAATTCAATTCATCTTCTATTTGTTCAAATAGAGCGAACTCCTTACCTCTATCCag                                                                             <  1:5084744/100‑2 (MQ=255)
                aaCCACGTTGATATGGACAACCCAAATCACTAAAATAAATACTACAATTCAATTCATCTTCTATTTGTTCAAATAGAGCGAACTCCTTACCTCTATCCaa                                                                             <  1:5400412/100‑2 (MQ=255)
                aaCCACGTTGATATGGACAACCCAAATCACTAAAATAAATACTACAATTCAATTCATCTTCTATTTGTTCAAATAGAGCGAACTCCTTACCTCTATCCAt                                                                             <  1:1176976/100‑1 (MQ=255)
                aaCCACGTTGATATGGACAACCCAAATCACTAAAATAAATACTACAATTCAATTCATCTTCTATTTGTTCAAATAGAGCGAACTCCTTACCTCTATCCAt                                                                             <  1:1369859/100‑1 (MQ=255)
                aaCCACGTTGATATGGACAACCCAAATCACTAAAATAAATACTACAATTCAATTCATCTTCTATTTGTTCAAATAGAGCGAACTCCTTACCTCTATCCAt                                                                             <  1:6374591/100‑1 (MQ=255)
                aaCCACGTTGATATGGACAACCCAAATCACTAAAATAAATACTACAATTCAATTCATCTTCTATTTGTTCAAATAGAGCGAACTCCTTACCTCTATCCAt                                                                             <  1:6368693/100‑1 (MQ=255)
                aaCCACGTTGATATGGACAACCCAAATCACTAAAATAAATACTACAATTCAATTCATCTTCTATTTGTTCAAATAGAGCGAACTCCTTACCTCTATCCAt                                                                             <  1:4659130/100‑1 (MQ=255)
                aaCCACGTTGATATGGACAACCCAAATCACTAAAATAAATACTACAATTCAATTCATCTTCTATTTGTTCAAATAGAGCGAACTCCTTACCTCTATCCAt                                                                             <  1:4273966/100‑1 (MQ=255)
                aaCCACGTTGATATGGACAACCCAAATCACTAAAATAAATACTACAATGCAATTCATCTTCTATTTGTTCAAATAGAGCGAACTCCTTACCTCTATCCAt                                                                             <  1:6112561/100‑1 (MQ=255)
                 aCCACGTTGATATGGACAACCCAAATCACTAAAATAAATACTACAATTCAATTCATCTTCTATTTGTTCAAATAGAGCGAACTCCTTACCTCTATCCAtg                                                                            <  1:307105/100‑2 (MQ=255)
                 aCCACGTTGATATGGACAACCCAAATCACTAAAATAAATACTACAATTCAATTCATCTTCTATTTGTTCAAATAGAGCGAACTCCTTACCTCTATCCAtc                                                                            <  1:5989745/100‑2 (MQ=255)
                 aCCACGTTGATATGGACAACCCAAATCACTAAAATAAATACTACAATTCAATTCATCTTCTATTTGTTCAAATAGAGCGAACTCCTTACCTCTATCCAtc                                                                            <  1:1640175/100‑2 (MQ=255)
                 aCCACGTTGATATGGACAACCCAAATCACTAAAATAAATACTACAATTCAATTCATCTTCTATTTGTTCAAATAGAGCGAACTCCTTACCTCTATCCAtc                                                                            <  1:2782363/100‑2 (MQ=255)
                 aCCACGTTGATATGGACAACCCAAATCACTAAAATAAATACTACAATTCAATTCATCTTCTATTTGTTCAAATAGAGCGAACTCCTTACCTCTATCCAtc                                                                            <  1:2649364/100‑2 (MQ=255)
                 aCCACGTTGATATGGACAACCCAAATCACTAAAATAAATACTACAATTCAATTCATCTTCTATTTGTTCAAATAGAGCGAACTCCTTACCTCTATCCATa                                                                            <  1:2566684/100‑1 (MQ=255)
                  ccACGTTGATATGGACAACCCAAATCACTAAAATAAATACTACAATTCAATTCTCCTTCTAGTTGTTCAAATAGAGCGAACTCCTTACCTCTATCCAtgt                                                                           <  1:5870505/100‑3 (MQ=255)
                  ccACGTTGATATGGACAACCCAAATCACTAAAATAAATACTACAATTCAATTCATCTTCTGTTTGTTCAAATAGAGCGAACTCCTTACCTCTATCCAtca                                                                           <  1:3600532/100‑3 (MQ=255)
                  ccACGTTGATATGGACAACCCAAATCACTAAAATAAATACTACAATTCAATTCATCTTCTATTTGTTCAAATAGAGCGAACTCCTTGCCTCTATCCAtgc                                                                           <  1:4550021/100‑3 (MQ=255)
                  ccACGTTGATATGGACAACCCAAATCACTAAAATAAATACTACAATTCAATTCATCTTCTATTTGTTCAAATAGAGCGAACTCCTTACCTCTATCCagct                                                                           <  1:2574585/100‑4 (MQ=255)
                  ccACGTTGATATGGACAACCCAAATCACTAAAATAAATACTACAATTCAATTCATCTTCTATTTGTTCAAATAGAGCGAACTCCTTACCTCTATCCAtgg                                                                           <  1:4366307/100‑3 (MQ=255)
                  ccACGTTGATATGGACAACCCAAATCACTAAAATAAATACTACAATTCAATTCATCTTCTATTTGTTCAAATAGAGCGAACTCCTTACCTCTATCCATAg                                                                           <  1:3275145/100‑1 (MQ=255)
                   cACGTTGATATGGGCAACCCAAATCACTAAAATAAATACTACAATTCAATTCATCTTCTATTTGTTCAAATAGAGCGAACTCCTTACCTCTATCCATAGt                                                                          <  1:2704057/100‑1 (MQ=255)
                   cACGTTGATATGGACAACCCAAATCACTAAAATAAATACTACAATTCAATTCATCTTCTATTTGTTCAAATAGAGCGAACTCCTTACCTCTATCCAtgat                                                                          <  1:86642/100‑4 (MQ=255)
                   cACGTTGATATGGACAACCCAAATCACTAAAATAAATACTACAATTCAATTCATCTTCTATTTGTTCAAATAGAGCGAACTCCTTACCTCTATCCATAGt                                                                          <  1:2129793/100‑1 (MQ=255)
                    aCGTTGATATGGACAACCCAAATCACTAAAATAAATACTACAATTCAATTCATCTTCTATTTGTTCAAATAGAGCGAACTCCTTACCTCTATCCATAGtt                                                                         <  1:944290/100‑1 (MQ=255)
                    aCGTTGATATGGACAACCCAAATCACTAAAATAAATACTACAATTCAATTCATCTTCTATTTGTTCAAATAGAGCGAACTCCTTACCTCTATCCATAGtt                                                                         <  1:4873377/100‑1 (MQ=255)
                    aCGTTGATATGGACAACCCAAATCACTAAAATAAATACTACAATTCAATTCATCTTCTATTTGTTCAAATAGAGCGAACTCCTTACCTCTATCCATAGtt                                                                         <  1:5888045/100‑1 (MQ=255)
                    aCGTTGATATGGACAACCCAAATCACTAAAATAAATACTACAATTCAATTCATCTTCTATTTGTTCAAATAGAGCGAACTCCTTACCTCTATCCATAGtt                                                                         <  1:6125835/100‑1 (MQ=255)
                    aCGTTGATATGGACAACCCAAATCACTAAAATAAATACTACAATTCAATTCATCTTCTATTTGTTCAAATAGAGCGAACTCCTTACCTCTATCCATAGtt                                                                         <  1:5741061/100‑1 (MQ=255)
                    aCGTTGATATGGACAACCCAAATCACTAAAATAAATACTACAATTCAATTCATCTTCTATTTGTTCAAATAGAGCGAACTCCTTACCTCTATCCATAGtt                                                                         <  1:141615/100‑1 (MQ=255)
                      gTTGATATGGACAACCCAAATCACTAAAATAAATACTACAATTCAATTCATCATCTATTTGTTCAAATAGAGCGAACTCCTTACCTCTATCCAtgggcac                                                                       <  1:3582260/100‑7 (MQ=255)
                        tGATATGGACAACCCAAATCACTAAAATAAATACTACAATTCAATTCATCTTCTATTTGTTCAAATAGAGCGAACTCCTTGCCTCTATCCATAGTtatgt                                                                     <  1:4757977/100‑3 (MQ=255)
                        tGATATGGACAACCCAAATCACTAAAATAAATACTACAATTCAATTCATCTTCTATTTGTTCAAATAGAGCGAACTCCTTACCTCTATCCATAGTtagat                                                                     <  1:5272517/100‑4 (MQ=255)
                         gATATGGACAACCCAAATCACTAAAATAAATACTACAATTCAATTCATCTTCTATTTGTTCAAATAGAGCGAACTCCTTACCTCTATCCATAGTTATGCt                                                                    <  1:5382323/100‑1 (MQ=255)
                         gATATGGACAACCCAAATCACTAAAATAAATACTACAATTCAATTCATCTTCTATTTGTTCAAATAGAGCGAACTCCTTACCTCTATCCATAGTTATGCt                                                                    <  1:4567181/100‑1 (MQ=255)
                          atatGGACAACCCAAATCACTAAAATAAATACTACAATTCAATTCATCTTCTATTTGTTCAAATAGAGCGAACTCCTTACCTCTATCCATAGTTATGCtt                                                                   <  1:5789476/100‑1 (MQ=255)
                          atatGGACAACCCAAATCACTAAAATAAATACTACAATTCAATTCATCTTCTATTTGTTCAAATAGAGCGAACTCCTTACCTCTATCCATAGTTATGCtt                                                                   <  1:2832346/100‑1 (MQ=255)
                          atatGGACAACCCAAATCACTAAAATAAATACTACAATTCAATTCATCTTCTATTTGTTCAAATAGAGCGAACTCCTTACCTCTATCCATAGTTATGCtt                                                                   <  1:6331595/100‑1 (MQ=255)
                           tatGGACAACCCAAATCACTAAAATAAATACTACAATTCAATTCATCTTCTATTTGTTCAAATAGAGCGAACTCCTTACCTCTATCCATAGTTATGCttt                                                                  <  1:1663822/100‑1 (MQ=255)
                           tatGGACAACCCAAATCACTAAAATAAATACTACAATTCAATTCATCTTCTATTTGTTCAAATAGAGCGAACTCCTTACCTCTATCCATAGTTATGCttg                                                                  <  1:4442829/100‑2 (MQ=255)
                            atGGACAACCCAAATCGCTAAAATAAATACTACAATTCAATTCATCTTCTATTTGTTCAAATAGAGCGAACTCCTTACCTCTATCCATAGTTATGCtttg                                                                 <  1:994071/100‑2 (MQ=255)
                                         aTCACTAAAATAAATACTACAATTCAATTCATCTTCTATTTGTTCAAATAGAGCGAACTCCTTACCTCTATCCATAGTTATGCTTTTTGGTAATAATGGt                                                    <  1:5810732/100‑1 (MQ=255)
                                           cACTAAAATAAATACTACAATTCAATTCATCTTCTATTTGTTCAAATAGAGCGAACTCCTTGCCTCTATCCATAGTTATGCTTTTTGGTAATAATGGttt                                                  <  1:6484515/100‑1 (MQ=255)
                                                 aataaataCTACAATTCAATTCATCTTCTATTTGTTCAAATAGAGCGAACTCCTTACCTCTATCCATAGTTATGCTTTTTGGTAATAATGGTTTCataat                                            <  1:2311478/100‑1 (MQ=255)
                                                                   attcaTCTTCTATTTGTTCAAATAGAGCGAACTCCTTACCTCTATCCATAGTTATGCTTTTTGGTAATAATGGTTTCATAATATTTATAATCTCATTACa                          <  1:6299408/100‑1 (MQ=255)
                                                                   attcaGCTTCTATTTGTGCAAATAGAGCGAACTCCTTTCCGCTACCCATATTTATGCTTTGTGGTAAGAATGGTTTCATAATATTTATAATCTCATTACa                          <  1:3063454/100‑1 (MQ=255)
                                                                       aTCTTCTATTTGTTCAAATAGAGCGAACTCCTTACCTTTATCCATAGTTATGCTTTTTGGTAATAATGGTTTCATAATATTTATAATCTCATTACAAACt                      <  1:363823/100‑1 (MQ=255)
                                                                           tctATTTGTTCAAATAGAGCGAACTCCTTACCTCTATCCATAGTTATGCTTTTTGGTAATAATGGTTTCATAATATTTATAATCTCATTACAAACTACag                  <  1:3705594/100‑2 (MQ=255)
                                                                            ctATTTGTTCAAATAGAGCGAACTCCTTACCTCTATCCATAGTTATGCTTTTTGGTAATAATGGTTTCATAATATTTATAATCTCATTACAAACTACAtt                 <  1:5732396/100‑1 (MQ=255)
                                                                              aTTTGTTCAAATAGAGCGAACTCCTTACCTCTATCCATAGTTATGCTTTTTGGTAATAATGGTTTCATAATATTTATAATCTCATTACAAACTACATTgg               <  1:3757010/100‑1 (MQ=255)
                                                                                        aTAGAGCGAACTCCTTACCTCTATCCATAGTTATGCTTTTTGGTAATAATGGTTTCATAATATTTATAATCTCATTACAAACTACATTGGATTTACgtgt     <  1:3023433/100‑1 (MQ=255)
                                                                                        aTAGAGCGAACTCCCTACCTCTATCCATAGTTATGCTTTTTGGTAATAATGGTTTCATAATATTTATAATCTCATTACAAACTACATTGGATTTACgtgg     <  1:4991132/100‑2 (MQ=255)
                                                                                           gagCGAACTCCTTACCTCTATCCATAGTTATGCTTTTTGGTAATAATGGTTTCATAATATTTATAATCTCATTACAAACTACATTGGATTTACGTGTAtt  <  1:6476554/100‑1 (MQ=255)
                                                                                               |                                                                                               
TTTGCATTTTCAATAGAACCACGTTGATATGGACAACCCAAATCACTAAAATAAATACTACAATTCAATTCATCTTCTATTTCTTCAAATAGAGCAAACTCCTTACCTCTATCCATAGTTATACTTTTGGGTAATAATGGTTTCATAATATTTATAATCTCATTACAAACTACATTGGATTTACGTGTATT  >  NC_007793/1631150‑1631340

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: ‑