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 231 total aligned reads displayed.

GGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGCGACACAATTGCGAAAGGTCATAATGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGGGCCATTGGCAACCA  >  NC_000913/4094942‑4095125
                                                                                           |                                                                                            
gtggggacggcttcgttggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGt                                                                                             <  1:3594623‑M2/2‑1 (MQ=255)
   gggacggcttcgttggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGAt                                                                                          >  1:2039085‑M2/89‑93 (MQ=255)
    ggacggcttcgttggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATc                                                                                         <  1:4386787‑M2/6‑1 (MQ=255)
     gacggcttcgttggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATcc                                                                                        <  1:1560280‑M2/7‑1 (MQ=255)
     gacggcttcgttggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATcc                                                                                        <  2:4171857‑M2/7‑1 (MQ=255)
     gacggcttcgttggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATcc                                                                                        >  1:3523710‑M2/87‑93 (MQ=255)
     gacggcttcgttggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagccgacatatgGTGATcc                                                                                        >  1:4721008‑M2/87‑93 (MQ=255)
     gacggcttcgctggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATcc                                                                                        >  2:3642070‑M2/87‑93 (MQ=255)
      acggcttcgttggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCct                                                                                       <  2:3646521‑M2/8‑1 (MQ=255)
         gcttcgttggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCctgct                                                                                    >  1:4590319‑M2/83‑93 (MQ=255)
          cttcgttggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCctgctg                                                                                   >  2:230739‑M2/82‑93 (MQ=255)
          cttcgttggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCctgctg                                                                                   >  2:2687814‑M2/82‑93 (MQ=255)
           ttcgttggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGa                                                                                  >  1:857602‑M2/81‑93 (MQ=255)
            tcgttggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGaa                                                                                 <  2:796678‑M2/14‑1 (MQ=255)
             tgttggtaaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAAt                                                                                <  2:3553010‑M2/15‑1 (MQ=255)
              gttggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAAtt                                                                               >  2:1936083‑M2/78‑93 (MQ=255)
               ttggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAAttt                                                                              >  1:2279882‑M2/77‑93 (MQ=255)
               ttggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAAttt                                                                              >  2:3384183‑M2/77‑93 (MQ=255)
                 ggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCa                                                                            >  1:1847881‑M2/75‑93 (MQ=255)
                 ggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCa                                                                            <  1:3672427‑M2/19‑1 (MQ=255)
                  gttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCAt                                                                           <  2:4515760‑M2/20‑1 (MQ=255)
                  gttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCAt                                                                           >  1:285743‑M2/74‑93 (MQ=255)
                  gttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCAt                                                                           <  1:3936351‑M2/20‑1 (MQ=255)
                   ttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCAtt                                                                          >  1:2000545‑M2/73‑93 (MQ=255)
                    taacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCattc                                                                         >  2:4577454‑M2/72‑92 (MQ=255)
                    taacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTa                                                                         >  1:4469141‑M2/72‑93 (MQ=255)
                     aacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTAc                                                                        >  2:4698011‑M2/71‑93 (MQ=255)
                      acccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACg                                                                       >  1:4532482‑M2/70‑93 (MQ=255)
                      acccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACg                                                                       >  2:3444695‑M2/70‑93 (MQ=255)
                      acccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACg                                                                       >  2:3988634‑M2/70‑93 (MQ=255)
                       cccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGa                                                                      >  2:2603791‑M2/69‑93 (MQ=255)
                        ccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGAc                                                                     >  1:1705445‑M2/68‑93 (MQ=255)
                        ccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGAc                                                                     >  1:2556457‑M2/68‑93 (MQ=255)
                        ccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGAc                                                                     >  1:2780419‑M2/68‑93 (MQ=255)
                         cccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGAcc                                                                    <  1:660149‑M2/27‑1 (MQ=255)
                         cccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGAcc                                                                    >  1:494203‑M2/67‑93 (MQ=255)
                         cccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGAcc                                                                    >  1:4271466‑M2/67‑93 (MQ=255)
                           caaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAg                                                                  >  2:4421470‑M2/65‑93 (MQ=255)
                           caaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAg                                                                  <  2:640392‑M2/29‑1 (MQ=255)
                           caaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAg                                                                  >  2:1871017‑M2/65‑93 (MQ=255)
                           caaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAg                                                                  <  1:2208194‑M2/29‑1 (MQ=255)
                            aaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGt                                                                 >  2:2740474‑M2/64‑93 (MQ=255)
                            aaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGt                                                                 >  2:841774‑M2/64‑93 (MQ=255)
                               gaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTa                                                              >  1:3153892‑M2/61‑93 (MQ=255)
                                aatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTaa                                                             <  2:1594651‑M2/34‑1 (MQ=255)
                                   ttggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTaaaaa                                                          <  1:3779975‑M2/37‑1 (MQ=255)
                                   gtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTaaaaa                                                          >  1:2480881‑M2/57‑93 (MQ=255)
                                    tggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAg                                                         >  2:1719706‑M2/56‑93 (MQ=255)
                                      gtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAgcg                                                       >  1:3216889‑M2/54‑93 (MQ=255)
                                        agcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGcc                                                     >  1:2941447‑M2/52‑93 (MQ=255)
                                        agcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGcc                                                     >  1:526980‑M2/52‑93 (MQ=255)
                                         gcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCt                                                    >  2:1149852‑M2/51‑93 (MQ=255)
                                           ccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTTATTACGACCAGTCTAAAAAGCGCCTGa                                                  >  1:1361916‑M2/49‑93 (MQ=255)
                                           ccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacataggGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGa                                                  >  1:1831973‑M2/49‑93 (MQ=255)
                                            cgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGaa                                                 <  1:925934‑M2/46‑1 (MQ=255)
                                            cgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGaa                                                 <  1:3295025‑M2/46‑1 (MQ=255)
                                            cgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGaa                                                 >  2:3852831‑M2/48‑93 (MQ=255)
                                            cgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGaa                                                 >  1:4502803‑M2/48‑93 (MQ=255)
                                             gcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAAt                                                >  2:2039963‑M2/47‑93 (MQ=255)
                                              cgccgtcgctgtcgacgcatgcggatcctctagagtcgagatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAACCGCCtgaaat                                               <  2:4403843‑M2/48‑3 (MQ=255)
                                              cgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAAtt                                               >  1:2106177‑M2/46‑93 (MQ=255)
                                               gccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTc                                              <  1:2351947‑M2/49‑1 (MQ=255)
                                                 cgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTcgc                                            >  1:1194001‑M2/43‑93 (MQ=255)
                                                 cgtcgctggcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTcgc                                            <  1:1292729‑M2/51‑1 (MQ=255)
                                                  gtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATTCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTcgcg                                           >  2:141451‑M2/42‑93 (MQ=255)
                                                   tcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGa                                          >  1:267264‑M2/41‑93 (MQ=255)
                                                       tgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCtt                                      >  2:754603‑M2/37‑93 (MQ=255)
                                                       tgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCtt                                      >  1:2923630‑M2/37‑93 (MQ=255)
                                                        gtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTtc                                     <  2:1074053‑M2/58‑1 (MQ=255)
                                                        gtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTtc                                     >  1:2559851‑M2/36‑93 (MQ=255)
                                                          cgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTtctc                                   <  1:1000848‑M2/60‑1 (MQ=255)
                                                          cgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTtctc                                   >  1:1624050‑M2/34‑93 (MQ=255)
                                                           gacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGTCTGAATTCGCGACCTTCTCg                                  >  1:4348447‑M2/33‑93 (MQ=255)
                                                           gacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCg                                  >  2:343267‑M2/33‑93 (MQ=255)
                                                               catgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTACAAAGCGCCTGAATTCGCGACCTTCTCGTTAc                              >  1:681852‑M2/29‑93 (MQ=255)
                                                               catgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTAc                              >  2:2206392‑M2/29‑93 (MQ=255)
                                                                 tgcggatcctctagagtcgacatatgGTGATCATGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTg                            >  1:1484033‑M2/27‑93 (MQ=255)
                                                                    ggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACa                         <  2:2255950‑M2/70‑1 (MQ=255)
                                                                    ggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACa                         >  1:3430332‑M2/24‑93 (MQ=255)
                                                                      atcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAgg                       >  2:2445613‑M2/22‑93 (MQ=255)
                                                                       tcctctatggtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGa                      <  2:905955‑M2/73‑1 (MQ=255)
                                                                       tcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGa                      >  2:40609‑M2/21‑93 (MQ=255)
                                                                        cctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGaa                     >  2:1112356‑M2/20‑93 (MQ=255)
                                                                          tctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGaaaa                   >  1:2308711‑M2/18‑93 (MQ=255)
                                                                             agagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAAtga                >  2:3063719‑M2/15‑92 (MQ=255)
                                                                             agagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAAtga                >  2:4030562‑M2/15‑92 (MQ=255)
                                                                                gtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGcc             >  1:1678140‑M2/12‑93 (MQ=255)
                                                                                gtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGcc             >  2:2102992‑M2/12‑93 (MQ=255)
                                                                                gtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGcc             >  2:1489695‑M2/12‑93 (MQ=255)
                                                                                gtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGcc             >  1:2426508‑M2/12‑93 (MQ=255)
                                                                                  cgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCAt           >  1:655650‑M2/10‑93 (MQ=255)
                                                                                   gacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCTCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCAtt          <  2:1699199‑M2/85‑1 (MQ=255)
                                                                                    acatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTg         >  1:3050933‑M2/8‑93 (MQ=255)
                                                                                      atatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGc       <  2:4138936‑M2/88‑1 (MQ=255)
                                                                                        atgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCaa     >  1:1101156‑M2/4‑93 (MQ=255)
                                                                                         tgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAAc    >  1:1949394‑M2/3‑93 (MQ=255)
                                                                                         tgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAAc    <  2:3144074‑M2/91‑1 (MQ=255)
                                                                                         tgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAAc    >  2:3317569‑M2/3‑93 (MQ=255)
                                                                                          gGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAAcc   >  2:2827598‑M2/2‑93 (MQ=255)
                                                                                           gtgATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCa  <  2:1242152/93‑1 (MQ=255)
                                                                                           |                                                                                            
GGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGCGACACAATTGCGAAAGGTCATAATGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGGGCCATTGGCAACCA  >  NC_000913/4094942‑4095125

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: