Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 167,401 Δ5,915 bp fhuA[fhuB] fhuA, fhuC, fhuD, [fhuB]

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 167401 173315 5915 60 [0] [0] 60 fhuA–[fhuB] fhuA,fhuC,fhuD,[fhuB]

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 1674000 (0.000)60 (0.920) 52/490 0.1 100% intergenic (+136/‑84) mrcB/fhuA peptidoglycan glycosyltransferase/peptidoglycan DD‑transpeptidase MrcB/ferrichrome outer membrane transporter/phage receptor
?NC_000913 173316 = 0 (0.000)coding (1855/1983 nt) fhuB iron(III) hydroxamate ABC transporter membrane subunit

GCGATTACCCCGCTGTGCCGCCGCTGGCTGACCATTTTACCGCTGGGTGGTGATACCGCCCGAGCCGTAGGAATGGCGCTGACGCCGACGCGAATTGCGCTGCTGCTGTTAGCGGCTTGCCTGACGGCGACCGCGACGATGACTATTGGACCGTTGAGTTTTGTTGGTTTAATGGCACCGCATATTGCGCGGATGATGGGCTTTCGACGGACGATGCCACACATCGTAATTTCGGCGCTGGTGGGTGGTTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCG  >  NC_000913/173067‑173562
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                      
ttgatcagtcttcccagccgcagcaacagccgcaacagcaacctgctcagcaagagcagaaagacagcgacggtgtagccggttggatcaaggatatgtttggtagtaattaacatctaagcgtgaaataccggatggcgagttgccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaacct                                                                                                                                                                                                                                                        <  2:394478‑M2/1‑1 (MQ=255)
          ttctcagccgcagcaacagccgcaacagcaacctgctcagcaagagcagaaagacagcgacggtgtagccggttggatcaaggatatgtttggtagtaattaacatctaagcgtgaaataccggatggcgagttgccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGgtg                                                                                                                                                                                                                                             >  1:185435‑M2/240‑251 (MQ=255)
           tctcagccgcagcaacagccgcaacagcaacctgctcagcaagagcagaaagacagcgacggtgtagccggttggatcaaggatatgtttggtagtaattaacatctaagcgtgaaataccggatggcgagttgccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGgtgt                                                                                                                                                                                                                                            >  2:492557‑M2/239‑251 (MQ=255)
                gccgcagcaacagccgcaacagcaacctgctcagcaagagcagaaagacagcgacggtgtagccggttggatcaaggatatgtttggtagtaattaacatctaagcgtgaaataccggatggcgagttgccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCt                                                                                                                                                                                                                                       <  1:521070‑M2/18‑1 (MQ=255)
                      gcaacagccgcaacagcaacctgctcagcaagagcagaaagacagcgacggtgtagccggttggatcaaggatatgtttggtagtaattaacatctaagcgtgaaataccggatggcgagttgccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTgg                                                                                                                                                                                                                                 >  1:131550‑M2/228‑251 (MQ=255)
                      gcaacagccgcaacagcaacctgctcagcaagagcagaaagacagcgacggtgtagccggttggatcaaggatatgtttggtagtaattaacatctaagcgtgaaataccggatggcgagttgccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTgg                                                                                                                                                                                                                                 >  1:89022‑M2/228‑251 (MQ=255)
                         ccagccgcaacagcaacctgctcagcaagagcagaaagacagcgacggtgtagccggttggatcaaggatatgtttggtagtaattaacatctaagcgtgaaataccggatggcgagttgccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGgtgt                                                                                                                                                                                                                              <  2:89022‑M2/27‑1 (MQ=255)
                             ccgcaacagcaacctgctcagcaagagcagaaagacagcgacggtgtagccggttggatcaaggatatgtttggtagtaattaacatctaagcgtgaaataccggatggcgagttgccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGCGGGc                                                                                                                                                                                                                          <  1:504225‑M2/31‑1 (MQ=255)
                                            gctcagcaagagcagaaagacagcgacggtgtagccggttggatcaaggatatgtttggtagtaattaacatctaagcgtgaaataccggatggcgagttgccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGtt                                                                                                                                                                                                             <  1:233923‑M2/44‑1 (MQ=255)
                                                 gcaagagcagaaagacagcgacggtgtagccggttggatcaaggatatgtttggtagtaattaacatctaagcgtgaaataccggatggcgagttgccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGttt                                                                                                                                                                                                            <  1:213115‑M2/45‑1 (MQ=255)
                                                 gcaagagcagaaagacagcgacggtgtagccggttggatcaaggatatgtttggtagtaattaacatctaagcgtgaaataccggatggcgagttgccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGttt                                                                                                                                                                                                            >  2:213115‑M2/201‑245 (MQ=255)
                                                    agagcagaaagacagcgacggtgtagccggttggatcaaggatatgtttggtagtaattaacatctaagcgtgaaataccggatggcgagttgccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAg                                                                                                                                                                                                   <  1:151490‑M2/54‑1 (MQ=255)
                                                        cagaaagacagcgacggtgtagccggttggatcaaggatgtgtttggtagtaattaacatctaagcgtgaaataccggatggcgagttgccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaagccttcaggaacgctcagattgcggaccgcttgcgaaccTTACTGATGGTGTTCGCTGACTGGTGTGGGAGGGTGGTGCTGTTTCCATTCCAGATcc                                                                                                                                                                                               >  1:195579‑M2/194‑251 (MQ=255)
                                                             agacagcgacggtgtagccggttggctcaaggatatgtttggtagtaattaacatctaagcgtgaaataccggatggcgagttgccatacggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaaccctccaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGTTGTTCGCTGACTGGTGTGGGCGGATGGTGATGTTTCCATTCCAGCTCCcggcg                                                                                                                                                                                          >  1:459916‑M2/189‑251 (MQ=255)
                                                                 agcgacggtgtagccggttggatcaaggctatgtttggtagtaattaacatctaagcgtgaaataccggatggcgagttgccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGAGGGCGGATGGTGCTGTTTCCCTTCCAGATCCCGGcggggg                                                                                                                                                                                      >  2:298139‑M2/185‑250 (MQ=255)
                                                                 agcgacggtgtagccggttggatcaaggatatgtttggtagtaattaacatctaagcgtgaaataccggatggcgagttgccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGgc                                                                                                                                                                                      >  1:533785‑M2/185‑251 (MQ=255)
                                                                        gtgtagccggttggatcaaggatatgtttggtagtaattaacatctaagcgtgaaataccggatggcgagttgccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTc                                                                                                                                                                               <  1:155726‑M2/74‑1 (MQ=255)
                                                                               cggttggatcaaggatatgtttggtagtaattaacatctaagcgtgaaataccggatggcgagttgccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCttt                                                                                                                                                                        <  1:17122‑M2/81‑1 (MQ=255)
                                                                                    ggatcaaggatatgtttggtagtaattaacatctaagcgtgaaataccggatggcgagttgccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCg                                                                                                                                                                    <  1:492557‑M2/85‑1 (MQ=255)
                                                                                       tcaaggatatgtttggtagtaattaacatctaagcgtgaaataccggatggcgagttgccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGgcgc                                                                                                                                                                >  2:300727‑M2/163‑251 (MQ=255)
                                                                                          aggatatgtttggtagtaattaacatctaagcgtgaaataccggatggcgagttgccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttggttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGcc                                                                                                                                                             <  1:66104‑M2/92‑1 (MQ=255)
                                                                                               atgtttggtagtaattaacatctaagagtgaaataccggatggcgagttgccatccggaaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGCGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTTGCCACCGTTATCGGCGCGTCATAtt                                                                                                                                                        >  2:344768‑M2/155‑251 (MQ=255)
                                                                                                 gttgggtagtaattaacatccaatcgtgaantgccggatggcgggttgccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttttttattaacccttcaggaacgctcagattgcgtacctcttgcgaaccTTACTGCTGGTGTGCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCTGCGCGCCATAttt                                                                                                                                                       <  2:239141‑M2/98‑1 (MQ=255)
                                                                                                   ttggtagtaattaacatctaagcgtgaaataccggatggcgagttgccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcgtgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTa                                                                                                                                                     <  1:49002‑M2/100‑1 (MQ=255)
                                                                                                            attaacatctaagcgtgaaataccggatggcgagttgccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTc                                                                                                                                                                               <  1:295352‑M2/74‑1 (MQ=255)
                                                                                                            attaacatctaagcgtgaaataccggatggcgagttgccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTc                                                                                                                                                                               >  2:295352‑M2/142‑215 (MQ=255)
                                                                                                                                   cggatggcgagttgccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTTGGAAAGAAGAGCCGGTAAtttt                                                                                                                    >  2:167659‑M2/119‑251 (MQ=255)
                                                                                                                                    ggatggcgagttgccatccggtaaaataacctcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAAttttt                                                                                                                   <  1:192346‑M2/134‑1 (MQ=255)
                                                                                                                                      atggcgagttgccatccgttaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTcc                                                                                                                 <  2:47134‑M2/136‑1 (MQ=255)
                                                                                                                                      atggcgagttgccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTcc                                                                                                                 <  1:140017‑M2/136‑1 (MQ=255)
                                                                                                                                                   atccggtaaaataacatcccatcaaatatatnaaccctttcttttcatctggttgtttattatcccttcaggaacgctcagattgcgtaccgctttcgaaccTTACTGCTTGTGTTCGCTGACTGGTGTTGGCGGATGGTTCTGTTTCCATTCCAGATCCCGGCGGGTCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTGTTTGTTGATAAAGCAGAGCCGTTAATTGTTCCATATAACGGCCCt                                                                                                    <  2:250500‑M2/149‑1 (MQ=255)
                                                                                                                                                        gtaaaataacatcccatctaagatattaaccctttcttttcatctggctgtttataaacccttaaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACATTTCTCCGCGCGCCCTATTTTATCTCTTTTATGAGGAAACAGAGCCGTTAATTTTTCCATATAACGGGCCCGTACg                                                                                               >  2:482301‑M2/98‑251 (MQ=255)
                                                                                                                                                          aaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTTTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGcc                                                                                             >  1:338234‑M2/96‑251 (MQ=255)
                                                                                                                                                           aaataccatcccatctaagatattaaccctttcttttaatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgcaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTttccattc                                                                                                                                                                                                      >  2:421616‑M2/95‑145 (MQ=255)
                                                                                                                                                           aaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTttccattc                                                                                                                                                                                                      <  1:421616‑M2/51‑1 (MQ=255)
                                                                                                                                                           aaaaaacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCt                                                                                            <  2:338234‑M2/157‑1 (MQ=255)
                                                                                                                                                             ataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGGGAAAGCAGAGCcgttaatttttccatataagatcgggagagcacacc                                                                                          >  1:350860‑M2/93‑234 (MQ=255)
                                                                                                                                                             ataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCAtata                                                                                                            <  2:350860‑M2/141‑1 (MQ=255)
                                                                                                                                                                   tcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGGTGCTGTCAACCTTTTACGGCGCGCCATATTTTAACTTTTTGTTGAGAAAGCAGAGGCGTTAATTTTTCCATATAACGGACCGGTACGCCTGGacgggt                                                                                    >  1:205411‑M2/87‑249 (MQ=255)
                                                                                                                                                                           aagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGtt                                                                            <  1:468254‑M2/173‑1 (MQ=255)
                                                                                                                                                                           aaaatatgaaccctttctattcatctggttgtttattaacccttcaggaacgctcagcttgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGCGCTGTTTCCATTCCAGATCCCGGCGGGGGTGCTGTCAACCTTTTTTGGCGCGCCATATTTTCTCTATTTGTTGAGAAAGCAGAGCCGTTAatttttccctataacggcactgttgacatggccggcgaacacatt                                                                            >  2:451345‑M2/79‑229 (MQ=255)
                                                                                                                                                                              atattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTAc                                                                         <  1:446505‑M2/176‑1 (MQ=255)
                                                                                                                                                                                 ttaaccctttcttttcatctggttgtttnttaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTAc                                                                         <  2:251066‑M2/176‑1 (MQ=255)
                                                                                                                                                                                 ttaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTAc                                                                         >  1:251066‑M2/73‑248 (MQ=255)
                                                                                                                                                                                      cctttcttttcatctgtttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAAc                                                                 <  1:193706‑M2/184‑1 (MQ=255)
                                                                                                                                                                                               tcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGccc                                                         <  1:330469‑M2/192‑1 (MQ=255)
                                                                                                                                                                                                 atctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGAt                                                      >  1:110519‑M2/57‑251 (MQ=255)
                                                                                                                                                                                                      gttgtttattaacccttcaggaacactcagattgcgtaccgcttgcgatccTTACTTCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAATCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGAc                                                 <  1:296436‑M2/200‑1 (MQ=255)
                                                                                                                                                                                                          tttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGGCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTAAGCCCGGACGGATAAAGCGTTTACGCGACAACAGGTGACcgatgcggagctg                                             >  1:4513‑M2/48‑245 (MQ=255)
                                                                                                                                                                                                             attaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGTCGCGCCATATTTTATCTATTTGTTGGGAAAGCAGAGACGGTAATTTTTCCATAAAACGGCCCTGTACGCCTGGACGGATAAAACGTTTACGCGACAACGGACGCCCGGTGCGACGCTGgcg                                          >  2:101515‑M2/45‑251 (MQ=255)
                                                                                                                                                                                                             attaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGgcg                                          >  1:225755‑M2/45‑251 (MQ=255)
                                                                                                                                                                                                                 acccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTCTCTATTTGTTGGGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGGCAACGGATGCCCGATGCGACGCTggcgcggc                                      >  1:216602‑M2/41‑249 (MQ=255)
                                                                                                                                                                                                                  cccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTTTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGACCGATGCGACGCTGGCGCGTCt                                     >  2:414095‑M2/40‑251 (MQ=255)
                                                                                                                                                                                                                  cccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAACAGAGCCGGTAATTTTTCCATATAACGGCCCTGTAAGCCTGGACGGATCAAGCGTTTACGCGACAACGGATACCCGAGGCGACGCAGGCGCGTCt                                     >  1:273643‑M2/40‑251 (MQ=255)
                                                                                                                                                                                                                                tcagattgggtaccgattgcgaaccTTACTGCTGGGGTTCGCTGATTGGTGTGGGCGGATGGTGCTGTTTCCATTCAAGATCCCGGCGGGGATGCTGTCAACCTTTATCGGCGAGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTTATTTTTACATATAAAGGCCCTGTACGCCTGGAAGGATAAAGCGTTTAAGCGAAAACGGATGCCCGATGCGACGCAGGAGCGTGTTATTGGGCctaaa                        >  2:180641‑M2/26‑248 (MQ=255)
                                                                                                                                                                                                                                 cagattgcgtaccgcttgagaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCCACCTTTATCGGCGCGCCATATTTTTTCTATTTGTTTAGAAAGCATCGCCGTTAATTTTTCCATCTAACGGGCCTGTAAGCCTGGCCGGGTAAAGCGTTTACCCGACCACGGATGCCCGATTCGACGCTGGCGCGTCCTTtcgggccgccca                       >  2:442186‑M2/25‑245 (MQ=255)
                                                                                                                                                                                                                                         gtacctcttgcgaatcGTACTGCGGGTGTTCGCTGACTGGGGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCTTACAAACGGCCCCg              <  2:349709‑M2/235‑1 (MQ=255)
                                                                                                                                                                                                                                                  gcgaaccTTTCTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAg      <  1:70508‑M2/243‑1 (MQ=255)
                                                                                                                                                                                                                                                    gacccTTCCTGCGGGTGTGCGCTGGCTGGTGGGGGCGGATGGTGCGGTTTCCATTCCAGATCCCGGCGGGGCTGCGGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTc   <  2:216602‑M2/246‑1 (MQ=255)
                                                                                                                                                                                                                                                     atccTTACTGCGGGTGTTCTCTGTTTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCATCCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCATAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCg  <  1:524877‑M2/247‑1 (MQ=255)
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                      
GCGATTACCCCGCTGTGCCGCCGCTGGCTGACCATTTTACCGCTGGGTGGTGATACCGCCCGAGCCGTAGGAATGGCGCTGACGCCGACGCGAATTGCGCTGCTGCTGTTAGCGGCTTGCCTGACGGCGACCGCGACGATGACTATTGGACCGTTGAGTTTTGTTGGTTTAATGGCACCGCATATTGCGCGGATGATGGGCTTTCGACGGACGATGCCACACATCGTAATTTCGGCGCTGGTGGGTGGTTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCG  >  NC_000913/173067‑173562

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