Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 4,139,771 10483 bp→TG [ptsA][eptC] [ptsA], frwC, frwB, pflD, pflC, frwD, yijO, [eptC]

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 4139771 4150253 10483 57 [0] [0] 56 [ptsA]–[eptC] [ptsA],frwC,frwB,pflD,pflC,frwD,yijO,[eptC]

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 41397700 (0.000)56 (0.960)
+TG
54/490 0.1 100% coding (2452/2502 nt) ptsA putative PTS multiphosphoryl transfer protein PtsA
?NC_000913 4150254 = 0 (0.000)coding (12/1734 nt) eptC phosphoethanolamine transferase EptC

GGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCGGCGTTAACCAGCGCGTTACGGAATTCTTCGTGCATGATGCGGCGAGCGAGGCGCGAGAAAATGCGCATATGCTGATCGCCCGCAGCGTGTTTGTTCAGGGTTAACATAATGATGAATTGCGCTTCATCATCGCCCCAGCGCACCGGCGCTTGCAGACGCGCCACGCTGATGGTGGATTGCTCAATGTGTTCTGATTTGCTGTGTGGAATGGCAAAACTAAAGCCCAGACCGGTAGAGAAAACGGCCT  >  NC_000913/4139522‑4140018
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                        
ggCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtg                                                                                                                                                                                                                                                        <  1:420450‑M1/251‑3 (MQ=255)
     aTGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtgg                                                                                                                                                                                                                                                   >  2:121981‑M1/1‑244 (MQ=255)
          cTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgc                                                                                                                                                                                                                                              >  2:178208‑M1/1‑239 (MQ=255)
               cAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagt                                                                                                                                                                                                                                         >  2:62690‑M1/1‑234 (MQ=255)
                 aaCAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagta                                                                                                                                                                                                                                        >  1:37049‑M1/1‑232 (MQ=255)
                     gAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaac                                                                                                                                                                                                                                   >  1:329394‑M1/1‑228 (MQ=255)
                                  cccaATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaa                                                                                                                                                                                                                      <  1:375649‑M1/250‑37 (MQ=255)
                                   ccAATCGCTTTTGGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGATGTCACCCCGGCAATGGGGAATATGCGTGCCAGCCGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtggggaatgcatagtaagaacctgtcttgcaaaa                                                                                                                                                                                                                      >  2:338124‑M1/1‑214 (MQ=255)
                                     aaTCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGGGTCCAGATACAGTTACATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGTCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctttcttgacaaaat                                                                                                                                                                                                                    >  2:271964‑M1/1‑212 (MQ=255)
                                     aaTCGCTTTTTGCAGACGAGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaata                                                                                                                                                                                                                   <  1:338124‑M1/251‑40 (MQ=255)
                                      aTCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgcaaaaatat                                                                                                                                                                                                                  >  1:64033‑M1/1‑211 (MQ=255)
                                           tttttGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccg                                                                                                                                                                                                             >  2:46407‑M1/1‑206 (MQ=255)
                                                 cAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgta                                                                                                                                                                                                       <  1:446771‑M1/251‑52 (MQ=255)
                                                        ggCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaat                                                                                                                                                                                               >  1:59775‑M1/1‑192 (MQ=255)
                                                              aCTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattcctt                                                                                                                                                                                          >  1:320455‑M1/1‑187 (MQ=255)
                                                                 tCCCATATGGACTCCTTGCTGGCAGCGATAATGCGCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaa                                                                                                                                                                                        <  1:368681‑M1/250‑67 (MQ=255)
                                                                        tGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatct                                                                                                                                                                                >  1:154022‑M1/1‑177 (MQ=255)
                                                                                 tGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattcctttaggatatctgaaggtat                                                                                                                                                                        <  1:84920‑M1/250‑83 (MQ=255)
                                                                                   cTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatatt                                                                                                                                                                     >  1:179465‑M1/1‑166 (MQ=255)
                                                                                              aaTGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGTGTGCCAGACGTTCGACTTCTGCGCCGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgactgcaacgacaattccttaaggatatctgaaggaatattcagaatttga                                                                                                                                                           >  2:368158‑M1/1‑155 (MQ=255)
                                                                                               attccctggttaattgTCACCCCGGCAATGGGGAATATTCGTTCCAGACGTTCGGCTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaatgtatattcagaatgtgaat                                                                                                                                                         <  2:382463‑M1/244‑98 (MQ=255)
                                                                                                   tCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaa                                                                                                                                                      >  2:194055‑M1/1‑150 (MQ=255)
                                                                                                                cACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaat                                                                                                                                        <  2:164311‑M1/251‑115 (MQ=255)
                                                                                                                 aCCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaat                                                                                                                                        >  1:12635‑M1/1‑136 (MQ=255)
                                                                                                                  ccccGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaatat                                                                                                                                      <  1:185695‑M1/251‑117 (MQ=255)
                                                                                                                         aaTGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaatatattgaaa                                                                                                                               <  1:429583‑M1/251‑124 (MQ=255)
                                                                                                                          aTGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaatatattgaaa                                                                                                                               >  1:169298‑M1/1‑127 (MQ=255)
                                                                                                                            ggggAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaatatattgaaaacg                                                                                                                            <  1:121981‑M1/251‑127 (MQ=255)
                                                                                                                                   atGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaatatattgaaaacgagggtg                                                                                                                      <  1:96157‑M1/250‑133 (MQ=255)
                                                                                                                                           cAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaatatattgaaaacgagggtgttagaacag                                                                                                             <  2:329394‑M1/251‑142 (MQ=255)
                                                                                                                                                  tCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaatatattgaaaacgagggtgttagaacagaagtat                                                                                                       <  1:368158‑M1/250‑148 (MQ=255)
                                                                                                                                                        tCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccataaggatatctgaaggtatattcagaagttgaataaaatgcagacagaaatatattgaaaacgagggtgttagaacagaagtagttcag                                                                                                  <  2:156726‑M1/249‑153 (MQ=255)
                                                                                                                                                               cGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaatatattgaaaacgagggtgttagaacagaagtatttcagaaaaccctc                                                                                         <  1:46407‑M1/251‑162 (MQ=255)
                                                                                                                                                               cGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaatatattgaaaacgagggtgttagaacagaagtatttcagaaaaccctc                                                                                         <  1:178208‑M1/251‑162 (MQ=255)
                                                                                                                                                                  tAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaatatattgaaaacgagggtgttagaacagaagtatttcagaaaaccctcgcg                                                                                      >  2:448145‑M1/1‑87 (MQ=255)
                                                                                                                                                                         gtCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaatatattgaaaacgagggtgttagaacagaagtatttcagaaaaccctcgcgcaaaagc                                                                               <  2:59775‑M1/251‑172 (MQ=255)
                                                                                                                                                                               aTACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaatatattgaaaacgagggtgttagaacagaagtatttcagaaaaccctcgcgcaaaagcacgagg                                                                         >  1:309846‑M1/1‑74 (MQ=255)
                                                                                                                                                                                    gTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaatatattgaaaacgagggtgttagaacagaagtatttcagaaaaccctcgcgcaaaagcacgagggtttg                                                                    <  1:212462‑M1/251‑183 (MQ=255)
                                                                                                                                                                                      tCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaatatattgaaaacgagggtgttagaacagaagtatttcagaaaaccctcgcgcaaaagcacgagggtttgc                                                                   <  1:396699‑M1/250‑184 (MQ=255)
                                                                                                                                                                                                         cagttccagttcATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaatatattgaaaacgagggtgttagaacagaagtatttcagaaaaccctcgcgcaaaagcacgagggtttgcagaagaggaagattagccgg                                               >  1:189466‑M1/1‑48 (MQ=255)
                                                                                                                                                                                                         cagttccagttcATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaatatattgaaaacgagggtgttagaacagaagtatttcagaaaaccctcgcgcaaaagcacgagggtttgcagaagaggaagattagccgg                                               >  2:15594‑M1/1‑48 (MQ=255)
                                                                                                                                                                                                           gttccagttcATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaatatattgaaaacgagggtgttagaacagaagtatttcagaaaaccctcgcgcaaaagcacgagggtttgcagaagaggaagat                                                      <  2:345852‑M1/242‑197 (MQ=255)
                                                                                                                                                                                                           gttccagttcATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaatatattgaaaacgagggtgttagaacagaagtatttcagaaaaccctcgcgcaaaagcacgagggtttgcagaagaggaagat                                                      >  1:345852‑M1/1‑46 (MQ=255)
                                                                                                                                                                                                               cagttcATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaatatattgaaaacgagggtgttagaacagaagtatttcagaaaaccctcgcgcaaaagcacgagggtttgcagaagaggaagattagccggtattac                                         <  1:24719‑M1/251‑210 (MQ=255)
                                                                                                                                                                                                                  ttcATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagccagaaatatattgaaaacgagggtgttagaacagaagtatttcagaaaaccctcgcgcaaacgcacgagggtttgcagaagaggaagattagccggtattacgc                                       >  2:438762‑M1/1‑39 (MQ=255)
                                                                                                                                                                                                                       gTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaatatattgaaaacgagggtgttagaacagaagtatttcagaaaaccctcgcgcaaaagcacgagggtttgcagaagaggaagattagccgg                                               <  1:379491‑M1/237‑204 (MQ=255)
                                                                                                                                                                                                                       gTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaatatattgaaaacgagggtgttagaacagaagtatttcagaaaaccctcgcgcaaaagcacgagggtttgcagaagaggaagattagccgg                                               >  2:379491‑M1/1‑34 (MQ=255)
                                                                                                                                                                                                                            agcagGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaatatattgaaaacgagggtgttagaacagaagtatttcagaaaaccctcgcgcaaaagcacgagggtttgcagaagaggaagattagccggtattacgcatccctgccgc                            <  1:40167‑M1/251‑223 (MQ=255)
                                                                                                                                                                                                                            agcagGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaatatattgaaaacgagggtgttagaacagaagtatttcagaaaaccctcgcgcaaaagcacgagggtttgcagaagaggaagattagccggtattacgcatacctgccgc                            >  1:76876‑M1/1‑29 (MQ=255)
                                                                                                                                                                                                                                        aTAGCGGCGGCAGAGGGtggtttgaaggcatagtaagaacctttctttaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaatatattgaaaacgagggtgttagaacagaagtatttcagaaaaccctcgcgcaaaagcacgagggtttgcagaagaggaagattagccggtattacgcatacctgcctcaatcccggcaat                <  1:370685‑M1/251‑235 (MQ=255)
                                                                                                                                                                                                                                          aGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaatatattgaaaacgagggtgttagaacagaagtatttcagaaaaccctcgcgcaaaagcacgagggattgcagaagaggaagattagccggtattacgcatacctgccgcaatcccggcaat                >  2:87219‑M1/1‑15 (MQ=255)
                                                                                                                                                                                                                                          aGCGTCGGCAGAGGCggtggggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattaagaatttgaataaaatgcagacagaaatatattgaaaacgagggtgttagaacagaagtatttcagaaacccctcgcgcaaaagcacgagggtttgcagaagaggaagattagccggtattacgcatccctgccgcaatcccggcaatag              <  2:385630‑M1/251‑237 (MQ=255)
                                                                                                                                                                                                                                               cGGCAGAGTCtgtgtggaatgcagattaagaccctgtcttgaaaaaatatcgccgaatgtaacgccaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaatatattgaaaacgagggtgttagaacagaagtatttcagaaaaccctcgcgaaaaagctcgagggttggcagaagaggaagattagccggtattccgcattcctgccgcaatcccggcaatagtgac          <  2:303107‑M1/250‑241 (MQ=255)
                                                                                                                                                                                                                                                ggCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaatatattgaaaacgagggtgttagaacagaagtatttcagaaaaccctcgcgcaaaagcacgagggtttgcagaagaggaagattagccggtattacgcatacctgccgcaatcccggcaatagtgacc         >  1:397985‑M1/1‑9 (MQ=255)
                                                                                                                                                                                                                                                 gCAGAGGCtgtttggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaatatattgaaaacgagggtgttagaacagaagtatttcagaaacccctcgcgcaaaagcacgagggtttgcagaagaggaagattagccggtattacgcatacctgccgcaatcccggcaatagtgaccat       <  1:385903‑M1/251‑244 (MQ=255)
                                                                                                                                                                                                                                                     agGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaatatattgaaaacgagggtgttagaacagaagtatttcagaaaaccctcgcgcaaaagcacgagggtttgcagaagaggaagattagccggtattacgcatacctgccgcaatcccggcaatagtgaccattaac   >  1:301153‑M1/1‑4 (MQ=255)
                                                                                                                                                                                                                                                      ggCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaatatattgaaaacgagggtgttagaacagaagtatttcagaaaaccctcgcgcaaaagcacgagggtttgcagaagaggaagattagccggtattacgcatacctgccgcaatcccggcaatagtgaccattaacg  >  1:98450‑M1/1‑3 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                        
GGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCGGCGTTAACCAGCGCGTTACGGAATTCTTCGTGCATGATGCGGCGAGCGAGGCGCGAGAAAATGCGCATATGCTGATCGCCCGCAGCGTGTTTGTTCAGGGTTAACATAATGATGAATTGCGCTTCATCATCGCCCCAGCGCACCGGCGCTTGCAGACGCGCCACGCTGATGGTGGATTGCTCAATGTGTTCTGATTTGCTGTGTGGAATGGCAAAACTAAAGCCCAGACCGGTAGAGAAAACGGCCT  >  NC_000913/4139522‑4140018

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