Predicted mutation
evidence seq id position mutation annotation gene description
JC NC_000913 3,822,651 (CAGGAT)1→2 coding (252/2109 nt) spoT → bifunctional (p)ppGpp synthase/hydrolase SpoT

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 3822646 =0 (0.000)42 (0.930) 39/490 0.2 97.7% coding (247/2109 nt) spoT bifunctional (p)ppGpp synthase/hydrolase SpoT
?NC_000913 = 3822651 2 (0.040)coding (252/2109 nt) spoT bifunctional (p)ppGpp synthase/hydrolase SpoT

CCTTGTATCTGTTTGAAAGCCTGAATCAACTGATTCAAACCTACCTGCCGGAAGACCAAATCAAGCGTCTGCGGCAGGCGTATCTCGTTGCACGTGATGCTCACGAGGGGCAAACACGTTCAAGCGGTGAACCCTATATCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCGCTGCTGCATGACGTGATTGAAGATACTCCCGCCACCTACCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGGGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATAAGAAAGAGGCGCAGGCCGAAAACTTTCGCAAGATGATTATGGCGATGGTGCAGGATATCCGCGTCATCCTCATCAAACTTGCCGACCGTACCCACAACATGCGCACGCTGGGCTCACTTCGCCCGGACAAACGTCGCCGCATCGCCCGTGAAACTCTCGAAATTT  >  NC_000913/3822398‑3822895
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                         
atctgtttgaaagcctgaatcaactgattcaaacctacctgccggaagaccaaatcaagcgtctgcggcaggcgtatctcgttgcacgtgatgctcacgaggggcaaacacgttcaagcggtgaaccctatatcacgcacccggtagcggttgcctgcattctggccgagatgaaactcgactatgaaacgctgatggcggcgctgctgcatgacgtgattgaagatactcccgccacctaccaggatca                                                                                                                                                                                                                                                          >  2:145631‑M1/249‑250 (MQ=255)
      aTCTGTTTGAAAGCCTGAATCAACTGATTCAAACCTACCTGCCGGAAGACCAAATCAAGCGTCTGCGGCAGGCGTATCTCGTTGCACGTGATGCTCACGAGGGGCAAACACGTTCAAGCGGTGAACCCTATATCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCGCTGCTGCATGACGTGATTGAAGATACTCCCGCCACCTACCAGGATca                                                                                                                                                                                                                                                    >  2:145631‑M2/1‑248 (MQ=255)
           agcctgaatcaactgattcaaacctacctgccggaagaccaaatcaagcgtctgcggcaggcggatctcgttgcacgtgatgctcacgaggggcaaactcgttcaagcggtgaaccctatatcacgcacccggtagcggttgcctgcattctggccgagatgaaactcgactatgaaacgctgatggcggcgctgctgcatgacgtgattgaagatactcccgccacctaccaggatCAGGATATGGAACa                                                                                                                                                                                                                                              >  2:335142‑M1/238‑251 (MQ=255)
                 aGCCTGAATCAACTGATTCAAACCTACCTGCCGGAAGACCAAATCAAGCGTCTGCGGCAGGCGGATCTCGTTGCACGTGATGCTCACGAGGGGCAAACTCGTTCAAGCGGTGAACCCTATATCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCGCTGCTGCATGACGTGATTGAAGATACTCCCGCCACCTACCAGGATcaggatatggaaca                                                                                                                                                                                                                                        >  2:335142‑M2/1‑237 (MQ=255)
                     aactgattcaaacctgcctgactgaagaccaaatcaagcgtctgcggcaggcgtatctcggtgcacgtgatgctcacgaggggcaaacacgttcaagcggtgaaccctatatcacgcacccggtagcggttgcctgcattctggccgagatgaaactcgactatgaaacgctgatggcggcgctgctgcatgacgtgattgaagatactcccgccacctaccaggatCAGGATATGGAACAGCTTTTTGGt                                                                                                                                                                                                                                    >  2:429273‑M1/228‑251 (MQ=255)
                           aaCTGATTCAAACCTGCCTGACTGAAGACCAAATCAAGCGTCTGCGGCAGGCGTATCTCGGTGCACGTGATGCTCACGAGGGGCAAACACGTTCAAGCGGTGAACCCTATATCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCGCTGCTGCATGACGTGATTGAAGATACTCCCGCCACCTACCAGGATcaggatatggaacagctttttggt                                                                                                                                                                                                                              >  2:429273‑M2/1‑227 (MQ=255)
                              aaaccgacctgccggaagaccaaatcaagcgtctgcggcaggcgtatctcgttgcacgtgatgcgcacgaggggcaaacacgttcaagcggtgaaccctatatcacgcacccggtagcggttgcctgcattctggccgagatgaaactcgactatgaaacgctgatggcggcgctgctgcatgacgtgattgaagatactcccgccacctaccaggatCAGGATATGGAACAGCTTTTTGGTAAAAGCGTc                                                                                                                                                                                                                           >  2:117998‑M1/219‑251 (MQ=255)
                               aacctacctgccggaagaccaaatcaagcgtctgcggcaggcgtatctcgttgcacgtgatgctcacgaggggcaaacacgttcaagcggtgaaccctatatcacgcacccggtagcggttgcctgcattctggccgagatgaaactcgactatgaaacgctgatggcggcgctgctgcatgacgtgattgaagatactcccgccacctaccaggatCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCg                                                                                                                                                                                                                          <  2:339311‑M1/34‑1 (MQ=255)
                                    aaaCCGACCTGCCGGAAGACCAAATCAAGCGTCTGCGGCAGGCGTATCTCGTTGCACGTGATGCGCACGAGGGGCAAACACGTTCAAGCGGTGAACCCTATATCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCGCTGCTGCATGACGTGATTGAAGATACTCCCGCCACCTACCAGGATcaggatatggaacagctttttggtaaaagcgtc                                                                                                                                                                                                                     >  2:117998‑M2/1‑218 (MQ=255)
                                     cctgccggaagaccaaatcaagcgtctgcggcaggcgtatctcgttgcacgtgatgctcacgaggggcaaacacgttcaagcggtgaaccctatatcacgcacccggtagcggttgcctgcattctggccgagatgaaactcgactatgaaacgctgatggcggcgctgctgcctgacgtgattgaagatactcccgccacctaccaggaaCAGGATATGGAACAGCTTTTTGGAAAAAGCGTCGCCGAGc                                                                                                                                                                                                                    <  2:159431‑M1/40‑1 (MQ=255)
                                     aaCCTACCTGCCGGAAGACCAAATCAAGCGTCTGCGGCAGGCGTATCTCGTTGCACGTGATGCTCACGAGGGGCAAACACGTTCAAGCGGTGAACCCTATATCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCGCTGCTGCATGACGTGATTGAAGATACTCCCGCCACCTACCAGGATcaggatatggaacagctttttggtaaaagcgtcg                                                                                                                                                                                                                    <  2:339311‑M2/251‑35 (MQ=255)
                                           cctGCCGGAAGACCAAATCAAGCGTCTGCGGCAGGCGTATCTCGTTGCACGTGATGCTCACGAGGGGCAAACACGTTCAAGCGGTGAACCCTATATCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCGCTGCTGCCTGACGTGATTGAAGATACTCCCGCCACCTACCAGGAAcaggatatggaacagctttttggaaaaagcgtcgccgagc                                                                                                                                                                                                              <  2:159431‑M2/251‑41 (MQ=255)
                                                                     aggcgtatctcgttgcacgtgatgctcacgaggggcaaacacgttcaagcggtgaaccctatatcacgcacccggtagcggttgcctgcattctggccgagatgaaactcgactatgaaacgctgatagcggcgctgctgcatgacgtgattgaagatactcccgccacctaccaggatCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGGGGGTGTCGAAACTTGATAAACTc                                                                                                                                                                                    <  1:95965‑M1/72‑1 (MQ=255)
                                                                           aGGCGTATCTCGTTGCACGTGATGCTCACGAGGGGCAAACACGTTCAAGCGGTGAACCCTATATCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATAGCGGCGCTGCTGCATGACGTGATTGAAGATACTCCCGCCACCTACCAGGATcaggatatggaacagctttttggtaaaagcgtcgccgagctggtagagggggtgtcgaaacttgataaactc                                                                                                                                                                              <  1:95965‑M2/251‑73 (MQ=255)
                                                                               cgttgcacgtgatgctcacgcggggcaaacacgttcaagcggtgaacactatctcacgcacccggtagcggttgcctgcattctggccgagatgaaactcgactatgaaacgatgctggcggcgctgctgcatgacgtgattgaagatactcccgcaacctaccaggatCAGGATATGGAACAGCTTTTTGGTACAAGCGTCGCCGAGCTGGTAGAGGGGGTGTCGAAACTTGATAAACTCAAGTTCcgcg                                                                                                                                                                          <  1:83071‑M1/82‑1 (MQ=255)
                                                                               cgttgcacgtgatgctcacgaggggcaaacacgttcaagcggtgaaccctatatcacgcacccggtagcggttgcctgcattctggccgagatgaaactcgactatgaaacgctgatggcggcgctgctgcatgacgtgattgaagatactcccgccacctaccaggatCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGGGGGTGTCGAAACTTGATAAACTCAAGTTCcgcg                                                                                                                                                                          >  2:145344‑M1/170‑251 (MQ=255)
                                                                                gttgcacgtgatgctcacgaggggcaaacacgctcacgcggtgaaccctatatcacgcacccggtagcggttgcctgcattctggccgagatgaaactcgactatgaaacgctgatggcggcgctgctgcatgacgtgattgaagatactcccgccacctaccaggatCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGACGCGCTGGTAGAGGGGGTGTCGAAACTTGATAAACTCAAGTTCcgcg                                                                                                                                                                          <  2:89325‑M1/82‑1 (MQ=255)
                                                                                     ccgtgatgctcacgaggggcaaacacgttcaagcggtgaaccctataccacgcacccggtagcggttgcctgcattctggccgagatgaaactcgactatgaaacgctgatggcggcgctgctgcatgacgtgattgaagatactcccgccacctaccaggatCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGGGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATaaga                                                                                                                                                                    >  1:374099‑M1/164‑251 (MQ=255)
                                                                                     cGTTGCACGTGATGCTCACGCGGGGCAAACACGTTCAAGCGGTGAACACTATCTCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGATGCTGGCGGCGCTGCTGCATGACGTGATTGAAGATACTCCCGCAACCTACCAGGATcaggatatggaacagctttttggtacaagcgtcgccgagctggtagagggggtgtcgaaacttgataaactcaagttccgcg                                                                                                                                                                    <  1:83071‑M2/251‑83 (MQ=255)
                                                                                     cGTTGCACGTGATGCTCACGAGGGGCAAACACGTTCAAGCGGTGAACCCTATATCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCGCTGCTGCATGACGTGATTGAAGATACTCCCGCCACCTACCAGGATcaggatatggaacagctttttggtaaaagcgtcgccgagctggtagagggggtgtcgaaacttgataaactcaagttccgcg                                                                                                                                                                    >  2:145344‑M2/1‑169 (MQ=255)
                                                                                      gTTGCACGTGATGCTCACGAGGGGCAAACACGCTCACGCGGTGAACCCTATATCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCGCTGCTGCATGACGTGATTGAAGATACTCCCGCCACCTACCAGGATcaggatatggaacagctttttggtaaaagcgtcgacgcgctggtagagggggtgtcgaaacttgataaactcaagttccgcg                                                                                                                                                                    <  2:89325‑M2/250‑83 (MQ=255)
                                                                                      cgggatgctcacgaggggcaaacacgttcaagcggtgaaccctatatcacgcacccggtagcggttgcctgcattctggccgagatgaaactcgactatgaaacgctgatggcggcgctgctgcatgacgtgattgaagatactcccgccacctaccaggatCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGGGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATaagaa                                                                                                                                                                   >  2:166297‑M1/163‑251 (MQ=255)
                                                                                          atgctcacgagggtcaaacacgttcaagcggtgagccctatatcacgcacccggtagcggttgcctgcattctggccgagatgaaactcgactatgaaacgctgatggcggcgctgctgcatgacgtgattgaagatactccctccacctaccaggatCAGGATATGGAACAGCTTTTTGGTAAAAGCTTCGCCGAGCTGGTAGAGGGGGTGTCGAAACTTGATAAACCTAAGTTCCGCGATAAGAAagag                                                                                                                                                               >  2:422596‑M1/159‑251 (MQ=255)
                                                                                          atgctcacgaggggcaaacacgttcaagcggtgaaccctatatcacgcacccggtagcggttgcctgcattctggccgagatgaaactcgactatgaaacgctgatggcggcgctgctgcatgacgtgattgaagatactcccgccacctaccaggatCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGGGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATAAGAAagag                                                                                                                                                               >  2:33857‑M1/159‑251 (MQ=255)
                                                                                           ccgTGATGCTCACGAGGGGCAAACACGTTCAAGCGGTGAACCCTATACCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCGCTGCTGCATGACGTGATTGAAGATACTCCCGCCACCTACCAGGATcaggatatggaacagctttttggtaaaagcgtcgccgagctggtagagggggtgtcgaaacttgataaactcaagttccgcgataaga                                                                                                                                                              >  1:374099‑M2/2‑163 (MQ=255)
                                                                                            cgggatgCTCACGAGGGGCAAACACGTTCAAGCGGTGAACCCTATATCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCGCTGCTGCATGACGTGATTGAAGATACTCCCGCCACCTACCAGGATcaggatatggaacagctttttggtaaaagcgtcgccgagctggtagagggggtgtcgaaacttgataaactcaagttccgcgataagaa                                                                                                                                                             >  2:166297‑M2/4‑162 (MQ=255)
                                                                                                aTGCTCACGAGGGTCAAACACGTTCAAGCGGTGAGCCCTATATCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCGCTGCTGCATGACGTGATTGAAGATACTCCCTCCACCTACCAGGATcaggatatggaacagctttttggtaaaagcttcgccgagctggtagagggggtgtcgaaacttgataaacctaagttccgcgataagaaagag                                                                                                                                                         >  2:422596‑M2/1‑158 (MQ=255)
                                                                                                aTGCTCACGAGGGGCAAACACGTTCAAGCGGTGAACCCTATATCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCGCTGCTGCATGACGTGATTGAAGATACTCCCGCCACCTACCAGGATcaggatatggaacagctttttggtaaaagcgtcgccgagctggtagagggggtgtcgaaacttgataaactcaagttccgcgataagaaagag                                                                                                                                                         >  2:33857‑M2/1‑158 (MQ=255)
                                                                                                   aggggcaaacacggtcaagcggtgaaccctatatcccgcacccggtagcggttgcctgcattctggccgagatgaaactcgactatgaaacgctgatggcggcgctgctgcatgacgtgattgaagatactcccgccacctaccaggatCAGGATATGGAACATCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGGGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATACGAAAGAGGCGCAGGcc                                                                                                                                                      >  2:213790‑M1/150‑251 (MQ=255)
                                                                                                         aGGGGCAAACACGGTCAAGCGGTGAACCCTATATCCCGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCGCTGCTGCATGACGTGATTGAAGATACTCCCGCCACCTACCAGGATcaggatatggaacatctttttggtaaaagcgtcgccgagctggtagagggggtgtcgaaacttgataaactcaagttccgcgatacgaaagaggcgcaggcc                                                                                                                                                >  2:213790‑M2/1‑149 (MQ=255)
                                                                                                              cgtgcaagcggtgaacccgatttcacgcacccggtagcggttgcctgcattctggccgagatgaaactcgactatgaaacgctgatggcggcgctgctgcatgacgtgattgaagatactcccgccacctaccaggatCAGGATATGGATCAGGTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGGGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATAAGAAAGAGGCGCAGGGCGAAAACTTTCg                                                                                                                                           >  2:11135‑M1/139‑251 (MQ=255)
                                                                                                               gttcaagcggtgaaccctatatcacgcacccggtagcggttgcctgcattctggccgagatgaaactcgactatgaaacgctgatggcggcgctgctgcatgacgtgattgaagatactcccgccacctaccaggatCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGGGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATAAGAAAGAGGCGCAGGCCGAAAACTTTCGc                                                                                                                                          >  2:136692‑M1/138‑251 (MQ=255)
                                                                                                               gttcaagcggtgaaccctagatcacgcacccggtagcggttgcctgcattctggccgagatgaaactcgactatgaaacgctgatggcggcgctgctgcatgacgtgattgaagatactcccgccacctaccaggatCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGGGGGTGTCGAAACTTGATAAACGCAAGTTCCGCGATAAGAAAGAGGCGCAGGCCGAAAACTTTCGc                                                                                                                                          >  2:193525‑M1/138‑251 (MQ=255)
                                                                                                                ttcaagcggtgaaccctatatcacgcacccggtagcggttgcctgcattctggccgagatgaaactcgactatgaaacgctgatggcggcgctgctgcatgacgtgattgaagatactcccgccacctaccaggatCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGGGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATAAGAAAGAGGCGCAGGCCGAAAACTTTCGc                                                                                                                                          <  1:72592‑M1/114‑1 (MQ=255)
                                                                                                                    cGTGCAAGCGGTGAACCCGATTTCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCGCTGCTGCATGACGTGATTGAAGATACTCCCGCCACCTACCAGGATcaggatatggatcaggtttttggtaaaagcgtcgccgagctggtagagggggtgtcgaaacttgataaactcaagttccgcgataagaaagaggcgcagggcgaaaactttcg                                                                                                                                     >  2:11135‑M2/1‑138 (MQ=255)
                                                                                                                     gcggtgaaccctatatcacgctcccggtagcggttgcctgcattctggccgagatgaaactcgactatgaaacgctgatggcggcgctgctgcctgacgtgattgaagatactcccgcctcctaccaggatCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGGGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATAAGAAATAGGCGCAGGCCGAAAACTTTCGCAAgatg                                                                                                                                    >  1:130195‑M1/132‑251 (MQ=255)
                                                                                                                     gTTCAAGCGGTGAACCCTATATCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCGCTGCTGCATGACGTGATTGAAGATACTCCCGCCACCTACCAGGATcaggatatggaacagctttttggtaaaagcgtcgccgagctggtagagggggtgtcgaaacttgataaactcaagttccgcgataagaaagaggcgcaggccgaaaactttcgc                                                                                                                                    >  2:136692‑M2/1‑137 (MQ=255)
                                                                                                                     gTTCAAGCGGTGAACCCTAGATCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCGCTGCTGCATGACGTGATTGAAGATACTCCCGCCACCTACCAGGATcaggatatggaacagctttttggtaaaagcgtcgccgagctggtagagggggtgtcgaaacttgataaacgcaagttccgcgataagaaagaggcgcaggccgaaaactttcgc                                                                                                                                    >  2:193525‑M2/1‑137 (MQ=255)
                                                                                                                      ttCAAGCGGTGAACCCTATATCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCGCTGCTGCATGACGTGATTGAAGATACTCCCGCCACCTACCAGGATcaggatatggaacagctttttggtaaaagcgtcgccgagctggtagagggggtgtcgaaacttgataaactcaagttccgcgataagaaagaggcgcaggccgaaaactttcgc                                                                                                                                    <  1:72592‑M2/250‑115 (MQ=255)
                                                                                                                      cggtgaaccctatatcacgcacacggtagcggttgcctgcattctggccgagatgaaactcgactatgaaacgctgatggcggcgctgctgcatgacgagattgaagatactcccgccacctaccaggatCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGGGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATAAGAAAGAGGCGCAGGCCGAAAACTTTCGCAAgatg                                                                                                                                    <  1:383110‑M1/120‑1 (MQ=255)
                                                                                                                          gaaccctgtatcacgcacccggtagcggttgccggcattctggccgagaggaaactcgactatgaaacgctgatggcggcgctgctgcatgacgtgattgaagatactcccgccacctaccaggatCAGGATATGGAACAGCTTGTTGGTAAAAGCGTCGCCGAGCTGGTAGAGGGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATAAGAAAGAGGCGCAGTCCGAAAACTTTCGCAAGATGATTAt                                                                                                                               >  2:34743‑M1/127‑251 (MQ=255)
                                                                                                                           gCGGTGAACCCTATATCACGCTCCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCGCTGCTGCCTGACGTGATTGAAGATACTCCCGCCTCCTACCAGGATcaggatatggaacagctttttggtaaaagcgtcgccgagctggtagagggggtgtcgaaacttgataaactcaagttccgcgataagaaataggcgcaggccgaaaactttcgcaagatg                                                                                                                              >  1:130195‑M2/1‑131 (MQ=255)
                                                                                                                            cGGTGAACCCTATATCACGCACACGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCGCTGCTGCATGACGAGATTGAAGATACTCCCGCCACCTACCAGGATcaggatatggaacagctttttggtaaaagcgtcgccgagctggtagagggggtgtcgaaacttgataaactcaagttccgcgataagaaagaggcgcaggccgaaaactttcgcaagatg                                                                                                                              <  1:383110‑M2/250‑121 (MQ=255)
                                                                                                                            accctatatcacgcacccggtagcggttgcctgcattctggccgagatgaaactcgactatgaaacgctgatggcggcgctgctgcatgacgtgattgaagatactcccgccacctaccaggatCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGGGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATAAGAAAGAGGCGCAGGCCGAAAACTTTCGCAAGATGATTATgg                                                                                                                             >  1:182803‑M1/125‑251 (MQ=255)
                                                                                                                              cctatatcacgcacccggtagcggttgcctgcattctggccgagatgaaactcgactatgaaacgctgatggcggcgctgctgcatgacgtgattgaagatactcccgccacctaccaggatCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGGGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATAAGAAAGAGGCGCAGGCCGAACACTTTCGCAAGATGATTATGGCg                                                                                                                           <  1:249577‑M1/129‑1 (MQ=255)
                                                                                                                               ctatatcacgcacccggtagcggttgcctgcattctggccgagatgaaactcgactatgaaacgctgatggcggcgctgctgcatgacgtgattgaagatactcccgccacctaccaggatCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGGGGGTGTCGAAACtt                                                                                                                                                                                             <  1:158144‑M1/63‑1 (MQ=255)
                                                                                                                               ctatatcacgcacccggtagcggttgcctgcattctggccgagatgaaactcgactatgaaacgctgatggcggcgctgctgcatgacgtgattgaagatactcccgccacctaccaggatCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGGGGGTGTCGAAACtt                                                                                                                                                                                             >  2:158144‑M1/122‑184 (MQ=255)
                                                                                                                                gAACCCTGTATCACGCACCCGGTAGCGGTTGCCGGCATTCTGGCCGAGAGGAAACTCGACTATGAAACGCTGATGGCGGCGCTGCTGCATGACGTGATTGAAGATACTCCCGCCACCTACCAGGATcaggatatggaacagcttgttggtaaaagcgtcgccgagctggtagagggggtgtcgaaacttgataaactcaagttccgcgataagaaagaggcgcagtccgaaaactttcgcaagatgattat                                                                                                                         >  2:34743‑M2/1‑126 (MQ=255)
                                                                                                                                  aCCCTATATCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCGCTGCTGCATGACGTGATTGAAGATACTCCCGCCACCTACCAGGATcaggatatggaacagctttttggtaaaagcgtcgccgagctggtagagggggtgtcgaaacttgataaactcaagttccgcgataagaaagaggcgcaggccgaaaactttcgcaagatgattatgg                                                                                                                       >  1:182803‑M2/1‑124 (MQ=255)
                                                                                                                                   atcacgcacccggtagcggttgcctgcattctggccgagatgaaactcgactatgaaacgctgatggcggcgctgctgcatgacgtgattgaagatactcccgccacctaccaggatCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGGGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATAAGAAAGAGGCGCAGGCCGAAAACTTTCGCAAGATGATTATGGCGAtggg                                                                                                                      <  1:325967‑M1/134‑2 (MQ=255)
                                                                                                                                   atcacgcacccggtagcggttgcctgcattctggccgagatgaaactcgactatgaaacgctgatggcggcgctgctgcatgacgtgattgaagatactcccgccacctaccaggatCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGGGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATAAGAAAGAGGCGCAGGCCGAAAACTTTCGCAAGATGAGTATGGCGATgg                                                                                                                       >  2:30785‑M1/118‑250 (MQ=255)
                                                                                                                                    ccTATATCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCGCTGCTGCATGACGTGATTGAAGATACTCCCGCCACCTACCAGGATcaggatatggaacagctttttggtaaaagcgtcgccgagctggtagagggggtgtcgaaacttgataaactcaagttccgcgataagaaagaggcgcaggccgaacactttcgcaagatgattatggcg                                                                                                                     <  1:249577‑M2/251‑130 (MQ=255)
                                                                                                                                     cTATATCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCGCTGCTGCATGACGTGATTGAAGATACTCCCGCCACCTACCAGGATcaggatatggaacagctttttggtaaaagcgtcgccgagctggtagagggggtgtcgaaactt                                                                                                                                                                                       >  2:158144‑M2/1‑121 (MQ=255)
                                                                                                                                     cTATATCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCGCTGCTGCATGACGTGATTGAAGATACTCCCGCCACCTACCAGGATcaggatatggaacagctttttggtaaaagcgtcgccgagctggtagagggggtgtcgaaactt                                                                                                                                                                                       <  1:158144‑M2/184‑64 (MQ=255)
                                                                                                                                         atCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCGCTGCTGCATGACGTGATTGAAGATACTCCCGCCACCTACCAGGATcaggatatggaacagctttttggtaaaagcgtcgccgagctggtagagggggtgtcgaaacttgataaactcaagttccgcgataagaaagaggcgcaggccgaaaactttcgcaagatgattatggcgatggg                                                                                                                <  1:325967‑M2/251‑135 (MQ=255)
                                                                                                                                         atCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCGCTGCTGCATGACGTGATTGAAGATACTCCCGCCACCTACCAGGATcaggatatggaacagctttttggtaaaagcgtcgccgagctggtagagggggtgtcgaaacttgataaactcaagttccgcgataagaaagaggcgcaggccgaaaactttcgcaagatgagtatggcgatgg                                                                                                                 >  2:30785‑M2/1‑117 (MQ=255)
                                                                                                                                           cccggtagcggttgcctgcattctggccgagatgaaactcgactatgaaacgctgatggcggcgctgctgcatgacgtgattgaagatactcccgccacctaccaggatCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGGGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATAAGAAAGAGGCGCAGGCCGAAAACTTTCGCAAGATGATTATGGCGATGGTGCAGGata                                                                                                              >  1:296017‑M1/110‑251 (MQ=255)
                                                                                                                                              ggtagcggttgcctgcattctggccgagatgaaactcgactatgaaacgctgatggcggcgctgctgcatgacgtgattgaagatactcccgccacctaccaggatCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGGGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATAAGAAAGAGGCGCAGGCAGAAAACTTTCGCAAGATGATTATGGCGATGGTGCAGGATATcc                                                                                                           <  2:72433‑M1/145‑1 (MQ=255)
                                                                                                                                                 cccGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCGCTGCTGCATGACGTGATTGAAGATACTCCCGCCACCTACCAGGATcaggatatggaacagctttttggtaaaagcgtcgccgagctggtagagggggtgtcgaaacttgataaactcaagttccgcgataagaaagaggcgcaggccgaaaactttcgcaagatgattatggcgatggtgcaggata                                                                                                        >  1:296017‑M2/1‑109 (MQ=255)
                                                                                                                                                    ggTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCGCTGCTGCATGACGTGATTGAAGATACTCCCGCCACCTACCAGGATcaggatatggaacagctttttggtaaaagcgtcgccgagctggtagagggggtgtcgaaacttgataaactcaagttccgcgataagaaagaggcgcaggcagaaaactttcgcaagatgattatggcgatggtgcaggatatcc                                                                                                     <  2:72433‑M2/251‑146 (MQ=255)
                                                                                                                                                        gcctgcattctggccgagatgaaactcgactatgaaacgctgatggcggcgctgctgcatgacgtgattgaagatactcccgccacctaccaggatCAGGATATGGAACAGCTTTGTGGTAAAAGCGTCGCCGAGCTGGTAGAGGGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATAAGAAAGAGGCGCAGGCCGAAAACTTGCGCAAGATGATTAGGGCGATGGTGCAGGATATCCGCGTCATcc                                                                                                  >  2:396358‑M1/97‑250 (MQ=255)
                                                                                                                                                          ctgcattctggccgagatgaaactcgactatgaaacgctgatggcggcgctgctgcatgacgtgattgaagatactcccgccacctaccaggatCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGGGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATAAGAAAGAGGCGCAGGCCGAAAACTTTCGCAAGATGATTATGGCGATGGTGCAGGATATCCGCGTCATCCtc                                                                                                <  2:122153‑M1/156‑1 (MQ=255)
                                                                                                                                                              gCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCGCTGCTGCATGACGTGATTGAAGATACTCCCGCCACCTACCAGGATcaggatatggaacagctttgtggtaaaagcgtcgccgagctggtagagggggtgtcgaaacttgataaactcaagttccgcgataagaaagaggcgcaggccgaaaacttgcgcaagatgattagggcgatggtgcaggatatccgcgtcatcc                                                                                            >  2:396358‑M2/1‑96 (MQ=255)
                                                                                                                                                                cTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCGCTGCTGCATGACGTGATTGAAGATACTCCCGCCACCTACCAGGATcaggatatggaacagctttttggtaaaagcgtcgccgagctggtagagggggtgtcgaaacttgataaactcaagttccgcgataagaaagaggcgcaggccgaaaactttcgcaagatgattatggcgatggtgcaggatatccgcgtcatcctc                                                                                          <  2:122153‑M2/250‑157 (MQ=255)
                                                                                                                                                                                        tgaaacgctgatggcggcgctgctgcatgacgtgattgaagatactcccgccacctaccaggatCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGGGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATAAGAAAGAGGCGCAGGCCGAAAACTTTCGCAAGCTGATTATGGCGATGGTGCAGGATATCCGCGTCATCCTCATCAAACTTGCCGACCGAACCCACAACATgc                                                                 <  1:299833‑M1/187‑1 (MQ=255)
                                                                                                                                                                                           aacgctgatgtcggcgctgctgcatgacgtgattgaagatactcccgccacctaccaggatCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGGGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATAAGAAAGAGGCGCAGGCCGAAAACTTTCGCAAGATGATTATGGCGATGGTGCAGGATATCCGCGTCATCCTCATCAAACTTGCCGACCGTACCCACAACATgcgc                                                               <  1:243517‑M1/189‑1 (MQ=255)
                                                                                                                                                                                              tGAAACGCTGATGGCGGCGCTGCTGCATGACGTGATTGAAGATACTCCCGCCACCTACCAGGATcaggatatggaacagctttttggtaaaagcgtcgccgagctggtagagggggtgtcgaaacttgataaactcaagttccgcgataagaaagaggcgcaggccgaaaactttcgcaagctgattatggcgatggtgcaggatatccgcgtcatcctcatcaaacttgccgaccgaacccacaacatgc                                                           <  1:299833‑M2/251‑188 (MQ=255)
                                                                                                                                                                                                 aaCGCTGATGTCGGCGCTGCTGCATGACGTGATTGAAGATACTCCCGCCACCTACCAGGATcaggatatggaacagctttttggtaaaagcgtcgccgagctggtagagggggtgtcgaaacttgataaactcaagttccgcgataagaaagaggcgcaggccgaaaactttcgcaagatgattatggcgatggtgcaggatatccgcgtcatcctcatcaaacttgccgaccgtacccacaacatgcgc                                                         <  1:243517‑M2/250‑190 (MQ=255)
                                                                                                                                                                                                       ggcgctgctgcatgacgtgattgaagatactcccgccacctaccaggatCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGGGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATAAGAAAGAGGCGCAGGCCGAAAACTTTCGCAAGATGATTATGGCGATGGTGCAGGATATCCGCGTCATCCTCATCAAACTTGCCGACCGTACCCACAACATGCGCACGCTGGGCTCAc                                                  <  1:364634‑M1/202‑1 (MQ=255)
                                                                                                                                                                                                             ggcgCTGCTGCATGACGTGATTGAAGATACTCCCGCCACCTACCAGGATcaggatatggaacagctttttggtaaaagcgtcgccgagctggtagagggggtgtcgaaacttgataaactcaagttccgcgataagaaagaggcgcaggccgaaaactttcgcaagatgattatggcgatggtgcaggatatccgcgtcatcctcatcaaacttgccgaccgtacccacaacatgcgcacgctgggctcac                                            <  1:364634‑M2/251‑203 (MQ=255)
                                                                                                                                                                                                                gcatgacgtgattgaagatactcccgccacctaccaggatCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGGGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATAAGAAAGAGGCGCAGGCCGAAAACTTTCGCAAGATGATTATGGCGATGGTGCAGGATATCCGCGTCATCCTCATCAAACTTGCCGACCGTACCCACAACATGCGCACGCTGGGCTCACTTCGCCCgg                                         >  1:122153‑M1/41‑251 (MQ=255)
                                                                                                                                                                                                                gcatgaagtgattgaagatcctcccgccccctaccaggatCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGGGGGTGACGAAACTTGATAAACTCAAGTTCCGCGATAAGAAAGAGGCGCAGGCCGAAAACTTTCGCAAGAGGATTATGGCGATGGTGCATGATATCCGCGTCATCCTCATCAAACGTGCCGACCGTACCCTCAACATGCGCACGCTTGGCTCACTtcgcctgg                                         >  2:253616‑M1/41‑248 (MQ=255)
                                                                                                                                                                                                                 catgacgtgattgaagatactcccgccacctaccaggatCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGGGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATAAGAAAGAGGCGCAGGCCGAAAACTTTCGCAAGATGATTATGGCGATGGTGCAGGATATCCGCGTCATCCTCATCAAACTTGCCGACCGTACCCACAACATGCGCACGCTGGGCTCACTTCGCCCgg                                         <  1:290393‑M1/211‑1 (MQ=255)
                                                                                                                                                                                                                      gcATGACGTGATTGAAGATACTCCCGCCACCTACCAGGATcaggatatggaacagctttttggtaaaagcgtcgccgagctggtagagggggtgtcgaaacttgataaactcaagttccgcgataagaaagaggcgcaggccgaaaactttcgcaagatgattatggcgatggtgcaggatatccgcgtcatcctcatcaaacttgccgaccgtacccacaacatgcgcacgctgggctcacttcgcccgg                                   >  1:122153‑M2/1‑40 (MQ=255)
                                                                                                                                                                                                                      gcATGAAGTGATTGAAGATCCTCCCGCCCCCTACCAGGATcaggatatggaacagctttttggtaaaagcgtcgccgagctggtagagggggtgacgaaacttgataaactcaagttccgcgataagaaagaggcgcaggccgaaaactttcgcaagaggattatggcgatggtgcatgatatccgcgtcatcctcatcaaacgtgccgaccgtaccctcaacatgcgcacgcttggctcacttcgcctgg                                   >  2:253616‑M2/1‑40 (MQ=255)
                                                                                                                                                                                                                       cATGACGTGATTGAAGATACTCCCGCCACCTACCAGGATcaggatatggaacagctttttggtaaaagcgtcgccgagctggtagagggggtgtcgaaacttgataaactcaagttccgcgataagaaagaggcgcaggccgaaaactttcgcaagatgattatggcgatggtgcaggatatccgcgtcatcctcatcaaacttgccgaccgtacccacaacatgcgcacgctgggctcacttcgcccgg                                   <  1:290393‑M2/250‑212 (MQ=255)
                                                                                                                                                                                                                             gaagatactcccgccacctaccaggatCAGTATATGGAACAGCGTTTTGGTAAAAGCGTCGCCGAGCGGGTAGAGGGGGTGTCGAAACTTGATAAACGCAAGTTCCGCGATAAGAAAGAGGCGCAGGCCGAAAACTTTCGCAAGATGATTATGGCGATGGTGCAGGATATCCGCGTCATCCTCATCAAACGTGCCGACCTTACCCACAACATGCGCACGCTGGGCTCACTTCGCCCGGACAAACGTcgccg                            >  2:254256‑M1/28‑251 (MQ=255)
                                                                                                                                                                                                                             gaagatactcccgccacctaccaggatCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGGGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATAAGAAAGAGGCGCAGGCCGAAAACTTTCGCAAGATGATTATGGCGATGGTGCAGGATATCCGCGTCATCCTCATCAAACTTGCCGACCGTACCCACAACATGCGCACGCTGGGCTCACTTCGCCCGGACAAACGTcgccg                            <  1:313010‑M1/224‑1 (MQ=255)
                                                                                                                                                                                                                              aagatactcccgccacctaccaggatCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGGGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATAAGAAAGAGGCGCAGGCCGAAAACTTTCGCAAGATGATTATGGCGATGGTGCAGGATATCCGCGTCATCCTCATCAAACTTGCCGACCGTACCCACAACATGCGCACGCTGGGCTCACTTCGCCCGGACAAACGTcgccgc                           <  1:10004‑M1/225‑1 (MQ=255)
                                                                                                                                                                                                                                   gAAGATACTCCCGCCACCTACCAGGATcagtatatggaacagcgttttggtaaaagcgtcgccgagcgggtagagggggtgtcgaaacttgataaacgcaagttccgcgataagaaagaggcgcaggccgaaaactttcgcaagatgattatggcgatggtgcaggatatccgcgtcatcctcatcaaacgtgccgaccttacccacaacatgcgcacgctgggctcacttcgcccggacaaacgtcgccg                      >  2:254256‑M2/1‑27 (MQ=255)
                                                                                                                                                                                                                                   gAAGATACTCCCGCCACCTACCAGGATcaggatatggaacagctttttggtaaaagcgtcgccgagctggtagagggggtgtcgaaacttgataaactcaagttccgcgataagaaagaggcgcaggccgaaaactttcgcaagatgattatggcgatggtgcaggatatccgcgtcatcctcatcaaacttgccgaccgtacccacaacatgcgcacgctgggctcacttcgcccggacaaacgtcgccg                      <  1:313010‑M2/251‑225 (MQ=255)
                                                                                                                                                                                                                                    aaGATACTCCCGCCACCTACCAGGATcaggatatggaacagctttttggtaaaagcgtcgccgagctggtagagggggtgtcgaaacttgataaactcaagttccgcgataagaaagaggcgcaggccgaaaactttcgcaagatgattatggcgatggtgcaggatatccgcgtcatcctcatcaaacttgccgaccgtacccacaacatgcgcacgctgggctcacttcgcccggacaaacgtcgccgc                     <  1:10004‑M2/251‑226 (MQ=255)
                                                                                                                                                                                                                                      cccgccacctaccaggatCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGGGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATAAGAAAGAGGCGCAGGCCGAAAACTTTCGCAAGATGATTATGGCGATGGTGCAGGATATCCGCGTCATCCTCATCAAACTTGCCGACCGTACCCACAACATGCGCACGCTGGGCTCACTTCGCCCGGACAAACGTCGCCGCATCGCCCg                   <  2:75675‑M1/233‑1 (MQ=255)
                                                                                                                                                                                                                                            cccGCCACCTACCAGGATcaggatatggaacagctttttggtaaaagcgtcgccgagctggtagagggggtgtcgaaacttgataaactcaagttccgcgataagaaagaggcgcaggccgaaaactttcgcaagatgattatggcgatggtgcaggatatccgcgtcatcctcatcaaacttgccgaccgtacccacaacatgcgcacgctgggctcacttcgcccggacaaacgtcgccgcatcgcccg             <  2:75675‑M2/251‑234 (MQ=255)
                                                                                                                                                                                                                                                       gcAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGGGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATAAGAAAGAGGCGCAGGCCGAAAACTTTCGCAAGATGATTATGGCGATGGTGCAGGATATCCGCGTCATCCTCATCAAACTTGCCGACCGTACCCACAACATGCGCACGCTGGGCTCACTTCGCCAGGACAAACGTCGCCGCATCGCCCGTGAAACTCTCGAAAttt  <  2:10004/250‑1 (MQ=255)
                                                                                                                                                                                                                                                        cAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGGGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATAAGAAAGAGGCGCAGGCCGAAAACTTTCGCAAGATGATTATGGCGATGGTGCAGGATATCCGCGTCATCCTCATCAAACTTGCCGACCGTACCCACAACATGCGCACGCTGGGCTCACTTCGCCCGGACAAACGTCGCCGCATCGCCCGTGAAACTCTCGAAAttt  <  1:9654/250‑1 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                         
CCTTGTATCTGTTTGAAAGCCTGAATCAACTGATTCAAACCTACCTGCCGGAAGACCAAATCAAGCGTCTGCGGCAGGCGTATCTCGTTGCACGTGATGCTCACGAGGGGCAAACACGTTCAAGCGGTGAACCCTATATCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCGCTGCTGCATGACGTGATTGAAGATACTCCCGCCACCTACCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGGGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATAAGAAAGAGGCGCAGGCCGAAAACTTTCGCAAGATGATTATGGCGATGGTGCAGGATATCCGCGTCATCCTCATCAAACTTGCCGACCGTACCCACAACATGCGCACGCTGGGCTCACTTCGCCCGGACAAACGTCGCCGCATCGCCCGTGAAACTCTCGAAATTT  >  NC_000913/3822398‑3822895

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