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

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 1978503–1979270 1979278 9–776 173 [0] [1] 172 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)148 (1.200) 77/162 NT 100% intergenic (‑305/+16) flhD/insB1 flagellar class II regulon transcriptional activator, with FlhC/IS1 transposase B
?NC_000913 1979279 = 0 (0.000)intergenic (‑64/‑474) insA/uspC IS1 repressor TnpA/universal stress protein

Only 100 of 173 total aligned reads displayed.

CTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCG  >  NC_000913/1978403‑1978601
                                                                                                   |                                                                                                   
cTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATG                                                                                                     <  2:3710852‑M1/100‑1 (MQ=255)
cTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATG                                                                                                     >  2:971778‑M1/1‑100 (MQ=255)
  gTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGta                                                                                                   <  2:3367885‑M1/100‑3 (MQ=255)
       cGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTTAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataag                                                                                               >  1:1712770‑M1/1‑93 (MQ=255)
       cGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagt                                                                                              >  2:361370‑M1/1‑93 (MQ=255)
        gTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtc                                                                                             >  2:2142599‑M1/1‑92 (MQ=255)
          cACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcat                                                                                           <  1:1002680‑M1/100‑11 (MQ=255)
           aCTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcata                                                                                          <  1:3149343‑M1/100‑12 (MQ=255)
           aCTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcata                                                                                          >  2:2855473‑M1/1‑89 (MQ=255)
             tACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATTtataagtcataat                                                                                        <  2:2567463‑M1/100‑14 (MQ=255)
               cacGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttt                                                                                      <  2:1894805‑M1/100‑16 (MQ=255)
                 cGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttg                                                                                    <  2:1492887‑M1/100‑18 (MQ=255)
                   cacaTACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgtt                                                                                  >  1:960884‑M1/1‑81 (MQ=255)
                   cacaTACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgtt                                                                                  >  1:2255145‑M1/1‑81 (MQ=255)
                   cacaTACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgtt                                                                                  <  1:1981199‑M1/100‑20 (MQ=255)
                     caTACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgtttt                                                                                >  2:2430029‑M1/1‑79 (MQ=255)
                     caTACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgtttt                                                                                <  1:222055‑M1/100‑22 (MQ=255)
                      aTACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttg                                                                               <  2:1146768‑M1/100‑23 (MQ=255)
                       tACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgg                                                                              >  2:2388355‑M1/1‑77 (MQ=255)
                        aCAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttggg                                                                             <  2:3576029‑M1/100‑25 (MQ=255)
                        aCAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttggg                                                                             <  1:2407350‑M1/100‑25 (MQ=255)
                        aCAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttggg                                                                             <  2:3447509‑M1/100‑25 (MQ=255)
                        aCAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttggg                                                                             <  2:39523‑M1/100‑25 (MQ=255)
                        aCAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttggg                                                                             <  1:2695519‑M1/100‑25 (MQ=255)
                         cAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACCAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggt                                                                            >  1:294739‑M1/1‑75 (MQ=255)
                         cAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggt                                                                            >  2:558816‑M1/1‑75 (MQ=255)
                         cAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggt                                                                            >  1:801651‑M1/1‑75 (MQ=255)
                         cAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacgtttgttttgggt                                                                            <  2:2511828‑M1/100‑26 (MQ=255)
                               aGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtattt                                                                      <  2:2954532‑M1/100‑32 (MQ=255)
                                gggtggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttc                                                                     <  1:541453‑M1/100‑33 (MQ=255)
                                ggggggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttc                                                                     <  1:2125623‑M1/100‑33 (MQ=255)
                                 gggggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttca                                                                    <  2:3745702‑M1/100‑34 (MQ=255)
                                  ggggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaa                                                                   <  1:1620483‑M1/100‑35 (MQ=255)
                                    ggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaa                                                                   <  1:1075514‑M1/98‑35 (MQ=255)
                                     gCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatct                                                                >  2:44661‑M1/1‑63 (MQ=255)
                                     gCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatct                                                                >  2:889430‑M1/1‑63 (MQ=255)
                                     gCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaa                                                                   <  1:1180048‑M1/97‑35 (MQ=255)
                                       tGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgt                                                              >  2:2649642‑M1/1‑61 (MQ=255)
                                       tGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgt                                                              >  2:3255915‑M1/1‑61 (MQ=255)
                                       tGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgt                                                              >  2:104914‑M1/1‑61 (MQ=255)
                                       tGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgtattgggtgtatttcaatctgt                                                              >  1:940250‑M1/1‑61 (MQ=255)
                                        gCGATTTTCAATAATGCGTGATGCAGATCACACACAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgtt                                                             <  1:2007927‑M1/100‑41 (MQ=255)
                                        gCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgtt                                                             >  1:3772538‑M1/1‑60 (MQ=255)
                                        gCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgtt                                                             >  1:2557328‑M1/1‑60 (MQ=255)
                                         cGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgtta                                                            >  1:1241905‑M1/1‑59 (MQ=255)
                                         cGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgtta                                                            <  1:2941462‑M1/100‑42 (MQ=255)
                                         cGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgtta                                                            >  1:1166895‑M1/1‑59 (MQ=255)
                                          gATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaa                                                           >  2:287352‑M1/1‑58 (MQ=255)
                                            ttttCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaa                                                         >  1:2417875‑M1/1‑56 (MQ=255)
                                             tttCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaagat                                                        >  2:2235076‑M1/1‑55 (MQ=255)
                                              ttCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaag                                                       <  1:543337‑M1/100‑47 (MQ=255)
                                              ttCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaag                                                       >  2:1623680‑M1/1‑54 (MQ=255)
                                               tCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagt                                                      >  2:2491811‑M1/1‑53 (MQ=255)
                                                cAATAATGCGTGATGCAGATCCCACAAAACACTCAATTACTTAACATAAATGtataagtcatcattttgttttgggtgtattttaatttgttaaaaagt                                                      >  1:1017495‑M1/1‑52 (MQ=255)
                                                 aataatGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagttt                                                    >  2:2667658‑M1/1‑51 (MQ=255)
                                                    aatGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttgaatctgttaaaaagtttttc                                                 <  2:3590544‑M1/100‑53 (MQ=255)
                                                     atGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcg                                                >  2:2402057‑M1/1‑47 (MQ=255)
                                                     atGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatactctcgttttgggtctagttcaatctcttaaacagtctttcg                                                >  1:1982397‑M1/1‑47 (MQ=255)
                                                      tGCGTGATGCATATCACACAAAACACTCAATTACTTAACATAATTTtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgc                                               <  1:739159‑M1/100‑55 (MQ=255)
                                                      tGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgc                                               <  2:1401866‑M1/100‑55 (MQ=255)
                                                      tGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtattgcaatctgttaaaaagtttttcgc                                               <  1:2530274‑M1/100‑55 (MQ=255)
                                                        cGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgcta                                             <  1:391174‑M1/100‑57 (MQ=255)
                                                         gTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctac                                            >  1:1369710‑M1/1‑43 (MQ=255)
                                                           gATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgc                                          <  1:1415915‑M1/100‑60 (MQ=255)
                                                           gATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgc                                          >  2:2440770‑M1/1‑41 (MQ=255)
                                                             tgaagatCACACAAAACACTCAATTACTTAACATAAATGtataagttatacttttgttttgtgtgtatttcaatctgttaaatagtttttcgctacgct                                         <  1:2892740‑M1/96‑61 (MQ=255)
                                                             tGCATAGCACACAAAACACTTAATTACTTAACATAAATGtataagtcatacttttgttttcggtgtatttcaatctgttaaaaagtttttcgatacgc                                          <  1:3266983‑M1/98‑60 (MQ=255)
                                                             tGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgcta                                        >  1:1171669‑M1/1‑39 (MQ=255)
                                                               cAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctggtaaaacgtttttcgctacgctagc                                      >  2:1098960‑M1/1‑37 (MQ=255)
                                                                aGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagca                                     <  2:117346‑M1/100‑65 (MQ=255)
                                                                aGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagca                                     >  1:2388966‑M1/1‑36 (MQ=255)
                                                                  aTCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttctctacgctagcaag                                   >  1:2822743‑M1/1‑34 (MQ=255)
                                                                    cacacaAAACACTCAATTACTTATCATACATGtagaaggcatactcttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaag                                   <  1:1727256‑M1/98‑67 (MQ=255)
                                                                     ncacaaAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaa                                >  1:2899789‑M1/2‑31 (MQ=255)
                                                                      cacaAAACACTCTATTACTTAACATAAATGtataagtcatatttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaa                               <  1:361548‑M1/100‑71 (MQ=255)
                                                                        caAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaa                             >  1:3426294‑M1/1‑28 (MQ=255)
                                                                        caAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaa                             >  2:2476873‑M1/1‑28 (MQ=255)
                                                                         aaaaCACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaan                            <  2:2899803‑M1/100‑74 (MQ=255)
                                                                         aaaaCACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagatcggaa                            >  1:248061‑M1/1‑27 (MQ=255)
                                                                          aaaCACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatg                           <  2:2322370‑M1/100‑75 (MQ=255)
                                                                              acTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaac                       <  1:1765039‑M1/100‑79 (MQ=255)
                                                                              acTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaac                       <  1:997359‑M1/100‑79 (MQ=255)
                                                                              acTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaac                       <  2:3629972‑M1/100‑79 (MQ=255)
                                                                               cTCAATTACTTAACATAAATGtataagtcatacttttgtttttggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaac                       <  2:2737446‑M1/99‑79 (MQ=255)
                                                                               cTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaaca                      >  2:3144971‑M1/1‑21 (MQ=255)
                                                                                 cAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacagg                    <  2:2479649‑M1/100‑82 (MQ=255)
                                                                                 cAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacagg                    >  2:1488106‑M1/1‑19 (MQ=255)
                                                                                     tACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaata                <  1:261635‑M1/100‑86 (MQ=255)
                                                                                     tACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaata                >  2:2872651‑M1/1‑15 (MQ=255)
                                                                                         tAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcg            <  1:2833344‑M1/100‑90 (MQ=255)
                                                                                         tAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcg            >  2:1868699‑M1/1‑11 (MQ=255)
                                                                                              tAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatn       <  1:3496965‑M1/100‑95 (MQ=255)
                                                                                              tAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatg       >  2:1962375‑M1/1‑6 (MQ=255)
                                                                                              tAAATGtataagtcatacttttgttttgggggtatttcaatctggtaaaaagttgttcgctacgctagcaagaaaaaatgaaacaggaataatcgaaat        >  1:3750864‑M1/1‑6 (MQ=255)
                                                                                               aaaTGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgg      >  2:3079320‑M1/1‑5 (MQ=255)
                                                                                                aaTGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaacatgaaacaggaataatcgaaatggg     >  1:213295‑M1/1‑4 (MQ=255)
                                                                                                aaTGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatggg     <  1:1337500‑M1/100‑97 (MQ=255)
                                                                                                 aTGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatggga    <  2:2130029‑M1/100‑98 (MQ=255)
                                                                                                 aTGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatggga    >  2:581396‑M1/1‑3 (MQ=255)
                                                                                                   gtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatg  >  1:4346‑M1/1‑1 (MQ=255)
                                                                                                   |                                                                                                   
CTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCG  >  NC_000913/1978403‑1978601

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