Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,519,014 G→A pseudogene (294/504 nt) insB1 → IS1 transposase B,IS, phage, Tn, Transposon‑related functions, extrachromosomal, transposon related

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,519,0140GA100.0% 210.9 / NA 74pseudogene (294/504 nt)insB1IS1 transposase B,IS, phage, Tn, Transposon‑related functions, extrachromosomal, transposon related
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (37/37);  total (37/37)

TGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAGATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATAATCTGAATCTGAGACAACATCTGGCAAGGCTGGGACG  >  NC_000913/4518892‑4519142
                                                                                                                          |                                                                                                                                
tggctgTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGAtggctggc                                                                                                                  <  2:1357753/139‑1 (MQ=255)
  gctgTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCg                                                                                                                >  1:1138170/1‑139 (MQ=255)
            gcgTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGcc                                                                                                                 >  2:312016/1‑128 (MQ=255)
            gcgTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGcc                                                                                                                 <  1:312016/128‑1 (MQ=255)
            gcgTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAAt                                                                                                      >  2:679905/1‑139 (MQ=255)
                aTGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACg                                                                                                  <  2:105715/139‑1 (MQ=255)
                  gACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGcc                                                                                                >  2:229217/1‑139 (MQ=255)
                      ggATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGaa                                                                                            <  2:1147268/139‑1 (MQ=255)
                       gctccGGAAGACGGTTGTTGCGCACGTATTCGGTGAACGCACTATGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGTGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAgg                                                                                          >  2:1180577/5‑139 (MQ=11)
                         tccGGAAGACGGTTGTTGCGCACGTATTCGGTGAACGCACTATGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGTGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGa                                                                                        <  1:1255119/137‑1 (MQ=11)
                           cGGAAGACGGTTGTTGCGCACGTATTCGGTGAACGCACTATGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGTGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGaaa                                                                                      >  2:1039678/1‑139 (MQ=11)
                            ggaggaCGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGaaa                                                                                      >  1:4908/1‑139 (MQ=255)
                                   ggTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCAc                                                                               >  1:51169/1‑139 (MQ=255)
                                      tgtgGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGcc                                                                                                                 <  1:217181/102‑1 (MQ=255)
                                      tgtgGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGcc                                                                                                                 >  2:217181/1‑102 (MQ=255)
                                        ttgcgcACGTATTCGGTGAACGCACTATGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGTGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGcc                                                                                                                 <  1:773736/97‑1 (MQ=14)
                                        ttgcgcACGTATTCGGTGAACGCACTATGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGTGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGcc                                                                                                                 >  2:773736/3‑99 (MQ=14)
                                        tgGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTAt                                                                          >  2:753983/1‑139 (MQ=255)
                                          gcgcACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGAtggct                                                                                                                     >  1:729738/1‑94 (MQ=255)
                                          gcgcACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGAtggct                                                                                                                     <  2:729738/94‑1 (MQ=255)
                                            gcACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGc                                                                      <  2:4908/139‑1 (MQ=255)
                                            gcACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGc                                                                      <  1:337676/139‑1 (MQ=255)
                                            gcACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGc                                                                      <  2:825314/139‑1 (MQ=255)
                                               cgtattCGGTGAACGCACTATGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGTGGCCTTTGAGGTCGTGGTATAAATGACGGAtggct                                                                                                                     >  1:1408193/5‑88 (MQ=14)
                                               cGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGAtggctggc                                                                                                                  >  2:213627/1‑92 (MQ=255)
                                               cGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGAtggctggc                                                                                                                  <  1:213627/92‑1 (MQ=255)
                                               cGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGc                                                                                                               <  2:607593/95‑1 (MQ=255)
                                               cGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGc                                                                                                               >  1:607593/1‑95 (MQ=255)
                                               cGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAg                                                                   <  1:1115527/139‑1 (MQ=255)
                                               cGTATTCGGTGAACGCACTATGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGTGGCCTTTGAGGTCGTGGTATAAATGACGGAtggct                                                                                                                     <  2:1408193/88‑1 (MQ=17)
                                                    tCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTa                                                              <  1:679905/139‑1 (MQ=255)
                                                       gTGAACGCACTATGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGTGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCgatatacg                                                          >  1:1334811/1‑132 (MQ=14)
                                                            cGCACTATGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGTGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCg                                                                                                                >  2:218729/1‑80 (MQ=18)
                                                            cGCACTATGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGTGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCg                                                                                                                <  1:218729/80‑1 (MQ=18)
                                                            cGCACTATGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGTGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCt                                                                                                              >  1:1325116/1‑82 (MQ=17)
                                                            cGCACTATGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGTGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCt                                                                                                              <  2:1325116/82‑1 (MQ=17)
                                                               aCTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCa                                                   <  2:1173206/139‑1 (MQ=255)
                                                                    ggCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAg                                              >  2:203042/1‑139 (MQ=255)
                                                                     gcgacgcTGGAGCGTCTTCTGAGCCTGCTGTGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCg                                            >  1:1138092/7‑139 (MQ=11)
                                                                       cacaCTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGa                                           <  1:203042/139‑1 (MQ=255)
                                                                        acgcTGGAGCGTCTTCTGAGCCTGCTGTGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGa                                               >  1:866242/4‑133 (MQ=14)
                                                                        acgcTGGAGCGTCTTCTGAGCCTGCTGTGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGa                                               <  2:866242/130‑1 (MQ=11)
                                                                         cgcTGGAGCGTCTTCTGAGCCTGCTGTGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGa                                               <  1:825071/130‑1 (MQ=14)
                                                                         cgcTGGAGCGTCTTCTGAGCCTGCTGTGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGa                                               >  2:825071/3‑132 (MQ=14)
                                                                         cacTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACa                                         <  2:51169/139‑1 (MQ=255)
                                                                           cTGGAGCGTCTTCTGAGCCTGCTGTGGCCTTTGAGGTCGTGGTATAAATGACGGAtggctggc                                                                                                                  >  1:197708/1‑63 (MQ=37)
                                                                           cTGGAGCGTCTTCTGAGCCTGCTGTGGCCTTTGAGGTCGTGGTATAAATGACGGAtggctggc                                                                                                                  <  2:197708/63‑1 (MQ=37)
                                                                           cTGGAGCGTCTTCTGAGCCTGCTGTGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCggtataa                                      >  2:1069770/1‑133 (MQ=11)
                                                                               aGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATaatct                                   >  1:286504/1‑139 (MQ=255)
                                                                                   tcttctGAGCCTGCTGTGGCCTTTGAGGTCGTGGTATAAATGACGGAtggctggc                                                                                                                  >  2:1012920/1‑55 (MQ=25)
                                                                                   tcttctGAGCCTGCTGTGGCCTTTGAGGTCGTGGTATAAATGACGGAtggctggc                                                                                                                  <  1:1012920/55‑1 (MQ=37)
                                                                                   tcttctGAGCCTACTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGAtggct                                                                                                                     <  1:251449/53‑1 (MQ=255)
                                                                                   tcttctGAGCCTACTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGAtggct                                                                                                                     >  2:251449/1‑53 (MQ=255)
                                                                                     ttctGAGCCTGCTGTGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAg                                                                                     <  2:659924/82‑1 (MQ=255)
                                                                                     ttctGAGCCTGCTGTGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAg                                                                                     >  1:659924/1‑82 (MQ=255)
                                                                                        tGAGCCTGCTGTGGCCTTTGAGGTCGTGGTATAAATGACGGAtggct                                                                                                                     >  1:87451/1‑47 (MQ=255)
                                                                                        tGAGCCTGCTGTGGCCTTTGAGGTCGTGGTATAAATGACGGAgggct                                                                                                                     <  2:87451/47‑1 (MQ=255)
                                                                                             ctgctgTGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCa                                                                                                    <  2:472835/59‑1 (MQ=21)
                                                                                             ctgctgTGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCa                                                                                                    >  1:472835/1‑59 (MQ=21)
                                                                                              tgctgTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATAATCTGAATCTGAGacaaca                    >  2:1366836/1‑139 (MQ=255)
                                                                                               gctgTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATaatctgaat                               <  2:43486/127‑1 (MQ=255)
                                                                                               gctgTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATaatctgaat                               >  1:43486/1‑127 (MQ=255)
                                                                                               gctgTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATAATCTGAATCTGAGACAACAt                   >  2:340977/1‑139 (MQ=255)
                                                                                                ctgTCGGCCTTTGAGGTCGTGGTATAAATGACGGAtggct                                                                                                                     >  2:67230/1‑40 (MQ=255)
                                                                                                ctgTCGGCCTTTGAGGTCGTGGTATAAATGACGGAtggct                                                                                                                     <  1:67230/40‑1 (MQ=255)
                                                                                                      gCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTg                                                                                  <  2:1152638/69‑1 (MQ=255)
                                                                                                      gCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTg                                                                                  >  1:1152638/1‑69 (MQ=255)
                                                                                                      gCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGa                                               >  2:1114306/1‑104 (MQ=1)
                                                                                                      gCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGa                                               <  1:1114306/104‑1 (MQ=1)
                                                                                                           tGAGGTCGTGGTATAAATGACGGATGGCTGGTCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGc                                             >  2:268794/1‑101 (MQ=255)
                                                                                                           tGAGGTCGTGGTATAAATGACGGATGGCTGGTCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGc                                             <  1:268794/101‑1 (MQ=255)
                                                                                                             aGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATAATCTGAATCTGAGACAACATCTGGCAAGGCTggg     <  2:339629/139‑1 (MQ=255)
                                                                                                              ggTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATAATCTGAATCTGAGACAACATCTGGCAAGGCTGGGa    <  1:753983/139‑1 (MQ=255)
                                                                                                                tCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATAATCTGAATCTGAGACAACATCTGGCAAGGCTGGGACg  <  1:229217/139‑1 (MQ=255)
                                                                                                                          |                                                                                                                                
TGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAGATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATAATCTGAATCTGAGACAACATCTGGCAAGGCTGGGACG  >  NC_000913/4518892‑4519142

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: