Predicted mutation
evidence seq id position mutation freq annotation gene description
MC JC NC_000913 1,978,503 Δ776 bp 100% insB1insA insB1, insA

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 1978503–1979270 1979278 9–776 90 [0] [1] 89 insB1–insA insB1, insA

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 19785020 (0.000)77 (0.920) 53/162 NT 100% intergenic (‑305/+16) flhD/insB1 flagellar class II regulon transcriptional activator, with FlhC/IS1 transposase B
?NC_000913 1979279 = 0 (0.000)intergenic (‑64/‑474) insA/uspC IS1 repressor TnpA/universal stress protein

AACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGAT  >  NC_000913/1979182‑1979376
                                                                                                 |                                                                                                 
tactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgtat                                                                                                 <  2:2329753‑M2/3‑1 (MQ=255)
 actcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgtata                                                                                                >  1:1021346‑M2/97‑100 (MQ=255)
  ctcgtcactacacgcacatacaacggagggggcctgcgattttccataaagcgtgctgcagatcacaccaaacactcaattacttaacataaatgTATaa                                                                                               >  1:625419‑M2/96‑100 (MQ=255)
   tcgtcactacacgcacatacaacgtaggggcgctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAg                                                                                              <  1:187326‑M2/6‑1 (MQ=255)
   tcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAg                                                                                              <  2:2045316‑M2/6‑1 (MQ=255)
   tcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAg                                                                                              <  2:472476‑M2/6‑1 (MQ=255)
    cgtcactacacgcacatacaacggaggggggctgcgtttttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGt                                                                                             >  2:1825272‑M2/94‑100 (MQ=255)
    cgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGt                                                                                             >  1:2405002‑M2/94‑100 (MQ=255)
       cactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCAt                                                                                          >  2:2312430‑M2/91‑100 (MQ=255)
         ctacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATAc                                                                                        <  2:1395140‑M2/12‑1 (MQ=255)
         ctacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATAc                                                                                        >  1:1498634‑M2/89‑100 (MQ=255)
            cacgcacatacaacggaggggggctgagattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACtt                                                                                      >  2:1577039‑M2/86‑99 (MQ=255)
              ctcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttttg                                                                                   <  2:534997‑M2/17‑1 (MQ=255)
              cgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttttg                                                                                   >  1:872830‑M2/84‑100 (MQ=255)
              cgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttttg                                                                                   >  2:122751‑M2/84‑100 (MQ=255)
              cgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttttg                                                                                   <  1:2029099‑M2/17‑1 (MQ=255)
                cacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATacttttgag                                                                                 >  1:534995‑M2/82‑98 (MQ=255)
                cacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttttgtt                                                                                 >  1:1413119‑M2/82‑100 (MQ=255)
                  catacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttttgtttt                                                                               >  2:2008738‑M2/80‑100 (MQ=255)
                  catacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttttgtttt                                                                               >  2:2008165‑M2/80‑100 (MQ=255)
                  catacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttttgtttt                                                                               >  1:2560561‑M2/80‑100 (MQ=255)
                  catacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttttgtttt                                                                               <  2:1553101‑M2/21‑1 (MQ=255)
                   atacaacggaggggggctgcgattttcaataattcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttttcttctg                                                                              >  1:1222803‑M2/79‑100 (MQ=255)
                   atacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttttgttttg                                                                              >  1:1049370‑M2/79‑100 (MQ=255)
                      caacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTggg                                                                            <  1:2400304‑M2/24‑1 (MQ=255)
                      caacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGgt                                                                           >  1:571094‑M2/76‑100 (MQ=255)
                       aacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGgtg                                                                          >  1:57152‑M2/75‑100 (MQ=255)
                         cggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTa                                                                        >  2:426975‑M2/73‑100 (MQ=255)
                           gaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTAtt                                                                      <  2:1518292‑M2/30‑1 (MQ=255)
                            aggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTAttt                                                                     >  2:1497896‑M2/70‑100 (MQ=255)
                            aggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTAttt                                                                     <  2:1039336‑M2/31‑1 (MQ=255)
                            aggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTAttt                                                                     <  2:311179‑M2/31‑1 (MQ=255)
                             ggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTc                                                                    <  2:317378‑M2/32‑1 (MQ=255)
                              gggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTc                                                                    <  1:1942345‑M2/32‑1 (MQ=255)
                               ggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCaa                                                                  >  2:2532405‑M2/67‑100 (MQ=255)
                               ggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCaa                                                                  <  2:1944084‑M2/34‑1 (MQ=255)
                                 ggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATc                                                                <  1:2023443‑M2/36‑1 (MQ=255)
                                   ctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTg                                                              >  2:30451‑M2/63‑100 (MQ=255)
                                    tgcgattttcaataatgcgtgatgcagatcacacaaaacacttaagtgcttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGGATTTCAATCt                                                               <  1:930004‑M2/37‑1 (MQ=255)
                                    tgcaattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGt                                                             <  1:1973902‑M2/39‑1 (MQ=255)
                                          tttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTaaaaa                                                       <  2:1515705‑M2/45‑1 (MQ=255)
                                           ttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTCTATTTCAATCTGTTAAAAAg                                                      <  1:1577030‑M2/46‑1 (MQ=255)
                                            tcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGttaaaaatt                                                     >  1:1992144‑M2/54‑98 (MQ=255)
                                             caataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGtt                                                    >  2:638308‑M2/53‑100 (MQ=255)
                                               ataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGtttt                                                  <  2:1176620‑M2/50‑1 (MQ=255)
                                                taatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGttttt                                                 <  1:980214‑M2/51‑1 (MQ=255)
                                                 aatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTc                                                <  1:1990688‑M2/52‑1 (MQ=255)
                                                      gtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgctac                                           <  1:1114827‑M2/57‑1 (MQ=255)
                                                          tgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgctacgcta                                       >  2:1285842‑M2/40‑100 (MQ=255)
                                                          tgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgctacgcta                                       >  2:958259‑M2/40‑100 (MQ=255)
                                                          ngcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgctacgcta                                       >  1:671591‑M2/40‑100 (MQ=255)
                                                           gcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgctacgcta                                       <  1:1324617‑M2/61‑1 (MQ=255)
                                                           gcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTag                                      >  2:92999‑M2/39‑100 (MQ=255)
                                                           gcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTag                                      >  2:93655‑M2/39‑100 (MQ=255)
                                                             agatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTagca                                    >  2:2055578‑M2/37‑100 (MQ=255)
                                                             agatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTagca                                    <  2:1750497‑M2/64‑1 (MQ=255)
                                                             agatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTagca                                    <  2:86084‑M2/64‑1 (MQ=255)
                                                             agatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTagca                                    <  2:880609‑M2/64‑1 (MQ=255)
                                                             agatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTAAATCTGTTAAAAAGTTTTTCGCTACGCTagca                                    >  2:792009‑M2/37‑100 (MQ=255)
                                                              gatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTagcaa                                   >  2:1424186‑M2/36‑100 (MQ=255)
                                                                tcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTagcaagc                                 <  2:310006‑M2/67‑1 (MQ=255)
                                                                tcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTagcaag                                  >  1:432244‑M2/34‑99 (MQ=255)
                                                                  acacaaaacactcaattacttaacataaatgTATAATTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAgcaagcaa                               <  2:1815224‑M2/69‑1 (MQ=255)
                                                                  acacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAgcaagcaa                               <  2:167014‑M2/69‑1 (MQ=255)
                                                                   cacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCaaa                              >  1:1449957‑M2/31‑100 (MQ=255)
                                                                   cacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCaaa                              >  2:356385‑M2/31‑100 (MQ=255)
                                                                   cacaaaacactcaattacttaacataaatgTAGAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCaaa                              <  1:1028119‑M2/70‑1 (MQ=255)
                                                                      aaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAAt                           >  2:960464‑M2/28‑100 (MQ=255)
                                                                      aaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAAt                           >  2:2147245‑M2/28‑100 (MQ=255)
                                                                         acactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGaa                        <  1:2218407‑M2/76‑1 (MQ=255)
                                                                         acactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGaa                        <  2:40741‑M2/76‑1 (MQ=255)
                                                                          cactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGaaa                       <  2:1449965‑M2/77‑1 (MQ=255)
                                                                           actcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAAc                      <  2:191334‑M2/78‑1 (MQ=255)
                                                                            ctcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTTTCTATGTTAGCAAGCAAAAAtgaaaaa                     >  1:454780‑M2/22‑98 (MQ=255)
                                                                            ctcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACa                     >  2:1867700‑M2/22‑100 (MQ=255)
                                                                             tcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAg                    <  1:122750‑M2/80‑1 (MQ=255)
                                                                              caattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAgg                   <  2:1708501‑M2/81‑1 (MQ=255)
                                                                              caattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAgg                   <  2:1708913‑M2/81‑1 (MQ=255)
                                                                              caattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAgg                   <  1:1429797‑M2/81‑1 (MQ=255)
                                                                              caattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTAGCTACGCTAGCAAGCAAAAATGAAACAg                    >  1:1140750‑M2/20‑99 (MQ=255)
                                                                                       aacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGa          <  2:1054733‑M2/90‑1 (MQ=255)
                                                                                       aacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGa          <  1:1825262‑M2/90‑1 (MQ=255)
                                                                                        acataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGaa         >  2:268182‑M2/10‑100 (MQ=255)
                                                                                         cataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGaaa        <  2:264391‑M2/92‑1 (MQ=255)
                                                                                           taaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATg      >  2:1123702‑M2/7‑100 (MQ=255)
                                                                                             aatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTTTCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATggg    >  1:506894‑M2/5‑100 (MQ=255)
                                                                                             aatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATggg    >  1:15869‑M2/5‑100 (MQ=255)
                                                                                              aTGTATAAGTCATACTTTTGTTTTAAGTAAATAAAAAAATGTAAAAAAGTTTTTCGATAAAAAAAAAAAAAAAAATGAAACAGGAATAATCGAAATGGGa   <  1:521709/100‑1 (MQ=255)
                                                                                               tgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGAt  >  2:965726‑M2/3‑100 (MQ=255)
                                                                                                 |                                                                                                 
AACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGAT  >  NC_000913/1979182‑1979376

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