Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 2,765,456 (TATGGCAC)6→5 intergenic (‑303/+50) yfjL ← / ← yfjM CP4‑57 putative defective prophage, DUF4297/DUF1837 polymorphic toxin family protein/CP4‑57 prophage; uncharacterized protein

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 2765412 2765419 8 50 [0] [0] 49 yfjL/yfjM CP4‑57 putative defective prophage, DUF4297/DUF1837 polymorphic toxin family protein/CP4‑57 prophage; uncharacterized protein

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 2765420 =0 (0.000)35 (0.810) 33/210 0.5 100% intergenic (‑267/+93) yfjL/yfjM CP4‑57 putative defective prophage, DUF4297/DUF1837 polymorphic toxin family protein/CP4‑57 prophage; uncharacterized protein
?NC_000913 = 2765411 0 (0.000)intergenic (‑258/+102) yfjL/yfjM CP4‑57 putative defective prophage, DUF4297/DUF1837 polymorphic toxin family protein/CP4‑57 prophage; uncharacterized protein

TTTTCAGGATTCCACAAATATTGTAACCCTGACGTCTAATTTGGTTATATAACATACTGCCCCAACTAACGTAAGGTAGGCACAGCCTTAATCAAGGTTTTGGTTTTCAAATAGTTCAGAGCTGAGACCGCCAACCGCACTATGGCACTATGGCACTATGGCACTATGGCACTATGGCACTATGGCACCGCCAGCGATTGTAATCACACTTGATATTATAAAACACAGTTGCACGCATTATTTCCTGGTTGGTAGGGTCATATCTCGATGCTCTTTGAGCAATGTCAACATCGCG  >  NC_000913/2765276‑2765570
                                                                                                                                                |                                                                                                                                                      
tatcaagcttttcaggattccacaaatattgtaaccctgacgtctaatttggttatataacatactgccccaactaacgtaaggtaggcacagccttaatcaaggttttggttttcaaatagttcagagctgagaccgccaaccgcactat                                                                                                                                                  <  1:343054‑M1/7‑1 (MQ=255)
    aagcttttcaggattccacaaatattgtaaccctgacgtctaatttggttatataacatactgccccaactaacgtaaggtaggcacagccttaatcaaggttttggttttcaaatagttcagagctgagaccgccaaccgcactatggca                                                                                                                                              >  2:20616‑M1/141‑151 (MQ=255)
    aagcttttcaggattccacaaatattgtaaccctgacgtctaatttggttatataacatactgccccaactaacgtaaggtaggcacagccttaatcaaggttttggttttcaaatagttcagagctgagaccgccaaccgcactatggca                                                                                                                                              >  1:245153‑M1/141‑151 (MQ=255)
    aagcttttcaggattccacaaatattgtaaccctgacgtctaatttggttatataacatactgccccaactaacgtaaggtaggcacagccttaatcaaggttttggttttcaaatagttcagagctgagaccgccaaccgcactatgg                                                                                                                                                <  2:781025‑M1/9‑1 (MQ=255)
     agcttttcaggattccacaaatattgtaaccctgacgtctaatttggttatataacatactgccccaactaacgtaaggtaggcacagccttaatcaaggttttggttttcaaatagttcagagctgagaccgccaaccgcactatggcac                                                                                                                                             >  1:14011‑M1/140‑151 (MQ=255)
         tttcaggattccacaaatattgtaaccctgacgtctaatttggttatataacatactgccccaactaacgtaaggtaggcacagccttaatcaaggttttggttttcaaatagttcagagctgagaccgccaaccgcgctatggcacta                                                                                                                                           >  2:767864‑M1/136‑149 (MQ=255)
          ttcaggattccacaaatattgtaaccctgacgtctaatttggttatataacatactgccccaactaacgtaaggtaggcacagccttaatcaaggttttggttttcaaatagttcagagctgagaccgccaaccgcactatggcactatgg                                                                                                                                        <  1:133779‑M1/17‑1 (MQ=255)
            caggattccacaaatattgtaaccctgacgtctaatttggttatataacatactgccccaactaacgtaaggtaggcacagccttaatcaaggttttggttttcaaatagttcagagctgagaccgccaaccgcactatggcactatggc                                                                                                                                       >  2:96065‑M1/133‑150 (MQ=255)
                   ccacaaatattgtaaccctgacgtctaatttggttatataacatactgccccaactaacgtaaggtaggcacagccttaatcaaggttttggttttcaaatagttcagagctgagaccgccaaccgcactatggcactatggcactatggc                                                                                                                               <  2:730016‑M1/26‑1 (MQ=255)
                       aaatattgtaaccctgacgtctaatttggttatataacatactgccccaactaacgtaaggtaggcacagccttaatcaaggttttggttttcaaatagttcagagctgagaccgccaaccgcactatggcactatggcactatggcact                                                                                                                            <  1:417707‑M1/29‑1 (MQ=255)
                                  ccctgacgtctaatttggttatataacatactgccccaactaacgtaaggtaggcacagccttaatcaaggttttggttttcaaatagttcagagctgagaccgccaaccGcactatggcactatggcactatggcactatggcactatgg                                                                                                                >  2:492987‑M1/111‑151 (MQ=255)
                                    ctgacgtctaatttggttatataacatactgccccaactaacgtaaggtaggcacagccttaatcaaggttttggttttcaaatagttcagagctgagaccgccaaccGCactatggcactatggcactatggcactatggcactatggc                                                                                                               <  2:710947‑M1/42‑1 (MQ=255)
                                     tgacgtctaatttggttatataacatactgccccaactaacgtaaggtaggcacagccttaatcaaggttttggttttcaaatagttcagagctgagaccgccaaccGCACtatggcactatggcactatggcactatggcactatggcac                                                                                                             <  1:208210‑M1/44‑1 (MQ=255)
                                          tctaatttggttatataacatactgccccaactaacgtaaggtaggcacagccttaatcaaggttttggttttcaaatagttcagagctgagaccgccaaccGCACTATGGCACTATGGCACTATGGCACTATGGCACTATGGCACCGcc                                                                                                         <  2:674955‑M1/48‑1 (MQ=255)
                                                ttggttatataacatactgccccaactaacgtaaggtaggcacagccttaatcaaggttttggttttcaaatagttcagagctgagaccgccaaccGCACTATGGCACTATGGCACTATGGCACTATGGCACTATGGCACCGCCAGCGAtt                                                                                                  <  1:806116‑M1/55‑1 (MQ=255)
                                                ttggttatataacatactgccccaactaacgtaaggtaggcacagccttaatcaaggttttggttttcaaatagttcagagctgagaccgccaaccGCACTATGGCACTATGGCACTATGGCACTATGGCACTATGGCACCGCCAGCGAtt                                                                                                  <  1:927910‑M1/55‑1 (MQ=255)
                                                ttggttatataacatactgccccaactaacgtaaggtaggcacagccttaatcaaggttttggttttcaaatagttcagagctgagaccgccaaccGCACTATGGCACTATGGCACTATGGCACTATGGCACTATGGCACCGCCAGCGAt                                                                                                   >  2:896960‑M1/97‑150 (MQ=255)
                                                  ggttatataacatactgccccaactaacgtaaggtaggcacagccttaatcaaggttttggttttcaaatagttcagagctgagaccgccaaccGCACTATGGCACTATGGCACTATGGCACTATGGCACTATGGCACCGCCAGCGATTGt                                                                                                >  1:253360‑M1/95‑151 (MQ=255)
                                                        ataacatactgccccaactaacgtaaggtaggcacagccttaatcaaggttttggttttcaaatagttcagagctgagaccgccaaccGCACTATGGCACTATGGCACTATGGCACTATGGCACTATGGCACCGCCAGCGATTGTAATca                                                                                           >  1:401683‑M1/89‑150 (MQ=255)
                                                          aacatactgccccaactaacgtaaggtaggcacagccttaatcaaggttttggttttcaaatagttcagagctgagaccgccaaccGCACTATGGCACTATGGCACTATGGCACTATGGCACTATGGCACCGCCAGCGATTGTAATCacac                                                                                        >  1:312707‑M1/87‑151 (MQ=255)
                                                            catactgccccaactaacgtaaggtaggcacagccttaatcaaggttttggttttcaaatagttcagagctgagaccgccaaccGCACTATGGCACTATGGCACTATGGCACTATGGCACTATGGCACCGCCAGCGATTGTAATCACACtt                                                                                      >  1:335819‑M1/85‑151 (MQ=255)
                                                                ctgccccaactaacgtaaggtaggcacagccttaatcaaggctttggttttcaaatagttcagagctgagaccgccaaccGCACTATGGCACTATGGCACTATGGCACTATGGCACTATGGCACCGCCAGCGATTGTAATCACACTTGat                                                                                   <  2:665076‑M1/70‑1 (MQ=255)
                                                                     ccaactaacgtacggtaggcacagccttaatcaaggttttggttttcaaatagttcagagctgagaccgccaaccGCACTATGGCACTATGGCACTATGGCACTATGGCACTATGGCACCGCCAGCGATTGTAATCACACTTGAtatta                                                                               <  2:796311‑M1/74‑1 (MQ=255)
                                                                        actaacgtaaggtaggcacagccttaatcaaggttttggttttcaaatagttcagagctgagaccgccaaccGCACTATGGCACTATGGCACTATGGCACTATGGCACTATGGCACCGCCAGCGATTGTAATCACACTTGATATTATaaaa                                                                          >  1:880520‑M1/73‑151 (MQ=255)
                                                                           aacgtaaggtaggcacagccttaatcaaggttttggttttcaaatagttcagagctgagaccgccaaccGCACTATGGCACTATGGCACTATGGCACTATGGCACTATGGCACCGCCAGCGATTGTAATCACACTTGATATTATAAAac                                                                         >  1:179439‑M1/70‑149 (MQ=255)
                                                                            acgtaaggtaggcacagccttaatcaaggttttggttttcaaatagttcagagctgagaccgccaaccGCACTATGGCACTATGGCACTATGGCACTATGGCACTATGGCACCGCCAGCGATTGTAATCACACTTGATATTATAAAacac                                                                       <  1:783012‑M1/82‑1 (MQ=255)
                                                                             cgtaaggtaggcacagccttaatcaaggttttggttttcaaatagttcagagctgagaccgccaaccGCACTATGGCACTATGGCACTATGGCACTATGGCACTATGGCACCGCCAGCGATTGTAATCACACTTGATATTATAAAacac                                                                       <  2:375943‑M1/82‑1 (MQ=255)
                                                                                aaggtaggcccagccttaatcaaggttttggttttcaaatagttcagagctgagaccgccaaccGCACTATGGCACTATGGCACTATGGCACTATGGCACTATGGCACCGCCAGCGATTGTAATCACACTTGATATTATAAAACACCGTTg                                                                  <  2:169856‑M1/87‑1 (MQ=255)
                                                                                 aggtaggcacagccttaatcaaggttttggttttcaaatatttcagagctgagaccgccaaccGCTCTATGGCACTATGGCACTATGTCACTATGGAACTATGGCACCGCCAGCGATTGTAATCACACTTGATATTATAAACCACAGtt                                                                   >  2:495465‑M1/64‑149 (MQ=255)
                                                                                      gacacagccttaatcaaggttttggttttcaaatagttcagagctgagaccgccaaccGCATTATGGCACTATGGCACTATGGCACTATGGCACTATGGCACCGCCAGCGATTGTAATCACACTTGATATTATAAAACACAGTTGCACGCa                                                            >  2:759490‑M1/59‑151 (MQ=255)
                                                                                           agccttaatcaaggttttggttttcaaatagttcagagctgagaccgccaaccGCACTATGGCACTATGGCACTATGGCACTATGGCACTATGGCACCGCCAGCGATTGTAATCACACTTGATATTATAAAACACAGTTGCACGCattatt                                                       <  2:592422‑M1/98‑1 (MQ=255)
                                                                                             ccttaatcaaggttttggttttcaaatagttcagagctgagaccgccaaccAGATTATGGCACTATAGCACTATGGCATTATGGCACTATGGCACCGCCAGCGATTGTAATCACACTTGATATTATAAAACACAGTTGCACGCATTAttt                                                      >  2:512099‑M1/52‑150 (MQ=255)
                                                                                                    caaggttttggttttcaaatagttcagagctgagaccgccaaccGCACTATGGCACTATGGCACTATGGCACTATGGCACTATGGCACCGCCAGCGATTGTAATCACACTTGATATTATAAAACACAGTTGCACGCATTATTTCCtggttg                                              >  2:13120‑M1/45‑151 (MQ=255)
                                                                                                    caaggttttggttttcaaatagttcagagctgagaccgccaaccGCACTATGGCACTATGGCACTATGGCACTATGGCACTATGGCACCGCCAGCGATTGTAATCACACTTGATATTATAAAACACAGTTGCACGCATTATTTCCtggtt                                               >  1:207576‑M1/45‑150 (MQ=255)
                                                                                                        gttttggttttcaaatagttcagagctgagaccgccaaccGCACTATGGCACTATGGCACTATGGCACTATGGCACTATGGCACCGCCAGCGATTGTAATCACACGTGATATTATAAAACACAGTTGCACGCATTATTTCCTGGTTGGTAg                                          >  2:86053‑M1/41‑151 (MQ=255)
                                                                                                           ttggttttcaaatagttcagagctgagaccgccaaccGCACTATGGCACTATGGCACTATGGCACTATGGCACTATGGCACCGCCAGCGATTGTAATCACACTTGATATTATAAAACACAGTTGCACGCATTATTTCCTGGTTGGTAGGGt                                       <  2:336332‑M1/114‑1 (MQ=255)
                                                                                                              gttttcaaatagttcagagctgagaccgccaaccGCACTATGGCACTATGGCACTATGGCACTATGGCACTATGGCACCGCCAGCGATTGTAATCACACTTGATATTATAAAACACAGTTGCACGCATTATTTCCTGGTTGGTAGGGTCat                                    <  1:879998‑M1/117‑1 (MQ=255)
                                                                                                                    aaatagttcagagctgagaccgccaaccGCACTATGGCACTATGGCACTATGGCACTATGGCACTATGGCACCGCCAGCGATTGTAATCACACTTGATATTATAAAACACAGTTGCACGCATTATTTCCTGGTTGGTAGGGTCATATCTCg                              >  1:231999‑M1/29‑151 (MQ=255)
                                                                                                                        agttcagagctgagcccgccaaccGCACTATGGCACTATGGCACTATGGCACTATGGCACTATGGCACCGCCAGCGATTGTAATCACACTTGATATTATAAAACACAGTTGCACGCATTATTTCCTGGTTGGTAGGGTCATATCTCGATg                           >  2:821498‑M1/25‑150 (MQ=255)
                                                                                                                        agttcagagctgagaccgccaaccGCACTATGGCACTATGGCACTATGGCACTATGGCACTATGGCACCGCCAGCGATTGTAATCACACTTGATATTATAAAACACAGTTGCACGCATTATTTCCTGGTTGGTAGGGTCATATCTCGATg                           >  2:847194‑M1/25‑150 (MQ=255)
                                                                                                                          ttcagagctgagaccgccaaccGCACTATGGCACTATGGCACTATGGCACTATGGCACTATGGCACCGCCAGCGATTGTAATCACACTTGATATTATAAAACACAGTTGCACGCATTATTTCCTGGTTGGTAGGGTCATATCTCGATGct                         <  2:344‑M1/128‑1 (MQ=255)
                                                                                                                                gctgagaccgccaaccGCACTATGGCACTATGGCACTATGGCACTATGGCACTATGGCACCGCCAGCGATTGTAATCACACTTGATATTATAAAACACAGTTGCACGCATTATTTCCTGGTTGGTAGGGTCATATCTCGATGCTCTTTGa                   >  2:482894‑M1/17‑150 (MQ=255)
                                                                                                                                   gagaccgccaaccGCACTATGGCACTATGGCACTATGGCACTATGGCACTATGGCACCGCCAGCGATTGTAATCACACTTGATATTATAAAACACAGTTGCACGCATTATTTCCTGGTTGGTAGGGTCATATCTCGATGCTCTTTGAGCa                <  2:83870‑M1/137‑1 (MQ=255)
                                                                                                                                     gaccgccaaccGCACTATGGCACTATGGCACTATGGCACTATGGCACTATGGCACCGCCAGCGATTGTAATCACACTTGATATTATAAAACACAGTTGCACGCATTATTTCCTGGTTGGTAGGGTCATATCTCGATGCTCTTTGAGCAATg             >  2:300247‑M1/12‑151 (MQ=255)
                                                                                                                                      accgccaaccGCACTATGGCACTATGGCACTATGGCACTATGGCACTATGGCACCGCCAGCGATTGTAATCACACTTGATATTATAAAACACAGTTGCACGCATTATTTCCTGGTTGGTAGGGTCATATCTCGATGCTCTTTGAGCAATGt            <  2:810312‑M1/141‑1 (MQ=255)
                                                                                                                                       ccgccaaccGCACTATGGCACTATGGCACTATGGCACTATGGCACTATGGCACCGCTAGCGATTGTAATCACACTTGATATTATAAAACACAGTTGCACGCATTATTTCCTGGTTGGTAGGGTCATATCTCGATGCTCTTTGAGCAATGTc           >  2:227930‑M1/10‑151 (MQ=255)
                                                                                                                                          ccaaccGCACTATGGCACTATGGCACTATGGCACTATGGCACTATGGCACCGCCAGCGATTGTAATCACACTTGATATTATAAAACACAGTTGCACGCATTATTTCCTGGTTGGTAGGGTCATATCTCGATGCTCTTTGAGCAATGTCAAc        >  1:711355‑M1/7‑151 (MQ=255)
                                                                                                                                              ccGCACTATGGCACTATGGCACTATGGCACTATGGCACTATGGCACCGCCAGCGATTGTAATCACACTTGATATTATAAAACACAGTTGCACGCATTATTTCCTGGTTGGTAGGGTCATATCTCGATGCTCTTTGAGCAATGTCAACATcg    >  2:903544‑M1/3‑151 (MQ=255)
                                                                                                                                                GCACTATGGCACTATGGCACTATGGCACTATGGCACTATGGCACCGCCAGCGATTGTAATCACACTTGATATTATAAAACACAGTTGCACGCATTATTTCCTGGTTGGTAGGGTCATATCTCGATGCTCTTTGAGCAATGTCAACATcgcg  >  1:452608‑M1/1‑151 (MQ=255)
                                                                                                                                                |                                                                                                                                                      
TTTTCAGGATTCCACAAATATTGTAACCCTGACGTCTAATTTGGTTATATAACATACTGCCCCAACTAACGTAAGGTAGGCACAGCCTTAATCAAGGTTTTGGTTTTCAAATAGTTCAGAGCTGAGACCGCCAACCGCACTATGGCACTATGGCACTATGGCACTATGGCACTATGGCACTATGGCACCGCCAGCGATTGTAATCACACTTGATATTATAAAACACAGTTGCACGCATTATTTCCTGGTTGGTAGGGTCATATCTCGATGCTCTTTGAGCAATGTCAACATCGCG  >  NC_000913/2765276‑2765570

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