Predicted mutation
evidence seq id position mutation annotation gene description
JC JC NC_000913 1,431,062 IS2 (–) +5 bp coding (2014‑2018/3363 nt) stfR → putative prophage side tail fiber protein StfR

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 1431062 =0 (0.000)22 (1.250) 20/292 0.0 100% coding (2014/3363 nt) stfR putative prophage side tail fiber protein StfR
?NC_000913 = 1469240 NA (NA)noncoding (1/1331 nt) IS2 repeat region
* ? NC_000913 = 14310660 (0.000)20 (1.140) 17/292 0.1 100% coding (2018/3363 nt) stfR putative prophage side tail fiber protein StfR
?NC_000913 = 1651548 NA (NA)noncoding (706/706 nt) IS2 repeat region

GACTGACAGGGGGCGTTATGCTTATGGCATGTTCTGGCAATATCAAAACAATGAACGAGCCATTCACTCAATAATGATGAGTAATAAGGGCGATGATTTGCGCTCTGTGTTCTATGTTGATGGCGCTGCTTTCCCTGTTTTTGCGTTTATTGAAGATGGCCTGTCAATATCCGCACCTGGTGCTGATCTCGTTGTTAATGATACGACCTATAAGTTTGGGGCAACAAATCCGGCGACTGAATGTATCGCGGCGGACGTTATCCTTGATTTTAAGAGTGGGCGTGGTTTTTATGAGTCTCATTCGTTAATCGTTAACGATAACTTGTCGTGCAAAAAACTTTTTGCCACAGACGAAATTGTAGCGCGTGGTGGTAATCAGATTCGAATGATAGGTGGGGAGTATGGTGCATTATGGCGTAATGATGGCGCTAAAACTTACCTGCTGCTTACC  >  NC_000913/1430833‑1431283
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                         
gACTGCCAGGGGGCGTTATGCTTATGGCATGTTCTGGCAATATCAAAACAATGAACGAGCCATTCACTCAATAATGATGAGTAATAAGGGCGATGATTTGCGCTCTGTGTTCTATGTTGATGGCGCTGCTTTCCCTGTTTTTGCGTTTATTGAAGATGGCCTGTCAATATCCGCACCTGGTGCTGATCTCGTTGTTAATGATACGCCCTATAAGTTTGGGGGAAAAAATCCGGCtggatttgcccctatat                                                                                                                                                                                                          >  1:340789‑M1/1‑234 (MQ=255)
                  accatgttgccgggcaacgagggagaccgtcatccccggttcaaagctctgctgaacaattgcgatcttttcctgtgtggtacgccgtctgcgtttctccggccctaagacatcaatcatctgttctccaatgactagtctaaaaactagtattaagactatcacttatttaagtgatattggttgtctggagattcagggggccagtctaCCGGCGACTGAATGTATCGCGGCGGACGTTATCCTTGAtt                                                                                                                                                                                        <  1:372174‑M1/40‑1 (MQ=255)
                                       ggagaccgtcatccccggttcaaagctcttctgaacaattgcgatcttttcctgtgtggtacgccgtctgcgtttctccggccctaagacatcaatcatctgttctccaatgactagtataaaaactagtattaagactatcacttatttaagtgatattggttgtctggagagtcagggggccattctaCCGGCGACTGAATGTATCGCGTCGGACGTTATCCTTGATTTTAAGAGTGGGCGTGgttttn                                                                                                                                                                   >  1:47525‑M1/191‑250 (MQ=255)
                                       ggagaccgtcatccccggttcaaagctctactgaacaattgggatcttttcctgtgtggtacgccgtctgcgtttctccggccctaagacatcaatcatctgttctccaatgactagtctaaaaactagtattaagactatcacttatttaagtgatattggttgtctggagattcagggggccagtctaCCGGCGACTGAATGTATCGCGGCGGACGTTATCCTTGATtttaagagtgggggataggaa                                                                                                                                                                    <  1:1656‑M1/60‑10 (MQ=255)
                                               aCAATGAACGAGCCATTCACTCAATAATGATGAGTACTAAGGGCGATGATTTGCGCTCTGTGTTCTATGTTGATGGCGCTGCTTTCCCTGTTTTTGCGTTTATTGAAGATGGCCTGTCAATATCCGCACCTGGTGCTGATCTCGTTGTTAATGATACGACCTATAAGTTTGGGGCAACAAATCCGGCtggatttgcccctatatttccagacatctgttatcacttaacccattacaagcccgctgccgc                                                                                                                                                            >  1:360503‑M1/1‑187 (MQ=255)
                                                        gttcaaagctctgctgaaaaattgcgatcttttcctgtgtggtacgccgtctgcgtttctccggccctaagacatcaatcatatgttctccaatgactcgtctaaaaactagtattaagactatcacttatttaagtgatattggttgtctggagattcagggggccagtctaCCGGCGACTGAATGTATCGCGGCGGACGTTATCCTTGCTTTTAAGAGTGGGCGTGGTTTTTATGAGTCTCCTtcgtt                                                                                                                                                   <  1:192579‑M1/77‑1 (MQ=255)
                                                                  ctgctgaacaattgcgatcttttcctgtgtggtacgccgtctgcgtttctccggccctaagacatcaatcatctgttctccaatgactagtctaaaaactagtattaagactatcacttatttaagtgatattggttgtctggagattcagggggccagtctaCCGGCGACTGAATGTATCGCGGCGGACGTTATCCTTGATTTTAag                                                                                                                                                                                   <  1:371935‑M1/45‑1 (MQ=255)
                                                                                cgatcttttcctgtgtggtacgccgtctgcgtttctccggccctaagacatcaatcatctgttctcccatgactagtctaaaaactagtattaagactatcacttatttcagtgatattggttgtctggagattcagggggccagtctcCCGGCGACTGAATGTATCGCGGCGGACGTTATCCTTGATTTTAAGAGTGGGCGTGGTTTTTATGAGTCTCATTCGTTAATCGTTAACGATAACTTGTCGTGc                                                                                                                          <  1:241201‑M1/102‑1 (MQ=255)
                                                                                        tccggtgtggtacgccttctgcgtttctccggctctaagacatcaatcatctgttctccaatgactagtctaaaaactagtattaagactatcacttatttaagggatattggttgtctggagatccagggggccagtctaCCGGCGACTGAATGTATCGCGGCGGACGTTATCCTTGATTTTAAGAGGGGGCGCGGTTTTTATGAGTCTCATTCGTTAATAGTTAACGATAACTTGTCGTGCAAAAAAc                                                                                                                   >  1:239682‑M1/142‑250 (MQ=255)
                                                                                          cGATGATTTGCGCTCTGTGTTCTATGTTGATGCCGCTGCTTTCCCTGTTTTTGCGTTTATTGAAGCTGGCCTGTCAATATCCGCACCTGGGGCTGATCTCGTTGTTAATGATACCACCTATAAGTTTGGGGCAACAAATCCGGCtggatttgcccctatatttccagacatctgttatcacttaacccattacaagcccgctgccgcagatattcccgtggcgagc                                                                                                                                         >  1:97895‑M1/1‑144 (MQ=255)
                                                                                          cGATGATTTGCCCTCTGTGTTCTATGTTGGTGGCACTGCTTTCCCTGTTTTTGCGTTTATTGAAGATGGCCTGTCACTATCCGCACCTGGTGCTGTTCTCGTTGTTAATGATACGACCTATAAGTTTGGGGCAACAAATCCGGCtggatttgcccctatatttccagacatctgttatcacttaacccattacaagcccgctgccgcagatattcccgtggcgagcgataacccagcgcac                                                                                                                          >  1:146159‑M1/1‑144 (MQ=255)
                                                                                                ggtacgccgtctgcgtttctccgcccctaagacctcaatcatctgttctccaatgactagtctaaaaactagtattaagactatcacttatttaagtgatattggttgtctggcgattcagggggccagtctaCCGGCGACTGCATGTATCGCGGCGGACGTTATCCTTg                                                                                                                                                                                           <  1:135051‑M1/37‑1 (MQ=255)
                                                                                                             cgtttctccggccctaagacatcaatcatctgttctccaatgactagtctaaaaactagtattaagactatcacttatttaagtgatattggttgtctggagattcagggggccagtctaCCGGCGACTGAATGTATCGCGGCGGACGTTATCCTTGATTTTAAGAGTGGGCGTGGTTTTTAt                                                                                                                                                                 <  1:220624‑M1/63‑1 (MQ=255)
                                                                                                                ttctccggccctaagacatcaatcatctgttctccaatgactagtctaaaaactagtattaagactatcacttatttaagtgatattggttggctggagatacagggggccaggctaCCGGCGACGGAATGTATCGCGGCGGACGTTATCCTTGATTTTAAGAGTGGGCGTGGTTTTTATTAGTCTCATTCGTTAATCATTAACGATAACTTGTCGTGCAAAAAACTTTTTGCCACAGACGAAATTGTAg                                                                                           >  1:432872‑M1/118‑250 (MQ=255)
                                                                                                                      ggccctaagacatcaatcatctgttctccaatgactagactaacaactagtattaagactatcacatatttaagagatattggttgtctggagattcagggggccagtctaCCGGCGACTGAATGTATCGCGGCGGACGTTATCCTTGATTTTAag                                                                                                                                                                                   <  1:333547‑M1/45‑1 (MQ=255)
                                                                                                                           ngcntctntnncNNTTTNGGCGTTTATTGAAGATGGCCTNNCAATATCCGNACCTGGTGCTGATCTCTTTGGTAATGATACGACCTATAAGTTTGGGGCAACAANTCCGGCtggnnttncccctntnnntnnnganatctgttagcacttatcccnttnnaagcccnctgccgcagatattcccgtggcgagc                                                                                                                                         <  2:97895‑M1/188‑83 (MQ=255)
                                                                                                                           cgcTGCTTTCCCTGTTTTTGCGTTTATTGAAGATGGCCTGTCAATATCCGCTCCTGGTGATGATATAGTTGTTAATGATACGACCTCTAAGTTTGGGGCAACATATCAGGCtggatttgcccctatatttacagacatctgttatcacttaacccattacaagcccgctgcggcagatattcccgtggcgagcgataacccagcgcactatgcggatgccattcgttataatgctcgaacgcctctgca                                                                                 >  1:308042‑M1/1‑111 (MQ=255)
                                                                                                                               gctTTCCCTGTTTTTGCGTTTATTGAAGATGGCCTGTCAATATCCGCACCTGGTGCTGATCTCGTTGTTAATGATACGACCTATAAGTTTGGGGCAACAAATCCGGCtggatttgcccctatatttccagacatctgttatcacttaacccattacaagcccgctgccgcagatagtcccgtggcgagcgataacccagcgcactatgcggatgccattcgttataatgctcgaacgcctctgcaaggttc                                                                           <  1:333160‑M1/251‑145 (MQ=255)
                                                                                                                                    cccTGTTTTTGCGTTTATTGAAGATGGCCTGTCAATATCCGCACCTGGTGCTGATCTCGTTGTTAATGATACGACATATCAGTTTGGGGCAACAAATCCGGCtggatttgaccctatatttccagacatctgttatcacttaaccccttaca                                                                                                                                                                         >  1:318356‑M1/1‑102 (MQ=255)
                                                                                                                                    cccTGTTTTTGCGTTTATTGAAGATGGCCTGTCAATATCCGCACCTGGTGCTGATCTCGGTGTNAATNATACNCCNTATAAGTTTGGGGCANCAAATCCGGCtggatttgcccctatatttccagacatctgttatnncnnnanccattaca                                                                                                                                                                         <  2:318356‑M1/152‑51 (MQ=255)
                                                                                                                                         tttttGCGTTTATTGAAGATGGCCTGTCAATATCCGCACCTGGTGCTGATCTCGTTGTTAATGATACGACCTATAAGTTTGGGGCAACAAATCCGGCtggatttgcccctatatttccagacatctgttatcacttaacccattacaag                                                                                                                                                                       <  1:146285‑M1/149‑53 (MQ=255)
                                                                                                                                           tgttctccaatgactagtctaaaaactagtattaagactatcacttatttaagtgatattggttgtctggagattcagggggccagtctaCCGTCGACTGAATGTATCGCGGCGGACGTGATCCTTGATTTGAAGAGTGGGCGTGGGTTTTATGAGTATCATTCGTTAATCGTTAACGATAACTTGTCGTGCAAAAAACTTTTTGCCACAGAc                                                                                                     >  1:404483‑M1/91‑213 (MQ=255)
                                                                                                                                               cGTTTATTGAAGATGGCCTGTCAATATCCGCACCTGGTGCTGATCTCGTTGTTAATGATACGACCTATAAGTTTGGGTCAACAAATCCGGCtggatttgcccctatatttccagacatctgttatcacttaacccattacaagcccgctgccgcagatattcccgtggcgatcgataacccagcgcactatgcggatgccattcgttataatgctcgaacgcctctgcaaggttctttgctgccgttaacc                                                           <  1:284702‑M1/251‑161 (MQ=255)
                                                                                                                                                    aTTGAAGATGGCCTGTCAATATCCGCACCTGGTGCTGATCTCGTTGTTAATGATACGACCTATAAGTTTGGGGCAACAAATCCGGCtggatttgcccctatatttccagacatctgttatcacttaacccagtacaagcccgctgccgcagatattcccgtggcgagcgataacccagcgcactatgcggatgccattcgttataatgctcgaacgcctctgcaaggttctgtgctgccgttaacccgtct                                                      <  1:147228‑M1/251‑166 (MQ=255)
                                                                                                                                                      nactagtctaaaaactagtattaagactatcacttatttaagtgatattggttgtctggagattcagggggccagtctaCCGGCGACTGAATGTATCGCGGCGGACGTTATCCTTGATTTTAagag                                                                                                                                                                                 <  1:89240‑M1/47‑1 (MQ=255)
                                                                                                                                                      gactagtctaaaaactagtattaagactatcacttatttaagtgatattggttgtctggagattcatggggccaggctaCCGGCGACTGCATGTATCGCGGCGGACGTTATCCTTGAttt                                                                                                                                                                                       >  1:29721‑M1/80‑120 (MQ=255)
                                                                                                                                                      gactagtctaaaaactagtattaagacgatcacttatttaagtgatattggttgtcttgagagtcatggggccgtgctaCCGGCGACTGAATGTATCGCGGCGGACGTTATCCTTGATTTTAAGAGTGGGCGTGGTTTTTATGa                                                                                                                                                               >  1:215475‑M1/80‑144 (MQ=255)
                                                                                                                                                          agtctaaaaactagtattaagactatcacttatttaagtgatattggttgtctggagattcagggggccagtctaCCGGCGACTGAATGTATCGCGGCGGACGTTATCCTTGATTTTAAGAGTGGGCGTGGTTTTTATGAGTCTCa                                                                                                                                                         <  1:184796‑M1/71‑1 (MQ=255)
                                                                                                                                                                aaaactagtattaagactatcacttatttaagtgatattggtggtctggagattcagggggccagtctaCCGGCGACTGAAGGTATCGCGGCGGACGTTATCCTTGATTTTAAGAGTGGGc                                                                                                                                                                            >  1:392486‑M1/70‑121 (MQ=255)
                                                                                                                                                                    cAATATCCGCACCTGGTGCTGATCTCGTTGTTAATGATACGACCTATAAGATTGGGGCAACACATACGGCtggatatgcccctatatttccagccatctgttatcacttaacccattacaagcccgctgccgccgctattcccgtggcgagcgataacccagcgcactatgcggatgccattcgttataatgctcgaccgcctctgcaaggttctttgctgccgttaacccgtctggtttgggcatgatac                                      >  1:382140‑M1/1‑70 (MQ=255)
                                                                                                                                                                      agtattaagactatcacttatttaagtgatattggttgtctggagattcagggggccagtctaCCGGCGACTGAATGTATCGCGGCGGACGTTATCCTTGAGTTTAAGAGTGGGCGTGGTTTTTATGAGTCTCATTCTTTAATCGGTAACGATAACTTGTCGTGCAAAAAACTTTTTGCCACAAACGAAATTGTATCGCGTGGTGGTAATCAGATTCGAATGATAGGTGGGGAGTATGGTGCATTATGGc                                     >  1:62835‑M1/64‑250 (MQ=255)
                                                                                                                                                                                    tGCTGATCTCGTTGTTAATGATACGACCTATAAGTTTGGGGCAACAAATCCGGAtggatttgcccctatatttccagacatctgctatcacttaacccattacaagcccgctgccgcagatattcacgtggcgagcgataacccagcgcacaatgcggatgccattcgttataatgctcgaacg                                                                                         >  1:245923‑M1/1‑54 (MQ=255)
                                                                                                                                                                                      ntgATCTCGTTGTTAATGATACGACCTATAAGTTTGGGGCAACAAATCCGGCtggatttgcccctatatttccagacatctgttatcacttaacccattacaagccc                                                                                                                                                                    >  1:265062‑M1/2‑52 (MQ=255)
                                                                                                                                                                                      cTGATCTCGTTGTTAATGATACGACCTATAAGTTTGGGGCAACAAATCCGGCtggatttgcccctatatttccagacatctgttatcacttaacccattacaagccc                                                                                                                                                                    <  2:265062‑M1/107‑56 (MQ=255)
                                                                                                                                                                                        natCTCGTTGTTAATGATACGACCTATAAGTTTGGGGCAACAAATCCGGAtggatttgcccctatatttccagacatctgttatcacttaacccattacaagcccgctgccgcagatattcccgtggcgagcgataacccagcgcactatgcggatgccattcgttataatgctcgaacgcctctgcaaggttctttgctgccgttaacccgtctggtttgggcatgatactgatgtagtcacgctttatc                  >  1:90950‑M1/2‑50 (MQ=255)
                                                                                                                                                                                                 gttAATGATACGACCTATAAGTTTGGGGCAACAAATCCGGCtggatttgcccctatatttccagacatctgttatcacttaacccattacaagcccgctgccgcagatattcccgtggcgagcgataacccag                                                                                                                               <  1:153955‑M1/133‑93 (MQ=255)
                                                                                                                                                                                                                  gattcagggggccagtctaCCGGCGACTGAATGTATCGCGGCGGACGTTATCCTTGATTTTCAGAGTGGGCGTGGTTTTTATGAGTCTCATTCGTTAATCGTTAACGCTAACTTGTCGTGCAAAAAACTTTTTGCCACAGACGAAATTGTAGCGCGTGGAGGTAATCAGATTCGAATGATAGGTGGGGAGTATGGTGCATTATGGCGTAat                                <  1:394418‑M1/192‑1 (MQ=255)
                                                                                                                                                                                                                    aGTTTGGGGCAACAAATCCGGCtggatttgcccctatatttccagacatctgttatcacttaacccattacaagcccgctgccgcagatatgcccgtggcgagcgataacccagcgcacgatgcggatgccaatcgttagaatcctcgaacgcccctgcaaggggctttgctgccgttaacccgtctggtttgggcatgat                                        <  1:197812‑M1/201‑180 (MQ=255)
                                                                                                                                                                                                                          gggCAACAAATCCGGCtggatttgcccctatatttccagacatctgttatcacttaacccattacaagcccgctgccgcagctattcccgtggcgagcgataacccagcgcactatgcggatgccattcgttataa                                                                                                   >  1:173412‑M1/1‑16 (MQ=255)
                                                                                                                                                                                                                               agtctaCCGGCGACTGAATGTATCGCGGCGGACGTTATCCGTGATTTTAAGAGTGGGCGTGGTTTTTATGAGTCTCATTCGTTAATCGTTAACGATAACTTGTCGTGCAAAAAACTTTTTGCCACAGACGAAATTGTAGCGCGTGGTGGTAAGCAGATTCGAATGATAGGTGGGGAGTATGGTGCATTATGGTGTAATGATGGCGCTAAAACTTACCTGCTGCTTAcc  >  1:159616‑M1/7‑228 (MQ=255)
                                                                                                                                                                                                                                gtctaCCGGCGACTGAATGTATCGCGGCGGACGTTATCCTTGATTTTAAGAGTTGGCTGGGTTTTTATGAGTCTCATTCGTGAATCGTTAACGATAACTTGTTGTGCAATAACCGTTTTGCCACAGACGAAATTGTAGCGCGTGGTGGTAATCAGATTCGAATGATa                                                              >  1:196792‑M1/6‑167 (MQ=255)
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                         
GACTGACAGGGGGCGTTATGCTTATGGCATGTTCTGGCAATATCAAAACAATGAACGAGCCATTCACTCAATAATGATGAGTAATAAGGGCGATGATTTGCGCTCTGTGTTCTATGTTGATGGCGCTGCTTTCCCTGTTTTTGCGTTTATTGAAGATGGCCTGTCAATATCCGCACCTGGTGCTGATCTCGTTGTTAATGATACGACCTATAAGTTTGGGGCAACAAATCCGGCGACTGAATGTATCGCGGCGGACGTTATCCTTGATTTTAAGAGTGGGCGTGGTTTTTATGAGTCTCATTCGTTAATCGTTAACGATAACTTGTCGTGCAAAAAACTTTTTGCCACAGACGAAATTGTAGCGCGTGGTGGTAATCAGATTCGAATGATAGGTGGGGAGTATGGTGCATTATGGCGTAATGATGGCGCTAAAACTTACCTGCTGCTTACC  >  NC_000913/1430833‑1431283

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