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)14 (0.650) 14/488 0.5 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 = 11053800 (0.000)16 (0.740) 15/490 0.4 100% intergenic (+87/‑34) csgC/ymdA curlin chaperone/uncharacterized protein YmdA
?NC_000913 1467910 = NA (NA)noncoding (1331/1331 nt) IS2 repeat region

CAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTCATTTATCACAACAATGGCCGCCCTCTTCAGAAAAGTCTTAATTTGTTGAAATATCGAGCATAAGATGAATCTGGAGAGAATGGTCTGCTGCGAATCAGCCAACCTGAAAGTATGGATAACACAACCCTCAAGGATGACTAATCATTGAGGAAATAGAATAAATGTTCAGACCTTTTTTAAACTCTCTTATGCTCGGCAGTTTGTTTTTTCCTTTTATTGCCATTGCTGGAAGCACCGTGCAAGGGGGCGTGATCCATTTTTATGGCCAAATTGTGGAACCGGCATGTGACGTCAGCACCCAGTCATCACCCGTAGAAATGAACTGCCCACAGAATGGTTCTATTCCGGGCAAAACCTACTCCAGCAAAGCGTTA  >  NC_000913/1105132‑1105626
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                      
cAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTCATTTATCACAACAATGGCCGCCCTCTTCAGAAAAGTCTTAATTTGTTGAAATATCGAGCATAAGATGAATCTGGAGAGAATGGTCTGCTGCGAATCAGCCAACCTGAAAGTATGGATAACACAACCCTtg                                                                                                                                                                                                                                                      >  1:156477‑M1/1‑249 (MQ=255)
   acgccagagaaataactggctggctgctacaccatgttgccgggcaacgagggagaccgtcatccccggttcaaagctctgctgaacaattgcgatcttttcctgtgtggtacgccgtctgcgtttctccggccctaagacatcaatcatctgttctccaatgactagtctaaaaactagtattaagactatcacgtatttaagtgatattggttgtctggagattcagggggccagtctaACCCTCAAg                                                                                                                                                                                                                                                    <  2:179222‑M2/9‑1 (MQ=255)
     aaaGACCCGTTCATTGCCTGCTAATCAACCCATTTCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATTGACAGTCACTTCATTTATCACAACAATGGCCGCCCTCTTCAGAAAAGTCTTAATTTGTTGAAATATCGAGCATAAGATGAAGCTGGAGAGAATGGTCTGCTGCGAATCTGCCAACCTGAAAGTATGGATAACACAACCCTtggattt                                                                                                                                                                                                                                                 <  2:154862‑M1/251‑8 (MQ=255)
              ataactggctggctgctacaccatgttgccgggcaacgagggagaccgtcatccccggttcaaagctctgctgaacaattgcgatcttttcctgtgtggtacgccgtctgcgtttctccggccctaagacatcaatcatctgttctccaatgactagtctaaaaactagtattaagactatcacttatttaagtgatattggttgtctggagattcagggggccagtctaACCCTCAAGGATGACTAATc                                                                                                                                                                                                                                         >  2:78307‑M2/231‑250 (MQ=255)
                          ctgctacaccatgttgccgggcaacgagggagaccgtcatccccggttcaaagctctgctgaacaattgcgatcttttcctgtgtggtacgccgtctgcgtttctccggccctaagacatcaatcatctgttctccaatgactagtctaaaaactagtcttaaaactatcacttatttaagtgatattggttgtctggagattcagggggccagtctaACCCTCAAGGATGACTAATCATTGAGGAAAt                                                                                                                                                                                                                              >  1:83167‑M2/219‑249 (MQ=255)
                              cAACCCATTGCTGTGACGAAGTTGAGGTTAAATATTTCCCCGGACGATCGCGTGAAAGTAGTTGTTAGTGTTTCGGATAGACAGTCACTTCATTTATCACAACAATGGCCGCCCTCTTGAGAAAGGACTTAATTTGTTGAAATATAGAGCATAAGATGAATCTGGAGAGAATGGTCTGCTGCGAATCAGCCAACATGGAAGTATGGATAACACAACCCTtggaatagcccctatatttgcagacatctgat                                                                                                                                                                                                                        >  1:35316‑M1/1‑219 (MQ=255)
                                                   gagggagaccgcaatccccggatcaaagctctgctgaacaattgcgatcttttcctgtgtggtacgccgtctgcgtttctccggccctaagacatcaatcatctgttctccaatgactagtctaaaaactagtattaagactatcacttatttaagtgatattggttgtctggagattcagggggccagtctaACCCTCAAGGATGACTAATCATTGAGGCAATAGAATAAATGTTCAGACCTTTTTTaa                                                                                                                                                                                                    >  2:107710‑M2/194‑250 (MQ=255)
                                                                         tcaaagctctgctgaacaattgcgatcttttcctgtgtggtacgccgtctgcgtttctccggccctaagacatcaatcatctgttctccaatgactagtctaaaaactagtattaagactatcacttatttaagtgatattggttgtctggagattcagggggccagtctaACCCTCAAGGATGACTAATCATTGAGGAAATAGAATAAATGTTCAGACCTTTTTTAAACTCTCTTATGCTCGGCAgttt                                                                                                                                                                              <  1:131504‑M2/79‑1 (MQ=255)
                                                                                 ctgctgaacaattgcgatcttttcctgtgtggtacgccgtctgcgtttctccggccctaagacatcaatcatctgttctccaatgactagtctaaaaactagtattaagactatcacttatttaagtgatattggttgtctggagattcagggggccagtctaACCCTCAAGGATGACTAATCATTGAGGAAATAGAATAAATGTTCAGACCTTTTTTAAACTCTCTTATGCTCGGCAGTTTGTTTTTTcc                                                                                                                                                                     >  1:161691‑M2/164‑251 (MQ=255)
                                                                                           ttgttACTGCTTCTGATGGACAGGCACTACATTTATCACAACAATGCCTGCCCTCTTCAGAAAAGTCTTAATTTGTTGAAATATCGAGCATACGATGAATATGGAGAGAATGGTCTGCTGCGAATCAGCCAACCTGAAAGTATGGATAACACAACCCTgcgatttgcccctatatttccagacatctgttaccacttaacccattacaagcccgctgccgccgctatgcccgttgcgagcgataacccagc                                                                                                                                                           >  2:220064‑M1/1‑158 (MQ=255)
                                                                                                             gtggtacgccgtctgcgtttctccggccctaagacatcaatcatctgttctccaatgactagtctaaaaactagtattaagactatcacttatttaagtgatattggttgtctggagattcagggggccagtctaACCCTCAAGGATGACTAATCATTGAGGAAATAGAATAAATGTTCAGACCTTTTTTAAACTCTCTTATGCTCGGCAGTTTGTTTTTTCCTTTTATTGCCATTGCTGGAAGCACCGTg                                                                                                                                         <  1:78307‑M2/116‑1 (MQ=255)
                                                                                                                                        ggCCGACCTCTTCAGAAAAGACTTAATTTGTTGAAATATCGAGCATAAGATGAATCTGGAGAGAATGGTCTGCTGCGAATCAGCCAACCTGAAAGTATGGATAACACAACCCTtggatttgcccctatatttccagacatctgttatcacttaacccattaccagaccgctgccgcagatattcccgaggcgagcgataacccagcgcactatgcggatgccattcgttataatgctcgaacgcatctgc                                                                                                               >  2:56870‑M1/1‑113 (MQ=255)
                                                                                                                                               catcaatcatctgttctccaatgactagtctaaaaactagtattaagactatcacttatttaagtgatattggttgtctggagattcagggggccagtctaACCCTCAAGGATGACTAATCATTGAGGAAATAGAATAAATGTTCAGACCTTTTTTAAACTCTCTTATGCTCGGCAGTTTGTTTTTTCCTTTTATTGCAATTGCTGGAAGCACCGTGCAAGGGGGCGTGATCCATTTTTATGGCCAAATtg                                                                                                       >  2:65564‑M2/102‑251 (MQ=255)
                                                                                                                                                  caatcatctgttctccaatgactagtctaaaaactagtattaagactatcacttatttaagtgatattggttgtctggagattcagggggccagtctaACCCTCAAGGATGACTAATCATTGAGGAAATAGAATAAATGTTCAGACCTTTTTTAAACTCTCTTATGCTCGGCAGTTTGTTTTTTCCTTTTATTGCCATTGCTGGAAGCACCGTGCAAGGGGGCGTGATCCATTTTTATGGCCAAATtgtg                                                                                                     <  1:57119‑M2/152‑1 (MQ=255)
                                                                                                                                                   aatcatctgttctccaatgactagtctaaaaactagtattaagactatcacttatttaagtgatattggttgtctggagattcagggggccagtctaACCCTCAAGGATGACTAATCATTGAGGAAATAGAATAAATGTTCAGACCTTTTTTAAACTCTCTTATGCTCGGCAGTTTGTTTTTTCCTTTTATTGCCATTGCTGGAAGCACCGTGCAAGGGGGCGTGATCCATTTTTATGGCCAAATTGTgg                                                                                                    >  1:131030‑M2/98‑250 (MQ=255)
                                                                                                                                                         aaGTCTTAATTTGTTGAAATATCGAGCATAAGATGAATCTGGAGAGAATGGTCTGCTGCGAATCAGCCAACCTGAAAGTATGGATAACACAACCCTtggatttgcccctatatttccagacatctgttatcacttaacccattacaagcccgctgccgcagatattcccgtggcgagcgataacccagcgcactatgcggatgccattcgttataatgctcgaacgcctctgcaaggttctttgctgccgt                                                                                             >  2:69310‑M1/1‑96 (MQ=255)
                                                                                                                                                                                         aTGAATCTGGAGAGAATGGTCTGCTGCGAAGCAGCCAACCTGAAAGTATGGATAACACAACCCTtggatttgcccctatatttccagacatctgttatcacttaacccattacaagcccgctgccgcagatattcccgtggcgagcgataacccagcgcactatgcggatgccattcgttataatgctcgaacgcctctgcaaggttctttgctgccgttaacccgtctggtttgggcatgatactgatgt                                                             <  2:73454‑M1/251‑188 (MQ=255)
                                                                                                                                                                                              tCTGGAGAGAATGGTCTGCTGCGAATCAGCCAACCTGAAAGTATGGATAACACAACCCTtggatttgcccctatatttccagacatctgttatctcttaacccattacaagcccgctgccgcagatattcccgtggcgagcgataacccagcgcactatgcggatgccattcgttataatgctcgaacgcctctgcaaggttctttgctgccgttaacccgtctggtttgggcatgatactgatgtagtcc                                                        >  1:96849‑M1/1‑59 (MQ=255)
                                                                                                                                                                                                   agagaATGGTCTGCTGCGAATCAGCCAACCTGAAAGTATGGATAACACAACCCTtggatttgcccctatatttccagacatctgttatcacttaacccattacaagcccgctgccgcagatattcccgtggcgagcgatatcccagcgcactatgcggatgccattcgttataatgctcgaacgcctctgcaaggttctttgcttccgttaacccgtctggtttgggcatgatactgatgtagtcacgctt                                                   <  1:19278‑M1/251‑198 (MQ=255)
                                                                                                                                                                                                           gTCTGCTGCGAATCAGCCAACCTGAAAGTATGGATAACACAACCCTtggatttgcccctatatttccagacatctgttatcacttaacccattacaagcccgctgccgcagatattcccgtggcgagcgataacccagcgcactatgcggatgccattcgttataatgctcgaacgcctctgcaaggttctttgctgccgttaacccgtctggtttgggcatgatactgatgtagtcacgctttatcgttt                                           >  2:219109‑M1/1‑46 (MQ=255)
                                                                                                                                                                                                                      ggttgtctggagcttcagggggccagtctaACCCTCAAGGATGACTAATCATAGAGGAAATAGAATAAATGTTCAGACCTTTTTTAAACTCTCTTATGCTCGGCAGTTTGTTTTTTCCTTTTATTGCCATTGCTGGAAGCACCGTGCAAGGGGGCGTGATCCATTTTTATGGCCAAATTGTGGAACCGGCATGTGACGTCAGCACCCAGTCATCACCCGTAGAAATGAACTGCCAACAGAATGGTTCTAtt                                >  2:122941‑M2/31‑251 (MQ=255)
                                                                                                                                                                                                                            cAACCTGAAAGTATGGATAACACAACCCTttgatttgcccctatatttccagacatctgttatcacttaacccattacaagcccgctgccgcagatattcccgtggcgagcgatatcccagcgcactatgcggatgccattcgttataatgctcgaacgcctctgcaaggttctttgctgccgttaacccgtctggtttgggcatgattctaatgtagttacgctttatcgtagtcacgaagctctctgct                          <  1:61469‑M1/251‑223 (MQ=255)
                                                                                                                                                                                                                              ggagattcagggggccagtctaACCCTCAAGGATGACTAATCATTGAGGAAATAGAATAAATGTTCAGACCTTTTTTAAACTCTCTTATGCTCGGCAGTTTGTTTTTTCCTTTTATTGCCATTGCTGGAAGCACCGTGCAAGGGGGCGTGATCCATTTTTATGGCCAAATTGTGGAACCGGCATGTGACGTCAGCACCCAGTCATCACCCGTAGAAATGAACTGCCCACAGAATGGTTCTATTCCGGGc                          <  1:208807‑M2/227‑1 (MQ=255)
                                                                                                                                                                                                                                agattcagggggccagtctaACCCTCAAGGATGACTAATCATTGAGGAAATAGAATAAATGTTCAGACCTTTTTTAAACTCTCTTATGCTCGGCAGTTTGTTTTTTCCTTTTATTGCCATTGCTGGAAGCACCGTGCAAGGGGGCGTGATCCATTTTTATGGCCAAATTGTGGAACCGGCATGTGACGTCAGCACCCAGTCATCACCCGTAGAAATGAACTGCCCACAGAATGGTTCTATTCCGGGCaaaa                      <  2:89320‑M2/231‑1 (MQ=255)
                                                                                                                                                                                                                                  gAAAGTATGGATAACACAACCCTtggatttgcccctatatttccagacatctgttatcacttaacccattacaagcccgctgccgcagatattcccgtggtgagcgataacccagcgcactatgcggatgccattcgttataatgctcgaacgcctctgcaaggttctttgctgccgttaacccgtctggtttgggcatgatactgaggtagtcacgctttatcgttttcacgaagctctctgctattccg                    >  2:56652‑M1/1‑23 (MQ=255)
                                                                                                                                                                                                                                  gAAAGTATGGATAACACAACCCTtggatttgcccctatagttccagacatctgttatcacttaacccattacaagcccgctgccgcagatattcccgtggcgagcgataacccagcgcactatgcggatgccattcgttataatgctcgaacgcctctgcaaggttctttgctgccgttaacccgtctggtttgggcatgatactgatgtagtcacgcgttatcgttttcacgaagctctctgctattc                      <  2:156477‑M1/249‑227 (MQ=255)
                                                                                                                                                                                                                                              aaCACAACCCTtggatttgcccctatatttccagacatctgttatcacttaacccattacaagcccgctgccgcagatattcccgtggcgagcgataacccagcgcactatgcggatgccattcgttataatgctcgaacgcctctgcaaggttctttgctgccgttaacccgtctggtttgggcatgatactgatgtagtcacgctttatcgttttcacgaagctctctgctattccgttactctccgg         <  2:194063‑M1/250‑240 (MQ=255)
                                                                                                                                                                                                                                                tctaACCCTCAAGGATGACTAATCATTGAGGAAATAGAATAAATGTTCAGACCTTTTTTAAACTCTCTTATGCTCGGCAGTTTGTTTTTTCCTTTTATTGCCATTGCTGGAAGCACCGTGCAAGGGGGCGTGATCCATTTTTATGGCCAAATTGTGGAACCGGCATGTGACGTCAGCACCCAGTCATCACCCGTAGAAATGAACTGCCCACAGAATGGTTCTATTCCGGGCAAAACCTACTCCAGCAAAGc      >  1:38810‑M2/5‑251 (MQ=255)
                                                                                                                                                                                                                                                  caACCCTtggatttgcccctatatttccagacatctgttatcacttaacccattacaagcccgctgccgcagatattcccgtggcgagcgataacccagcgcactatgcggatgccattcgttataatgctcgaacgcctctgcaaggttctttgctgccgttaacccgtctggtttgggcatgatactgatgtagtcacgctttatcgttttcacgaagctctctgctattccgttactctccggactcc    >  1:64848‑M1/1‑7 (MQ=255)
                                                                                                                                                                                                                                                  caACCCTtggatttgcccctatatttccagacatctgttatcacttaacccattacaagcacgctgccgcagatattcccgtggcgagcgctaacccagcgcactatgcggatgcccttcgttataatgatcgaacgcctatgcaaggttctttgctgccgttaccccgtctggttcgggcatgatactgatgtagtcacgctttatcgttttcacgaagctctctgctattccgttactctccggactcc    >  1:218402‑M1/1‑7 (MQ=255)
                                                                                                                                                                                                                                                    aCCCTCAAGGATGACTAATCATTGAGGAAATAGAATAAATGTTCAGACCTTTTTTAAACTCTCTTATGCTCGGCAGTTTGTTTTTTCCTTTTATTGCCATTGCTGGAAGCACCGTGCAAGGGGGCGTGATCCATTTTTATGGCCAAATTGTGGAACCGGCATGTGACGTCAGCACCCAGTCATCACCCGTAGAAATGAACTGCCCACAGAGTGGTTCTAGTCCGGGCAAAACCTACTCCAGCAAAGCGTTa  <  1:75698/251‑1 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                      
CAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTCATTTATCACAACAATGGCCGCCCTCTTCAGAAAAGTCTTAATTTGTTGAAATATCGAGCATAAGATGAATCTGGAGAGAATGGTCTGCTGCGAATCAGCCAACCTGAAAGTATGGATAACACAACCCTCAAGGATGACTAATCATTGAGGAAATAGAATAAATGTTCAGACCTTTTTTAAACTCTCTTATGCTCGGCAGTTTGTTTTTTCCTTTTATTGCCATTGCTGGAAGCACCGTGCAAGGGGGCGTGATCCATTTTTATGGCCAAATTGTGGAACCGGCATGTGACGTCAGCACCCAGTCATCACCCGTAGAAATGAACTGCCCACAGAATGGTTCTATTCCGGGCAAAACCTACTCCAGCAAAGCGTTA  >  NC_000913/1105132‑1105626

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.