New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 4583412 =114 (0.610)5 (0.030) 5/226 NT 4.6% coding (3350/3513 nt) hsdR endonuclease R Type I restriction enzyme
?NC_000913 4583454 = 98 (0.550)coding (3308/3513 nt) hsdR endonuclease R Type I restriction enzyme
Rejected: Frequency below/above cutoff threshold.

Only 100 of 144 total aligned reads displayed.

CCAGATATAATCGCTGAATTTGCCCAGCAGGGTATCGAGATTATCGTCAAAGGTTCTTTGCAGCATCGCCTTCCCGCCGCGACGGTGGAAGTTGCCGGTTTTGAAGACATCGTCGTCGAGCACCACTTTCTCTTTCAGCGCCTGCGCTAAACGATCGAGCCAGCTTAATTGCTCGCTGCTCCAGTCGTTTTCGCCCTTAATGCGCGTCAGCGCGTGATCGACACGTTCCTCAAACGGTTTCAGCGCATCGCCCACCGCAGCGCGGCGAATATGACCAATCAGCCGGGCGGCGATATC  >  NC_000913/4583264‑4583560
                                                                                                                                                    |                                                                                                                                                    
ccAGATATAATCGCTGAATTTGCCCAGCAGGGTATCGAGATTATCGTCAAAGGTTCTTTGCAGCATCGCCTTCCCGCCGCGACGGTGGAAGTTGCCGGTTTTGAAGACATCGTCGTCGAGCACCACTTTCTCTTTCAGCGCCTGCGCTa                                                                                                                                                      >  4:1352806/1‑149 (MQ=255)
  aGATATAATCGCTGAATTTGCCCAGCAGGGTATCGAGATTATCGTCAAAGGTTCTTTGCAGCATCGCCTTCCCGCCGCGACGGTGGAAGTTGCCGGTTTTGAAGACATCGTCGTCGAGCACCACTTTCTCTTTCAGCGCCTGCGCTaaa                                                                                                                                                    <  3:1352806/149‑1 (MQ=255)
   gATATAATCGCTGAATTTGCCCAGCAGGGTATCGAGATTATCGTCAAAGGTTCTTTGCAGCATCGCCTTCCCGCCGCGACGGTGGAAGTTGCCGGTTTTGAAGACATCGTCGTCGAGCACCACTTTCTCTTTCAGCGCCTGCGCTAAAc                                                                                                                                                   >  2:417176/1‑149 (MQ=255)
    ttccctgcgcaaagcatggaaagagacgcgcaatgaagatatcgccgcccggctgattggtcatattcgccgcgctgcggtgggcgatgcgctgaaaccgtttgaggaacgtgtcgatcacgcgctgacgcgcattaagggcgaAAACg                                                                                                                                                  <  4:1166328‑M1/5‑1 (MQ=255)
    ttccctgcgcaaagcatggaaagagacgcgcaatgaagatatcgccgcccggctgattggtcatattcgccgcgctgcggtgggcgatgcgctgaaaccgtttgaggaacgtgtcgatcacgcgctgacgcgcattaagggcgaAAACg                                                                                                                                                  <  1:879150‑M1/5‑1 (MQ=255)
     tccctgcgcaaagcatggaaagagacgcgcaatgaagatatcgccgcccggctgattggtcatattcgccgcgctgcggtgggcgatgcgctgaaaccgtttgaggaacgtgtcgatcacgcgctgacgcgcattaagggcgaAAACGa                                                                                                                                                 >  2:879150‑M1/144‑149 (MQ=255)
     tataATCGCTGAATTTGCCCAGCAGGGTATCGAGATTATCGTCAAAGGTTCTTTGCAGCATCGCCTTCCCGCCGCGACGGTGGAAGTTGCCGGTTTTGAAGACATCGTCGTCGAGCACCACTTTCTCTTTCAGCGCCTGCGCTAAACGa                                                                                                                                                 <  2:457320/149‑1 (MQ=255)
     tataATCGCTGAATTTGCCCAGCAGGGTATCGAGATTATCGTCAAAGGTTCTTTGCAGCATCGCCTTCCCGCCGCGACGGTGGAAGTTGCCGGTTTTGAAGACATCGTCGTCGAGCACCACTTTCTCTTTCAGCGCCTGCGCTAAACGa                                                                                                                                                 >  4:361527/1‑149 (MQ=255)
      ataATCGCTGAATTTGCCCAGCAGGGTATCGAGATTATCGTCAAAGGTTCTTTGCAGCATCGCCTTCCCACCGCGACGGTGGAAGTTGCCGGTTTTGAAGACATCGTCGTCGAGCACCACTTTCTCTTTCCGCGCCTGCGCTAAACCAt                                                                                                                                                <  4:1266501/149‑1 (MQ=255)
         aTCGCTGAATTTGCCCAGCAGGGTATCGAGATTATCGTCAAAGGTTCTTTGCAGCATCGCCTTCCCGCCGCGACGGTGGAAGTTGCCGGTTTTGAAGACATCGTCGTCGAGCACCACTTTCTCTTTCAGCGCCTGCGCTAAACGATCGa                                                                                                                                             <  3:361527/149‑1 (MQ=255)
              tGAATTTGCCCAGCAGGGTATCGAGATTATCGTCAAAGGTTCTTTGCAGCATCGCCTTCCCGCCGCGACGGTGGAAGTTGCCGGTTTTGAAGACATCGTCGTCGAGCACCACTTTCTCTTTCAGCGCCTGCGCTAAACGATCGAGCCAg                                                                                                                                        <  2:1102058/149‑1 (MQ=255)
              tGAATTTGCCCAGCAGGGTATCGAGATTATCGTCAAAGGTTCTTTGCAGCATCGCCTTCCCGCCGCGACGGTGGAAGTTGCCGGTTTTGAAGACATCGTCGTCGAGCACCACTTTCTCTTTCAGCGCCTGCGCTAAACGATCGAGCCAg                                                                                                                                        >  1:1102058/1‑149 (MQ=255)
                  tttGCCCAGCAGGGTATCGAGATTATCGTCAAAGGTTCTTTGCAGCATCGCCTTCCCGCCGCGACGGTGGAAGTTGCCGGTTTTGAAGACATCGTCGTCGAGCACCACTTTCTCTTTCAGCGCCTGCGCTAAAc                                                                                                                                                   <  2:1333628/134‑1 (MQ=255)
                  tttGCCCAGCAGGGTATCGAGATTATCGTCAAAGGTTCTTTGCAGCATCGCCTTCCCGCCGCGACGGTGGAAGTTGCCGGTTTTGAAGACATCGTCGTCGAGCACCACTTTCTCTTTCAGCGCCTGCGCTAAAc                                                                                                                                                   >  1:1333628/1‑134 (MQ=255)
                  tttGCCCAGCAGGGTATCGAGATTATCGTCAAAGGTTCTTTGCAGCATCGCCTTCCCGCCGCGACGGTGGAAGTTGCCGGTTTTGAAGACATCGTCGTCGAGCACCACTTTCTCTTTCAGCGCCTGCGCTAAACGAt                                                                                                                                                <  2:1134015/137‑1 (MQ=255)
                  tttGCCCAGCAGGGTATCGAGATTATCGTCAAAGGTTCTTTGCAGCATCGCCTTCCCGCCGCGACGGTGGAAGTTGCCGGTTTTGAAGACATCGTCGTCGAGCACCACTTTCTCTTTCAGCGCCTGCGCTAAACGAt                                                                                                                                                >  1:1134015/1‑137 (MQ=255)
                  catggaaagagacgcgcaatgaagatatcgccgcccggctgattggtcatattcgccgcgctgcggtgggcgatgcgctgaaaccgtttgaggaacgtgtcgatcacgcgctgacgcgcattaagggcgaAAAc                                                                                                                                                   <  2:687419‑M1/4‑1 (MQ=255)
                    tGCCCAGCAGGGTATCGAGATTATCGTCAAAGGTTCTTTGCAGCATCGCCTTCCCGCCGCGACGGTGGAAGTTGCCGGTTTTGAAGACATCGTCGTCGAGCACCACTTTCTCTTTCAGCGCCTGCGCTAAACGATCGAGCCAGCTTAAt                                                                                                                                  <  2:145025/149‑1 (MQ=255)
                     gCCCAGCAGGGTATCGAGATTATCGTCAAAGGTTCTTTGCAGCATCGCCTTCCCGCCGCGACGGTGGAAGTTGCCGGTTTTGAAGACATCGTCGTCGAGCACCACTTTCTCTTTCAGCGCCTGCGCTAAACGATCGAGCCAGCTTAAtt                                                                                                                                 <  1:476972/149‑1 (MQ=255)
                     gCCCAGCAGGGTATCGAGATTATCGTCAAAGGTTCTTTGCAGCATCGCCTTCCCGCCGCGACGGTGGAAGTTGCCGGTTTTGAAGACATCGTCGTCGAGCACCACTTTCTCTTTCAGCGCCTGCGCTAAACGATCGAGCCAGCTTAAtt                                                                                                                                 >  3:1697814/1‑149 (MQ=255)
                     gCCCAGCAGGGTATCGAGATTATCGTCAAAGGTTCTTTGCAGCATCGCCTTCCCGCCGCGACGGTGGAAGTTGCCGGTTTTGAAGACATCGTCGTCGAGCACCACTTTCTCTTTCAGCGCCTGCGCTAAACGATCGAGCCAGCTTAAtt                                                                                                                                 <  1:417176/149‑1 (MQ=255)
                       ccAGCAGGGTATCGAGATTATCGTCAAAGGTTCTTTGCAGCATCGCCTTCCCGCCGCGACGGTGGAAGTTGCCGGTTTTGAAGACATCGTCGTCGAGCACCACTTTCTCTTTCAGCGCCTGCGCTAAACGATCGAGCCAGCTTAATTGc                                                                                                                               >  3:649703/1‑149 (MQ=255)
                       ccAGCAGGGTATCGAGATTATCGTCAAAGGTTCTTTGCAGCATCGCCTTCCCGCCGCGACGGTGGAAGTTGCCGGTTTTGAAGACATCGTCGTCGAGCACCACTTTCTCTTTCAGCGCCTGCGCTAAACGATCGAGCCAGCTTAATTGc                                                                                                                               >  3:1219175/1‑149 (MQ=255)
                       ccAGCAGGGTATCGAGATTATCGTCAAAGGTTCTTTGCAGCATCGCCTTCCCGCCGCGACGGTGGAAGTTGCCGGTTTTGAAGACATCGTCGTCGAGCACCACTTTCTCTTTCAGCGCCTGCGCTAAACGATCGAGCCAGCTTAATTGc                                                                                                                               >  1:1070860/1‑149 (MQ=255)
                       ccAGCAGGGTATCGAGATTATCGTCAAAGGTTCTTTGCAGCATCGCCTTCCCGCCGCGACGGTGGAAGTTGCCGGTTTTGAAGACATCGTCGTCGAGCACCACTTTCTCTTTCAGCGCCTGCGCTAAACGATCGAGCCAGCTTAATTGc                                                                                                                               >  3:996949/1‑149 (MQ=255)
                         agcagGGTATCGAGATTATCGTCAAAGGTTCTTTGCAGCATCGCCTTCCCGCCGCGACGGTGGAAGTTGCCGGTTTTGAAGACATCGTCGTCGAGCACCACTTTCTCTTTCAGCGCCTGCGCTAAACGATCGAGCCAGCTTAATTGCTc                                                                                                                             <  3:576743/149‑1 (MQ=255)
                         agcagGGTATCGAGATTATCGTCAAAGGTTCTTTGCAGCATCGCCTTCCCGCCGCGACGGTGGAAGTTGCCGGTTTTGAAGACATCGTCGTCGAGCACCACTTTCTCTTTCAGCGCCTGCGCTAAACGATCGAGCCAGCTTAATTGCTc                                                                                                                             <  3:68644/149‑1 (MQ=255)
                         agcagGGTATCGAGATTATCGTCAAAGGTTCTTTGCAGCATCGCCTTCCCGCCGCGACGGTGGAAGTTGCCGGTTTTGAAGACATCGTCGTCGAGCACCACTTTCTCTTTCAGCGCCTGCGCTAAACGATCGAGCCAGCTTAATTGCTc                                                                                                                             <  1:1164886/149‑1 (MQ=255)
                            agGGTATCGAGATTATCGTCAAAGGTTCTTTGCAGCATCGCCTTCCCGCCGCGACGGTGGAAGTTGCCGGTTTTGAAGACATCGTCGTCGAGCACCACTTTCTCTTTCAGCGCCTGCGCTaaa                                                                                                                                                    <  4:586096/123‑1 (MQ=255)
                            agGGTATCGAGATTATCGTCAAAGGTTCTTTGCAGCATCGCCTTCCCGCCGCGACGGTGGAAGTTGCCGGTTTTGAAGACATCGTCGTCGAGCACCACTTTCTCTTTCAGCGCCTGCGCTAAACGATCGAGCCAGCTTAATTGCTCgct                                                                                                                          <  1:417120/149‑1 (MQ=255)
                                tATCGAGATTATCGTCAAAGGTTCTTTGCAGCATCGCCTTCCCGCCGCGACGGTGGAAGTTGCCGGTTTTGAAGACATCGTCGTCGAGCACCACTTTCTCTTTCAGCGCCTGCGCTAAACGATCGAGCCAGCTTAATTGCTCGCTGCTc                                                                                                                      <  1:1754433/149‑1 (MQ=255)
                                     agaTTATCGTCAAAGGTTCTTTGCAGCATCGCCTTCCCGCCGCGACGGTGGAAGTTGCCGGTTTTGAAGACATCGTCGTCGAGCACCACTTTCTCTTTCAGCGCCTGCGCTAAACGATCGAGCCAGCTTAATTGCTCGCTGCTCCAGTc                                                                                                                 <  4:492393/149‑1 (MQ=255)
                                        ttATCGTCAAAGGTTCTTTGCAGCATCGCCTTCCCGCCGCGACGGTGGAAGTTGCCGGTTTTGAAGACATCGTCGTCGAGCACCACTTTCTCTTTCAGCGCCTGCGCTAAACGATCGAGCCAGCTTAATTGCTCGCTGCTCCAGTCGtt                                                                                                              <  4:1697814/149‑1 (MQ=255)
                                           tatcgccgcccggctgattggtcatattcgccgcgctgcggtgggcgatgcgctgaaaccgtttgaggaacgtgtcgatcacgcgctgacgcgcattaagggcgaAAACGATCGAGCCAGCTTAATTGCTCGCTGCTCCAGTCGTTTTc                                                                                                           >  4:1826651‑M1/106‑149 (MQ=255)
                                           tCGTCAAAGGTTCTTTGCAGCATCGCCTTCCCGCCGCGACGGTGGAAGTTGCCGGTTTTGAAGACATCGTCGTCGAGCACCACTTTCTCTTTCAGCGCCTGCGCTAAACGATCGAGCCAGCTTAATTGCTCGCTGCTCCAGTCGTTTTc                                                                                                           >  1:9640/1‑149 (MQ=255)
                                             gTCAAAGGTTCTTTGCAGCATCGCCTTCCCGCCGCGACGGTGGAAGTTGCCGGTTTTGAAGACATCGTCGTCGAGCACCACTTTCTCTTTCAGCGCCTGCGCTAAACGATc                                                                                                                                               <  2:1099269/111‑1 (MQ=255)
                                              tCAAAGGTTCTTTGCAGCATCGCCTTCCCGCCGCGACGGTGGAAGTTGCCGGTTTTGAAGACATCGTCGTCGAGCACCACTTTCTCTTTCAGCGCCTGCGCTAAACGATCGAGCCAGCTTAATTGCt                                                                                                                              <  3:1453230/127‑1 (MQ=255)
                                              tCAAAGGTTCTTTGCAGCATCGCCTTCCCGCCGCGACGGTGGAAGTTGCCGGTTTTGAAGACATCGTCGTCGAGCACCACTTTCTCTTTCAGCGCCTGCGCTAAACGATCGAGCCAGCTTAATTGCt                                                                                                                              >  4:1453230/1‑127 (MQ=255)
                                                ccgcccggctgattggtcatattcgccgcgctgcggtgggcgatgcgctgaaaccgtttgaggaacgtgtcgatcacgcgctgacgcgcattaagggcgaAAACGa                                                                                                                                                 <  3:1871538‑M1/6‑1 (MQ=255)
                                                ccgcccggctgattggtcatattcgccgcgctgcggtgggcgatgcgctgaaaccgtttgaggaacgtgtcgatcacgcgctgacgcgcattaagggcgaAAACGa                                                                                                                                                 >  4:1871538‑M1/101‑106 (MQ=255)
                                                      tCTTTGCAGCATCGCCTTCCCGCCGCGACGGTGGAAGTTGCCGGTTTTGAAGACATCGTCGTCGAGCACCACTTTCTCTTTCAGCGCCTGCGCTAAACGATCGAGCCAGCTTAATTGCTCGCTGCTCCAGTCGTTTTCGCCCTTAATgc                                                                                                <  4:649703/149‑1 (MQ=255)
                                                         ttGCAGCATCGCCTTCCCGCCGCGACGGTGGAAGTTGCCGGTTTTGAAGACATCGTCGTCGAGCACCACTTTCTCTTTCAGCGCCTGCGCTAAACGATCGAGCCAGCTTAATTGCTCGCTGCTCCAGTCGTTTTCGCCCTTAATGcgcg                                                                                             >  1:1047409/1‑149 (MQ=255)
                                                            cagcaTCGCCTTCCCGCCGCGACGGTGGAAGTTGCCGGTTTTGAAGACATCGTCGTCGAGCACCACTTTCTCTTTCAGCGCCTGCGCTAAACGATCGAGCCAGCTTAAt                                                                                                                                  >  3:434322/1‑109 (MQ=255)
                                                            cagcaTCGCCTTCCCGCCGCGACGGTGGAAGTTGCCGGTTTTGAAGACATCGTCGTCGAGCACCACTTTCTCTTTCAGCGCCTGCGCTAAACGATCGAGCCAGCTTAAt                                                                                                                                  <  4:434322/109‑1 (MQ=255)
                                                               gtcatattcgccgcgctgcggtgggcgatgcgctgaaaccgtttgaggaacgtgtcgatcacgcgctgacgcgcattaagggcgaAAACGATCGAGCCAGCTTAATTGCTCGCTGCTCCAGTCGt                                                                                                               <  3:1738648‑M1/40‑1 (MQ=255)
                                                                aTCGCCTTCCCGCCGCGACGGTGGAAGTTGCCGGTTTTGAAGACATCGTCGTCGAGCACCACTTTCTCTTTCAGCGCCTGCGCTAAACGATCGAGCCAGCTTAATTGCTCGCTGCTCCAGTCGTTTTCGCCCTTAATgc                                                                                                <  1:1580595/139‑1 (MQ=255)
                                                                aTCGCCTTCCCGCCGCGACGGTGGAAGTTGCCGGTTTTGAAGACATCGTCGTCGAGCACCACTTTCTCTTTCAGCGCCTGCGCTAAACGATCGAGCCAGCTTAATTGCTCGCTGCTCCAGTCGTTTTCGCCCTTAATGCGCGTCAgcgc                                                                                      >  4:85196/1‑149 (MQ=255)
                                                                  cGCCTTCCCGCCGCGACGGTGGAAGTTGCCGGTTTTGAAGACATCGTCGTCGAGCACCACTTTCTCTTTCAGCGCCTGCGCTAAACGATCGAGCCAGCTTAATTGCTCGCTGCTCCAGTCGTTTTCGCCCTTAATGCGCGTCAGCGCGt                                                                                    <  1:103362/149‑1 (MQ=255)
                                                                   gCCTTCCCGCCGCGACGGTGGAAGTTGCCGGTTTTGAAGACATCGTCGTCGAGCACCACTTTCTCTTTCAGCGCCTGCGCTAAACGATCGAGCCAGCTTAATTGCTCGCTGCTCCAGTCGttt                                                                                                             <  3:422245/123‑1 (MQ=255)
                                                                   gCCTTCCCGCCGCGACGGTGGAAGTTGCCGGTTTTGAAGACATCGTCGTCGAGCACCACTTTCTCTTTCAGCGCCTGCGCTAAACGATCGAGCCAGCTTAATTGCTCGCTGCTCCAGTCGttt                                                                                                             <  1:546019/123‑1 (MQ=255)
                                                                    ccTTCCCGCCGCGACGGTGGAAGTTGCCGGTTTTGAAGACATCGTCGTCGAGCACCACTTTCTCTTTCAGCGCCTGCGCTAAACGAt                                                                                                                                                <  3:1085994/87‑1 (MQ=255)
                                                                      ttCCCGCCGCGACGGTGGAAGTTGCCGGTTTTGAAGACATCGTCGTCGAGCACCAGTTTCTCTTTCAGCGCCTGCGCTAAACGATCGAGCCAGCTTAATTGCTCGCTGCTCCAGTCGTTTTCGCCCTTAATGCGCGTCAGCGCGTGATc                                                                                >  2:1613538/1‑149 (MQ=255)
                                                                           gccgcGACGGTGGAAGTTGCCGGTTTTGAAGACATCGTCGTCGAGCACCACTTTCTCTTTCAGCGCCTGCGCTAAACGATCGAGCCAGCTTAATTGCTCGCTGCTCCAGTCGTTTTCGccc                                                                                                       <  3:862339/121‑1 (MQ=255)
                                                                           gccgcGACGGTGGAAGTTGCCGGTTTTGAAGACATCGTCGTCGAGCACCACTTTCTCTTTCAGCGCCTGCGCTAAACGATCGAGCCAGCTTAATTGCTCGCTGCTCCAGTCGTTTTCGccc                                                                                                       >  4:862339/1‑121 (MQ=255)
                                                                                gACGGTGGAAGTTGCCGGTTTTGAAGACATCGTCGTCGAGCACCACTTTCTCTTTCAGCGCCTGCGCTAAACGATCGAGCCAGCTTAATTGCTCGCTGCTCCAGTCGTTTTCGCCCTTAATGCGCGTCAGCGCGTGATCGACACGTTcc                                                                      <  4:957391/149‑1 (MQ=255)
                                                                                  ggtgggcgatgcgctgaaaccgtttgaggaacgtgtcgatcacgcgctgacgcgcattaagggcgaAAACGATCGAGCCAGCTTAATTGCTCGCTGCTCCAGTCGtt                                                                                                              >  1:517462‑M1/67‑107 (MQ=255)
                                                                                  ggtgggcgatgcgctgaaaccgtttgaggaacgtgtcgatcacgcgctgacgcgcattaagggcgaAAACGATCGAGCCAGCTTAATTGCTCGCTGCTCCAGTCGtt                                                                                                              <  2:517462‑M1/41‑1 (MQ=255)
                                                                                    gTGGAAGTTGCCGGTTTTGAAGACATCGTCGTCGAGCACCACTTTCTCTTTCAGCGCCTGCGCTAAACGATCGAGCCAGCTTAATTGCTCGCTGCTCCAGTCGTTTTCGCCCTTAATGCGCGTCAGCGCGTGATc                                                                                >  4:744708/1‑135 (MQ=255)
                                                                                    gTGGAAGTTGCCGGTTTTGAAGACATCGTCGTCGAGCACCACTTTCTCTTTCAGCGCCTGCGCTAAACGATCGAGCCAGCTTAATTGCTCGCTGCTCCAGTCGTTTTCGCCCTTAATGCGCGTCAGCGCGTGATc                                                                                <  3:744708/135‑1 (MQ=255)
                                                                                    gTGGAAGTTGCCGGTTTTGAAGACATCGTCGTCGAGCACCACTTTCTCTTTCAGCGCCTGCGCTAAACGATCGAGCCAGCTTAATTGCTCGCTGCTCCAGTCGTTTTCGCCCTTAATGCGCGTCAGCGCGTGATCGACACGTTCCTCaa                                                                  >  1:1352653/1‑149 (MQ=255)
                                                                                            tGCCGGTTTTGAAGACATCGTCGTCGAGCACCACTTTCTCTTTCAGCGCCTGCGCTAAACGATCGAGCCAGCTTAATTGCTCGCTGCTCCAGTCGTTTTCGCCCTTAATGCGCGTCAGCGCGTGATCGACACGTTCCTCAAACGGTTTc                                                          <  4:1219175/149‑1 (MQ=255)
                                                                                              ccGGTTTTGAAGACATCGTCGTCGAGCACCACTTTCTCTTTCAGCGCCTGCGCTAAACGATCGAGCCAGCTTAATTGCTCGCTGCTCCAGTCGTTTTCGCCCTTAATGCGCGTCAGCGCGTGATCGACACGTTCCTCAAACGGTTTCAg                                                        >  1:143806/1‑149 (MQ=255)
                                                                                                       aaGACATCGTCGTCGAGCACCACTTTCTCTTTCAGCGCCTGCGCTAAACGATCGAGCCAGCTTAATTGCTCGCTGCTCCAGTCGTTTTCGCCCTTAATGCGCGTCAGCGCGTGATCGACACGt                                                                         >  4:377162/1‑123 (MQ=255)
                                                                                                       aaGACATCGTCGTCGAGCACCACTTTCTCTTTCAGCGCCTGCGCTAAACGATCGAGCCAGCTTAATTGCTCGCTGCTCCAGTCGTTTTCGCCCTTAATGCGCGTCAGCGCGTGATCGACACGTTCCTCAAACGGttt                                                           >  4:1786409/1‑137 (MQ=255)
                                                                                                       aaGACATCGTCGTCGAGCACCACTTTCTCTTTCAGCGCCTGCGCTAAACGATCGAGCCAGCTTAATTGCTCGCTGCTCCAGTCGTTTTCGCCCTTAATGCGCGTCAGCGCGTGATCGACACGTTCCTCAAACGGttt                                                           <  1:91580/137‑1 (MQ=255)
                                                                                                       aaGACATCGTCGTCGAGCACCACTTTCTCTTTCAGCGCCTGCGCTAAACGATCGAGCCAGCTTAATTGCTCGCTGCTCCAGTCGTTTTCGCCCTTAATGCGCGTCAGCGCGTGATCGACACGTTCCTCAAACGGttt                                                           >  2:91580/1‑137 (MQ=255)
                                                                                                       aaGACATCGTCGTCGAGCACCACTTTCTCTTTCAGCGCCTGCGCTAAACGATCGAGCCAGCTTAATTGCTCGCTGCTCCAGTCGTTTTCGCCCTTAATGCGCGTCAGCGCGTGATCGACACGTTCCTCAAACGGTTTCAGCGCATCGcc                                               >  3:659259/1‑149 (MQ=255)
                                                                                                        aGACATCGTCGTCGAGCACCACTTTCTCTTTCAGCGCCTGCGCTAAACGATCGAGCCAGCTTAATTGCTCGCTGCTCCAGTCGTTTTCGCCCTTAATGCGCGTCAGCGCGTGATCGACACGTTCCTCAAACGGTTTCAGCGCATCGccc                                              >  4:152859/1‑149 (MQ=255)
                                                                                                            aTCGTCGTCGAGCACCACTTTCTCTTTCAGCGCCTGCGCTAAACGATCGAGCCAGCTTAATTGCTCGCTGCTCCAGTCGTTTTCGCCCTTAATGCGCGTCa                                                                                          >  3:822928/1‑101 (MQ=255)
                                                                                                              cgtcgtcgAGCACCACTTTCTCTTTCAGCGCCTGCGCTAAACGATCGAGCCAGCTTAATTGCTCGCTGCTCCAGTCGTTTTCGCCCTTAATGCGCGTCAGCGCGTGATCGACACGTTCCTCAAACGGTTTCAGCGCATCGCCCACCGCa                                        <  2:9640/149‑1 (MQ=255)
                                                                                                              cgtcgtcgAGCACCACTTTCTCTTTCAGCGCCTGCGCTAAACGATCGAGCCAGCTTAATTGCTCGCTGCTCCAGTCGTTTTCGCCCTTAATGCGCGTCAGCGCGTGATCGACACGTTCCCCAAACGGTTTCAGCGCATCGCCCACCGCa                                        <  4:204079/149‑1 (MQ=255)
                                                                                                               gtcgtcgAGCACCACTTTCTCTTTCAGCGCCTGCGCTAAACGATCGAGCCAGCTTAATTGCTCGCTGCTCCAGTCGTTTTCGCCCTTAATGCGCGTCAGCGCGTGATCGACACGTTCCTCAAACGGTTTCAGCGCATCGCCCACCGCAg                                       >  4:690676/1‑149 (MQ=255)
                                                                                                                     gAGCACCACTTTCTCTTTCAGCGCCTGCGCTAAACGATCGAGCCAGCTTAATTGCTCGCTGCTCCAGTCGTTTTCGCCCTTAATGCGCGTCAGCGCGTGATCGACACGTTCCTCAAACGGTTTCAGCGCATCGCCCACCGCAGCgcggc                                 <  2:1760788/149‑1 (MQ=255)
                                                                                                                     gAGCACCACTTTCTCTTTCAGCGCCTGCGCTAAACGATCGAGCCAGCTTAATTGCTCGCTGCTCCAGTCGTTTTCGCCCTTAATGCGCGTCAGCGCGTGATCGACACGTTCCTCAAACGGTTTCAGCGCATCGCCCACCGCAGCgcggc                                 >  2:266377/1‑149 (MQ=255)
                                                                                                                        caccacTTTCTCTTTCAGCGCCTGCGCTAAACGATCGAGCCAGCTTAATTGCTCGCTGCTCCAGTCGTTTTCGCCCTTAATGCGCGTCAGCGCGTGATCGACACGTTCCTCAAACGGTTTCAgc                                                       >  4:1312964/1‑124 (MQ=255)
                                                                                                                        caccacTTTCTCTTTCAGCGCCTGCGCTAAACGATCGAGCCAGCTTAATTGCTCGCTGCTCCAGTCGTTTTCGCCCTTAATGCGCGTCAGCGCGTGATCGACACGTTCCTCAAACGGTTTCAgc                                                       <  3:1312964/124‑1 (MQ=255)
                                                                                                                         accacTTTCTCTTTCAGCGCCTGCGCTAAACGATCGa                                                                                                                                             >  2:144384/1‑37 (MQ=255)
                                                                                                                         accacTTTCTCTTTCAGCGCCTGCGCTAAACGATCGa                                                                                                                                             <  1:144384/37‑1 (MQ=255)
                                                                                                                          ccacTTTCTCTTTCAGCGCCTGCGCTAAACGATCGAGCCAGCTTAATTGCTCGCTGCTCCAGTCGTTTTCGCCCTTAATGCGCGTCAGCGCGt                                                                                    <  3:1877027/93‑1 (MQ=255)
                                                                                                                            acTTTCTCTTTCAGCGCCTGCGCTAAACGATCGAGCCAGCTTAATTGCTCGCTGCTCCAGTCGTTTTCGCCCTTAATGCGCGTCAGCGCGTGAt                                                                                 <  4:1023853/94‑1 (MQ=255)
                                                                                                                            acTTTCTCTTTCAGCGCCTGCGCTAAACGATCGAGCCAGCTTAATTGCTCGCTGCTCCAGTCGTTTTCGCCCTTAATGCGCGTCAGCGCGTGAt                                                                                 >  3:1023853/1‑94 (MQ=255)
                                                                                                                            acTTTCTCTTTCAGCGCCTGCGCTAAACGATCGAGCCAGCTTAATTGCTCGCTGCTCCAGTCGTTTTCGCCCTTAATGCGCGTCAGCGCGTGATCGACACGTTCCTCAAACGGTTTCAGCGCATCGCCCACCGCAGCGCGGCGAATATg                          >  1:1847883/1‑149 (MQ=255)
                                                                                                                             cTTTCTCTTTCAGCGCCTGCGCTAAACGATCGAGCCAGCTTAATTGCTCGCTGCTCCAGTCGTTTTCGCCCTTAATGCGCGTCAGCGCGTGAt                                                                                 <  2:134004/93‑1 (MQ=255)
                                                                                                                             cTTTCTCTTTCAGCGCCTGCGCTAAACGATCGAGCCAGCTTAATTGCTCGCTGCTCCAGTCGTTTTCGCCCTTAATGCGCGTCAGCGCGTGATCGACACGTTCCTCAAACGGTTTCAGCGCATCGCCCACCGCAGCGCGGCGAATATGa                         <  3:152859/149‑1 (MQ=255)
                                                                                                                              tttCTCTTTCAGCGCCTGCGCTAAACGATCGAGCCAGCTTAATTGCTCGCTGCTCCAGTCGTTTTCGCCCTTAATGCGCGTCAGCGCGTGATCGACACGTTc                                                                       <  1:203672/102‑1 (MQ=255)
                                                                                                                              tttCTCTTTCAGCGCCTGCGCTAAACGATCGAGCCAGCTTAATTGCTCGCTGCTCCAGTCGTTTTCGCCCTTAATGCGCGTCAGCGCGTGATCGACACGTTc                                                                       >  2:203672/1‑102 (MQ=255)
                                                                                                                              tttCTCTTTCAGCGCCTGCGCTAAACGATCGAGCCAGCTTAATTGCTCGCTGCTCCAGTCGTTTTCGCCCTTAATGCGCGTCAGCGCGTGATCGACACGTTCCTCAAACGGTTTCAGCGCATCGCCCACCGCAGCGCGGCGAat                             <  2:482223/144‑1 (MQ=255)
                                                                                                                                tctcTTTCAGCGCCTGCGCTAAACGATCGAGCCAGCTTAATTGCTCGCTGCTCCAGTCGTTTTCGCCCTTAATGCGCGTCAGCGCGTGATCGACACGTTCCTCAAACGGTTTCAGCGCATCGCCCACCGCAGCGCGGCGAATATGa                         <  3:335341/146‑1 (MQ=255)
                                                                                                                                tctcTTTCAGCGCCTGCGCTAAACGATCGAGCCAGCTTAATTGCTCGCTGCTCCAGTCGTTTTCGCCCTTAATGCGCGTCAGCGCGTGATCGACACGTTCCTCAAACGGTTTCAGCGCATCGCCCACCGCAGCGCGGCGAATATGa                         >  4:335341/1‑146 (MQ=255)
                                                                                                                                 ctctTTCAGCGCCTGCGCTAAACGATCGAGCCAGCTTAATTGCTCGCTGCTCCAGTCGTTTTCGCCCTTAATGCGCGTCAGCGCGTGATCGACACGTTCCTCAAACGGTTTCAGCGCATCGCCCa                                             <  4:1608719/125‑1 (MQ=255)
                                                                                                                                 ctctTTCAGCGCCTGCGCTAAACGATCGAGCCAGCTTAATTGCTCGCTGCTCCAGTCGTTTTCGCCCTTAATGCGCGTCAGCGCGTGATCGACACGTTCCTCAAACGGTTTCAGCGCATCGCCCa                                             >  3:1608719/1‑125 (MQ=255)
                                                                                                                                 ctctTTCAGCGCCTGCGCTAAACGATCGAGCCAGCTTAATTGCTCGCTGCTCCAGTCGTTTTCGCCCTTAATGCGCGTCAGCGCGTGATCGACACGTTCCTCAAACGGTTTCAGCGCATCGCCCACCGCAGCGCGGCGAATATGACCaa                     <  2:1047409/149‑1 (MQ=255)
                                                                                                                                   ctTTCAGCGCCTGCGCTAAACGATCGAGCCAGCTTAATTGCTCGCTGCTCCAGTCGTTTTCGCCCTTAATGCGCGTCAGCGCGTGATCGACACGTTCCTCAAACGGTTTCAGCGCATCGCCCACCGCAGCGCGGCGAATATGACCAATc                   <  3:355680/149‑1 (MQ=255)
                                                                                                                                        aGCGCCTGCGCTAAACGATCGAGCCAGCTTAATTGCTCGCTGCTCCAGTCGTTTTCGCCCTTAATGCGCGTCAGCGCGTGATCGACACGTTCCTCAAACGGTTTCAGCGCATCGCCCACCGCAGCGCGGCGAATATGACCAATCAGCCg              >  1:333006/1‑149 (MQ=255)
                                                                                                                                              tGCGCTAAACGATCGAGCCAGCTTAATTGCTCGCTGCTCCAGTCGTTTTCGCCCTTAATGCGCGTCAGCGCGTGATCGACACGTTCCTCAAACGGTTTCAGCGCATCGCCCACCGCAGCGCGGCGAATATGACCAATCAGCCGggcggc        >  1:1765961/1‑149 (MQ=255)
                                                                                                                                                 gcTAAACGATCGAGCCAGCTTAATTGCTCGCTGCTCCAGTCGTTTTCGCCCTTAATGCGCGTCAGCGCGTGATCGACACGTTc                                                                       >  3:908033/1‑83 (MQ=255)
                                                                                                                                                 gcTAAACGATCGAGCCAGCTTAATTGCTCGCTGCTCCAGTCGTTTTCGCCCTTAATGCGCGTCAGCGCGTGATCGACACGTTc                                                                       <  4:908033/83‑1 (MQ=255)
                                                                                                                                                    AAACGATCGAGCCAGCTTAATTGCTCGCTGCTCCAGTCGTTTTCGCCCTTAATGCGCGTCAGCGCGTGATCGACACGTTCCTCAAACGGTTTCAGCGCATCGCCCACCGCAGCGCGGCGAATATGACCAATCAGCCGGGCGGCGATATc  >  4:659259‑M1/1‑149 (MQ=255)
                                                                                                                                                    AAACGATCGAGCCAGCTTAATTGCTCGCTGCTCCAGTCGTTTTCGCCCTTAATGCGCGTCAGCGCGTGATCGACACGTTCCTCAAACGGTTTCAGCGCATCGCCCACCGCAGCGCGGCGAATATGACCAATCAGCCGGGCGGCGATATc  <  3:457238‑M1/149‑1 (MQ=255)
                                                                                                                                                    AAACGATCGAGCCAGCTTAATTGCTCGCTGCTCCAGTCGTTTTCGCCCTTAATGCGCGTCAGCGCGTGATCGACACGTTCCTCAAACGGTTTCAGCGCATCGCCCACCGCAGCGCGGCGAATATGACCAATCAGCCGGGCGGCGATATc  >  2:1765961‑M1/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CCAGATATAATCGCTGAATTTGCCCAGCAGGGTATCGAGATTATCGTCAAAGGTTCTTTGCAGCATCGCCTTCCCGCCGCGACGGTGGAAGTTGCCGGTTTTGAAGACATCGTCGTCGAGCACCACTTTCTCTTTCAGCGCCTGCGCTAAACGATCGAGCCAGCTTAATTGCTCGCTGCTCCAGTCGTTTTCGCCCTTAATGCGCGTCAGCGCGTGATCGACACGTTCCTCAAACGGTTTCAGCGCATCGCCCACCGCAGCGCGGCGAATATGACCAATCAGCCGGGCGGCGATATC  >  NC_000913/4583264‑4583560

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑
Reads not counted as support for junction
read_name Not counted due to insufficient overlap past the breakpoint.
read_name Not counted due to not crossing MOB target site duplication.