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

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

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

Only 100 of 162 total aligned reads displayed.

ACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGT  >  NC_000913_3_bme_pgi/1979183‑1979378
                                                                                                |                                                                                                   
actcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgtat                                                                                                   >  1:2074011‑M2/97‑99 (MQ=255)
actcgtcactacacgcacatacaacggaggggggccgcgaatttcaataatgcgtgatgcagatcacacaaaaaactcaattacttaacataaatgtat                                                                                                   >  2:1430217‑M2/97‑99 (MQ=255)
 ctcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgtata                                                                                                  <  2:4533230‑M2/4‑1 (MQ=255)
    gtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAg                                                                                                <  1:438623‑M2/6‑1 (MQ=255)
    gtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTc                                                                                              <  2:4272367‑M2/8‑1 (MQ=255)
       actacacgcacatacaacggaggggggctgcgattttcaataatgcgttatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATa                                                                                           <  2:2853152‑M2/11‑1 (MQ=255)
       actacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCAt                                                                                            <  2:1141615‑M2/10‑1 (MQ=255)
       actacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATa                                                                                           >  2:3788658‑M2/90‑100 (MQ=255)
        ctacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATAc                                                                                          <  1:368381‑M2/12‑1 (MQ=255)
        ctacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATAc                                                                                          >  1:1239921‑M2/89‑100 (MQ=255)
          acacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACtt                                                                                        >  1:186833‑M2/87‑100 (MQ=255)
           cacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttt                                                                                       <  2:4575450‑M2/15‑1 (MQ=255)
             cgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttttg                                                                                     <  1:2013301‑M2/17‑1 (MQ=255)
               cacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttttgtt                                                                                   >  2:4646980‑M2/82‑100 (MQ=255)
               cacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttttgtt                                                                                   >  1:4082107‑M2/82‑100 (MQ=255)
                acatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttttgtt                                                                                   <  2:1114020‑M2/19‑1 (MQ=255)
                 catacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttttgttt                                                                                  >  2:1538396‑M2/80‑99 (MQ=255)
                 catacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttttgttt                                                                                  <  2:3116764‑M2/20‑1 (MQ=255)
                 catacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttttgttt                                                                                  >  2:3347942‑M2/80‑99 (MQ=255)
                  atacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttttgttttg                                                                                <  2:1615387‑M2/22‑1 (MQ=255)
                  atacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttttgttttg                                                                                <  2:4467816‑M2/22‑1 (MQ=255)
                   tacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttttgttttg                                                                                >  2:1863143‑M2/78‑99 (MQ=255)
                   tacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttttgttt                                                                                  >  2:895671‑M2/78‑97 (MQ=255)
                   tacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTgg                                                                               >  2:3791307‑M2/78‑100 (MQ=255)
                     caacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTggg                                                                              >  2:3245082‑M2/76‑99 (MQ=255)
                     caacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTggg                                                                              <  1:2986275‑M2/24‑1 (MQ=255)
                     caacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGgt                                                                             >  1:8333‑M2/76‑100 (MQ=255)
                     caacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGgt                                                                             <  2:2713023‑M2/25‑1 (MQ=255)
                     caacggaggggggctgcgattttaaataatgcgtgatgcaggtcacacaaaacactcaattacttaacataaatgTATAAGTCATACttttgtt                                                                                   >  2:78935‑M2/76‑94 (MQ=255)
                       acggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGgtg                                                                            >  2:757508‑M2/74‑99 (MQ=255)
                        cggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTa                                                                          <  1:3061220‑M2/28‑1 (MQ=255)
                         ggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTAt                                                                         >  1:4539738‑M2/72‑100 (MQ=255)
                            ggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTAttt                                                                       >  1:4745564‑M2/69‑99 (MQ=255)
                             gggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTc                                                                      <  1:4393379‑M2/32‑1 (MQ=255)
                              ggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCaa                                                                    <  1:1579282‑M2/34‑1 (MQ=255)
                                ggctgcgattttcaataatgcgtgatgcagaacacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAAt                                                                   >  2:3581781‑M2/65‑99 (MQ=255)
                                  ctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTg                                                                >  2:81016‑M2/63‑100 (MQ=255)
                                  ctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTg                                                                >  1:4672332‑M2/63‑100 (MQ=255)
                                   tgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGt                                                               >  2:2350040‑M2/62‑100 (MQ=255)
                                     cgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGtt                                                              >  1:2776157‑M2/60‑99 (MQ=255)
                                     cgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGt                                                               <  2:2496952‑M2/39‑1 (MQ=255)
                                     cgattttcaataatgcgtgatgcagatcacacaaaacactcaactacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGtt                                                              >  1:4701459‑M2/60‑99 (MQ=255)
                                         tttcaataatgagtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTCTTGTTCTGGGTGTCTTACAATCTgtaaaagc                                                         >  1:3973750‑M2/56‑98 (MQ=255)
                                          ttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTaaaaa                                                         <  2:2165282‑M2/45‑1 (MQ=255)
                                             aataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGttt                                                     >  1:3139925‑M2/52‑100 (MQ=255)
                                                 atgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTc                                                  <  1:2841621‑M2/52‑1 (MQ=255)
                                                  tgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgc                                                <  2:3507594‑M2/54‑1 (MQ=255)
                                                   gcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgc                                                >  1:2294891‑M2/46‑99 (MQ=255)
                                                     gtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgctac                                             <  2:1923520‑M2/57‑1 (MQ=255)
                                                       gatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgctacgc                                           >  2:3310969‑M2/42‑100 (MQ=255)
                                                       gatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgctacgc                                           >  2:4146968‑M2/42‑100 (MQ=255)
                                                       gatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgctacgc                                           >  2:2531177‑M2/42‑100 (MQ=255)
                                                        atgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgctacgc                                           >  2:126157‑M2/41‑99 (MQ=255)
                                                         tgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgctacgcta                                         <  1:2646082‑M2/61‑1 (MQ=255)
                                                         tgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgctacgc                                           >  2:695080‑M2/40‑98 (MQ=255)
                                                         tgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgctacgc                                           >  2:1997548‑M2/40‑98 (MQ=255)
                                                           cagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTagc                                       <  1:1769494‑M2/63‑1 (MQ=255)
                                                           cagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTag                                        <  2:2630602‑M2/62‑1 (MQ=255)
                                                           cagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTag                                        <  1:145296‑M2/62‑1 (MQ=255)
                                                            agatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTagc                                       >  2:2350401‑M2/37‑99 (MQ=255)
                                                            agatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTagc                                       >  2:669359‑M2/37‑99 (MQ=255)
                                                            agatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTagc                                       >  1:449343‑M2/37‑99 (MQ=255)
                                                             gatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTTCGCTagc                                       >  1:2239240‑M2/36‑98 (MQ=255)
                                                             gatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATAAAAAAAAAAAAAAAAGTTTTTCGCTACGCTagca                                      <  2:557344‑M2/64‑1 (MQ=255)
                                                              atcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTagc                                       >  2:2858918‑M2/35‑97 (MQ=255)
                                                               tcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTagcaagc                                   <  2:4225183‑M2/67‑1 (MQ=255)
                                                               tcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTagcaag                                    >  1:379431‑M2/34‑99 (MQ=255)
                                                                  cacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAgcaagcaa                                 >  1:107740‑M2/31‑99 (MQ=255)
                                                                   acaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCaaa                                >  1:3297167‑M2/30‑99 (MQ=255)
                                                                     aaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAAt                             <  1:4226486‑M2/73‑1 (MQ=255)
                                                                     aaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAAt                             >  1:1385241‑M2/28‑100 (MQ=255)
                                                                     aaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAAt                             >  1:1385812‑M2/28‑100 (MQ=255)
                                                                     aaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAAt                             <  1:1076777‑M2/73‑1 (MQ=255)
                                                                     aaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAAt                             <  1:3026801‑M2/73‑1 (MQ=255)
                                                                        acactcaattacttaacataaatgTATAAGTAATACTTTTGTTTTGGGGGTATTTCAAACTGTTAAAAAGTTTTTCGCTACGCTAGAAAGAAAAAATGaa                          <  2:3744725‑M2/76‑1 (MQ=255)
                                                                         cactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGaaa                         <  1:1869887‑M2/77‑1 (MQ=255)
                                                                         cactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGaaa                         <  1:1949787‑M2/77‑1 (MQ=255)
                                                                          actcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAAc                        <  1:1012995‑M2/78‑1 (MQ=255)
                                                                          actcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAAc                        <  1:1012481‑M2/78‑1 (MQ=255)
                                                                           ctcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGaaa                         >  1:2395865‑M2/22‑98 (MQ=255)
                                                                            tcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAg                      >  1:3222337‑M2/21‑100 (MQ=255)
                                                                             caattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTTAATCTGTTAAAAAGTTTTTCGCTACGCTAGCTAGCAAAAATGaaa                         >  2:4225543‑M2/20‑96 (MQ=255)
                                                                             caattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAg                      >  2:1032528‑M2/20‑99 (MQ=255)
                                                                              aattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAgg                     <  1:957309‑M2/81‑1 (MQ=255)
                                                                              aattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTGTTCGCTACGCTAGCAAGCAAAAATGAAACAgg                     <  2:3622715‑M2/81‑1 (MQ=255)
                                                                                 nacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGaat                  >  2:2710308‑M2/16‑99 (MQ=255)
                                                                                  acttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGaata                 <  1:3907996‑M2/85‑1 (MQ=255)
                                                                                   cttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGaataat               >  2:11814‑M2/14‑100 (MQ=255)
                                                                                     taacataaatgTATAAGTCATACTTTTGTTTTTGGTGTATTTCAATCTGTTAAAAAGTTTTTCGTTACGTTAGCAAGCAAAAATGAAACAGGAATAATc              >  2:570342‑M2/12‑99 (MQ=255)
                                                                                     taacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATc              >  1:684017‑M2/12‑99 (MQ=255)
                                                                                     taacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATc              >  1:4235199‑M2/12‑99 (MQ=255)
                                                                                      aacgtaaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATc              <  2:4334916‑M2/88‑1 (MQ=255)
                                                                                       acataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGaa           <  2:1121124‑M2/91‑1 (MQ=255)
                                                                                        cataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGaaa          >  1:4503626‑M2/9‑100 (MQ=255)
                                                                                         ataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAAt         >  2:2035597‑M2/8‑100 (MQ=255)
                                                                                         ataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAAt         <  1:2582684‑M2/93‑1 (MQ=255)
                                                                                          tAAATGTATAAGTCAGACTTAAGCTTTGGGTGTATTTCTATCTGTTAAAAAGTTTTTCGTTACGCTAGCAAGTAAAATTGAAACAGGAATAATCGGAATg        >  2:1768480/1‑100 (MQ=255)
                                                                                            aatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATgg       <  2:306830‑M2/95‑1 (MQ=255)
                                                                                             atgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGa     >  2:1037792‑M2/4‑100 (MQ=255)
                                                                                                tataAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGt  >  1:3879182/1‑100 (MQ=255)
                                                                                                |                                                                                                   
ACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGT  >  NC_000913_3_bme_pgi/1979183‑1979378

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