Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 4,583,270 3543 bp→TGTAGGCTGGAGCTGCTTCGAAGTTCCTATACTTTCTAGAGAATAGGAACTTCGAACTGCAGGTCGACGGATCCCCGGAAT [hsdR] [hsdR]

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 4583270 4586812 3543 81 [0] [0] 84 [hsdR] [hsdR]

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 45832690 (0.000)64 (0.960)
+TGTAGGCTGGAGCTGCTTCGAAGTTCCTATACTTTCTAGAGAATAGGAACTTCGAACTGCAGGTCGACGGATCCCCGGAAT
34/424 0.6 100% coding (3493/3513 nt) hsdR type I restriction enzyme EcoKI endonuclease subunit
?NC_000913 4586813 = 0 (0.000)intergenic (‑52/‑136) hsdR/mrr type I restriction enzyme EcoKI endonuclease subunit/Type IV methyl‑directed restriction enzyme EcoKMrr

AGTTTTGATAGGAAACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTCATCGAGCCGCAGAAGAAGAAATTGCCGAGGGTAATATACACAAAATCATTCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATATAATCGCTGAATTTGCCCAGCAGGGTATCGAGATTATCGTCAAAGGTTCTTTGCAGCATCGCCTTCCCGCCGCGACGGTGGAAGTTGCCGGTTTTGAAGACATCGTCGTCGAGCACCACTTTCTCTTTCAGCGCCTGCGCTAAACGATCGAGCCAGCTTAATTGCTCGCTGCTCCAGTCGTTTTCGCCCTTAATGCGCGTCAGCGCGTGATCGACACGTTCCTCAAACGGTTTCAGCGCATCGCCCACCGCAGCGCGGCGAATATGACCAATCAGCCGGGCGGCGATATCTTCAT  >  NC_000913/4582970‑4583565
                                                                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                                                        
aGTTTTGATAGGAAACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTAATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTCATCGAGCCGCAGAAGAAGAAATTGCCGAGGGTAATATACACAAAATCATTCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGat                                                                                                                                                                                                                                                                                                          <  1:128616/300‑1 (MQ=255)
             aaCGCCGCCATCGCGCANGTTGTCGCACAGCTTCCANAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTCATCGAGCCGCAGAAGAAGAAATTGCCGAGGGTAATATACACAAAATCATTCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCNGATtgtaggctggagc                                                                                                                                                                                                                                                                                             <  1:629788‑M1/300‑14 (MQ=255)
             aaCGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTCATCGAGCCGCAGAAGAAGAAATTGCCGAGGGTAATATACACAAAATCATTCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATtgtaggctggagc                                                                                                                                                                                                                                                                                             <  1:710906‑M1/300‑14 (MQ=255)
                   gccATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTCATCGAGCCGCAGAAGAAGAAATTGCCGAGGGTAATATACACAAAATCATTCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATtgtaggctggagctgctt                                                                                                                                                                                                                                                                                        <  1:583611‑M1/299‑19 (MQ=255)
                      aTCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTCATCGAGCCGCAGAAGAAGAAATTGCCGAGGGTAATATACACAAAATCATTCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCAGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATtgtaggctggagctgcttcgaa                                                                                                                                                                                                                                                                                    >  1:198247‑M1/1‑278 (MQ=255)
                              ggTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTCATCGAGCCGCAGAAGAAGAAATTGCCGAGGGTAATATACACAAAATCATTCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCACGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCAAGATtgtaggctggagctgcttcgaagttcctat                                                                                                                                                                                                                                                                            >  2:459454‑M1/1‑270 (MQ=255)
                                ttGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTCATCGAGCCNCAGAAGAAGAAATTGCCGAGGGTAATATACACAAAATCATTCAGGTTGCATCAATGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATtgtaggctggagctgcttcgaagttcctatac                                                                                                                                                                                                                                                                          <  2:566393‑M1/300‑33 (MQ=255)
                                ttGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTCATCGAGCCGCATAAGAAGAAATTGCCGAGGGTAATATACACAAAATCATTCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATtgtaggctggagctgcttcgaagttcctatac                                                                                                                                                                                                                                                                          <  2:479932‑M1/300‑33 (MQ=255)
                                         agcttccacagcttcGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTCATCGAGCCGCAGAAGAAGAAATTGCCGAGGGTAATATACACAAAATCATTCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCAAGATtgtaggctggagctgcttcgaagttcctatactttctaga                                                                                                                                                                                                                                                                  >  1:282520‑M1/1‑259 (MQ=255)
                                          gcttccacagcttcGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTCATCGAGCCGCAGAAGAAGAAATTGCAGAGGGTAATATACACAAAATCATTCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCAAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTACAAGATtgtaggctggagctgcttcgaagttcctatactttctaga                                                                                                                                                                                                                                                                  >  2:367633‑M1/1‑258 (MQ=255)
                                           cttccacagcttcGCGACCAGATCTTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTCATCGAGCCGCAGAAGAAGAAATTGCCGAGGGTAATAGACACAAAATCATCCAGGTTGCNTCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGATCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATtgtaggctggagctgcttcgaagttcctatactttctag                                                                                                                                                                                                                                                                   <  2:230918‑M1/296‑40 (MQ=255)
                                                  atcttcgcGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTCATCGATCCGCAGAAGAAGAAATTGCCGAGGGTAATATACACAAAATCATTCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATtgtaggctggagctgcttcgaagttcctatactttctagagaataggaac                                                                                                                                                                                                                                                        <  2:281228‑M1/298‑51 (MQ=255)
                                                              aGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTCATCGAGCCGCAGAAGAAGAAATTGCCGAGGGTAATATACACAAAATCATTCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATtgtaggctggagctgcttcgaagttcctatactttctagagaataggaacttcgaactgcag                                                                                                                                                                                                                                            <  2:54733‑M1/300‑63 (MQ=255)
                                                                 tCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTCATCGAGCCGCAGAAGAAGAAATTGCCGAGGGTAATATACACAAAATCATTCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATtgtaggctggagctgcttcgaagttcctatactttctagagaataggaacttcgaactgcaggtc                                                                                                                                                                                                                                         <  1:250097‑M1/300‑66 (MQ=255)
                                                                 tCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTCATCGAGCCGCAGAAGAAGAAATTGCCGAGGGTAATATACACAAAATCATTCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCAAGATtgtaggctggagctgcttcgaagttcctatactttctagagaataggaacttcgaactgcaggtc                                                                                                                                                                                                                                         >  1:382575‑M1/1‑235 (MQ=255)
                                                                    ttATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTCATCGAGCCGCAGAAGAAGAAATTGCCGAGGGTAATATACACAAAATCATTCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATtgtaggctggagctgcttcgaagttcctatactttctagagaataggaacttcgaactgcaggtcgac                                                                                                                                                                                                                                      >  1:569488‑M1/1‑232 (MQ=255)
                                                                              ttGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTCATCGAGCCGCAGAAGAAGAAATTGCCGAGGGTAATATACACAAAATCATTCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATtgtaggctggagctgcttcgaagttcctatactttctagagaataggaacttcgaactgcaggtcgacggatccccgg                                                                                                                                                                                                                            <  2:180296‑M1/300‑79 (MQ=255)
                                                                                       ccGTAAATTAAGCAGCGGCCCAAATTCCTCGCGCCTCCGCAGAAGAACTTGCCGAGGGTCATATACACAAAATCATTCAGGTTGCATCACGACGGCAAGTGAGTGAATCCACGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATtgtaggctggagctgcttcgaagttcatatactttatagagaataggaacttcgaactgcaggtcgacggatccccggcatcatctc                                                                                                                                                                                                                   >  1:507042‑M1/1‑213 (MQ=255)
                                                                                       ccGTAAATTAAGCAGCGGCCCAAATTCATCGAGCCGCAGAAGAAGAAATTGCCGAGGGTAATATACACAAAATCATTCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATtgtaggctggagctgcttcgaagttcctatactttctagagaataggaacttcgaactgcaggtcgacggatccccggaatcatct                                                                                                                                                                                                                    >  1:507323‑M1/1‑213 (MQ=255)
                                                                                             aTTAAGCAGCGGCCCAAATTCATCGAGCCGCAGAAGAAGAAATTGCCGAGGGTAATATACACAAAATCATTCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATtgtaggctggagctgcttcgaagttcctatactttctagagaataggaacttcgaactgcaggtcgacggatccccggaatcatctcatctta                                                                                                                                                                                                             >  1:216485‑M1/1‑207 (MQ=255)
                                                                                                       ggCCCAAATTCATCGAGCCGCAGAAGAAGAAATTGCCGAGGGTAATATACACAAAATCATTCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATtgtaggctggagctgcttcgaagttcctatactttctagagaataggaacttcgaactgcaggtcgacggatccccggaatcatctcatcttagctttctgta                                                                                                                                                                                                   >  1:487540‑M1/1‑197 (MQ=255)
                                                                                                        gCCCAAATTCATCGAGCCGCAGAAGAAGAAATTGCCGAGGGTAATATACACAAAATCATTCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATtgtaggctggagctgcttcgaagttcctatactttctagagaataggaacttcgaactgcaggtcgacggatccccggaatcatctcatcttagctttctgtac                                                                                                                                                                                                  <  1:17860‑M1/300‑105 (MQ=255)
                                                                                                            aaaTTCATCGAGCCGCAGAAGAAGAAATTGCCGAGGGTAATATACACAAAATCATTCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATtgtaggctggagctgcttcgaagttcctatactttctagagaataggaacttcgaactgcaggtcgacggatccccggaatcatctcatcttagctttctgtacc                                                                                                                                                                                                 <  1:382753‑M1/297‑106 (MQ=255)
                                                                                                             aaTTCATCGAGCCGCAGAAGAAGAAATTGCCGAGGGTAATATACACAAAATCATTCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATtgtaggctggagctgcttcgaagttcctatactttctagagaataggaacttcgaactgcaggtcgacggatccccggaatcatctcatcttagctttctgtacctttc                                                                                                                                                                                             <  1:735649‑M1/300‑110 (MQ=255)
                                                                                                                  aTCGAGCCGCAGAAGAAGAAATTGCCGAGGGTAATATACACAAAATCATTCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATtgtaggctggagctgcttcgaagttcctatactttctagagaataggaacttcgaactgcaggtcgacggatccccggaatcatctcatcttagctttctgtacctttccgggc                                                                                                                                                                                        <  2:717660‑M1/300‑115 (MQ=255)
                                                                                                                              aagaagaaATTGCCGAGGGTAATATACACAAAATCATGCAGGTTGCATCAATGCGGCAAGTGAGTGAATCCCCTGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTAGACGTGTCAGGCCAGCNCGTCCCAGATtgtaggctggagctgcttcgaagttcctatactttctagagaataggaacttcgaactgcaggtcgacggatcccctgaatcatctcatcttagctttctgtacctttccgggcaatgaccacgg                                                                                                                                                                             <  2:520698‑M1/299‑126 (MQ=255)
                                                                                                                                    aaaTTGCCGAGGGTAATATACACAAAATCATTCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATtgtaggctggagctgcttcgaagttcctatactttctagagaataggaacttcgaactgcaggtcgacggatccccggaatcatctaatcttagttttttgtacctttccgggcaatgaccactgttacagc                                                                                                                                                                      >  2:407112‑M1/1‑168 (MQ=255)
                                                                                                                                                  aaTATACACAAAATCATTCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATtgtaggctggagctgcttcgaagttcctatactttctagagaataggaacttcgaactgcaggtcgacggatccccggaatcatctcatcttagctttctgtacctttccgggcaatgaccacggtcacagcaactgactcatttc                                                                                                                                                        >  2:438153‑M1/1‑154 (MQ=255)
                                                                                                                                                       acacAAAATCATTCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATtgtaggctggagctgcttcgaagttcctatactttctagagaataggaacttcgaactgcaggtcgacggatccccggaatcatctcaccttagctttctgtacctttccgggcaatgaccacggtcacagcaactgactcatttctaacg                                                                                                                                                   <  2:279276‑M1/300‑152 (MQ=255)
                                                                                                                                                         acaAAATCACTCAGGTTGCATCAAGGCGGCAAACGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGATCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAATGTATACGTTTCAGGCCAGCTCGTCCCAAATtttaggctggagctgcttcgaagttcctatacttgctagagaataggaacttcgaactgcaggtcgacggatccccggaatcatctcatcttagctttctgtacctttccgggcaatgaccacggtcacagcaactgactcatttctaccc                                                                                                                                                   <  2:271358‑M1/298‑152 (MQ=255)
                                                                                                                                                                     gggtTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATtgtaggctggagctgcttcgaagttcctatactttctagagaataggaacttcgaactgcaggtcgacggatccccggaatcatctcatcttagctttctgtacctttccgggcaatgaccacggtcacagcaactgactcatttctaacgtgttcgtctatttt                                                                                                                                     >  1:53103‑M1/2‑135 (MQ=255)
                                                                                                                                                                     aGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATtgtaggctggagctgcttcgaagttcctatactttctagagaataggaacttcgaactgcaggtcgacggatccccggaatcatctcatcttagctttctgtacctttccgggcaatgaccacggtcacagcaactgactcatttctaacgtgttcgtctatttt                                                                                                                                     >  1:51935‑M1/1‑135 (MQ=255)
                                                                                                                                                                     aGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCAGGGAGCGTACCGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATtgtaggctggagctgcttcgaagttcctatactttctagagaataggaacttcgaactgcaggtcgacggatccccggaatcatctaatcttagctttctgtacctttccgggcaatgaccacggtcacagcaactgactcatttctcaagtgttcgtctctttt                                                                                                                                     >  1:51924‑M1/1‑135 (MQ=255)
                                                                                                                                                                           cATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATtgtaggctggagctgcttcgaagttcctatactttctagagaataggaacttcgaactgcaggtcgacggatccccggaancatctcatcttagctttctgtacctttccgggcaatgaccacggtcacagcaactgactcatttctaacgtgttcgtctatttttgtagt                                                                                                                               >  2:567547‑M1/1‑129 (MQ=255)
                                                                                                                                                                                aGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATtgtaggctggagctgcttcgaagttcctatactttctagagaataggaacttcgaactgcaggtcgacggatccccggaatcatctcatcttagctttctgtacctttccgggcaatgaccacggtcacagcaactgactcatttctaacgtgttcgtctatttttgtagtgctat                                                                                                                          <  1:557215‑M1/300‑177 (MQ=255)
                                                                                                                                                                                  gcggcAAGTGAGTGAATCCCCGTGAGCGTACAGAGGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATtgtaggctggagcagcttcgaagttcctatactttctagagaataggaacttcgaactgcaggtcgacggatccccggaatcatctcatcttagctttctgtacctttccgggcaatgaccacggtcacagcaactgactcatttctaacgtgttcgtctatttttgtagtgctatag                                                                                                                        <  2:601047‑M1/300‑179 (MQ=255)
                                                                                                                                                                                  gcggcAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATtgtaggctggagctgcttcgaagttcctatactttctagagaataggaacttcgaactgcaggtcgacggatccccggaatcatctcatcttagctttctgtacctttccgggcaatgaccacggtcacagcaactgactcatttctaacgtgttcgtctatttttgtagtgctatag                                                                                                                        <  1:250102‑M1/300‑179 (MQ=255)
                                                                                                                                                                                  gcggcAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATtgtaggctggagctgcttcgaagttcctatactttctagagaataggaacgtcgaactgcaggtcgacggatccccggaatcatctcatcttagctttctgtacctttccgggcaatgaccacggtcacagcaactgactcatttctaacgtgttcgtctatttttgtagtgctatag                                                                                                                        <  1:250919‑M1/300‑179 (MQ=255)
                                                                                                                                                                                  gcggcAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGCGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTTGTCCCAGGTtggaggctggagctgcttcgaagttcctatacttgctagagaataggaacttcgaactgcaggtcgacggatccccggaatcatctcatcttagctttctgtacctttccgggcaatgaccacggtcacagcaactgactcatttctaacgtgtgcgtctatttttgtagtgctatag                                                                                                                        <  2:512552‑M1/300‑179 (MQ=255)
                                                                                                                                                                                    ggcAAGTGAGTGAATCCCCGGGAGCGTACATAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATGCGTGTCAGGCCAGCTCGTCCCAGATtgtaggcttgagctgcttctaagttcctatactttctagagaataggaacttcgaactgcaggtcgacggatccccggaatcatctcatcttagctttctttacctttccgtgcaatgaccacggtcacagcaactgactcatttctaacgtgttcgtctatttttgtagtgctatag                                                                                                                        <  2:76055‑M1/298‑179 (MQ=255)
                                                                                                                                                                                     gcAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATtgtaggctggagctgcttcgaagttcctatactttctagagaataggaacttcgaactgcagttcgacggatccccggaatcatctcatcttagctttctgtacctttccgggcaatgaccacggtcacagcaactgactcatttctaacgtgttcgtctatttttgtagtgctatag                                                                                                                        <  2:435476‑M1/297‑179 (MQ=255)
                                                                                                                                                                                        agtgagtgAATCCCCTGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGGCCCAGATtgtaggctggagctgcttcgaagttcctatactttctagagaataggaacttcgaactgcaggtcgacggatccccggaatcatctcatcttagctttctgtacctttccgggcaatgaccacggtcacagcaactgactcatttctaacgtgttcgtctatttttgtagtgctatagtagcc                                                                                                                   <  1:204250‑M1/299‑184 (MQ=255)
                                                                                                                                                                                         gtgagtgaATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATtgtaggctggagctgcttcgaagttcctatactttctagagaataggaacttcgaactgcaggtcgacggatccccggaatcatctcatcttagctttctgtacctttccgggcaatgaccacggtcacagcaactgactcatttctaacgtgttcgtctatttttgtagtgctatagtagccga                                                                                                                 >  1:448876‑M1/1‑115 (MQ=255)
                                                                                                                                                                                                             gTACAGAAGTACGTGACCGGGGTGAACGAGCGCGGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATtgtaggctggagctgcttcgaagttcctatactttctagagaataggaacttcgaactgcaggtcgacggatccccggaatcatctcatcttagctttctgtacctttccgggcaatgaccacggtcacagcaactgactcatttctcacgtgttcgtctatttttgtagtgctcgtagtagccgaaaaacatctacctgattct                                                                                             >  2:608811‑M1/1‑95 (MQ=255)
                                                                                                                                                                                                             gTACAGAAGTACGTGACCGGGGTGAACGAGCGCGGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATtgtaggctggagctgcttcgaagttcctatactttctagagaataggaacttcgaactgcaggtcgacggatccccggaatcatctcatcttagctttctgtacctttccgggcaatgaccacggtcacagcaactgactcatttctaacgtgttcgtctctttttgtagtgctcgtagtagcagaaaaacatctacctgattct                                                                                             >  2:609124‑M1/1‑95 (MQ=255)
                                                                                                                                                                                                             gTACAGAAGTACGTGACCGGGGTGAACGAGCGCGGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATtgtaggctggagctgcttcgaagttcctatactttctagagaataggaacttcgaactgcaggtcgacggatccccggaatcatctaatcttagctttctgtacctttccgggcaatgaccacggtcacagcaactgactcatttctaacgtgttcgtctatttttgtagtgctcgtagtagccgaaaaacatctacctgattct                                                                                             >  2:609125‑M1/1‑95 (MQ=255)
                                                                                                                                                                                                             gTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATtgtaggctggagctgcttcgaagttcctatactttctagagaataggaacttcgaactgcaggtcgacggatccccggaatcatctcatcttagctttctgtacctttccgggcaatgaccacggtcacagcaactgactcatttctaacgtgttcgtctatttttgtagtgctatagtagccgaaaaacatctacctgattct                                                                                              >  1:683277‑M1/1‑95 (MQ=255)
                                                                                                                                                                                                                         gttaccgGGGTGAACAAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCATGCCAGCTCGTCCCAGATtgtaggctggagctgcttcgaagttcctatactttctagagaataggaacttcgaactgcaggtcgacggatccccggaatcatctcatcttagctttctgtacctttccgggcaatgaccacggtcacagcaactgactcatttctaacgtgttcgtctatttttgtagtgctatagtagccgaaaaacatctacctgattctgcaaggatgtact                                                                                 <  2:180233‑M1/297‑218 (MQ=255)
                                                                                                                                                                                                                          tGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATtgtaggctggagctgcttcgaagttcctatactttctagagaataggaacttcgaactgcaggtcgacggatccccggaatcatctcatcttagctttctgtacctttcctggcaatgaccacggtcacagcaactgactcatttctaacgtgttcgtctatttttgtagtgctatagtagccgaaaaacatctacctgattctgcaaggatgtacta                                                                                <  2:727912‑M1/300‑219 (MQ=255)
                                                                                                                                                                                                                          tGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATtgtaggctggagctgcttcgaagttcctatactttctagagaataggaacttcgaactgcaggtcgacggatccccggaatcatctcatcttagctttctgtacctttccgggcaatgaccacggtcacagcaactgactcatttctaacgtgttcgtctatttttgtagtgctatagtagccgaaaaacatctacctgattctgcaaggatgtacta                                                                                <  1:601276‑M1/300‑219 (MQ=255)
                                                                                                                                                                                                                           gACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATtgtaggctggagctgcttcgaagttcctatactttctagagaataggaacttcgaactgcaggtcgacggatccccggaatcatctaatcttagctttctgtacctttccgggcaatgaccacggtcacagcaactgactnatttctaacgtgttcgtctatttttgtagtgctatagtagccgaaaaacatctacctgattctgcaaggatgtactat                                                                               >  1:141351‑M1/1‑81 (MQ=255)
                                                                                                                                                                                                                           gACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTAAGGCCAGCTCGTCCCAGATtgtaggctggagctgcttcgaagttcctatactttctagagaataggaacttcgaactgcaggtcgacggatccccggaatcatctcatcttagctttctgtacctttcctggcaatgaccacggtcacagcaactgactcatttctaacgtgttcgtctatttttgtagtgctatagtagccgaaaaacatctacctgattctgcaaggatgtacta                                                                                <  2:151263‑M1/299‑219 (MQ=255)
                                                                                                                                                                                                                               ggggTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATtgtaggctggagctgcttcgaagttcctatactttctagagaataggaacttcgaactgcaggtcgacggatccccggaatcatctcatcttagctttctgtacctttccgggcaatgaccacggtcacagcaactgactcatttctaacgtgttcgtctatttttgtagtgctatagtagccgaaaaacatctacctgattctgcaaggatgtactatga                                                                             <  2:593273‑M1/298‑222 (MQ=255)
                                                                                                                                                                                                                               ggggTGAACAAACGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTGTACGTGTCATGCCAGCTCGTCCCAGATtgtaggctggagctgcttcgaagttcctagactttctagagaataggaacttcgaactgcaggtcgacggatccccggaatcatctcatcttagctttcggtacccttccgggcaatgaccacggtcacagcaactgactcatttctaacgtgttcgtctatttttgtagtgctatagtagccgaaaaacatctacctgattctgcaaggatgtactatga                                                                             <  2:593510‑M1/298‑222 (MQ=255)
                                                                                                                                                                                                                                gggTGAACGAGCGCAGCCAACGCATAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATtgtaggctggagctgcttcgaagttcctatactttctagagaataggaacttcgaactgcaggtcgacggatccccggaatcatctcatcttagctttctgtacctttccgggcaatgaccacggtcacagcaactgactcatttctaacgtgttcgtctatttttgtagtgctatagtagccgaaaaacatctacctgattctgcaaggatgtactatgacgg                                                                          <  1:190444‑M1/300‑225 (MQ=255)
                                                                                                                                                                                                                                gggTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATtgtaggctggagctgcttcgaagttcctatactttctagagaataggaacttcgaactgcaggtcgacggatccccggaatcatctcatcttagctttctgtacctttccgggcaatgaccacggtcacagcaactgactcatttctaacgtgttcgtctatttttgtagtgctatagtagccgaaaaacatctacctgattctgcaaggatgtactatgac                                                                            >  1:25669‑M1/1‑76 (MQ=255)
                                                                                                                                                                                                                                        gAGCGCAGCCAACGCATAGGCAGCCTGAAGGATGAAGTGTATACGTGTCATGCCAGCTCGTCCCAGATtgtaggctggagctgcttcgaagttcctatactttctagagaataggaacttcgaactgcaggtcgacggatccccggaatcatctcatcttagctttctttacctttccgggcaatgaccacggtcacagcaactgacttatttctaacgtgttcgtctatttttgtagtgctatattagcctaaaaacatctacctgattctgcaaggatgtactatgacggttcctacc                                                                  <  2:303015‑M1/300‑233 (MQ=255)
                                                                                                                                                                                                                                        gAGCGCAGCCAACGCAGATGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATtgtaggctggagctgcttcgaagttcctatactttctagagaataggaacttcgaactgcaggtctacggatccccggaatcatctcatcttagctttctgtacctttccgggcaatgaccacggtcacagcaactgactcatttctaacgtgttcgtctatttttgtagtgctatagtagccgaaaaacatctacctgattctgcaaggatgtactatgacggttcctacc                                                                  <  2:161908‑M1/300‑233 (MQ=255)
                                                                                                                                                                                                                                                                 tGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATtgtaggctggagctgcttcgaagttcctatactttctagagaataggaacttcgaactgcaggtcgacggatccccggaatcatctcatcttagctttctgtacctttccgggcaatgaccacggtcacagcaactgactcatttctaacgtgttcgtctatttttgtagtgctatagtagccgaaaaacatctacctgattctgcaaggatgtactatgacggttcctacctatgacaaatttattgaacctgttc                                         <  1:678101‑M1/300‑258 (MQ=255)
                                                                                                                                                                                                                                                                     ggATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATtgtaggctggagctgcttcgaagttcctatactttctagagaataggaacttcgaactgcaggtcgacggatccccggaatcatctcatcttagctttctgtacctttccgggcaatgaccacggtcacagcaactgactcatttctaacgtgttcgtctatttttgtagtgctatagtagccgaaaaacatctacctgattctgcaaggatgtactatgacggttcctacctatgacaaatttattgaacctgttctgc                                      >  1:448344‑M1/1‑39 (MQ=255)
                                                                                                                                                                                                                                                                      gATGAAGTGTATACGTTTCATGCCAACTCGTCCCAGATgttaggctggagatgcttcgaatttccggtactttctagagaagaggaactttgaactgcagttcgacggatcccctgaatcatctcatcttagctttctgtacctttccgggcaatgaccacggtcacagcaactgactcatgtctaacgtgttcgtctatttttgtagtgctatagtagccgaaaaacatctacctgattctgcaaggatgtactatgacggttcctacctatgacaaatttattgaacctgttctgcgt                                    <  2:396179‑M1/300‑263 (MQ=255)
                                                                                                                                                                                                                                                                      gATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATtgtaggctggagctgcttcgaagttcctatactttctagagaataggaacttcgaactgcaggtcgacggatccccggaatcatctcatcttagctttctgtacctttccgggcaatgaccacggtcacagcaactgactcatttctaacgtgttcgtctatttttgtagtgctatagtagccgaaaaacatctacctgattctgcaaggatgtactatgacggttcctacctatgacaaatttattgaacctgttctgcgt                                    <  2:177340‑M1/300‑263 (MQ=255)
                                                                                                                                                                                                                                                                      gATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATtgtaggctggagctgcttcgaagttcctatactttctagagaataggaacttcgaactgcaggtcgacggatccccggaatcatctcatcttagctttctgtacctttccgggcaatgaccacggtcacagcaactgactcatttctaacgtgttcgtctatttttgtagtgctatagtagccgaaaaacatctacctgattctgcaaggatgtactatgacggttcctacctatgacaaatttattgaacctgttctgcgt                                    <  1:414135‑M1/300‑263 (MQ=255)
                                                                                                                                                                                                                                                                      gATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATtgtaggctggagctgcttcgaagttcctatactttctagagaataggaacttcgaacggcaggtcgacggatccccggaatcatctcatcttagctttctgtacctttccgggcaatgaccacggtcacagcaactgactcatttctaacgtgttcgtctatttttgtagtgctatagtagccgaaaaacatctacctgattctgcaaggatgtactatgacggttcctacctatgacaaatttattgaacctgttctgcgt                                    <  2:282520‑M1/300‑263 (MQ=255)
                                                                                                                                                                                                                                                                           aGTGTATACGTTTCAGGCCAGCTCGTCCCAGATtgtaggctggagctgcttcgaagttcctatactttctagagaataggaacttcgaactgcaggtcgacggatcccctgaatcatctcatcttagctttctgtacctttccgggcaatgaccacggtcacagcaaccgactcatttctaacgtgttcgtctatttttgtagtgctatagtagccgaaaaacatctacctgattctgcaaggatgtactatgacggttcctacctatgacaaatttattgaacctgttctgcgttatct                               <  2:198247‑M1/300‑268 (MQ=255)
                                                                                                                                                                                                                                                                           aGTGTATACGTGTCAGGCCAGCTCGTCCCAGATtgtaggcttgagctgcttcggagttcctatactttctagagaataggaacttcgaactgcaggtcgacggatccccggaatcatctcatcttagctttctgtacctttccgggcaatgaccacggtcacagcaactgactcatttctaacgtgttcgtctatttttgtagtgctatagtagccgaaaaacatctacctgattctgcaaggatgtactatgacggttcctacctatgacaaatttattgaacctgttctgcgttatct                               <  2:515760‑M1/300‑268 (MQ=255)
                                                                                                                                                                                                                                                                           aGTGTATACGTGTCAGGCCAGCTCGTCCCAGATtgtaggctggagctgcttcgaagttcctatactttctagagaataggaacttcgaactgcaggtcgacggatccccggaatcatctcatcttagctttctgtacctttccgggcaatgaccacggtcacagcaactgactcatttctaacgtgttcgtctatttttgtagtgctatagtagccgaaaaacatctacctgattctgcaaggatgtactatgacggttcctacctatgacaaatttattgaacctgttctgcgttatct                               <  1:563530‑M1/300‑268 (MQ=255)
                                                                                                                                                                                                                                                                            gtgtATACGTGTCAGGCCAGCTCGTCCCAGATtgtaggctggagctgcttcgaagttcctatactttctaggtaacaggaacttcgaactgcaggtcgacggatccccggaatcatctcatcttagatttctgtacctttccgggcaatgaacacggtcacagcaaatgactcatttctaacgtgtttgtctatttttgtagtgctatagtagccgcaaaacatctacctgcttctgcacggatgtactat                                                                               >  2:23989‑M1/1‑32 (MQ=255)
                                                                                                                                                                                                                                                                            gtgtATACGTGTCAGGCCAGCTCGTCCCAGATtgtaggctggagctgcttcgaagttcctatactttctagagaataggaacttcgaactgcaggtcgacggatccccggaatcatctcatcttagctttctgtacctttccgggcaatgaccacggtcacagcaactgactcatttctaacgtgttcgtctatttttgtagtgctatagtagccgaaaaacatctacctgattctgcaaggatgtactatgacggttcctacctatgacaaatttattgaacctgttctgcgttatct                               >  1:23841‑M1/1‑32 (MQ=255)
                                                                                                                                                                                                                                                                            gtgtATACGTGTCAGGCCAGCTCGTCCCAGATtgtaggctggagctgcttcgaagttcctatactttctagagaataggaacttcgaactgcaggtcgacggatccccggaatcatctcatcttagctttctgtacctttacgggcaatgaccacggtcacagcaactgactcatttctaacgtgttcgtctatttttgtagtgctatagtagccgaaaaacatctacctgattctgcaaggatgtactatgacggttcctacctatgacaaatttattgaacctgttctgcgttatct                               >  2:23683‑M1/1‑32 (MQ=255)
                                                                                                                                                                                                                                                                            gtgtATACGTGTCAGGCCAGCTCGTCCAAGATtgtaggctggagctgcttcgaagttcctatactttctcgagaataggaacttcgaactgcaggtcgacggatccccggaatcatctcatcttagctttctgtacctttccgggcactgaccccggtcacagcaactgactcatttctaacgtgttcgtctatttttgtagtgctatagtagccgaaaaacatctacctgattctgcaaggatgtcctatgacggttcatacctatgacaaatttattggacctgttctgcggtatct                               >  2:25258‑M1/1‑32 (MQ=255)
                                                                                                                                                                                                                                                                            gtgtATACGTGTCAGGCCAGCTCGTCCAAGATtgtaggctggagctgcttcgaagttcctatactttctagagaataggaacttcgaactgcaggtcgacggatccccggaatcatctcatcttagctttctgtacctttccgggcaatgaccacggtcacagcaactgactcatttctaacgtgttcgtctatttttgtagtgctatagtagccgaaaaacatctacctgattctgcaaggatgtactatgacggttcctacctatgacaaatttattgaacctgttctgcgttatct                               >  1:363860‑M1/1‑32 (MQ=255)
                                                                                                                                                                                                                                                                              gtATACGTGTCAGGCCAGCTCGTCCCAGATtgtaggctggagctgcttcgaagttcctatactttctagagaataggaacttcgaactgcaggtcgacggatccccggaatcatctaatcttagctttctgtacctttccgggcaatgaacacggtcacagcaactgactcatttctaacgtgttcgtctatttttgtagtgctatagtagccgaaaaacatctacctgattctgcaaagatgtactatgacggttcctacctatgacaaatttattgaacctgttctgcgttatctggc                            >  2:120538‑M1/1‑30 (MQ=255)
                                                                                                                                                                                                                                                                                          ggCCAGCTCGTCCCAGATtgtaggctggagctgcttcgaagttcctatactttctagagaataggaacttcgaactgcaggtcgacggatccccggaatcatctcatcttagctttctgtacctttccgggcaatgaccacggtcacagcaactgactcatttctaacgtgttcgtctatttttgtagtgctatagtagccgaaaaacatctacctgattctgcaaggatgtactatgacggttcctacctatgacaaatttattgaacctgttctgcgttatctggcaacaaaaccgga                >  1:270193‑M1/1‑18 (MQ=255)
                                                                                                                                                                                                                                                                                          ggCCAGCTCGTCCCAGATtgtaggctggagctgcttcgaagttcctatactttctagagaataggaacttcgaactgcaggtcgacggatccccggaatcatctcatcttagctttctgtacctttccgggcaatgaccacggtcacagcaactgactcatttctaacgtgttcgtctatttttgtagtgctatagtagccgaaaaacatctacctgattctgcaaggatgtactatgacggttcctaactatgacaaatttcttgaacctgttctgcgttatctggcaacaaaaccgga                >  1:649501‑M1/1‑18 (MQ=255)
                                                                                                                                                                                                                                                                                          ggCCAGCTCGTCCCAGATtgtaggctggagctgcttcgaagttcctatactttctagagaataggaacttcgaactgcaggtcgacggatccccggaatcatctcatcttagctttctgtacctttccgggcaatgaccacggtcaaagcaactgactcatttctaacgtgttcgtctatttttgtagtgctatagtagccgaaaaacatctacctgattctgcaaggatgtactatgacggttcctacctatgacaaatttcttgaacctgttctgcgttatctggcaacaaaaccgga                >  1:271684‑M1/1‑18 (MQ=255)
                                                                                                                                                                                                                                                                                          ggCCAGCTCGTCCCAGATtgtaggctggagctgcttcgaagttcctatactttctagagaataggaacttcgaactgcaggtcgacggatccccggaatcatctaatcttagctttctgtacctttccgggcaatgaccacggtcacagcaactgactcatttctaacgtgttcgtctatttttgtagtgctatagtagccgaaaaacatctacctgattctgcaaggatgtactatgacggttcctacctatgacaaatttattgaacctgttctgcgttatctggcaacaaaaccgga                >  1:270358‑M1/1‑18 (MQ=255)
                                                                                                                                                                                                                                                                                          ggCCAGCTCGTCCCAGATtgtaggctggagctgcttcgaagttcctatactttctagagaataggaacttcgaactgcaggtcgacggatccccggaatcatctaatcttagctttctgtacctttccgggcaatgaccacggtcacagcaactgactcatttctaacgtgttcgtctatttttgtagtgctatagtagccgaaaaacatctacctgattctgcaaggatgtactatgacggttcctaactatgacaaatttattgaacctgttctgcgttatctggcaacaaaaccgga                >  2:340260‑M1/1‑18 (MQ=255)
                                                                                                                                                                                                                                                                                                        aGATtgtaggctggagctgcttcgaagttcctatactttctagagaataggaacttcgaactgcaggtcgacggatccccggaatcatctcatcttagctttctgtacctttccgggcaatgaccacggtcacagcaactgactcatttctaacgtgttcgtctatttttgtagtgctatagtagccgaaaaacatctacctgattctgcaaggatgtactatgacggttcctacctatgacaaatttattgaacctgttctgcgttatctggcaacaaaaccggaaggtgcagccgcgc  >  2:425939‑M1/1‑4 (MQ=255)
                                                                                                                                                                                                                                                                                                        aGATtgtaggctggagctgcttcgaagttcctatactttctagagaataggaacttcgaactgcaggtcgacggatccccggaatcatctcatcttagctttctgtacctttccgggcaatgaccacggtcacagcaactgactcatttctaacgtgttcgtctatttttgtagtgctatagtagccgaaaaacatctacctgattctgcaaggatgtactatgacggttcctacctatgacaaatttattgaacctgttctgcgttatctggcaacaaaaccggaaggtgcagccg     <  1:677858‑M1/297‑294 (MQ=255)
                                                                                                                                                                                                                                                                                                        aGATtgtaggctggagctgcttcgaagttcctatactttctagagaataggaacttcgaactgcaggtcgacggatccccggaatcatctcatcttagctttctgtacctttccgggcaatgaccacggtcacagcaactgactcatttctaacgtgttcgtctatttttgtagtgctatagtagccgaaaaacatctacctgattctgcaaggatgtactatgacggttcctacctatgacaaatttattgaacctgttctgcgttatctggcaacaaaaccggaaggtgcagccg     <  2:448344‑M1/297‑294 (MQ=255)
                                                                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                                                        
AGTTTTGATAGGAAACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTCATCGAGCCGCAGAAGAAGAAATTGCCGAGGGTAATATACACAAAATCATTCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATATAATCGCTGAATTTGCCCAGCAGGGTATCGAGATTATCGTCAAAGGTTCTTTGCAGCATCGCCTTCCCGCCGCGACGGTGGAAGTTGCCGGTTTTGAAGACATCGTCGTCGAGCACCACTTTCTCTTTCAGCGCCTGCGCTAAACGATCGAGCCAGCTTAATTGCTCGCTGCTCCAGTCGTTTTCGCCCTTAATGCGCGTCAGCGCGTGATCGACACGTTCCTCAAACGGTTTCAGCGCATCGCCCACCGCAGCGCGGCGAATATGACCAATCAGCCGGGCGGCGATATCTTCAT  >  NC_000913/4582970‑4583565

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 24 ≤ 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.