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 210 [0] [0] 242 [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)143 (0.860)
+GCATGCGGATCCTCTAGAGTCGACATATG
80/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 242 total aligned reads displayed.

AGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGCGACACAATTGCGAAAGGTCATAATGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGGGCCATTGGCAACCA  >  NC_000913/4094941‑4095125
                                                                                            |                                                                                            
agtggggacggcttcgttggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgg                                                                                              >  2:2896954‑M2/93‑93 (MQ=255)
    gggacggcttcgttggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGAt                                                                                          <  1:177991‑M2/5‑1 (MQ=255)
      gacggcttcgttggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATcc                                                                                        >  1:3059679‑M2/87‑93 (MQ=255)
       acggcttcgttggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCct                                                                                       >  1:1095956‑M2/86‑93 (MQ=255)
       acggcttcgttggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCct                                                                                       >  1:3692807‑M2/86‑93 (MQ=255)
         ggcttcgttggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCctgc                                                                                     >  1:6069370‑M2/84‑93 (MQ=255)
          gcttcgttggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCctgct                                                                                    <  1:4234350‑M2/11‑1 (MQ=255)
             tcgttggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGaa                                                                                 >  1:384774‑M2/80‑93 (MQ=255)
               gttggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAAtt                                                                               >  1:5795332‑M2/78‑93 (MQ=255)
                ttggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAAttt                                                                              <  1:2702509‑M2/17‑1 (MQ=255)
                 tggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacctatgGTGATCCTTCTGAattta                                                                             >  2:1721046‑M2/76‑92 (MQ=255)
                 tggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTc                                                                             <  2:1312079‑M2/18‑1 (MQ=255)
                  ggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCa                                                                            <  1:3457413‑M2/19‑1 (MQ=255)
                   gttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCAt                                                                           <  2:4310418‑M2/20‑1 (MQ=255)
                     taacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTa                                                                         <  2:878222‑M2/22‑1 (MQ=255)
                      aacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctatagtcgacatatgGTGATCCTGCTGAATTTcatttc                                                                        >  1:4449347‑M2/71‑91 (MQ=255)
                      aacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTAc                                                                        >  2:3168332‑M2/71‑93 (MQ=255)
                      aacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTAc                                                                        <  1:4114409‑M2/23‑1 (MQ=255)
                        cccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGa                                                                      <  2:4539391‑M2/25‑1 (MQ=255)
                         ccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGAc                                                                     >  1:4739330‑M2/68‑93 (MQ=255)
                         ccccaaagaatggggtagcccgcgccggcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGAc                                                                     <  2:571142‑M2/26‑1 (MQ=255)
                           ccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCa                                                                   <  1:319348‑M2/28‑1 (MQ=255)
                            caaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAg                                                                  >  2:4054237‑M2/65‑93 (MQ=255)
                              aagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTc                                                                <  2:4770140‑M2/31‑1 (MQ=255)
                                gaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTa                                                              >  1:5233625‑M2/61‑93 (MQ=255)
                                gaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTa                                                              <  1:1950746‑M2/33‑1 (MQ=255)
                                 aatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGGCCAGTCTaa                                                             <  2:4760530‑M2/34‑1 (MQ=255)
                                 aatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTaa                                                             >  2:4230401‑M2/60‑93 (MQ=255)
                                  atgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTaaa                                                            >  2:2122623‑M2/59‑93 (MQ=255)
                                   tgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTaaaa                                                           >  1:693172‑M2/58‑93 (MQ=255)
                                   tgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTaaaa                                                           <  2:799861‑M2/36‑1 (MQ=255)
                                     tggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAg                                                         <  1:5505077‑M2/38‑1 (MQ=255)
                                      ggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAgc                                                        <  2:4911120‑M2/39‑1 (MQ=255)
                                      ggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAgc                                                        <  2:5780639‑M2/39‑1 (MQ=255)
                                       gtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAgcg                                                       >  1:4822308‑M2/54‑93 (MQ=255)
                                         agcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGcc                                                     <  2:2350156‑M2/42‑1 (MQ=255)
                                         accccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGcc                                                     <  1:1607866‑M2/42‑1 (MQ=255)
                                          gcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCt                                                    <  2:3474714‑M2/43‑1 (MQ=255)
                                             cgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGaa                                                 >  1:6259221‑M2/48‑93 (MQ=255)
                                              gcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAAt                                                <  2:5123796‑M2/47‑1 (MQ=255)
                                               cgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAAtt                                               >  2:3431844‑M2/46‑93 (MQ=255)
                                                gccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGGTCCTGCTGAATTTCATTACGACCACTCTAAAAAGCGCCTGAATTc                                              >  1:2497294‑M2/45‑93 (MQ=255)
                                                gccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTc                                              <  1:494421‑M2/49‑1 (MQ=255)
                                                 ccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTcg                                             >  1:2820080‑M2/44‑93 (MQ=255)
                                                 ccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTcg                                             >  1:4405007‑M2/44‑93 (MQ=255)
                                                  cgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTcgc                                            <  2:2488098‑M2/51‑1 (MQ=255)
                                                  cgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTcgc                                            <  1:3488958‑M2/51‑1 (MQ=255)
                                                   gtcgctgtcgacgtatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTcgcg                                           >  2:1410443‑M2/42‑93 (MQ=255)
                                                    tcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATtcgcgc                                          >  2:2469224‑M2/41‑92 (MQ=255)
                                                    tcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGa                                          >  2:1594860‑M2/41‑93 (MQ=255)
                                                     cgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGAc                                         >  2:1273913‑M2/40‑93 (MQ=255)
                                                     cgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGAc                                         <  2:2043128‑M2/54‑1 (MQ=255)
                                                       ctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGCATTCGCGACCt                                       >  1:2758633‑M2/38‑93 (MQ=255)
                                                         gtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTtc                                     >  1:3004191‑M2/36‑93 (MQ=255)
                                                         gtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTtc                                     >  2:4424040‑M2/36‑93 (MQ=255)
                                                         gtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTtc                                     >  1:3774602‑M2/36‑93 (MQ=255)
                                                          tcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTtct                                    <  2:3448301‑M2/59‑1 (MQ=255)
                                                          tcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTtct                                    <  2:2812319‑M2/59‑1 (MQ=255)
                                                            gacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCg                                  <  2:3880052‑M2/61‑1 (MQ=255)
                                                              cgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGtt                                <  1:3649843‑M2/63‑1 (MQ=255)
                                                               gcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTGGCGACCTTCTCgttc                               >  1:1216513‑M2/30‑92 (MQ=255)
                                                                catgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTAc                              <  2:2983590‑M2/65‑1 (MQ=255)
                                                                catgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTAc                              <  1:93526‑M2/65‑1 (MQ=255)
                                                                catgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTAc                              <  1:3099464‑M2/65‑1 (MQ=255)
                                                                catgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTAc                              <  2:4023045‑M2/65‑1 (MQ=255)
                                                                  tgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTg                            <  2:5906851‑M2/67‑1 (MQ=255)
                                                                   gcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGa                           <  2:177962‑M2/68‑1 (MQ=255)
                                                                   gcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGa                           <  2:1117607‑M2/68‑1 (MQ=255)
                                                                     ggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTTGCGACCTTCTCGTTACTGACa                         >  2:5395560‑M2/24‑93 (MQ=255)
                                                                         cctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGaa                     >  2:2209181‑M2/20‑93 (MQ=255)
                                                                         cctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGaa                     >  1:4873360‑M2/20‑93 (MQ=255)
                                                                          ctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGaaa                    >  2:4210895‑M2/19‑93 (MQ=255)
                                                                          ctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGaaa                    >  2:6198949‑M2/19‑93 (MQ=255)
                                                                          ctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGACTTAAAAGAGCCTGAATTCGCGAACTTCTCGTTACTGACAGGaaa                    >  2:3524707‑M2/19‑93 (MQ=255)
                                                                           tctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGaaaa                   <  1:1519033‑M2/76‑1 (MQ=255)
                                                                            ctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAAt                  >  1:1046474‑M2/17‑93 (MQ=255)
                                                                            ctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAAt                  >  1:5068768‑M2/17‑93 (MQ=255)
                                                                             tagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATg                 >  2:2595286‑M2/16‑93 (MQ=255)
                                                                             tagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATg                 >  2:6230081‑M2/16‑93 (MQ=255)
                                                                              agagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAAtga                <  2:5282546‑M2/79‑2 (MQ=255)
                                                                              agagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAAtga                >  2:3043492‑M2/15‑92 (MQ=255)
                                                                              agagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAAtga                <  1:1148377‑M2/79‑2 (MQ=255)
                                                                               gagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAatgag               >  1:1296218‑M2/14‑91 (MQ=255)
                                                                               gagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAatgag               >  1:5638235‑M2/14‑91 (MQ=255)
                                                                                  tcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCa            >  1:3560674‑M2/11‑93 (MQ=255)
                                                                                  tcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCa            <  1:2124929‑M2/83‑1 (MQ=255)
                                                                                   cgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCAt           <  2:1830819‑M2/84‑1 (MQ=255)
                                                                                    gacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCAtt          <  1:5831635‑M2/85‑1 (MQ=255)
                                                                                    gacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCAtt          >  1:3989625‑M2/9‑93 (MQ=255)
                                                                                    gacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCAtt          <  1:887462‑M2/85‑1 (MQ=255)
                                                                                     acatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGTCCTTCTCGTTACTGACAGGAAAATGAGCCATTg         <  2:2541741‑M2/86‑1 (MQ=255)
                                                                                      catatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTgg        <  2:2058859‑M2/87‑1 (MQ=255)
                                                                                      catatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTgg        <  1:1421036‑M2/87‑1 (MQ=255)
                                                                                        tatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCCTTTGCa      >  2:6071338‑M2/5‑93 (MQ=255)
                                                                                        tatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCa      >  2:6001189‑M2/5‑93 (MQ=255)
                                                                                        tatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCa      >  1:3335549‑M2/5‑93 (MQ=255)
                                                                                          tgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAAc    >  2:5337698‑M2/3‑93 (MQ=255)
                                                                                           gGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAAcc   <  2:1963178‑M2/92‑1 (MQ=255)
                                                                                           gGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAAcc   <  1:5893878‑M2/92‑1 (MQ=255)
                                                                                            gtgATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCa  >  1:4555178/1‑93 (MQ=255)
                                                                                            |                                                                                            
AGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGCGACACAATTGCGAAAGGTCATAATGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGGGCCATTGGCAACCA  >  NC_000913/4094941‑4095125

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