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)44 (1.260) 41/490 0.1 96.7% coding (247/2109 nt) spoT bifunctional (p)ppGpp synthase/hydrolase SpoT
?NC_000913 = 3822651 3 (0.090)coding (252/2109 nt) spoT bifunctional (p)ppGpp synthase/hydrolase SpoT

TGTTTGAAAGCCTGAATCAACTGATTCAAACCTACCTGCCGGAAGACCAAATCAAGCGTCTGCGGCAGGCGTATCTCGTTGCACGTGATGCTCACGAGGGGCAAACACGTTCAAGCGGTGAACCCTATATCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCGCTGCTGCATGACGTGATTGAAGATACTCCCGCCACCTACCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGGGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATAAGAAAGAGGCGCAGGCCGAAAACTTTCGCAAGATGATTATGGCGATGGTGCAGGATATCCGCGTCATCCTCATCAAACTTGCCGACCGTACCCACAACATGCGCACGCTGGGCTCACTTCGCCCGGACAAACGTCGCCGCATCGCCCGTGAAACTCTCGAA  >  NC_000913/3822407‑3822891
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                     
aaagcctgaatcaactgtttcaaacctacctgccggaagaccaaatcaagcgtctgcggcaggcgtatctcgttgcacgtgatgctcacgaggggcaaacacgttcaagcggtgaaccctatatcacgcacccggtagcggttgcctgcattctggccgagatgaaactcgactatgaaacgctgatggcggcgctgctgcatgacgtgattgaagatactcccgccacctaccaggatCAGGATATGGa                                                                                                                                                                                                                                             <  1:299642‑M1/11‑1 (MQ=255)
      aaaGCCTGAATCAACTGTTTCAAACCTACCTGCCGGAAGACCAAATCAAGCGTCTGCGGCAGGCGTATCTCGTTGCACGTGATGCTCACGAGGGGCAAACACGTTCAAGCGGTGAACCCTATATCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCGCTGCTGCATGACGTGATTGAAGATACTCCCGCCACCTACCAGGATcaggatatgga                                                                                                                                                                                                                                       <  1:299642‑M2/250‑12 (MQ=255)
       gaatcaactgattcaaacctacctgccggaagaccaaatcaagcgtctgcggcaggcgtatctcgttgcacgtgatgctcacgaggggcaaacacgttcaagcggtgaaccctatatcacgcacccggtagcggttgcctgcattctggccgagatgaaactcgactatgaaacgctgatggcggcgctgctgcatgacgtgattgaagatactcccgccacctaccaggatCAGGATATGGAACAGCtt                                                                                                                                                                                                                                      >  1:190150‑M1/233‑250 (MQ=255)
       gaatcaactgattcaaacctacctgccggaagaccaaatcaagcgtctgcggcaggcgtatctcgttgcacgtgatgctcacgaggggcaaacacgttcaagcggtgaaccctatatcacgcacccggtagcggttgcctgcattctggccgagatgaaactcgactatgaaacgctgatggcggcgctgctgcaggacgtgattgaagatactcccgccacctaccaggatCAGGATATGGAACAGCttt                                                                                                                                                                                                                                     >  2:231496‑M1/233‑251 (MQ=255)
       gaatcaactgattcaaaccaacctgccggaagaccaaatcaagcgtctgcggcaggcgtatctcgttgcacgtgatgctcacgaggggcaaacacgttcaagcggtgaaccctatatcacgcacccggtagcggttgcctgcattctggccgagatgaaactctactatgaaacgctgatggcggcgctgctgcatgacgtgattgaagatactcccgccacctaccagtatCAGGATATGGAACAGCttt                                                                                                                                                                                                                                     >  2:231498‑M1/233‑251 (MQ=255)
           caactgattcaaccctacctgccggaagaccaaatcaagcgtctgcggcaggcgtatctcgttgcacgtgatgctcacgaggggcaaacacgttcaagcggtgaaccctatatcacgcacccggtagcggttgcctgcattctggccgagatgaaactcgactatgaaacgctgatggcggcgctgctgcatgacgtgattgaagatactcccgccacctaccaggatCAGGATATGGAACAGCTTTTTgg                                                                                                                                                                                                                                 >  2:631‑M1/229‑251 (MQ=255)
           caactgattcaaacctacctgccggaagaccaaatcaagcgtctgcggcaggcgtatctcgttgcacgtgatgctcacgaggggcaaacacgttcaagcggtgaaccctatatcacgcacccggtagcggttgcctgcattctggccgagatgaaactcgactatgaaacgctgatggcggcgctgctgcatgacgtgattgaagatactcccgccacctaccaggatCAGGATATGGAACAGCTTTTTgg                                                                                                                                                                                                                                 >  2:140527‑M1/229‑251 (MQ=255)
             gAATCAACTGATTCAAACCTACCTGCCGGAAGACCAAATCAAGCGTCTGCGGCAGGCGTATCTCGTTGCACGTGATGCTCACGAGGGGCAAACACGTTCAAGCGGTGAACCCTATATCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCGCTGCTGCATGACGTGATTGAAGATACTCCCGCCACCTACCAGGATcaggatatggaacagctt                                                                                                                                                                                                                                >  1:190150‑M2/1‑232 (MQ=255)
             gAATCAACTGATTCAAACCTACCTGCCGGAAGACCAAATCAAGCGTCTGCGGCAGGCGTATCTCGTTGCACGTGATGCTCACGAGGGGCAAACACGTTCAAGCGGTGAACCCTATATCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCGCTGCTGCAGGACGTGATTGAAGATACTCCCGCCACCTACCAGGATcaggatatggaacagcttt                                                                                                                                                                                                                               >  2:231496‑M2/1‑232 (MQ=255)
             gAATCAACTGATTCAAACCAACCTGCCGGAAGACCAAATCAAGCGTCTGCGGCAGGCGTATCTCGTTGCACGTGATGCTCACGAGGGGCAAACACGTTCAAGCGGTGAACCCTATATCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCTACTATGAAACGCTGATGGCGGCGCTGCTGCATGACGTGATTGAAGATACTCCCGCCACCTACCAGTATcaggatatggaacagcttt                                                                                                                                                                                                                               >  2:231498‑M2/1‑232 (MQ=255)
             actgattcaaacctacctgccggaagaccaaatcaagcgtctgcggcaggcgtatctcgttgcacgtgatgctcacgaggggcaaacacgttccagcggtgaaccctatatcacgcacccggtagcggttgcctgcattctggccgagatgaaactcgactatgaaacgctgatggcggcgctgctgcatgacgtgattgacgatactcccgccacctaccaggatCAGGATATGGAACAGCTTTTTGGt                                                                                                                                                                                                                                <  2:326627‑M1/24‑1 (MQ=255)
                 cAACTGATTCAACCCTACCTGCCGGAAGACCAAATCAAGCGTCTGCGGCAGGCGTATCTCGTTGCACGTGATGCTCACGAGGGGCAAACACGTTCAAGCGGTGAACCCTATATCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCGCTGCTGCATGACGTGATTGAAGATACTCCCGCCACCTACCAGGATcaggatatggaacagctttttgg                                                                                                                                                                                                                           >  2:631‑M2/1‑228 (MQ=255)
                 cAACTGATTCAAACCTACCTGCCGGAAGACCAAATCAAGCGTCTGCGGCAGGCGTATCTCGTTGCACGTGATGCTCACGAGGGGCAAACACGTTCAAGCGGTGAACCCTATATCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCGCTGCTGCATGACGTGATTGAAGATACTCCCGCCACCTACCAGGATcaggatatggaacagctttttgg                                                                                                                                                                                                                           >  2:140527‑M2/1‑228 (MQ=255)
                   aCTGATTCAAACCTACCTGCCGGAAGACCAAATCAAGCGTCTGCGGCAGGCGTATCTCGTTGCACGTGATGCTCACGAGGGGCAAACACGTTCCAGCGGTGAACCCTATATCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCGCTGCTGCATGACGTGATTGACGATACTCCCGCCACCTACCAGGATcaggatatggaacagctttttggt                                                                                                                                                                                                                          <  2:326627‑M2/250‑25 (MQ=255)
                    caaacctacctgccggaagaccaaatcaagcgtctgcggcaggcgtatctcgttgcacgtgatgctcacgaggggcaaacacgttcaagtggtgaaccctatatcacgcacccggtagcggttgcctgcattctggccgagatgaaactcgactatgaaacgctgatggcggcgctgctgcatgacgggattgaagatactcccgccacctaccaggatCAGGATATGGAACAGCTTTTTGGTAAAAGCGt                                                                                                                                                                                                                        >  2:62322‑M1/220‑251 (MQ=255)
                     aaacctacctgccggaagaccaaatcaagcgtctgcggcaggcgtatctcgttgcacgtgatgctcacgaggggcaaacacgttcaagcggtgaaccctatatcacgcacccggtagcggttgcctgcattctggccgagatgaaactcgactatgaaacgctgatggcggcgctgctgcatgacgtgattgaagatactcccgccacctaccaggatCAGGATATGGAACAGCTTTTTGGTAAAAGCGTc                                                                                                                                                                                                                       >  1:137951‑M1/219‑251 (MQ=255)
                          cAAACCTACCTGCCGGAAGACCAAATCAAGCGTCTGCGGCAGGCGTATCTCGTTGCACGTGATGCTCACGAGGGGCAAACACGTTCAAGTGGTGAACCCTATATCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCGCTGCTGCATGACGGGATTGAAGATACTCCCGCCACCTACCAGGATcaggatatggaacagctttttggtaaaagcgt                                                                                                                                                                                                                  >  2:62322‑M2/1‑219 (MQ=255)
                           aaaCCTACCTGCCGGAAGACCAAATCAAGCGTCTGCGGCAGGCGTATCTCGTTGCACGTGATGCTCACGAGGGGCAAACACGTTCAAGCGGTGAACCCTATATCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCGCTGCTGCATGACGTGATTGAAGATACTCCCGCCACCTACCAGGATcaggatatggaacagctttttggtaaaagcgtc                                                                                                                                                                                                                 >  1:137951‑M2/1‑218 (MQ=255)
                             ctgccggaagaccaaatcaagcgtctgcggcaggcgtatctcgttgcacgtgatgctcacgaggggcaaacacgttcaagcggtgaaccctatatcacgcacccggtagcggttgcctgcattctggccgagatgaaactcgactatgaaacgctgatggcggcgctgctgcatgacgtgattgaagatactcccgccacctaccaggatCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCt                                                                                                                                                                                                               <  1:23872‑M1/41‑1 (MQ=255)
                                   ctGCCGGAAGACCAAATCAAGCGTCTGCGGCAGGCGTATCTCGTTGCACGTGATGCTCACGAGGGGCAAACACGTTCAAGCGGTGAACCCTATATCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCGCTGCTGCATGACGTGATTGAAGATACTCCCGCCACCTACCAGGATcaggatatggaacagctttttggtaaaagcgtcgccgagct                                                                                                                                                                                                         <  1:23872‑M2/251‑42 (MQ=255)
                                    aagaccaaatcaagcgtctgcggcaggcgtatctcgttgcacgtgatgctcacgaggggcaaacacgttcaagcggtgaaccctatatcacgcacccggtagcggttgcctgcattctggccgagatgaaactcgactatgaaacgctgatggcggcgctgctgcatgacgtgattgaagatactcccgccacctaccaggatCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTagag                                                                                                                                                                                                        >  2:259708‑M1/204‑251 (MQ=255)
                                     agaccaaatcaagcgtctgcggcaggcgtatctcgttgcacgtgatgctcacgaggggcaaacacgttcaagcggtgaaccctatatcacgcacccggtagcggttgcctgcattctggccgagatgaaactcgactatgaaacgctgatggcggcgctgctgcatgacgtgattgaagatactcccgccacctaccaggatTAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAgg                                                                                                                                                                                                       <  1:114918‑M1/49‑1 (MQ=255)
                                      gaccaaatcaagcgtctgcggcaggcgtatctcgttgcacgtgatgctcacgaggggcaaacacgttcaagcggtgaaccctatatcacgcacccggtagcggttgcctgcattctggccgagatgaaactcgactatgaaacgctgatggcggcgctgctgcatgacgtgattgaagatactcccgccacctaccaggatCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAggg                                                                                                                                                                                                      <  1:245949‑M1/50‑1 (MQ=255)
                                      gaccaaatcaagcgtctgcggcaggcgtatctcgttgcacgtgatgctcacgaggggcaaacacgttcaagcggtgaaccctatatcacgcacccggtagcggttgcctgcattctggccgagatgaaactcgactatgaaacgctgatggcggcgctgctgcatgacgtgattgaagatactcccgccacctaccaggatCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAggg                                                                                                                                                                                                      >  1:55775‑M1/202‑251 (MQ=255)
                                          aaatcaagcgtctgcggcaggcgtatctcgttgcacgtgatgctcacgaggggcaaacacgttcaagcggtgaaccctatatcacgcacccggtagcggttgcctgcattctggccgagatgaaactcgactatgaaacgctgatggcggcgctgctgcatgacgtgattgaagatactcccgccacctaccaggatCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGGGGgtg                                                                                                                                                                                                  >  2:184192‑M1/198‑251 (MQ=255)
                                          aaGACCAAATCAAGCGTCTGCGGCAGGCGTATCTCGTTGCACGTGATGCTCACGAGGGGCAAACACGTTCAAGCGGTGAACCCTATATCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCGCTGCTGCATGACGTGATTGAAGATACTCCCGCCACCTACCAGGATcaggatatggaacagctttttggtaaaagcgtcgccgagctggtagag                                                                                                                                                                                                  >  2:259708‑M2/1‑203 (MQ=255)
                                           aGACCAAATCAAGCGTCTGCGGCAGGCGTATCTCGTTGCACGTGATGCTCACGAGGGGCAAACACGTTCAAGCGGTGAACCCTATATCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCGCTGCTGCATGACGTGATTGAAGATACTCCCGCCACCTACCAGGATtaggatatggaacagctttttggtaaaagcgtcgccgagctggtagagg                                                                                                                                                                                                 <  1:114918‑M2/251‑50 (MQ=255)
                                            gACCAAATCAAGCGTCTGCGGCAGGCGTATCTCGTTGCACGTGATGCTCACGAGGGGCAAACACGTTCAAGCGGTGAACCCTATATCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCGCTGCTGCATGACGTGATTGAAGATACTCCCGCCACCTACCAGGATcaggatatggaacagctttttggtaaaagcgtcgccgagctggtagaggg                                                                                                                                                                                                <  1:245949‑M2/251‑51 (MQ=255)
                                            gACCAAATCAAGCGTCTGCGGCAGGCGTATCTCGTTGCACGTGATGCTCACGAGGGGCAAACACGTTCAAGCGGTGAACCCTATATCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCGCTGCTGCATGACGTGATTGAAGATACTCCCGCCACCTACCAGGATcaggatatggaacagctttttggtaaaagcgtcgccgagctggtagaggg                                                                                                                                                                                                >  1:55775‑M2/1‑201 (MQ=255)
                                               aagcgtctgcggcaggcgtatctcgttgcacgtgatgctcacgaggggcaaacacgttcaagcggtgaaccctatatcacgcacccggtagcggttgcctgcattctggccgagatgaaactcgactatgaaacgctgatggcggcgctgctgcatgacgtgattgaagatactcccgccacctaccaggatCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGGGGGTGTCGa                                                                                                                                                                                              <  1:70129‑M1/58‑1 (MQ=255)
                                                aaaTCAAGCGTCTGCGGCAGGCGTATCTCGTTGCACGTGATGCTCACGAGGGGCAAACACGTTCAAGCGGTGAACCCTATATCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCGCTGCTGCATGACGTGATTGAAGATACTCCCGCCACCTACCAGGATcaggatatggaacagctttttggtaaaagcgtcgccgagctggtagagggggtg                                                                                                                                                                                            >  2:184192‑M2/1‑197 (MQ=255)
                                                  cgtctgcggcaggcgtatctcgttgcacgtgatgctcacgaggggcaaacacgttcaagcggtgaaccctatatcacgcacccggtagcggttgcctgcattctggccgagatgaaactcgactatgaaacgctgatggcggcgctgctgcatgacgtgattgaagatactcccgccacctaccaggatCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGGGGGTGTCGAAACt                                                                                                                                                                                          <  1:325108‑M1/62‑1 (MQ=255)
                                                     aaGCGTCTGCGGCAGGCGTATCTCGTTGCACGTGATGCTCACGAGGGGCAAACACGTTCAAGCGGTGAACCCTATATCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCGCTGCTGCATGACGTGATTGAAGATACTCCCGCCACCTACCAGGATcaggatatggaacagctttttggtaaaagcgtcgccgagctggtagagggggtgtcga                                                                                                                                                                                        <  1:70129‑M2/250‑59 (MQ=255)
                                                        cggcaggcgtatctcgttgcacgtgatgctcacgaggggcaaacacgttcaagcggtgaaccctatatcacgcacccggtagcggttgcctgcattctggccgagatgaaactcgactatgaaacgctgatggcggcgctgctgcatgacgtgattgaagatactcccgccacctaccaggatCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGGGGGTGTCGAAACTTGATaa                                                                                                                                                                                    >  2:71778‑M1/184‑251 (MQ=255)
                                                        cGTCTGCGGCAGGCGTATCTCGTTGCACGTGATGCTCACGAGGGGCAAACACGTTCAAGCGGTGAACCCTATATCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCGCTGCTGCATGACGTGATTGAAGATACTCCCGCCACCTACCAGGATcaggatatggaacagctttttggtaaaagcgtcgccgagctggtagagggggtgtcgaaact                                                                                                                                                                                    <  1:325108‑M2/251‑63 (MQ=255)
                                                         ggcaggcgtatctcgttgcacgtgatgctcacgaggggcaaacacgttcaagcggtgaaccctatatcacgcacccggtagcggttgcctgcattctggccgagatgaaactcgactatgaaacgctgatggcggcgctgctgcatgacgtgattgaagatactcccgccacctaccaggatCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGGGGGTGTCGAAACTTGATaa                                                                                                                                                                                    <  1:262944‑M1/68‑1 (MQ=255)
                                                              cGGCAGGCGTATCTCGTTGCACGTGATGCTCACGAGGGGCAAACACGTTCAAGCGGTGAACCCTATATCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCGCTGCTGCATGACGTGATTGAAGATACTCCCGCCACCTACCAGGATcaggatatggaacagctttttggtaaaagcgtcgccgagctggtagagggggtgtcgaaacttgataa                                                                                                                                                                              >  2:71778‑M2/1‑183 (MQ=255)
                                                               ggCAGGCGTATCTCGTTGCACGTGATGCTCACGAGGGGCAAACACGTTCAAGCGGTGAACCCTATATCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCGCTGCTGCATGACGTGATTGAAGATACTCCCGCCACCTACCAGGATcaggatatggaacagctttttggtaaaagcgtcgccgagctggtagagggggtgtcgaaacttgataa                                                                                                                                                                              <  1:262944‑M2/250‑69 (MQ=255)
                                                                            acgtgatgctcacgaggggcaaacacgttcaagcggtgaaccctatatcacgcacccggtagcggttgcctgcattctggccgagatgaaactcgactatgaaacgctgatggcggcgctgctgcatgacgtgattgaagatactcccgccacctaccaggatCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGGGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATaaga                                                                                                                                                                >  1:153865‑M1/164‑251 (MQ=255)
                                                                                  aCGTGATGCTCACGAGGGGCAAACACGTTCAAGCGGTGAACCCTATATCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCGCTGCTGCATGACGTGATTGAAGATACTCCCGCCACCTACCAGGATcaggatatggaacagctttttggtaaaagcgtcgccgagctggtagagggggtgtcgaaacttgataaactcaagttccgcgataaga                                                                                                                                                          >  1:153865‑M2/1‑163 (MQ=255)
                                                                                   gctcacgaggggcaaacacgttcaagcggtgaaccctatatcacgcacccggtagcggttgcctgcattctggccgagatgaaactcgactatgaaacgctgatggcggcgctgctgcatgacgtgattgaagatactcccgccacctaccaggatCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGGGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATAAGAAAGAGgc                                                                                                                                                         >  2:220310‑M1/157‑251 (MQ=255)
                                                                                         gCTCACGAGGGGCAAACACGTTCAAGCGGTGAACCCTATATCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCGCTGCTGCATGACGTGATTGAAGATACTCCCGCCACCTACCAGGATcaggatatggaacagctttttggtaaaagcgtcgccgagctggtagagggggtgtcgaaacttgataaactcaagttccgcgataagaaagaggc                                                                                                                                                   >  2:220310‑M2/1‑156 (MQ=255)
                                                                                           ggggcaaacacgttcaagcggtgaaccctatatcacgcacccggtagcggttgcctgcattctggccgagatgaaactcgactatgaaacgctgatggcggcgctgctgcatgacgtgattgaagatactcccgccacctaccaggatCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGGGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATAAGAAAGAGGCGCAGGCCg                                                                                                                                                 >  1:168383‑M1/149‑251 (MQ=255)
                                                                                                 ggggCAAACACGTTCAAGCGGTGAACCCTATATCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCGCTGCTGCATGACGTGATTGAAGATACTCCCGCCACCTACCAGGATcaggatatggaacagctttttggtaaaagcgtcgccgagctggtagagggggtgtcgaaacttgataaactcaagttccgcgataagaaagaggcgcaggccg                                                                                                                                           >  1:168383‑M2/1‑148 (MQ=255)
                                                                                                 aacacgttcaagcggtgaaccctatatcacgcacccggtagcggttgcctgcattctggccgagatgaaactcgactatgaaacgctgatggcggcgctgctgcatgacgtgattgaagatactcccgccacctaccaggatCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGGGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATAAGAAAGAGGCGCAGGCCGAAAACt                                                                                                                                           <  1:52715‑M1/109‑1 (MQ=255)
                                                                                                    acgttcaagcggtgaaccctatatcacgcacccggtagcggttgcctgcattctggccgagatgaaactcgactatgaaacgctgatggcggcgcttctgcatgacgtgattgaagatactcccgccacctaccaggatCAGGATATGGAACAGCTTTTTGGTAACAGCGTCTCCGAGCTGGTAGAGGGGGTGTCGAAACTTGATAAACTCAATTTCCGCGATAAGAAAGAGTCGCCGGCCGAAAACTTTc                                                                                                                                        <  2:227494‑M1/112‑1 (MQ=255)
                                                                                                    acgttcaagcggtgaaccctatatcacgcacccggtagcggttgcctgcattctggccgagatgaaactcgactatgaaacgctgatggcggcgctgctgcatgacgtgattgaagatactcccgccacctaccaggatCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGGGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATAAGAAAGAGGCGCAGGCCGAAAACTTTc                                                                                                                                        <  1:209453‑M1/112‑1 (MQ=255)
                                                                                                       aaCACGTTCAAGCGGTGAACCCTATATCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCGCTGCTGCATGACGTGATTGAAGATACTCCCGCCACCTACCAGGATcaggatatggaacagctttttggtaaaagcgtcgccgagctggtagagggggtgtcgaaacttgataaactcaagttccgcgataagaaagaggcgcaggccgaaaact                                                                                                                                     <  1:52715‑M2/251‑110 (MQ=255)
                                                                                                         caagcggtgaaccctatatcacgcacccggtagcggttgcctgcattctggccgagatgaaactcgactatgaaacgctgatggcggcgctgctgcatgacgtgattgaagatactcccgccacctaccaggatCAGGATATGGAACAGCTGTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGGGGGTGTCGAAACTTGATAAACTCAAGTGCCGCGATAAGGAAGAGTCGCAGGGCGAAAACTTTCGCAAg                                                                                                                                   >  2:361235‑M1/135‑251 (MQ=255)
                                                                                                          acGTTCAAGCGGTGAACCCTATATCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCGCTTCTGCATGACGTGATTGAAGATACTCCCGCCACCTACCAGGATcaggatatggaacagctttttggtaacagcgtctccgagctggtagagggggtgtcgaaacttgataaactcaatttccgcgataagaaagagtcgccggccgaaaactttc                                                                                                                                  <  2:227494‑M2/251‑113 (MQ=255)
                                                                                                          acGTTCAAGCGGTGAACCCTATATCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCGCTGCTGCATGACGTGATTGAAGATACTCCCGCCACCTACCAGGATcaggatatggaacagctttttggtaaaagcgtcgccgagctggtagagggggtgtcgaaacttgataaactcaagttccgcgataagaaagaggcgcaggccgaaaactttc                                                                                                                                  <  1:209453‑M2/251‑113 (MQ=255)
                                                                                                           agcggtgaaccctatatcacgcacccggtagcggttgcctgcattctggccgagatgaaactcgactatgaaacgctgatggcggctcagctgcatgacgtgattgaagatactcccgccacctaccaggatCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCGGGTAGAGGGAGTGTCGAAACTTGATAAACTCAAGTTCCGCGATAAGAAAGAGGCGCAGGCCGAAAACTTTCGCAAgat                                                                                                                                 >  1:338844‑M1/133‑251 (MQ=255)
                                                                                                              ggtgaaccctatatcacgcacccggtagcggttgcctgcattctggccgagatgaaactcgactatgaaacgctgatggcggcgctgctgcatgacgtgattgaagatactcccgccacctaccaggatCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGGGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATAAGAAAGAGGCGCAGGCCGAAAACTTTCGCAAgatgat                                                                                                                              <  1:267905‑M1/122‑1 (MQ=255)
                                                                                                               gtgaaccctatatcacgcacccggtagcggttgcctgcattctggccgagatgaaactcgactatgaaacgctgatggcggcgctgctgcatgacgtgattgaagatactcccgccacctaccaggatCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGGGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATAAGAAAGAGGCGCAGGCCGAAAACTTTCGCAAGATGAtt                                                                                                                             <  1:164125‑M1/123‑1 (MQ=255)
                                                                                                               cAAGCGGTGAACCCTATATCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCGCTGCTGCATGACGTGATTGAAGATACTCCCGCCACCTACCAGGATcaggatatggaacagctgtttggtaaaagcgtcgccgagctggtagagggggtgtcgaaacttgataaactcaagtgccgcgataaggaagagtcgcagggcgaaaactttcgcaag                                                                                                                             >  2:361235‑M2/1‑134 (MQ=255)
                                                                                                                 aGCGGTGAACCCTATATCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCTCAGCTGCATGACGTGATTGAAGATACTCCCGCCACCTACCAGGATcaggatatggaacagctttttggtaaaagcgtcgccgagcgggtagagggagtgtcgaaacttgataaactcaagttccgcgataagaaagaggcgcaggccgaaaactttcgcaagat                                                                                                                           >  1:338844‑M2/1‑132 (MQ=255)
                                                                                                                    ggTGAACCCTATATCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCGCTGCTGCATGACGTGATTGAAGATACTCCCGCCACCTACCAGGATcaggatatggaacagctttttggtaaaagcgtcgccgagctggtagagggggtgtcgaaacttgataaactcaagttccgcgataagaaagaggcgcaggccgaaaactttcgcaagatgat                                                                                                                        <  1:267905‑M2/251‑123 (MQ=255)
                                                                                                                     gTGAACCCTATATCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCGCTGCTGCATGACGTGATTGAAGATACTCCCGCCACCTACCAGGATcaggatatggaacagctttttggtaaaagcgtcgccgagctggtagagggggtgtcgaaacttgataaactcaagttccgcgataagaaagaggcgcaggccgaaaactttcgcaagatgatt                                                                                                                       <  1:164125‑M2/251‑124 (MQ=255)
                                                                                                                        atatcacgcacccggtagcggttgcctgcattctggccgagatgaaactcgactatgaaacgctgatggcggcgctgctgcatgacgtgattgaagatactcccgccacctaccaggatCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGGGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATAAGAAAGAGGCGCAGGCCGAAAACTTTCGCAAGATGATTATGGCGATg                                                                                                                    >  1:18206‑M1/120‑251 (MQ=255)
                                                                                                                          atcacgcacccggtagcggttgcctgcattctggccgagatgaaactcgactatgaaacgctgatggcggcgctgctgcatgacgtgattgaagatactcccgccacctaccaggatCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGGGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATAAGAAAGAGGCGCAGGCCGAAAACTTTCGCAAGATGATTATGGCGATgg                                                                                                                   >  1:326627‑M1/118‑250 (MQ=255)
                                                                                                                              atatCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCGCTGCTGCATGACGTGATTGAAGATACTCCCGCCACCTACCAGGATcaggatatggaacagctttttggtaaaagcgtcgccgagctggtagagggggtgtcgaaacttgataaactcaagttccgcgataagaaagaggcgcaggccgaaaactttcgcaagatgattatggcgatg                                                                                                              >  1:18206‑M2/1‑119 (MQ=255)
                                                                                                                                cacccggtagcggttgcctgcattctggccgagatgaaactcgactatgaaacgctgatggcggcgctgctgcatgacgtgattgaagatactcccgccacctaccaggatCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGGGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATAAGAAAGAGGCGCAGGCCGAAAACTTTCGCAAGATGATTATGGCGATGGTGCAGGa                                                                                                            >  1:40366‑M1/112‑251 (MQ=255)
                                                                                                                                atCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCGCTGCTGCATGACGTGATTGAAGATACTCCCGCCACCTACCAGGATcaggatatggaacagctttttggtaaaagcgtcgccgagctggtagagggggtgtcgaaacttgataaactcaagttccgcgataagaaagaggcgcaggccgaaaactttcgcaagatgattatggcgatgg                                                                                                             >  1:326627‑M2/1‑117 (MQ=255)
                                                                                                                                      cACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCGCTGCTGCATGACGTGATTGAAGATACTCCCGCCACCTACCAGGATcaggatatggaacagctttttggtaaaagcgtcgccgagctggtagagggggtgtcgaaacttgataaactcaagttccgcgataagaaagaggcgcaggccgaaaactttcgcaagatgattatggcgatggtgcagga                                                                                                      >  1:40366‑M2/1‑111 (MQ=255)
                                                                                                                                       tagcggttgcctgcattctggccgagatgaaactcgactatgaaacgctgatggcggcgctgctgcatgacgtgattgaagatactcccgccacctaccaggatCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGGGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATAAGAAAGAGGCGCAGGCCGAAAACTTTCGCAAGATGATTATGGCGATGGTGCAGGATATCcgc                                                                                                     <  1:14401‑M1/147‑1 (MQ=255)
                                                                                                                                             tAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCGCTGCTGCATGACGTGATTGAAGATACTCCCGCCACCTACCAGGATcaggatatggaacagctttttggtaaaagcgtcgccgagctggtagagggggtgtcgaaacttgataaactcaagttccgcgataagaaagaggcgcaggccgaaaactttcgcaagatgattatggcgatggtgcaggatatccgc                                                                                               <  1:14401‑M2/251‑148 (MQ=255)
                                                                                                                                              tgcctgcattctggccgagatgaaactcgactatgaaacgctgatggcggcgctgctgcatgacgtgattgaagatactcccgccacctaccaggatCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGGGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATAAGAAAGAGGCGCAGGCCGAAAACTTTCGCAAGATGATTATGGCGATGGTGCAGGATATCCGCGTCATc                                                                                               >  1:58891‑M1/98‑250 (MQ=255)
                                                                                                                                              ggcctgcattctggccgagatgaaactcgactatgaaacgctgatggcggcgctgctgcatgacgtgattgaagatactcccgccacctaccaggatCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGGGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATAAGAAAGAGGCGCAGGCCGAAAACTTTCGCAAGATGATTATGGCGATGGTGCAGGATATCCGCGTCATcc                                                                                              >  2:106296‑M1/98‑251 (MQ=255)
                                                                                                                                                    tGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCGCTGCTGCATGACGTGATTGAAGATACTCCCGCCACCTACCAGGATcaggatatggaacagctttttggtaaaagcgtcgccgagctggtagagggggtgtcgaaacttgataaactcaagttccgcgataagaaagaggcgcaggccgaaaactttcgcaagatgattatggcgatggtgcaggatatccgcgtcatc                                                                                         >  1:58891‑M2/1‑97 (MQ=255)
                                                                                                                                                    ggcCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCGCTGCTGCATGACGTGATTGAAGATACTCCCGCCACCTACCAGGATcaggatatggaacagctttttggtaaaagcgtcgccgagctggtagagggggtgtcgaaacttgataaactcaagttccgcgataagaaagaggcgcaggccgaaaactttcgcaagatgattatggcgatggtgcaggatatccgcgtcatcc                                                                                        >  2:106296‑M2/2‑97 (MQ=255)
                                                                                                                                                                gatgaaactcgactatgaaacgctgatggcggcgctgctgcatgacgtgattgaagatactcccgccacctaccaggatCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGGGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATAAGAAAGAGGCGCAGGCCGAAAACTTTCGCAAGATGATTATGGCGATGGTGCAGGATATCCGCGTCATCCTCATCAAACTTGCCGAcc                                                                            <  1:83841‑M1/172‑1 (MQ=255)
                                                                                                                                                                      gaTGAAACTCGACTATGAAACGCTGATGGCGGCGCTGCTGCATGACGTGATTGAAGATACTCCCGCCACCTACCAGGATcaggatatggaacagctttttggtaaaagcgtcgccgagctggtagagggggtgtcgaaacttgataaactcaagttccgcgataagaaagaggcgcaggccgaaaactttcgcaagatgattatggcgatggtgcaggatatccgcgtcatcctcatcaaacttgccgacc                                                                      <  1:83841‑M2/251‑173 (MQ=255)
                                                                                                                                                                             tatgaaacgctgatggcggcgctgctgcatgacgtgattgaagatactcccgccacctaccaggatCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGGGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATAAGAAAGAGGCGCAGGCCGAAAACTTTCGCAAGATGATTATGGCGATGGTGCAGGATATCCGCGTCATCCTCATCAAACTTGCCGACCGTACCCACAACAt                                                               <  2:115121‑M1/185‑1 (MQ=255)
                                                                                                                                                                              atgaaacgctgatggcggcgctgctgcatgacgtgattgaagatactcccgccacctaccaggatCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGGGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATAAGAAAGAGGCGCAGGCCGAAAACTTTCGCAAGATGATTATGGCGATGGTGCAGGATATCCGCGTCATCCTCATCAAACTTGCCGACCGTACCCACAACAt                                                               >  1:305011‑M1/66‑250 (MQ=255)
                                                                                                                                                                                   tATGAAACGCTGATGGCGGCGCTGCTGCATGACGTGATTGAAGATACTCCCGCCACCTACCAGGATcaggatatggaacagctttttggtaaaagcgtcgccgagctggtagagggggtgtcgaaacttgataaactcaagttccgcgataagaaagaggcgcaggccgaaaactttcgcaagatgattatggcgatggtgcaggatatccgcgtcatcctcatcaaacttgccgaccgtacccacaacat                                                         <  2:115121‑M2/251‑186 (MQ=255)
                                                                                                                                                                                    aTGAAACGCTGATGGCGGCGCTGCTGCATGACGTGATTGAAGATACTCCCGCCACCTACCAGGATcaggatatggaacagctttttggtaaaagcgtcgccgagctggtagagggggtgtcgaaacttgataaactcaagttccgcgataagaaagaggcgcaggccgaaaactttcgcaagatgattatggcgatggtgcaggatatccgcgtcatcctcatcaaacttgccgaccgtacccacaacat                                                         >  1:305011‑M2/1‑65 (MQ=255)
                                                                                                                                                                                      gtgatggcggcgctgctgcatgacgtgatggaagatactcccgccacctaccaggatCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGGGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATAAGAAAGAGGCGCAGGCCCAAAACTTTCGCAAGATGAGTATGGCGATGGTGCAGGATATCCGCGTCATCCTCATCAAACTTTCCGACCGTACGCACAGCATGCGCACGCt                                                      >  2:358763‑M1/58‑251 (MQ=255)
                                                                                                                                                                                            gtgATGGCGGCGCTGCTGCATGACGTGATGGAAGATACTCCCGCCACCTACCAGGATcaggatatggaacagctttttggtaaaagcgtcgccgagctggtagagggggtgtcgaaacttgataaactcaagttccgcgataagaaagaggcgcaggcccaaaactttcgcaagatgagtatggcgatggtgcaggatatccgcgtcatcctcatcaaactttccgaccgtacgcacagcatgcgcacgct                                                >  2:358763‑M2/2‑57 (MQ=255)
                                                                                                                                                                                                                   tgaagatactcccgccacctaccaggatCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGGGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATAAGAAAGAGGCGCAGGCCGAAAACTTTCGCAAGATGATTATGGCGATGGTGCAGGATATCCGCGTCATCCTCATCAAACTTGCCGACCGTACCCACAACATGCGCACGCTGGGCTCACTTCGCCCGGACGAACGTcgcc                         <  2:121676‑M1/223‑1 (MQ=255)
                                                                                                                                                                                                                      agatactcccgccacctaccaggatCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGGGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATAAGAAAGAGGCGCAGGCCGAAAACTTTCGCAAGATGATTATGGCGATGGTGCAGGATATCCGCGTCATCCTCATCAAACTTGCCGACCGTACCCACAACATGCGCACGCTGGGCTCACTTCGCCCGGACAAACGTCGCCGCa                      >  2:325108‑M1/26‑251 (MQ=255)
                                                                                                                                                                                                                         tactcccgccacctaccaggatCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGCGCTGGTAGAGGGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATAAGAAAGAGGCGCAGGCCGAAAACTTTCGCAAGATGATTATGGCGATGGTGCAGGATATCCGCGTCATCCTCATCAAACTTGCCGACCGTACCCACAACATGCGCACGCTGGGCTCACTTCGCCCGGACAAACGTCGCCGCATCg                   <  2:39877‑M1/229‑1 (MQ=255)
                                                                                                                                                                                                                         tGAAGATACTCCCGCCACCTACCAGGATcaggatatggaacagctttttggtaaaagcgtcgccgagctggtagagggggtgtcgaaacttgataaactcaagttccgcgataagaaagaggcgcaggccgaaaactttcgcaagatgattatggcgatggtgcaggatatccgcgtcatcctcatcaaacttgccgaccgtacccacaacatgcgcacgctgggctcacttcgcccggacgaacgtcgcc                   <  2:121676‑M2/251‑224 (MQ=255)
                                                                                                                                                                                                                            aGATACTCCCGCCACCTACCAGGATcaggatatggaacagctttttggtaaaagcgtcgccgagctggtagagggggtgtcgaaacttgataaactcaagttccgcgataagaaagaggcgcaggccgaaaactttcgcaagatgattatggcgatggtgcaggatatccgcgtcatcctcatcaaacttgccgaccgtacccacaacatgcgcacgctgggctcacttcgcccggacaaacgtcgccgca                >  2:325108‑M2/1‑25 (MQ=255)
                                                                                                                                                                                                                             cccgccacctaccaggatCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGGGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATAAGAAAGAGGCGCAGGCCGAAAACTTTCGCAAGATGATTATGGCAATGGTGCAGGATATCCGCGTCATCCTCATCAAACTTGCCGACCGTACCCACAACATGCGCACGCTGGGCTCACTTCGCCCGGACAAACGTCGCCGCATCGccc                >  1:272696‑M1/19‑250 (MQ=255)
                                                                                                                                                                                                                               tACTCCCGCCACCTACCAGGATcaggatatggaacagctttttggtaaaagcgtcgccgcgctggtagagggggtgtcgaaacttgataaactcaagttccgcgataagaaagaggcgcaggccgaaaactttcgcaagatgattatggcgatggtgcaggatatccgcgtcatcctcatcaaacttgccgaccgtacccacaacatgcgcacgctgggctcacttcgcccggacaaacgtcgccgcatcg             <  2:39877‑M2/251‑230 (MQ=255)
                                                                                                                                                                                                                                   cccGCCACCTACCAGGATcaggatatggaacagctttttggtaaaagcgtcgccgagctggtagagggggtgtcgaaacttgataaactcaagttccgcgataagaaagaggcgcaggccgaaaactttcgcaagatgattatggcaatggtgcaggatatccgcgtcatcctcatcaaacttgccgaccgtacccacaacatgcgcacgctgggctcacttcgcccggacaaacgtcgccgcatcgccc          >  1:272696‑M2/1‑18 (MQ=255)
                                                                                                                                                                                                                                          aggatCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGGGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATAAGAAAGAGGCGCAGGCCGAAAACTTTCGCAAGATGATTATGGCGATGGTGCAGGATATCCGCGTCATCCTCATCAAACTTGCCGACCGTACCCACAACATGCGCACGCTGGGCTCACTTCGCCCGGACAAACGTCGCCGCATCGCCCGTGAAACTCTCGaa  >  1:360554‑M1/6‑251 (MQ=255)
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                     
TGTTTGAAAGCCTGAATCAACTGATTCAAACCTACCTGCCGGAAGACCAAATCAAGCGTCTGCGGCAGGCGTATCTCGTTGCACGTGATGCTCACGAGGGGCAAACACGTTCAAGCGGTGAACCCTATATCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCGCTGCTGCATGACGTGATTGAAGATACTCCCGCCACCTACCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGGGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATAAGAAAGAGGCGCAGGCCGAAAACTTTCGCAAGATGATTATGGCGATGGTGCAGGATATCCGCGTCATCCTCATCAAACTTGCCGACCGTACCCACAACATGCGCACGCTGGGCTCACTTCGCCCGGACAAACGTCGCCGCATCGCCCGTGAAACTCTCGAA  >  NC_000913/3822407‑3822891

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