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

CTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGCGACACAATTGCGAAAGGTCATAATGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGGGCCATTGGCAACCA  >  NC_000913/4094944‑4095125
                                                                                         |                                                                                            
ggggacggcttcgttggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGa                                                                                           >  2:2382759‑M2/90‑93 (MQ=255)
 gggacggcttcgttggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGAt                                                                                          >  1:5712665‑M2/89‑93 (MQ=255)
 gggacggcttcgttggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagacgacatatgGTGAt                                                                                          >  2:365551‑M2/89‑93 (MQ=255)
   gacggcttcgttggttaccccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATcc                                                                                        <  1:6399819‑M2/7‑1 (MQ=255)
   gacggcttcgttggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATcc                                                                                        >  1:198085‑M2/87‑93 (MQ=255)
   gacggcttcgttggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATcc                                                                                        <  1:4674239‑M2/7‑1 (MQ=255)
   gacggcttcgttggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATcc                                                                                        >  1:6556237‑M2/87‑93 (MQ=255)
   gacggcttcgttggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATcc                                                                                        <  2:6139194‑M2/7‑1 (MQ=255)
      ggcttcgttggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCctgc                                                                                     >  1:2627900‑M2/84‑93 (MQ=255)
      ggcttcgttggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCctgc                                                                                     <  2:2941550‑M2/10‑1 (MQ=255)
      cgcttcgttggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCctgc                                                                                     <  1:280624‑M2/10‑1 (MQ=255)
       gcttcgttggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCctgct                                                                                    >  2:3950731‑M2/83‑93 (MQ=255)
       gcttcgttggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCctgct                                                                                    <  1:3557919‑M2/11‑1 (MQ=255)
       gcttcgttggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCctgct                                                                                    >  1:4397612‑M2/83‑93 (MQ=255)
         ttcgttggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGa                                                                                  >  1:6384343‑M2/81‑93 (MQ=255)
          tcgttggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGaa                                                                                 >  1:3326569‑M2/80‑93 (MQ=255)
           cgttggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAAt                                                                                <  2:5432911‑M2/15‑1 (MQ=255)
            gttggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAAtt                                                                               >  1:293903‑M2/78‑93 (MQ=255)
            gttggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAAtt                                                                               <  2:5025035‑M2/16‑1 (MQ=255)
            gttggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAAtt                                                                               <  1:560562‑M2/16‑1 (MQ=255)
               ggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCa                                                                            <  1:5554111‑M2/19‑1 (MQ=255)
               ggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCa                                                                            >  1:4658988‑M2/75‑93 (MQ=255)
               ggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCa                                                                            <  2:1559581‑M2/19‑1 (MQ=255)
                gttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCAt                                                                           >  2:4969961‑M2/74‑93 (MQ=255)
                gttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCAt                                                                           >  2:1204717‑M2/74‑93 (MQ=255)
                 ttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCAtt                                                                          <  1:5565283‑M2/21‑1 (MQ=255)
                   aacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTAc                                                                        <  2:4064005‑M2/23‑1 (MQ=255)
                   aacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTAc                                                                        >  2:6403203‑M2/71‑93 (MQ=255)
                     cccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGa                                                                      <  2:4759064‑M2/25‑1 (MQ=255)
                     cccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGa                                                                      >  1:4069704‑M2/69‑93 (MQ=255)
                     cccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGa                                                                      >  2:1065066‑M2/69‑93 (MQ=255)
                        ccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCa                                                                   >  2:1110477‑M2/66‑93 (MQ=255)
                           aagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTc                                                                <  2:4491289‑M2/31‑1 (MQ=255)
                           aagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTc                                                                >  1:3749461‑M2/63‑93 (MQ=255)
                           aagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTc                                                                <  2:3615868‑M2/31‑1 (MQ=255)
                           aagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTc                                                                <  2:320414‑M2/31‑1 (MQ=255)
                            agaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCt                                                               <  1:976706‑M2/32‑1 (MQ=255)
                            agaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCt                                                               >  2:847037‑M2/62‑93 (MQ=255)
                                 gtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTaaaaa                                                          >  1:5082288‑M2/57‑93 (MQ=255)
                                 gtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTaaaaa                                                          <  1:5930133‑M2/37‑1 (MQ=255)
                                 gtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTaaaaa                                                          <  2:5917188‑M2/37‑1 (MQ=255)
                                    gtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAgcg                                                       <  1:2521557‑M2/40‑1 (MQ=255)
                                           gcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAAt                                                <  2:3061228‑M2/47‑1 (MQ=255)
                                           gcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAAt                                                >  1:2071776‑M2/47‑93 (MQ=255)
                                           gcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAAt                                                <  1:2184522‑M2/47‑1 (MQ=255)
                                            cgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAAtt                                               <  1:3190666‑M2/48‑1 (MQ=255)
                                             gccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTc                                              >  2:1388088‑M2/45‑93 (MQ=255)
                                             gccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTc                                              >  1:1411196‑M2/45‑93 (MQ=255)
                                             gccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTc                                              <  2:750127‑M2/49‑1 (MQ=255)
                                             gccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTc                                              <  2:6715064‑M2/49‑1 (MQ=255)
                                             cccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTc                                              <  1:4245019‑M2/49‑1 (MQ=255)
                                              ccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTcg                                             >  1:6124443‑M2/44‑93 (MQ=255)
                                              ccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTcg                                             >  2:3538930‑M2/44‑93 (MQ=255)
                                                gtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTcgcg                                           >  1:1229637‑M2/42‑93 (MQ=255)
                                                gtcgctgtcgacgcatgcggatcctctagaggcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTcgcg                                           >  2:3334017‑M2/42‑93 (MQ=255)
                                                  cgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGAc                                         >  1:860266‑M2/40‑93 (MQ=255)
                                                  cgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGAc                                         >  1:4829656‑M2/40‑93 (MQ=255)
                                                   gctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGAcc                                        >  1:1698520‑M2/39‑93 (MQ=255)
                                                     tgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCtt                                      >  1:2403971‑M2/37‑93 (MQ=255)
                                                     tgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCtt                                      <  1:6083139‑M2/57‑1 (MQ=255)
                                                     tgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCtt                                      >  1:690284‑M2/37‑93 (MQ=255)
                                                      gtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTtc                                     >  2:5269206‑M2/36‑93 (MQ=255)
                                                       tcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTtct                                    >  2:1842225‑M2/35‑93 (MQ=255)
                                                       tcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTtct                                    <  1:3957737‑M2/59‑1 (MQ=255)
                                                       tcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTtct                                    <  2:4538045‑M2/59‑1 (MQ=255)
                                                        cgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTtctc                                   >  1:2746091‑M2/34‑93 (MQ=255)
                                                         gacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCg                                  <  2:1741141‑M2/61‑1 (MQ=255)
                                                          acgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGt                                 <  1:6033012‑M2/62‑1 (MQ=255)
                                                           cgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGtt                                >  2:5683421‑M2/31‑93 (MQ=255)
                                                             catgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTAc                              <  2:4170818‑M2/65‑1 (MQ=255)
                                                               tgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTg                            <  2:1974525‑M2/67‑1 (MQ=255)
                                                                  ggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACa                         <  2:5073598‑M2/70‑1 (MQ=255)
                                                                  ggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACa                         <  2:2975734‑M2/70‑1 (MQ=255)
                                                                  ggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACa                         >  2:4120309‑M2/24‑93 (MQ=255)
                                                                      cctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGaa                     >  1:4431773‑M2/20‑93 (MQ=255)
                                                                       ctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGaaa                    >  1:5441623‑M2/19‑93 (MQ=255)
                                                                       ctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGaaa                    >  1:5396991‑M2/19‑93 (MQ=255)
                                                                        tctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGaaaa                   <  2:2489326‑M2/76‑1 (MQ=255)
                                                                        tctagagtcgacatatgGTGATCCTGCAGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGaaaa                   <  2:3422589‑M2/76‑1 (MQ=255)
                                                                           agagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAAtga                <  1:4172830‑M2/79‑2 (MQ=255)
                                                                           agagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAAtga                >  1:4058334‑M2/15‑92 (MQ=255)
                                                                             agtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCATGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGc              >  2:782106‑M2/13‑93 (MQ=255)
                                                                              gtcgacatatgGTGATCCTGTTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGcc             <  2:3934818‑M2/82‑1 (MQ=255)
                                                                              gtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGcc             <  1:4614080‑M2/82‑1 (MQ=255)
                                                                              gtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGcc             >  2:6481628‑M2/12‑93 (MQ=255)
                                                                                cgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCAt           <  2:4059661‑M2/84‑1 (MQ=255)
                                                                                cgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCAt           >  1:2215216‑M2/10‑93 (MQ=255)
                                                                                cgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCAt           <  2:2306207‑M2/84‑1 (MQ=255)
                                                                                cgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCAt           <  2:2274232‑M2/84‑1 (MQ=255)
                                                                                cgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCAt           <  2:6210603‑M2/84‑1 (MQ=255)
                                                                                 gacatatgGTGATCGTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCAtt          <  1:4805965‑M2/85‑1 (MQ=255)
                                                                                 gacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCAtt          <  2:4843237‑M2/85‑1 (MQ=255)
                                                                                     tatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCa      <  1:1236144‑M2/89‑1 (MQ=255)
                                                                                       tgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAAc    >  1:5157600‑M2/3‑93 (MQ=255)
                                                                                         gtgATCGTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCa  <  1:241663/93‑1 (MQ=255)
                                                                                         gtgATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCa  <  1:6212003/93‑1 (MQ=255)
                                                                                         gtgATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCa  <  2:6606537/93‑1 (MQ=255)
                                                                                         gtgATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCa  >  1:2869787/1‑93 (MQ=255)
                                                                                         gtgATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCa  <  1:3063193/93‑1 (MQ=255)
                                                                                         gtgATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCa  >  2:4300672/1‑93 (MQ=255)
                                                                                         |                                                                                            
CTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGCGACACAATTGCGAAAGGTCATAATGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGGGCCATTGGCAACCA  >  NC_000913/4094944‑4095125

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