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 229 [0] [0] 231 [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)162 (1.310)
+GCATGCGGATCCTCTAGAGTCGACATATG
85/126 0.0 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 229 total aligned reads displayed.

GTGGGAAAGGAAGTGAGCCGGAAGTTCGGAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACTTTTACGATGCCTAAGTTAGCCGCAGGATCAAGCTGGACGTTACGGAAGAATTTGCCCGAGCCGGTGACAATAAACGGTGTATTTGCCAG  >  NC_000913/4091479‑4091661
                                                                                            |                                                                                          
gTGGGAAAGGAAGTGAGCCGGAAGTTCGGAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGAc                                                                                            >  2:2367802/1‑93 (MQ=255)
 tGGGAAAGGAAGTGAGCCGGAAGTTCGGAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACg                                                                                           >  1:2485177‑M1/1‑92 (MQ=255)
  gggAAAGGAAGTGAGCCGGAAGTTCGGAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGGCGACgc                                                                                          >  1:3077814‑M1/1‑91 (MQ=255)
    gAAAGGAAGTGAGCCGGAAGTTCGGAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcat                                                                                        >  2:1131999‑M1/1‑89 (MQ=255)
     aaaGGAAGTGAGCCGGAAGTTCGGAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatg                                                                                       <  2:2244196‑M1/93‑6 (MQ=255)
       aGGAAGTGAGCCGGAAGTTCGGAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcg                                                                                     <  1:3421380‑M1/93‑8 (MQ=255)
        ggAAGTGAGCCGGAAGTTCGGAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcgg                                                                                    >  2:4604895‑M1/1‑85 (MQ=255)
           aGTGAGCCGGAAGTTCGGAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatc                                                                                 <  1:2586797‑M1/93‑12 (MQ=255)
              gAGCCGGAAGTTCGGAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctc                                                                              <  1:3723502‑M1/93‑15 (MQ=255)
              gAGCCGGAAGTTCGGAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctc                                                                              <  1:3447241‑M1/93‑15 (MQ=255)
                gCCGGAAGTTCGGAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctcta                                                                            <  1:311669‑M1/93‑17 (MQ=255)
                gCCGGAAGTTCGGAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGAAgcatgcggatcctcta                                                                            >  2:804426‑M1/1‑77 (MQ=255)
                gCCGGAAGTGCGGAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctcta                                                                            <  1:93040‑M1/93‑17 (MQ=255)
                  cGGAAGTTCGGAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctaga                                                                          >  1:2945983‑M1/1‑75 (MQ=255)
                  cGGAAGTTCGGAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctaga                                                                          <  2:2663172‑M1/93‑19 (MQ=255)
                  cGGAAGTTCGGAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctaga                                                                          >  1:2476982‑M1/1‑75 (MQ=255)
                   ggAAGTTCGGAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagag                                                                         >  2:2422720‑M1/1‑74 (MQ=255)
                   ggAAGTTCGGAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagag                                                                         >  1:326328‑M1/1‑74 (MQ=255)
                    gAAGTTCGGAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagt                                                                        >  1:2081438‑M1/1‑73 (MQ=255)
                    gAAGTTCGGAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagc                                                                        >  2:3335225‑M1/1‑73 (MQ=255)
                       gTTCGGAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcga                                                                     <  2:4130446‑M1/93‑24 (MQ=255)
                         tCGGAAGTGGGGACGGCTTCGTTGGTTAACCCGCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgaca                                                                   <  2:3211361‑M1/93‑26 (MQ=255)
                         tCGGAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgaca                                                                   <  1:2564446‑M1/93‑26 (MQ=255)
                         tCGGAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgaca                                                                   <  1:4529733‑M1/93‑26 (MQ=255)
                           ggAAGTGGGGACGGCTTCGTTTGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacata                                                                 <  2:1585462‑M1/93‑28 (MQ=255)
                           ggAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacata                                                                 <  1:89697‑M1/93‑28 (MQ=255)
                            gAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGAGgcatgcggatcctctagaggcgagatat                                                                >  1:1239954‑M1/1‑65 (MQ=255)
                            gAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatat                                                                >  1:4576181‑M1/1‑65 (MQ=255)
                            gAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatat                                                                <  2:1482311‑M1/93‑29 (MQ=255)
                            gAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatat                                                                <  2:187393‑M1/93‑29 (MQ=255)
                            gAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatat                                                                <  1:1768651‑M1/93‑29 (MQ=255)
                                  gggACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgat                                                          >  1:2039085‑M1/1‑59 (MQ=255)
                                    gACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcc                                                        <  1:1560280‑M1/93‑37 (MQ=255)
                                    gACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcc                                                        <  2:4171857‑M1/93‑37 (MQ=255)
                                    gACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagccgacatatggtgatcc                                                        >  1:4721008‑M1/1‑57 (MQ=255)
                                    gACGGCTTCGCTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcc                                                        >  2:3642070‑M1/1‑57 (MQ=255)
                                     aCGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggaccctctagagcctacatatggtgatcct                                                       >  1:4222571‑M1/1‑56 (MQ=255)
                                      cGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctg                                                      >  1:1598324‑M1/1‑55 (MQ=255)
                                       ggCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgc                                                     >  2:93964‑M1/1‑54 (MQ=255)
                                        gCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcgggtcctctagagtcgacatatggtgatcctgcg                                                    >  1:118200‑M1/1‑53 (MQ=255)
                                        gCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgct                                                    <  1:3546251‑M1/93‑41 (MQ=255)
                                        gCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgct                                                    <  1:2777844‑M1/93‑41 (MQ=255)
                                           tCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaa                                                 >  1:2559185‑M1/1‑50 (MQ=255)
                                             gTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatt                                               >  1:3752809‑M1/1‑48 (MQ=255)
                                             gTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatt                                               >  2:1936083‑M1/1‑48 (MQ=255)
                                              ttGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaattt                                              >  2:3384183‑M1/1‑47 (MQ=255)
                                              ttGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaattt                                              >  1:2279882‑M1/1‑47 (MQ=255)
                                                ggTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttca                                            <  1:3672427‑M1/93‑49 (MQ=255)
                                                ggTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttca                                            <  1:1408142‑M1/93‑49 (MQ=255)
                                                ggTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttca                                            >  1:1847881‑M1/1‑45 (MQ=255)
                                                 gTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcat                                           <  1:3936351‑M1/93‑50 (MQ=255)
                                                 gTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcat                                           <  2:4515760‑M1/93‑50 (MQ=255)
                                                  ttAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcatt                                          <  1:2318907‑M1/93‑51 (MQ=255)
                                                   tAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattc                                         >  2:4577454‑M1/1‑42 (MQ=255)
                                                   tAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcatta                                         >  1:4469141‑M1/1‑42 (MQ=255)
                                                   tAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcatta                                         >  1:22308‑M1/1‑42 (MQ=255)
                                                    aaCCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattac                                        >  2:4698011‑M1/1‑41 (MQ=255)
                                                     aCCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacg                                       >  1:2866056‑M1/1‑40 (MQ=255)
                                                     aCCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacg                                       >  2:1708619‑M1/1‑40 (MQ=255)
                                                     aCCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacg                                       >  2:3988634‑M1/1‑40 (MQ=255)
                                                      cccccAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacga                                      >  2:2603791‑M1/1‑39 (MQ=255)
                                                       ccccAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgac                                     >  1:1705445‑M1/1‑38 (MQ=255)
                                                       ccccAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgac                                     <  2:778913‑M1/93‑56 (MQ=255)
                                                       ccccAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgac                                     >  1:2556457‑M1/1‑38 (MQ=255)
                                                        cccAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgacc                                    >  1:494203‑M1/1‑37 (MQ=255)
                                                        cccAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgacc                                    >  1:4271466‑M1/1‑37 (MQ=255)
                                                        cccAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgacc                                    <  2:3666486‑M1/93‑57 (MQ=255)
                                                          cAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgaccag                                  <  2:640392‑M1/93‑59 (MQ=255)
                                                          cAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgaccag                                  >  2:4421470‑M1/1‑35 (MQ=255)
                                                           aaaGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgaccagt                                 >  2:2740474‑M1/1‑34 (MQ=255)
                                                           aaaGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgaccagt                                 >  2:841774‑M1/1‑34 (MQ=255)
                                                            aaGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgaccagtc                                <  2:1270885‑M1/93‑61 (MQ=255)
                                                                aTGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgaccagtctaaa                            <  1:238544‑M1/93‑65 (MQ=255)
                                                                   tgGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgaccagtctaaaaag                         >  2:1719706‑M1/1‑26 (MQ=255)
                                                                    ggTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgaccagtctaaaaagc                        <  1:2191712‑M1/93‑69 (MQ=255)
                                                                       aGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgaccagtctaaaaagcgcc                     >  1:526980‑M1/1‑22 (MQ=255)
                                                                       aGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgaccagtctaaaaagcgcc                     >  1:657059‑M1/1‑22 (MQ=255)
                                                                       aGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgaccagtctaaaaagcgcc                     >  1:2941447‑M1/1‑22 (MQ=255)
                                                                          ccGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatagggtgatcctgctgaatttcattacgaccagtctaaaaagcgcctga                  >  1:1831973‑M1/1‑19 (MQ=255)
                                                                           cgcgCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgaccagtctaaaaagcgcctgaa                 >  2:3852831‑M1/1‑18 (MQ=255)
                                                                            gcgcCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgaccagtctaaaaagcgcctgaat                >  2:2039963‑M1/1‑17 (MQ=255)
                                                                             cgcCGTCGCTGTCGACgcatgcggatcctctagagtcgagatatggtgatcctgctgaatttcattacgaccagtctaaaaaccgcctgaaat               <  2:4403843‑M1/93‑78 (MQ=255)
                                                                              gcCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgaccagtctaaaaagcgcctgaattc              >  1:3599784‑M1/1‑15 (MQ=255)
                                                                               ccGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgaccagtctaaaaagcgcctgaattcg             >  1:4543540‑M1/1‑14 (MQ=255)
                                                                                cGTCGCTGGCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgaccagtctaaaaagcgcctgaattcgc            <  1:1292729‑M1/93‑81 (MQ=255)
                                                                                 gTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgattctgctgaatttcattacgaccagtctaaaaagcgcctgaattcgcg           >  2:141451‑M1/1‑12 (MQ=255)
                                                                                 gTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgaccagtctaaaaagcgcctgaattcgcg           <  1:3361368‑M1/93‑82 (MQ=255)
                                                                                  tCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgaccagtctaaaaagcgcctgaattcgcga          >  2:2932332‑M1/1‑11 (MQ=255)
                                                                                  tCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgaccagtctaaaaagcgcctgaattcgcga          >  1:2795281‑M1/1‑11 (MQ=255)
                                                                                    gCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgaccagtctaaaaagcgcctgaattcgcgacc        >  1:823205‑M1/1‑9 (MQ=255)
                                                                                    gCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgaccagtctaaaaagcgcctgaattcgcgacc        >  2:2487188‑M1/1‑9 (MQ=255)
                                                                                    gCTGTCGACgcatgcggatcctctacagtcgacatatggtgatcctgctgaatttcattacgaccagtctaaaaagcgcctgaatttgcgact        >  2:4059531‑M1/1‑9 (MQ=255)
                                                                                     cTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgaccagtctaaaaagcgcctgaattcgcgacct       <  1:677173‑M1/93‑86 (MQ=255)
                                                                                       gTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgaccagtctaaaaagcgcctgaattcgcgaccttc     >  2:3594565‑M1/1‑6 (MQ=255)
                                                                                       gTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgaccagtctaaaaagcgcctgaattcgcgaccttc     <  2:1732671‑M1/93‑88 (MQ=255)
                                                                                       gTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgaccagtctaaaaagcgcctgaattcgcgaccttc     <  2:1074053‑M1/93‑88 (MQ=255)
                                                                                        tCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgaccagtctaaaaagcgcctgaattcgcgaccttct    <  1:1983661‑M1/93‑89 (MQ=255)
                                                                                         cGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgaccagtctaaaaagcgcctgaattcgcgaccttctc   >  2:3451867‑M1/1‑4 (MQ=255)
                                                                                         cGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgaccagtctaaaaagcgcctgaattcgcgaccttctc   <  1:1000848‑M1/93‑90 (MQ=255)
                                                                                          gACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgaccagtctaaaaagcgtctgaattcgcgaccttctcg  >  1:4348447‑M1/1‑3 (MQ=255)
                                                                                            |                                                                                          
GTGGGAAAGGAAGTGAGCCGGAAGTTCGGAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACTTTTACGATGCCTAAGTTAGCCGCAGGATCAAGCTGGACGTTACGGAAGAATTTGCCCGAGCCGGTGACAATAAACGGTGTATTTGCCAG  >  NC_000913/4091479‑4091661

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: