Predicted mutation
evidence seq id position mutation annotation gene description
JC JC NC_000913 1,105,376 IS2 (–) +5 bp intergenic (+83/‑34) csgC → / → ymdA curlin chaperone/uncharacterized protein YmdA

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 381260 =NA (NA)12 (0.960) 8/308 0.3 100% noncoding (1/1331 nt) IS2 repeat region
?NC_000913 1105376 = 0 (0.000)intergenic (+83/‑38) csgC/ymdA curlin chaperone/uncharacterized protein YmdA
* ? NC_000913 = 11053801 (0.080)12 (0.950) 12/310 0.1 92.1% intergenic (+87/‑34) csgC/ymdA curlin chaperone/uncharacterized protein YmdA
?NC_000913 1467910 = NA (NA)noncoding (1331/1331 nt) IS2 repeat region

AACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTCATTTATCACAACAATGGCCGCCCTCTTCAGAAAAGTCTTAATTTGTTGAAATATCGAGCATAAGATGAATCTGGAGAGAATGGTCTGCTGCGAATCAGCCAACCTGAAAGTATGGATAACACAACCCTCAAGGATGACTAATCATTGAGGAAATAGAATAAATGTTCAGACCTTTTTTAAACTCTCTTATGCTCGGCAGTTTGTTTTTTCCTTTTATTGCCATTGCTGGAAGCACCGTGCAAGGGGGCGTGATCCATTTTTATGGCCAAATTGTGGAACCGGCATGTGACGTCAGCACCCAGTCATCACCCGTAGAAATGAACTGCCCACAGAATGGTTCTATTCCGGGCAAAACCTACTCCAGCA  >  NC_000913/1105163‑1105618
                                                                                                                                                                                                                         |                                                                                                                                                                                                                                              
ncaccatgttgccgggcaacgagggagaccgtcatccccggttcaaagctctgctgaacaattgcgatcttttcctgtgtggtacgccgtctgcgtttctccggccctaagacatcaatcatctgttctccaatgactagtctaaaaactagtattaagactatcacttatttaagtgatattggttgtctggagattcagggggccagtctaACCCTCAAGGATGACTAATCATTGAGGAAATAGAATa                                                                                                                                                                                                                >  1:174576‑M2/214‑250 (MQ=255)
                                   ccccggttcaaagctctgctgaacaattgcgatcttttcctgtgtggtaagccgtctgcgtgtctccggccttaagacatcaataatctgttctccaatgactagtctaaaaactagtattaagactatcacttatttaagtgatattggttgtctggagattcagggggccagtctaACCCTCAAGGATGACTAATCATTGAGGAAATAGAATAAATGTTCAGACCTTTTTTAAACTCTCTTATGCTCgg                                                                                                                                                                            >  1:191457‑M2/179‑251 (MQ=255)
                                     cccGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTCATTTATCACAACAATGGCCGCCCTCTTCAGAAAAGTCTTAATTTGTTGAAATATCGAGCATAAGATGAATCTGGAGAGAATGGTCTGCTGCGAATCAGCCAACCTGAAAGTATGGATAACACAACCCTtggatttgcccctatatttccagacatctgttatcacttaacccattacaagcccgctgccgcagatatt                                                                                                                                                                          >  1:92103‑M1/1‑181 (MQ=255)
                                              tCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGACACCTCATTTATCACAAAAATGGCAGCCCTCTCCAGAAAAGTCTTAATTTGTTGAAATATCGAGCATAAGATGAATCTGGAGAGAATGGTCTGCTGCGAATAAGCCAACCTGAAAGTATGGATAACACAACCCTtggatttgcccctatatttccagacatctgttatcacttaacccattacaaacccgctgccgcagatattaccgtggcg                                                                                                                                                                 >  1:259057‑M1/1‑172 (MQ=255)
                                                  ctgctgaacaattgcgatcgtttcctgtgtggtacgccgtctgcgtttctccggccctaagacatcaatcatctgttctccaatgactagtctaaaaactagtattaagactatcacttatttaagttatattggttgtctggagattcagggggccagtctaACCCTCAAGGATGACTAATCATTGAGCAAATAGAATAAATGTTCAGATCTTTTTTAAACTCTCGTAGGCTCGGCAGTTTGtttttt                                                                                                                                                               <  1:122815‑M2/86‑1 (MQ=255)
                                                       gaacccttgcgctctttttctgtgtggtacgccgtttgcgtttctccggccctaagccatcaatcatctgttctccaatgactagtctaaaaactagtattaagactatcacttatttaagtgatattggttgtctggagattcagggggccagtctaACCCTCAAGGATGACTAATCATTGCGGAAATAGAATAAATGTTCAGACCTTTTTTAAACTCTCTTATGCTCGGCAGTTTGTTTTTTCcttttt                                                                                                                                                        >  1:100740‑M2/159‑250 (MQ=255)
                                                                           aTGGACAGTCACTTCATTTATCACAACAATGGCCGCCCTCTTCAGAAAAGTCTTAATTTGTTGAAATATCGAGCATAAGATGAATCTGGAGAGAATGGTCTGCTGCGAATCAGCCAACCTGAAAGTATGGATAACACAACCCTtggatttgcccctatatttccagacatctgttagcacgtaacccattgcaatcccgctgccgcagatattcccgtggcgagcgataacccagcgcactatgcggatgc                                                                                                                                    <  1:278994‑M1/251‑109 (MQ=255)
                                                                              gACAGTCACTTCATTTATCACAACAATGGCCGCCCTCTTCAGAAAAGTCTTAATTTGTTGAAATATCGAGCATAAGATGAATCTGGAGAGAATGGTCTGCTGCGAATCAGCCAACCTGAAAGTATGGATAACACAACCCTtggatttgcccctatatttccagacatctgttatcacttaacc                                                                                                                                                                                                     >  1:43564‑M1/1‑140 (MQ=255)
                                                                                             cgtttctccggccctaagacatcaatcatctgttctccaatgactagtctaaaaactagtattaagactatcacttatttaagtgatattggttgtctggagattcagggggccagtctaACCCTCAAGGATGACTAATCATTGAGGAAATAGAATAAATGTTc                                                                                                                                                                                                         >  1:8745‑M2/121‑164 (MQ=255)
                                                                                                         cctaagacatcaatcatctgttctccaatgactagtctaaaaactagtattaagactatcacttatttaagtgatattggttgtctggagattcagggggccagtctaACCCTCAAGGATGACTAATCATTGAGGAAATAGAATAAATGTTCAGACCTTTTTTAAACTCTCTTATGCTCGGCAGTTTGTTTTTTCCTTTTATTGCCATTGCTGGAAGCACCGTGCAAGGGGGCGTGATCCATTTTTATGGc                                                                                                      <  1:4789‑M2/143‑1 (MQ=255)
                                                                                                               ccTCTTCAGAAAAGTCTTAATTTGTTGAAATATCGAGCATAAGATGAATCTGGAGAGAATGGTCTGCTGCGAATCGGCCAACCTGAAAGTATGGATAACACAACCCTtggatttgcccctatatttccagacatctgttatcacttaatccattacaagcctgctgccgcagatattcccgtggcgagcgataacccagcgcactatgcggatgtcattcgttataatgctcgaacgcctctgcaaggttc                                                                                                >  1:285514‑M1/1‑107 (MQ=255)
                                                                                                                      aGAAAAGTCTTAATTTGTTGAAATATCGAGCATAAGATGAATCTGGAGAGAATGGTCTGCTGCGAATCAGCCAACCTGAAAGTATGGATAACACAACCCTtggatttgcccctatatttccagacatctgttatcacttaacccattacaagcccgctgccgcagatagtcccgtggcgagcgataacccagcgcactatgcggatgccattcgttataatgctcgaacgcctctgcaaggttctttgctg                                                                                         <  1:38893‑M1/251‑152 (MQ=255)
                                                                                                                                           gtctaaaaactagtattaagactatcacttatttaagtgatattggttgtctggagattcagggggccagtctaACCCTCAAGGATGACTAATCATTGAGGAAATAGAATAAATGTTCAGACCTTTTTTAAACTCTCTTATGCTCGGCAGTTTGTTTTTTcc                                                                                                                                                             <  1:265373‑M2/88‑1 (MQ=255)
                                                                                                                                             ctaaaaactagtattaagactatcacttatttaagtgatattggttgtctggagattcagggggccagtctaACCCTCAAGGATGACTAATCATTGAGGAAATAGAATAAATGTTCAGACCTTTTTTAAACTCTCTTATGCTCGGCAGTTTGTTTTTTCCTTTTATTGCCATTGCTGGAAGCACCGTGCAAGGGGGCGTGATCCATTTTTATGGCCAAATTGTGGAACCGGCATGTGACGTCAGCAccc                                                                    >  1:117138‑M2/73‑249 (MQ=255)
                                                                                                                                                       nnnnnnagnnnnncacntatttaagtgatattggttgtctggagattcagggggccagtctaACCCTCAAGGATGACTAATCATTGAGGAAATAGAATAAATGTTCAGACCTTTTTTAAACTCTCTTATGCTCGGCAGTTTGTTTTTTcc                                                                                                                                                             <  2:174576‑M2/88‑1 (MQ=255)
                                                                                                                                                            gAATCTGGAGAGAATGGTCTGCTGCGAATCAGCAAACATGCCAGTATGGATAACACAACCCTtagatatgcccctatctttacagacatctgaaataacttaacccattccaagcccgctgccgaagatattcccgtggcgagcgataacccagcgcactatgaggatgccattcgttataatgctcgaacgcctctgcaaggatctttgctgccgttaacccgtctggtttgggcatgatactgatgt                                                     >  1:191616‑M1/1‑62 (MQ=255)
                                                                                                                                                                   gagagaATGGTCTGCTGCGAATCAGCCAACCTGAAAGTATGGATAACACAACCCTtggatttgcccctatatttccagccatctgttatcacttaacccattacaagcccgctgccgcagatattcccgtggcgagcgataacccagcgcactatgcggatgccattcgttataatgctcgaa                                                                                                                >  1:135819‑M1/1‑55 (MQ=255)
                                                                                                                                                                     acttatttaagtgatattggttgtctggagattcagggggccagtctaACCCTCAAGGATGACTAATCATTGAGGAAATAGAATAAATGTTCAGACCTTTTTTAAACTCTCTTATGCTCGGCAGTTTGGTTTTTCCTTTTATTGCCATTGCTGGAAGCACCGTGCAAGGGGGCGTGATCCATTTTTATgg                                                                                                       <  1:168745‑M2/142‑1 (MQ=255)
                                                                                                                                                                                          aGCCAACCTGAAAGTATGGATAACACAACCCTtggatttgaccctatatttccagacatctgttatcacttaacccattacaagcccgctgccgcagatattcccgtggcgagcgataacccagcgccctatgcggatgccattcgttataatgctcgaacgcctctgcaaggttctttgctgccgttaccccgtctggtttgggcatgatcctgctgtcgtcacgctttatcgttttcacgaagctctct                     >  1:182706‑M1/1‑32 (MQ=255)
                                                                                                                                                                                           gCCAACCTGAAAGTATGGATAACACAACCCTtggatttgcccctatatttccagacatctgttatcacttaacccatt                                                                                                                                                                                                 <  2:2841‑M1/78‑48 (MQ=255)
                                                                                                                                                                                           gCCAACCTGAAAGTATGGATAACACAACCCTtggatttgcccctatatttccagacatctgttatcacttaacccatt                                                                                                                                                                                                 >  1:2841‑M1/1‑31 (MQ=255)
                                                                                                                                                                                               ggagattcagggtgcccgtctaaCCCTCAAGGATGACTAATCATTGATGAAATAGAATAAATGTTCAGGCCTTTTTTAAACTCTGTTATGCTTGGCAGGTTGTTTTTTCCTGTTATGGCCTTTGCTGGAAGTACCGTGCAAGGGGGCGTGGTGCTTGTTTATGGCCAAATGGTGGAACCGGAATGTGTCGTCAGCACCCACTCATCAGGCGTAGAAATGACCTGCCCACAGAATGGTTCTATTCCGGGc                  <  1:278093/228‑1 (MQ=255)
                                                                                                                                                                                                ccTGAAAGTATGGATAACACAACCCTtggatttgcccctatatttccagacatctgttatcacttaacccattacaagcccgctgccgcagatattcccgtggcgagc                                                                                                                                                              <  1:78330‑M1/108‑83 (MQ=255)
                                                                                                                                                                                                              cagtctaACCCTCACGGATGACTAATCATTGAGGAAATAGAATAAATGTTCAGACCTTTTTTAAACTCTCTTATGCTCGGCAGTTTGTTTTTTCCTTTTATTGCCATTGCTGGAAGCACCGTGCAAGGGGGCGTGATCCATTTTTATGGCCAAATTGTGGAACCGGCATGTGACGTCAGCACCCAGTCATCACCCGTAGAAATGAACTGCCCACAGAATGGTTCTATTCCGGGCAAAACCTACTCCAGCa  >  1:231975‑M2/8‑250 (MQ=255)
                                                                                                                                                                                                                         |                                                                                                                                                                                                                                              
AACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTCATTTATCACAACAATGGCCGCCCTCTTCAGAAAAGTCTTAATTTGTTGAAATATCGAGCATAAGATGAATCTGGAGAGAATGGTCTGCTGCGAATCAGCCAACCTGAAAGTATGGATAACACAACCCTCAAGGATGACTAATCATTGAGGAAATAGAATAAATGTTCAGACCTTTTTTAAACTCTCTTATGCTCGGCAGTTTGTTTTTTCCTTTTATTGCCATTGCTGGAAGCACCGTGCAAGGGGGCGTGATCCATTTTTATGGCCAAATTGTGGAACCGGCATGTGACGTCAGCACCCAGTCATCACCCGTAGAAATGAACTGCCCACAGAATGGTTCTATTCCGGGCAAAACCTACTCCAGCA  >  NC_000913/1105163‑1105618

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