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

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

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

TGATGTCCGGCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTATAAATAA  >  NC_000913_3_pae_pgi/1979133‑1979427
                                                                                                                                                  |                                                                                                                                                    
aaatgggtcgttacacttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgtat                                                                                                                                                    <  2:651119‑M2/3‑1 (MQ=255)
  atgggtctttacatgtatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATaa                                                                                                                                                  <  2:1181568‑M2/5‑1 (MQ=255)
          ttacacttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACt                                                                                                                                          >  1:775646‑M2/137‑149 (MQ=255)
           tacacttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATAAtt                                                                                                                                         >  2:1015709‑M2/136‑149 (MQ=255)
               cttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttttgt                                                                                                                                     >  1:1014214‑M2/132‑149 (MQ=255)
                 tatctaagacttttcctaagtcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttttgttt                                                                                                                                   <  2:1073352‑M2/20‑1 (MQ=255)
                   tctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttttgttttg                                                                                                                                 <  1:735798‑M2/22‑1 (MQ=255)
                   tctaagatttctcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcatttacttaacataaatgTATAAGTCATACttttgttttg                                                                                                                                 <  2:1316828‑M2/22‑1 (MQ=255)
                    tgaatatttttcataaaacgacgcaactgtactcgttactacacgcacataaaatggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaataTATAAGTCATACCTTTGTTTTgg                                                                                                                                <  2:840688‑M2/23‑1 (MQ=255)
                     taagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAg                                                                                                                                                 >  1:164534‑M2/126‑131 (MQ=255)
                     taagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAg                                                                                                                                                 <  2:164534‑M2/6‑1 (MQ=255)
                      aagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGgt                                                                                                                              >  1:521532‑M2/125‑149 (MQ=255)
                                ctaaatctacgcaactgtacccgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAAt                                                                                                                    <  1:274489‑M2/35‑1 (MQ=255)
                                ctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaacccaaccagcttcgtaacataaaggTATAAGTCATACTTTTGTTTTGGGTGTGTTTCAAt                                                                                                                    >  1:1101919‑M2/115‑149 (MQ=255)
                                       gacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgatattcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTaa                                                                                                             >  1:480957‑M2/108‑149 (MQ=255)
                                        acgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgatggtcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTaaa                                                                                                            >  2:7029‑M2/107‑149 (MQ=255)
                                         cgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTaaaa                                                                                                           >  1:288240‑M2/106‑149 (MQ=255)
                                          gcaactgtactcgtcactacacggacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCAGTTaaaaa                                                                                                          <  2:222257‑M2/45‑1 (MQ=255)
                                             actgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGtt                                                                                                       <  1:792694‑M2/48‑1 (MQ=255)
                                              ctgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGttt                                                                                                      >  2:217103‑M2/101‑149 (MQ=255)
                                                   gtcgtcactaaacacgcatacaacggagggaggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaattTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgc                                                                                                 <  2:1399041‑M2/54‑1 (MQ=255)
                                                    tcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgct                                                                                                >  2:774690‑M2/95‑149 (MQ=255)
                                                    tcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgct                                                                                                >  2:436653‑M2/95‑149 (MQ=255)
                                                      gtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgctac                                                                                              >  2:231831‑M2/93‑149 (MQ=255)
                                                              acgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTagcaa                                                                                      <  1:900647‑M2/65‑1 (MQ=255)
                                                               cgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTGTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTagcaag                                                                                     <  1:604466‑M2/66‑1 (MQ=255)
                                                                    atacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTgg                                                                                                                                <  1:1195508‑M2/23‑1 (MQ=255)
                                                                    atacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTgg                                                                                                                                >  2:1195508‑M2/79‑101 (MQ=255)
                                                                     tacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCaaaaa                                                                               >  2:380557‑M2/78‑149 (MQ=255)
                                                                      acaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAAt                                                                              >  1:670137‑M2/77‑149 (MQ=255)
                                                                      acaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAAt                                                                              <  1:231831‑M2/73‑1 (MQ=255)
                                                                        aacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAAGGa                                                                            >  1:320518‑M2/75‑149 (MQ=255)
                                                                           ggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAAc                                                                         <  1:774690‑M2/78‑1 (MQ=255)
                                                                             aggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAg                                                                       <  1:1015709‑M2/80‑1 (MQ=255)
                                                                               gggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGa                                                                     <  1:1181267‑M2/82‑1 (MQ=255)
                                                                                     tgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgctac                                                                                              >  2:765650‑M2/62‑118 (MQ=255)
                                                                                     tgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgctac                                                                                              <  1:765650‑M2/57‑1 (MQ=255)
                                                                                      gcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCg                                                              >  1:1143671‑M2/61‑149 (MQ=255)
                                                                                       cgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGa                                                             <  1:361756‑M2/90‑1 (MQ=255)
                                                                                            tttaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATgg                                                        <  2:1262534‑M2/95‑1 (MQ=255)
                                                                                            ttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATgg                                                        >  1:965978‑M2/55‑149 (MQ=255)
                                                                                            ttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGCTCTTGTGTTCTCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATgg                                                        >  1:819528‑M2/55‑149 (MQ=255)
                                                                                                 taatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGt                                                   >  2:332622‑M2/50‑149 (MQ=255)
                                                                                                   gtgcgtgatgcatatcacacaaaacacgcaatcacttatcttaaatgTATGAGTCAGACTTTTGTTTTGGGTTTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGTGATGTTg                                                 <  2:704583‑M2/102‑1 (MQ=255)
                                                                                                    tgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTgc                                                >  1:1084053‑M2/47‑149 (MQ=255)
                                                                                                     cctttttgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTgcg                                               >  2:533160‑M2/46‑149 (MQ=255)
                                                                                                      cgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTgc                                                >  2:638221‑M2/45‑147 (MQ=255)
                                                                                                      cgtgatgcagatcacacaaaacaatcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTgc                                                <  1:638221‑M2/103‑1 (MQ=255)
                                                                                                       gtgatgcagatcacacaaaacactcaattacttaacataaatgTCTAAGTCATACTTTTGTTTTGGATGTATTTCAATCTGTTAAAAAGTTTTTCGTTACGCTTAGAAGAAAAAATGAAACAAGAATATTCAAAAAAggacgttagtca                                             >  2:1173052‑M2/44‑144 (MQ=255)
                                                                                                         gatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGcaca                                           <  1:348300‑M2/108‑1 (MQ=255)
                                                                                                          atgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCCCAg                                          >  1:113527‑M2/41‑149 (MQ=255)
                                                                                                                   acacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATaa                                 <  1:436653‑M2/118‑1 (MQ=255)
                                                                                                                    cacaaaacactcaattacttaacataaatgTATAAGTCATATTTTTGTTTTGGGTGTATTTCAATGTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCGAAAATGAAACAGAAATAATCGAAATGGGATGTTGCTCACCGTCAAAATAAc                                >  2:1333755‑M2/31‑149 (MQ=255)
                                                                                                                     acaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTTTTTTTTTTTTTAAGAAAAAATGAAACAGGAATTATCGAAATGGGATGTTGCGCACAGTAAATATAACt                               >  2:789421‑M2/30‑149 (MQ=255)
                                                                                                                     acaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCACccagtcaaataaaca                               >  2:290627‑M2/30‑143 (MQ=255)
                                                                                                                                   tacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACGGTCTAAATAACTCACCGTAAATAatc                 <  2:1014214‑M2/134‑1 (MQ=255)
                                                                                                                                       taacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACCGTCAAAATAACTACACGTAAAAAATCATCt             >  2:392193‑M2/12‑149 (MQ=255)
                                                                                                                                         acataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGc           >  2:1234236‑M2/10‑149 (MQ=255)
                                                                                                                                          cataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCt          >  1:305514‑M2/9‑149 (MQ=255)
                                                                                                                                              aatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTataa      >  2:886923‑M2/5‑149 (MQ=255)
                                                                                                                                               atgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTataaa     <  1:1238564‑M2/146‑1 (MQ=255)
                                                                                                                                                 gTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTTTTTTTTTTTTCAAGCAAAAATGAGACAGGCATAATCGAAATTGGATGTTGCGGACAGTCAAGATACCTCACCTTACTTAATCATCTCAtatatcta   >  2:1094108/1‑145 (MQ=255)
                                                                                                                                                 gTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTataaata   >  2:99508‑M2/2‑149 (MQ=255)
                                                                                                                                                  tataAGTCATACTTCTGTTTTGGGTGGATTCCAATCTGTTAAAAAGGTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTataaataa  <  1:392193/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
TGATGTCCGGCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTATAAATAA  >  NC_000913_3_pae_pgi/1979133‑1979427

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