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

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 1978503–1979268 1979278 11–776 74 [0] [1] 76 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)69 (1.000) 52/266 0.2 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

CTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGC  >  NC_000913/1978365‑1978644
                                                                                                                                          |                                                                                                                                              
ctgtacaCTTATCTAAGATGTTTCCTAAACCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataa                                                                                                                                           <  1:782856‑M1/140‑6 (MQ=255)
cTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTACGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataa                                                                                                                                           <  2:369199‑M1/143‑6 (MQ=255)
 tttACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtaaaag                                                                                                                                          >  2:778033‑M1/1‑137 (MQ=255)
    acacTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGCGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtca                                                                                                                                       >  1:1269706‑M1/1‑134 (MQ=255)
    acacTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTAACAAGGCCAGGGTGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtca                                                                                                                                       >  1:606547‑M1/1‑134 (MQ=255)
      acTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcata                                                                                                                                     >  1:1118452‑M1/1‑132 (MQ=255)
              aaGATTTTTCCTAAATCGACGCAACTGTGCTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgtt                                                                                                                             >  1:874072‑M1/1‑124 (MQ=255)
                 aTTTTTTCTATATCGACGTACATTTACTCGTTACTACACGCACATACAACGGAGGGGGGCTGCTATTTTCAATAATGCGTGATGCAGATCACACAATACACTCAATTACTTAACATAAATGtataagtcatacttttgttttg                                                                                                                          <  2:319001‑M1/143‑23 (MQ=255)
                  ttgctccgaAATCGACGTGACTGTACGCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAAGTACTTAACATAAATGtataagtcatacttttgttttgg                                                                                                                         <  2:924814‑M1/139‑24 (MQ=255)
                        cTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACGCTCAATTACTTAACATAAATGtataagtcatacctttgtgttgggtgtat                                                                                                                   >  1:609930‑M1/1‑114 (MQ=255)
                        cTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGACTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttggttgtat                                                                                                                   >  2:112204‑M1/1‑114 (MQ=255)
                         tAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGCTTTTCAATACTGAGTGATGCAGACCACACAAAACACTCAATTACTTAACATACCTGtatacgtcatactattgctctgggtgtatt                                                                                                                  >  2:302010‑M1/1‑113 (MQ=255)
                         tAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatt                                                                                                                  >  2:116027‑M1/1‑113 (MQ=255)
                         tAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatactttagttttggatgtatt                                                                                                                  >  2:22034‑M1/1‑113 (MQ=255)
                          aaaTCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtattt                                                                                                                 <  2:61579‑M1/143‑32 (MQ=255)
                            aTCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCGCACAAAACACCCAATTACTTAACATAAATGtataagtcatacgtttgttttgggtgtatttca                                                                                                               >  2:193200‑M1/1‑110 (MQ=255)
                              cGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATAtataagtcatacttttgtattgggtgtatatcaat                                                                                                             >  2:410136‑M1/1‑108 (MQ=255)
                              cGACGCAACGGTACTCGTCTCTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaat                                                                                                             <  2:1272329‑M1/143‑36 (MQ=255)
                               gACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatc                                                                                                            >  2:726242‑M1/1‑107 (MQ=255)
                                aCGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATCCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcaaacttttgtttttgttgtatataaatct                                                                                                           >  2:33321‑M1/1‑106 (MQ=255)
                                      cTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaa                                                                                                     <  1:6077‑M1/143‑44 (MQ=255)
                                       tGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaa                                                                                                    >  1:346078‑M1/1‑99 (MQ=255)
                                        gTACTCGTCTCTACACGCACGTACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACAAAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaa                                                                                                   <  2:1261094‑M1/143‑46 (MQ=255)
                                         tACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttttttgggtgtatttccatctgttaaaaagt                                                                                                  >  1:759304‑M1/1‑97 (MQ=255)
                                         tACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaa                                                                                                   >  1:1210928‑M1/1‑97 (MQ=255)
                                         tACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaa                                                                                                   <  2:1210928‑M1/142‑46 (MQ=255)
                                         tACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctcttaaaaag                                                                                                  >  2:500654‑M1/1‑97 (MQ=255)
                                            tCGTCAGTACACGCACTTACAACGGAGGTGCGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACTCTCAATTACTTAACTTAAACGtataagtcatacttttgttttgggtgtatttcaatctgtgaaaaagttt                                                                                               <  2:371520‑M1/143‑50 (MQ=255)
                                            tCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgaatttcaatctgttaaaaagttt                                                                                               >  1:448446‑M1/1‑94 (MQ=255)
                                              gTCACTACACGCACCTGCAACGGAGGCGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTGAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagttttt                                                                                             <  2:520701‑M1/143‑52 (MQ=255)
                                                  cTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgct                                                                                         >  1:477862‑M1/1‑88 (MQ=255)
                                                        gCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgcta                                                                                   >  2:114463‑M1/1‑82 (MQ=255)
                                                          acaTACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagttgttcgcttcgctagc                                                                                 >  2:852493‑M1/1‑80 (MQ=255)
                                                             tACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaag                                                                              <  2:247179‑M1/143‑67 (MQ=255)
                                                               cAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagca                                                                            >  1:263929‑M1/1‑75 (MQ=255)
                                                                   ggAGGTGCGCTTCGATTTTCAATAATGCGTGATGCATATCACATAAAGCACTCAATTACTTAACATAAATGtataagtcatacttttgtcttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaa                                                                        <  2:520324‑M1/143‑73 (MQ=255)
                                                                   ggAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaa                                                                        <  1:1047981‑M1/143‑73 (MQ=255)
                                                                   ggAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaa                                                                        <  1:25332‑M1/143‑73 (MQ=255)
                                                                   ggAGGGGGGCTGCACTTTTCAACAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgcgtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaa                                                                        <  2:621368‑M1/143‑73 (MQ=255)
                                                                    gAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGgataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaat                                                                       <  1:138081‑M1/143‑74 (MQ=255)
                                                                        ggggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaa                                                                   <  1:1293057‑M1/143‑78 (MQ=255)
                                                                          ggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaaca                                                                 <  1:1215019‑M1/143‑80 (MQ=255)
                                                                           gCTGCGATTTTCAATAATGCGTCATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacag                                                                <  1:980149‑M1/143‑81 (MQ=255)
                                                                            cTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatatgttaaaaagtttttcgctacgctagcaagccaaaatgaaatagg                                                               >  2:263353‑M1/1‑62 (MQ=255)
                                                                                   tttCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaccgtttatgtgcacgcgagcaagcaaaaatgaaacaggaataatc                                                        >  2:271345‑M1/1‑55 (MQ=255)
                                                                                     tCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagcttctgtcgacgcgagcaagcaaaaatgaaacaggaataatcga                                                      >  2:272223‑M1/1‑53 (MQ=255)
                                                                                         taatGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggactcatcgaaatg                                                  >  1:326789‑M1/1‑49 (MQ=255)
                                                                                         taagGCGTTATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttgtgttttgggtgtatttcaatctggtaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatg                                                  <  1:39227‑M1/143‑95 (MQ=255)
                                                                                          aatGCGTCATGCAGATCACGTAGAACACTCAAGTACTTAACGTAAATGtataagtcatacttttgttttgggtgtgtttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgg                                                 <  2:661478‑M1/143‑96 (MQ=255)
                                                                                                   aTGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgc                                         >  1:976131‑M1/1‑40 (MQ=255)
                                                                                                     gCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttctctacgctagcaagcaatcatgaaacaggaatcatcgaaatgggatactgcgc                                       >  2:1020145‑M1/1‑38 (MQ=255)
                                                                                                     gCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgc                                       >  2:882346‑M1/1‑38 (MQ=255)
                                                                                                       aGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcac                                     >  2:502281‑M1/1‑36 (MQ=255)
                                                                                                        gATTACACAAAAGACTCAATTACTTAACATAAATGtataagtcatacttttgttttggggggatttcaatctgttaaaaagtttttcgctacgcgagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcaca                                    <  2:554637‑M1/143‑109 (MQ=255)
                                                                                                        gATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcaca                                    <  2:706081‑M1/143‑109 (MQ=255)
                                                                                                        gATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcaca                                    <  1:778033‑M1/143‑109 (MQ=255)
                                                                                                          tCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagt                                  <  1:344713‑M1/143‑111 (MQ=255)
                                                                                                           cacacaAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtc                                 <  1:922488‑M1/143‑112 (MQ=255)
                                                                                                              acaAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaa                              >  1:1201659‑M1/1‑29 (MQ=255)
                                                                                                                aaaaCACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctaagctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcaaagtcaaaat                            >  1:976328‑M1/1‑27 (MQ=255)
                                                                                                                 aaaCACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaata                           >  1:726642‑M1/1‑26 (MQ=255)
                                                                                                                         aaTTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccg                   >  2:862537‑M1/1‑18 (MQ=255)
                                                                                                                           ttACGTAACATAAATGtataagtcatacttttgttttgggtgtagttcaatctgttaaaaagtgtttcgctacgctagcaagcaaaaatgaaccaagaatgatcgaaatgggatgttgcgcacagtcaaaataactcaccgta                 <  2:1076361‑M1/143‑128 (MQ=255)
                                                                                                                            tACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaacagttcttcgctacgctagcatgcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaa                >  1:971041‑M1/1‑15 (MQ=255)
                                                                                                                            tACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtccaaataactcaccgtaa                >  1:763756‑M1/1‑15 (MQ=255)
                                                                                                                              cTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatctaaatgggatgttgcgcacagtcaaaataactcaccataaat              >  1:590472‑M1/1‑13 (MQ=255)
                                                                                                                              cTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaat              >  1:1215701‑M1/1‑13 (MQ=255)
                                                                                                                               ttAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaat                                                                       <  1:176850‑M1/85‑74 (MQ=255)
                                                                                                                               ttAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaat                                                                       >  2:176850‑M1/1‑12 (MQ=255)
                                                                                                                                    aTAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagaaacaaacaacaagtgagaagcgaaatgggatgttgcgcacagtcaaaataactcaccgtaactaatcat        >  1:609107‑M1/1‑7 (MQ=255)
                                                                                                                                       aaTGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctg     >  2:849203‑M1/1‑4 (MQ=255)
                                                                                                                                          gtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctatcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgcta  >  1:1266816‑M1/1‑1 (MQ=255)
                                                                                                                                          gtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccctaaacaagcatccgctg  >  1:155765‑M1/1‑1 (MQ=255)
                                                                                                                                          gtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaaaaggaataatcgaaatgggatgttgcggacagacaaaataactcgcggtaagtaatcatatggta  >  1:742535‑M1/1‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                              
CTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGC  >  NC_000913/1978365‑1978644

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