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 232 [0] [0] 225 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)209 (1.150) 101/180 0.1 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 225 total aligned reads displayed.

CAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATG  >  NC_000913_3_bme_pgi/1979184‑1979377
                                                                                               |                                                                                                  
ctcgtcactacacgcacatacaacggaggggggctgcgatttttaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgtata                                                                                                 <  2:4675442‑M2/4‑1 (MQ=255)
ctcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgtata                                                                                                 <  1:1259441‑M2/4‑1 (MQ=255)
ctcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgtat                                                                                                  >  2:2843929‑M2/96‑98 (MQ=255)
ctcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgtat                                                                                                  >  2:1136055‑M2/96‑98 (MQ=255)
  cgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGt                                                                                              >  1:2410123‑M2/94‑100 (MQ=255)
   gtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTc                                                                                             <  2:3273166‑M2/8‑1 (MQ=255)
   gtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTc                                                                                             <  1:557769‑M2/8‑1 (MQ=255)
    tcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCa                                                                                            >  1:4004499‑M2/92‑100 (MQ=255)
      actacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCAt                                                                                           <  1:69528‑M2/10‑1 (MQ=255)
       ctacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatggagatcacacaaaacactcaattagttaacataaatgTATAAGTCATAc                                                                                         >  2:4491127‑M2/89‑100 (MQ=255)
         acacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACtt                                                                                       >  1:1045125‑M2/87‑100 (MQ=255)
          cacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATACGTCATACttt                                                                                      >  2:4847378‑M2/86‑100 (MQ=255)
            cgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttttg                                                                                    <  2:4209551‑M2/17‑1 (MQ=255)
            cgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttttg                                                                                    >  1:134793‑M2/84‑100 (MQ=255)
              cacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttttgtt                                                                                  <  1:2981551‑M2/19‑1 (MQ=255)
              cacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttttgtt                                                                                  >  2:2680585‑M2/82‑100 (MQ=255)
               acatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttttgttt                                                                                 >  2:1320598‑M2/81‑100 (MQ=255)
                catacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttttgttt                                                                                 <  1:1690908‑M2/20‑1 (MQ=255)
                  tacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTgg                                                                              <  2:4602143‑M2/23‑1 (MQ=255)
                  tacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTgg                                                                              >  1:1893550‑M2/78‑100 (MQ=255)
                  tacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTgg                                                                              >  2:1296732‑M2/78‑100 (MQ=255)
                  tacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTgg                                                                              >  2:2163495‑M2/78‑100 (MQ=255)
                    caacggaggggggctgcgattttcaataatgtgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGgt                                                                            >  1:3153812‑M2/76‑100 (MQ=255)
                    caacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaatcactcaattacttaacataaaagTATAAGTCATACTTTTGTTTTgg                                                                              <  1:4544785‑M2/23‑1 (MQ=255)
                    caacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTggg                                                                             <  1:4537266‑M2/24‑1 (MQ=255)
                    caacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTggg                                                                             <  2:213608‑M2/24‑1 (MQ=255)
                    caacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTgg                                                                              >  1:2678281‑M2/76‑98 (MQ=255)
                    caacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGgt                                                                            >  2:3991469‑M2/76‑100 (MQ=255)
                     aacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGgtg                                                                           >  2:2711637‑M2/75‑100 (MQ=255)
                      acggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTCTTGGgtg                                                                           >  1:2326326‑M2/74‑99 (MQ=255)
                         gagggtggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTggg                                                                             <  2:3227135‑M2/24‑1 (MQ=255)
                          aggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTAttt                                                                      <  1:2890001‑M2/31‑1 (MQ=255)
                           ggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTAttt                                                                      >  2:651360‑M2/69‑99 (MQ=255)
                               ggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATc                                                                 >  2:1735710‑M2/65‑100 (MQ=255)
                               ggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATc                                                                 >  2:907481‑M2/65‑100 (MQ=255)
                               ggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATc                                                                 >  1:2489869‑M2/65‑100 (MQ=255)
                                gctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCt                                                                <  2:1081702‑M2/37‑1 (MQ=255)
                                 ctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTg                                                               >  1:2050252‑M2/63‑100 (MQ=255)
                                 ctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTATGGGTGTATTTCaa                                                                   >  1:2644223‑M2/63‑96 (MQ=255)
                                     gattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTa                                                            >  2:3739468‑M2/59‑99 (MQ=255)
                                        tttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTaaa                                                          >  2:1938435‑M2/56‑98 (MQ=255)
                                           caataatgcgtgatgcagatcacacaatacactcaattacttaacataaatgTATAAGTCATACTATAGTTTTGGGTGTATTTCAATCTGTTAAAAAGtt                                                     >  1:1250747‑M2/53‑100 (MQ=255)
                                           caataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGtt                                                     >  2:2405783‑M2/53‑100 (MQ=255)
                                           caataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGt                                                      <  1:4104766‑M2/47‑1 (MQ=255)
                                             ataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGtttt                                                   <  2:3896585‑M2/50‑1 (MQ=255)
                                              taatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGttttt                                                  <  1:4289497‑M2/51‑1 (MQ=255)
                                               aatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGttttt                                                  >  2:414282‑M2/49‑99 (MQ=255)
                                                   cgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgcta                                             >  2:631198‑M2/45‑100 (MQ=255)
                                                   cgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgcta                                             <  1:506027‑M2/56‑1 (MQ=255)
                                                   cgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgct                                              <  1:4615232‑M2/55‑1 (MQ=255)
                                                    gtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgctac                                            >  2:3586499‑M2/44‑100 (MQ=255)
                                                     tgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgctac                                            >  1:4804125‑M2/43‑99 (MQ=255)
                                                      gatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgctacgc                                          >  1:812635‑M2/42‑100 (MQ=255)
                                                       atgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgctacgct                                         <  1:2721019‑M2/60‑1 (MQ=255)
                                                       atgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgctacgct                                         <  1:33790‑M2/60‑1 (MQ=255)
                                                       atgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgctacgct                                         <  2:225638‑M2/60‑1 (MQ=255)
                                                        tgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgctacgc                                          >  2:2322907‑M2/40‑98 (MQ=255)
                                                        tgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTGTcgctacgc                                          >  2:3115411‑M2/40‑98 (MQ=255)
                                                          cagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTag                                       <  2:25717‑M2/62‑1 (MQ=255)
                                                          cagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGATTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTag                                       >  2:4709375‑M2/38‑99 (MQ=255)
                                                           agatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTagca                                     <  2:3935984‑M2/64‑1 (MQ=255)
                                                           agatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTagc                                      >  1:4245608‑M2/37‑99 (MQ=255)
                                                            gatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTagc                                      >  2:3832625‑M2/36‑98 (MQ=255)
                                                            gatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTagc                                      >  2:3019271‑M2/36‑98 (MQ=255)
                                                             atcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTagc                                      >  1:3321733‑M2/35‑97 (MQ=255)
                                                              tcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTagcaag                                   <  1:3883101‑M2/66‑1 (MQ=255)
                                                               cacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTagcaagc                                  >  2:1007217‑M2/33‑99 (MQ=255)
                                                                acacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTagcaagca                                 >  2:1607637‑M2/32‑99 (MQ=255)
                                                                acacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTagcaagca                                 <  1:3522512‑M2/68‑1 (MQ=255)
                                                                acacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTagcaagca                                 >  2:2453405‑M2/32‑99 (MQ=255)
                                                                 cacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAgcaagcaa                                >  2:849036‑M2/31‑99 (MQ=255)
                                                                 cacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAgcaagcaa                                >  2:767980‑M2/31‑99 (MQ=255)
                                                                 cacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAgcaagcaa                                >  2:1821922‑M2/31‑99 (MQ=255)
                                                                   caaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCaaaaa                             <  1:3433372‑M2/72‑1 (MQ=255)
                                                                    aaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAAt                            <  1:3920953‑M2/73‑1 (MQ=255)
                                                                     aaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAAt                            <  2:1346109‑M2/73‑1 (MQ=255)
                                                                        cactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGaaa                        <  2:3402799‑M2/77‑1 (MQ=255)
                                                                        cactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGaaa                        >  1:4198542‑M2/24‑100 (MQ=255)
                                                                        cactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGaaa                        >  1:2526570‑M2/24‑100 (MQ=255)
                                                                         actcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGaaa                        >  2:3557121‑M2/23‑99 (MQ=255)
                                                                            caattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAg                     >  1:2879408‑M2/20‑99 (MQ=255)
                                                                            caattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAg                     >  2:4239726‑M2/20‑99 (MQ=255)
                                                                            caattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAg                     >  2:4309322‑M2/20‑99 (MQ=255)
                                                                            caattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTACTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAg                     >  1:3719597‑M2/20‑99 (MQ=255)
                                                                             aattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGa                   <  2:4021796‑M2/82‑1 (MQ=255)
                                                                             aattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCAGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGa                   <  1:3452937‑M2/82‑1 (MQ=255)
                                                                             aattacttaacataaatgTATAAGTCATACTTTTGTTTTGCGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAgg                    <  1:4733772‑M2/81‑1 (MQ=255)
                                                                               ttacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGa                   <  2:4350342‑M2/82‑1 (MQ=255)
                                                                                 acttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCATTCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGaataa               >  1:2202042‑M2/15‑100 (MQ=255)
                                                                                 acttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGaataa               >  2:2274248‑M2/15‑100 (MQ=255)
                                                                                 acttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGaataa               >  2:299091‑M2/15‑100 (MQ=255)
                                                                                  cttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGaataat              >  2:2489374‑M2/14‑100 (MQ=255)
                                                                                    taacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATc             >  2:1652282‑M2/12‑99 (MQ=255)
                                                                                     aacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCg            >  1:2082476‑M2/11‑99 (MQ=255)
                                                                                     aacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCg            <  1:3025816‑M2/89‑1 (MQ=255)
                                                                                       cataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGaaa         >  2:2048465‑M2/9‑100 (MQ=255)
                                                                                       cataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTATAGATACGCTAGCAAGCAAAAATGAAACAGGAATAATCGa           >  1:284554‑M2/9‑98 (MQ=255)
                                                                                        ataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAAt        >  1:3595768‑M2/8‑100 (MQ=255)
                                                                                            atgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGa    <  2:715667‑M2/97‑1 (MQ=255)
                                                                                              gTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATg  >  1:2422967‑M2/2‑100 (MQ=255)
                                                                                               |                                                                                                  
CAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATG  >  NC_000913_3_bme_pgi/1979184‑1979377

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