Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 4,091,572 3461 bp→GCATGCGGATCCTCTAGAGTCGACATATG [rhaD]rhaB [rhaD], rhaA, rhaB

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 4091572 4095032 3461 279 [0] [0] 255 [rhaD]–rhaB [rhaD],rhaA,rhaB

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 40915710 (0.000)176 (0.980)
+GCATGCGGATCCTCTAGAGTCGACATATG
81/126 0.6 100% coding (283/825 nt) rhaD rhamnulose‑phosphate aldolase
?NC_000913 4095033 = 0 (0.000)intergenic (‑4/‑284) rhaB/rhaS rhamnulokinase/RhaS

Only 100 of 279 total aligned reads displayed.

GTGGGAAAGGAAGTGAGCCGGAAGTTCGGAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACTTTTACGATGCCTAAGTTAGCCGCAGGATCAAGCTGGACGTTACGGAAGAATTTGCCCGAGCCGGTGACAATAAACGGTGTATTTGCCAGTA  >  NC_000913/4091479‑4091663
                                                                                            |                                                                                            
gTGGGAAAGGAAGTGAGCCGGAAGTTCGGAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGAc                                                                                              <  1:1829225/93‑1 (MQ=255)
  gggAAAGGAAGTGAGCCGGAAGTTCGGAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgc                                                                                            <  2:4788998‑M1/93‑3 (MQ=255)
  gggAAAGGAAGTGAGCCGGAAGTTCGGAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgc                                                                                            >  2:5196634‑M1/1‑91 (MQ=255)
   ggAAAGGAAGTGAGCCGGAAGTTCGGAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgca                                                                                           >  1:3823541‑M1/1‑90 (MQ=255)
    gAAAGGAAGTGAGCCGGAAGTTCGGAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcat                                                                                          <  1:6325069‑M1/93‑5 (MQ=255)
      aaGGAAGTGAGCCGGAAGTTCGGAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgc                                                                                        >  2:4318986‑M1/1‑87 (MQ=255)
      aaGGAAGTGAGCCGGAAGTTCGGAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgc                                                                                        >  1:1563539‑M1/1‑87 (MQ=255)
      aaGGAAGTGAGCCGGAAGTTCGGAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgc                                                                                        <  1:342305‑M1/93‑7 (MQ=255)
      aaGGAAGTGAGCCGCAAGTTCGGAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgc                                                                                        <  1:1251299‑M1/93‑7 (MQ=255)
       aGGAAGTGAGCCGGAAGTTCGGAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcg                                                                                       >  1:1078112‑M1/1‑86 (MQ=255)
       aGGAAGTGAGCCGGAAGTTCGGAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcg                                                                                       >  2:170996‑M1/1‑86 (MQ=255)
        ggAAGTGAGCCGGAAGTTCGGAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcgg                                                                                      >  1:1295023‑M1/1‑85 (MQ=255)
         gAAGTGAGCCGGAAGTTCGGAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcgga                                                                                     >  1:140618‑M1/1‑84 (MQ=255)
          aaGTGAGCCGGAAGTTCGGAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggat                                                                                    >  2:6389989‑M1/1‑83 (MQ=255)
           aGTGAGCCGGAAGTTCGGAAGTGGGGGCGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatc                                                                                   <  1:3171520‑M1/93‑12 (MQ=255)
           aGTGAGCCGGAAGTTCGGAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatc                                                                                   >  2:4936140‑M1/1‑82 (MQ=255)
           aGTGAGCCGGAAGTTCGGAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatc                                                                                   <  2:2081913‑M1/93‑12 (MQ=255)
           aGTGAGCCGGAAGTTCGGAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatc                                                                                   >  2:3252742‑M1/1‑82 (MQ=255)
            gTGAGCCGGAAGTTCGGAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcc                                                                                  <  1:2126168‑M1/93‑13 (MQ=255)
             tGAGCCGGAAGTTCGGAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcct                                                                                 <  2:6654462‑M1/93‑14 (MQ=255)
             tGAGCCGGAAGTTCGGAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcct                                                                                 <  1:1721783‑M1/93‑14 (MQ=255)
              gAGCCGGAAGTTCGGAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctc                                                                                <  1:4629342‑M1/93‑15 (MQ=255)
               aGCCGGAAGTTCGGAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctct                                                                               <  1:6027512‑M1/93‑16 (MQ=255)
               aGCCGGAAGTTCGGAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctct                                                                               >  1:5006880‑M1/1‑78 (MQ=255)
                gCCGGAAGTTCGGAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctcta                                                                              <  2:6140202‑M1/93‑17 (MQ=255)
                   ggAAGTTCGGAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagag                                                                           >  2:3882157‑M1/1‑74 (MQ=255)
                    gAAGTTCGGAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagt                                                                          >  2:2601732‑M1/1‑73 (MQ=255)
                    gAAGTTCGGAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagt                                                                          >  1:6191356‑M1/1‑73 (MQ=255)
                    gAAGTTCGGAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagt                                                                          <  2:1981708‑M1/93‑21 (MQ=255)
                     aaGTTCGGAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtc                                                                         <  1:3970784‑M1/93‑22 (MQ=255)
                      aGTTCGGAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcg                                                                        >  2:4568086‑M1/1‑71 (MQ=255)
                       gTTCGGAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgc                                                                       >  1:3290601‑M1/1‑70 (MQ=255)
                         tCGGAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgaca                                                                     >  1:2849737‑M1/1‑68 (MQ=255)
                          cGGAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGACTGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacat                                                                    >  2:2263547‑M1/1‑67 (MQ=255)
                          cGGAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacat                                                                    <  1:3641548‑M1/93‑27 (MQ=255)
                           ggAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacata                                                                   <  2:5020479‑M1/93‑28 (MQ=255)
                           ggAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacata                                                                   <  1:3342517‑M1/93‑28 (MQ=255)
                            gAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatat                                                                  >  2:6415746‑M1/1‑65 (MQ=255)
                            gAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCACCGTCGCCgcatgccgatcctctagagtcgacatat                                                                  >  2:607445‑M1/1‑65 (MQ=255)
                              aGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatgg                                                                >  2:4119324‑M1/1‑63 (MQ=255)
                               gTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGCGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatctcctagaggtcacatatgtt                                                               >  1:2312959‑M1/1‑62 (MQ=255)
                                  gggACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgat                                                            >  1:5712665‑M1/1‑59 (MQ=255)
                                  gggACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagacgacatatggtgat                                                            >  2:365551‑M1/1‑59 (MQ=255)
                                    gACGGCTTCGTTGGTTACCCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcc                                                          <  1:6399819‑M1/93‑37 (MQ=255)
                                    gACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatgctctagaagcgacatatggtgatcc                                                          >  1:4402900‑M1/1‑57 (MQ=255)
                                    gACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcc                                                          <  2:3772132‑M1/93‑37 (MQ=255)
                                    gACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcc                                                          >  1:6556237‑M1/1‑57 (MQ=255)
                                    gACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcc                                                          >  1:198085‑M1/1‑57 (MQ=255)
                                     aCGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcct                                                         <  1:3506131‑M1/93‑38 (MQ=255)
                                      cGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctg                                                        >  2:1253308‑M1/1‑55 (MQ=255)
                                       ggCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgc                                                       >  1:2627900‑M1/1‑54 (MQ=255)
                                        gCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgct                                                      >  2:3950731‑M1/1‑53 (MQ=255)
                                        gCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgct                                                      >  2:799521‑M1/1‑53 (MQ=255)
                                          ttCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctga                                                    >  1:2208799‑M1/1‑51 (MQ=255)
                                            cGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaat                                                  <  1:6145754‑M1/93‑45 (MQ=255)
                                             gTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatt                                                 >  2:6569138‑M1/1‑48 (MQ=255)
                                             gTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatt                                                 >  1:293903‑M1/1‑48 (MQ=255)
                                             gTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatt                                                 <  1:560562‑M1/93‑46 (MQ=255)
                                              ttGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaattt                                                <  2:2402133‑M1/93‑47 (MQ=255)
                                                  ttAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcatt                                            <  1:5565283‑M1/93‑51 (MQ=255)
                                                      cccccAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacga                                        >  2:1065066‑M1/1‑39 (MQ=255)
                                                      cccccAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacga                                        >  1:4069704‑M1/1‑39 (MQ=255)
                                                         ccAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgacca                                     >  2:1110477‑M1/1‑36 (MQ=255)
                                                         ccAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgacca                                     >  2:5051445‑M1/1‑36 (MQ=255)
                                                         ccAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgacca                                     <  1:4672138‑M1/93‑58 (MQ=255)
                                                           aaaGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgaccagt                                   <  1:5460367‑M1/93‑60 (MQ=255)
                                                            aaGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgaccagtc                                  <  2:3615868‑M1/93‑61 (MQ=255)
                                                             aGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgaccagtct                                 <  1:976706‑M1/93‑62 (MQ=255)
                                                             aGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgaccagtct                                 <  2:338648‑M1/93‑62 (MQ=255)
                                                              gAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgaccagtcta                                >  1:4453583‑M1/1‑31 (MQ=255)
                                                              gAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgaccagtcta                                >  2:5849016‑M1/1‑31 (MQ=255)
                                                                  gtgGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgaccagtctaaaaa                            <  1:5930133‑M1/93‑67 (MQ=255)
                                                                  gtgGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgaccagtctaaaaa                            <  2:5917188‑M1/93‑67 (MQ=255)
                                                                  gtgGTAGCCCGCGCCGTAGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgaccagtctaaaaa                            <  2:3545068‑M1/93‑67 (MQ=255)
                                                                    ggTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgaccagtctaaaaagc                          >  2:1583648‑M1/1‑25 (MQ=255)
                                                                     gTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgaccagtctaaaaagcg                         <  1:2521557‑M1/93‑70 (MQ=255)
                                                                       aGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgaccagtctaaaaagcgcc                       >  2:4062443‑M1/1‑22 (MQ=255)
                                                                        gCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgaccagtctaaaaagcgcct                      <  2:3517192‑M1/93‑73 (MQ=255)
                                                                        gCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgaccagtctaaaaagcgcct                      >  2:2966970‑M1/1‑21 (MQ=255)
                                                                          ccGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgaccagtctaaaaagcgcctga                    >  2:2928267‑M1/1‑19 (MQ=255)
                                                                          ccGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgaccagtctaaaaagcgcctga                    >  1:3726007‑M1/1‑19 (MQ=255)
                                                                            gcgcCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgaccagtctaaaaagcgcctgaat                  >  1:2071776‑M1/1‑17 (MQ=255)
                                                                            gcgcCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgaccagtctaaaaagcgcctgaat                  <  2:3061228‑M1/93‑77 (MQ=255)
                                                                              gcCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgaccagtctaaaaagcgcctgaattc                <  2:750127‑M1/93‑79 (MQ=255)
                                                                              gcCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgaccagtctaaaaagcgcctgaattc                >  2:1388088‑M1/1‑15 (MQ=255)
                                                                              cccgTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgaccagtctaaaaagcgcctgaattc                <  1:4245019‑M1/92‑79 (MQ=255)
                                                                                   cGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgaccagtctaaaaagcgcctgaattcgcgac           >  1:6113063‑M1/1‑10 (MQ=255)
                                                                                   cGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgaccagtctaaaaagcgcctgaattcgcgac           >  1:860266‑M1/1‑10 (MQ=255)
                                                                                   cGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgaccagtctaaaaagcgcctgaattcgcgac           >  1:4829656‑M1/1‑10 (MQ=255)
                                                                                      tGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgaccagtctaaaaagcgcctgaattcgcgacctt        <  1:6466549‑M1/93‑87 (MQ=255)
                                                                                      tGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgaccagtctaaaaagcgcctgaattcgcgacctt        >  1:2403971‑M1/1‑7 (MQ=255)
                                                                                      tGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgaccagtctaaaaagcgcctgaattcgcgacctt        <  1:6083139‑M1/93‑87 (MQ=255)
                                                                                      tGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgaccagtctaaaaagcgcctgaattcgcgacctt        >  1:690284‑M1/1‑7 (MQ=255)
                                                                                      tGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgaccagtctaaaaagcgcctgaattcgcgacctt        >  1:5208029‑M1/1‑7 (MQ=255)
                                                                                       gTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgaccagtctaaaaagcgcctgaattcgcgaccttc       >  2:5269206‑M1/1‑6 (MQ=255)
                                                                                        tCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgaccagtctaaaaagcgcctgaattcgcgaccttct      >  2:1842225‑M1/1‑5 (MQ=255)
                                                                                        tCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgaccagtctaaaaagcgcctgaattcgcgaccttct      <  2:4425574‑M1/93‑89 (MQ=255)
                                                                                         cGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgaccagtctaaaaagcgcctgaattcgcgaccttctc     <  2:2361172‑M1/93‑90 (MQ=255)
                                                                                          gACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgaccagtctaaaaagcgcctgaattcgcgaccttctcg    <  1:2750667‑M1/93‑91 (MQ=255)
                                                                                            cgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgaccagtctaaaaagcgcctgaattcgcgaccttctcgtt  >  2:5683421‑M1/1‑1 (MQ=255)
                                                                                            |                                                                                            
GTGGGAAAGGAAGTGAGCCGGAAGTTCGGAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACTTTTACGATGCCTAAGTTAGCCGCAGGATCAAGCTGGACGTTACGGAAGAATTTGCCCGAGCCGGTGACAATAAACGGTGTATTTGCCAGTA  >  NC_000913/4091479‑4091663

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