Predicted mutation
evidence seq id position mutation annotation gene description
JC JC NC_000913 4,541,967 IS5 (–) +4 bp intergenic (+408/‑67) fimB → / → fimE Type 1 fimbriae regulatory protein FimB/regulator for fimA

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 3652036 =NA (NA)23 (0.730) 21/486 0.4 100% noncoding (1195/1195 nt) IS5 repeat region
?NC_000913 = 4541970 0 (0.000)intergenic (+411/‑67) fimB/fimE Type 1 fimbriae regulatory protein FimB/regulator for fimA
* ? NC_000913 = 3653230NA (NA)29 (0.920) 26/486 0.3 100% noncoding (1/1195 nt) IS5 repeat region
?NC_000913 4541967 = 0 (0.000)intergenic (+408/‑70) fimB/fimE Type 1 fimbriae regulatory protein FimB/regulator for fimA

ACTAATTTTATTGCTGAACAACGAGCTGACTTACCAACTAATATAAAAGATATTCAAAGCCTGTATCAAAAAATGACGAAGTTATATATTGAGCACAGTGAAAATAAAAACAGGATGAAAGTCTTTGCCGGAACCAACTTTATTGATTTCAATATGACAGGACAGAACCTGTCCGGGTTCGTTCTAACCTTGTCAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCACTGACGTATCATTTGGTCCGCCCGAAACAGGTTGGCCAGCGTGAATAACATCGCCA  >  NC_000913/3651789‑3652255
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                           
aCTAATTTTATTGCTGAACAACGAGCTGACTTACCAACTAATATAAAAGATATTCAAAGCCTGTATCAAAAAATGACGAAGTTATATATTGAGCACAGTGAAAATAAAAACAGGATGAAAGTCTTTGCCGGAACCAACTTTATTGATTTCAATATGACAGGACAGAACCTGTCCGGGTTCGTTCTAACCTTGTCAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGaa                                                                                                                                                                                                                          <  2:183520/251‑1 (MQ=255)
                                              aaGATATTCAAAGCCTGTATCAAAAAATGACGAAGTTATATATTGAGCACAGTGAAAATAAAAACAGGATGAAAGTCTTTGCCGGAACCAACTTTATTGATTTCAATATGACAGGACAGAACCTGTCCGGGTTCGTTCTAACCTTGTCAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGtcatttc                                                                                                                                                                            <  1:328112/251‑1 (MQ=255)
                                                gATATTCAAAGCCTGTATCAAAAAATGACGAAGTTATATATTGAGAACAGTGAAAATAAAAACAGGATGAAAGTCTTTGCCGGAACCAACTTTATTGATTTCAATATGACAGGACAGAACCTGTCCGGGTTCGTTCTAACCTTGTCAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGtcatttact                                                                                                                                                                          >  2:36242/1‑248 (MQ=255)
                                                                    aaaaaTGACGAAGTTATATATTGAGCACAGTGAAAATAAAAACAGGATGAAAGTCTTTGCCGGAACCAACTTTATTGATTTCAATATGACAGGACAGAACCTGTCCGGGTTCGTTCTAACCTTGTCAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGAt                                                                                                                                                      <  1:291285/251‑1 (MQ=255)
                                                                                  tatataTTGAGCACAGTGAAAATAAAAACAGGATGAAAGTCTTTGCCGGAACCAACTTTATTGATTTCAATATGACAGGACAGAACCTGTCCGGGTTCGTTCTAACCTTGTCAAGGTTTTATTTCGAAGATCTACTTAATATACATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTa                                                                                                                                         >  1:40145/1‑250 (MQ=255)
                                                                                                                 gATGAAAGTCTTTGCCGGAACCAACTTTATTGATTTCAATATGACAGGACAGAACCTGTCCGGGTTCGTTCTAACCTTGTCAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTttc                                                                                                         >  1:102199/1‑251 (MQ=255)
                                                                                                                    gAAAGTCTTTGCCGGAACCAACTTTATTGATTTCAATATGACAGGACAGAACCTGTCCGGGTTCGTTCTAACCTTGTCAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTttcttc                                                                                                      <  1:223911/251‑1 (MQ=255)
                                                                                                                                   aaCCAACTTTATTGATTTCAATATGACAGGACAGAACCTGTCCGGGTTCGTTCTAACCTTGTCAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGtt                                                                                       >  2:237755/1‑251 (MQ=255)
                                                                                                                                    aCCAACTTTATTGATTTCAATATGACAGGACAGAACCTGTCCGGGTTCGTTCTAACCTTGTCAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTa                                                                                      <  1:237755/251‑1 (MQ=255)
                                                                                                                                        aCTTTATTGATTTCAATATGACAGGACAGAACCTGTCCGGGTTCGTTCTAACCTTGTCAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGTGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATccc                                                                                  <  2:319134/251‑1 (MQ=255)
                                                                                                                                           ttATTGATTTCAATATGACAGGACAGAACCTGTCCGGGTTCGTTCTAACCTTGTCAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAATGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCAGTTCATTTCTTCATGTGTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTGAGACCGTTTCGTCGCCATTTAAGGCGTTATCCCCa                                                                                <  2:125489/250‑1 (MQ=255)
                                                                                                                                           ttATTGATTTCAATATGACAGGACAGAACATGTCCGGGTTCGTTCTCACATTGTCAAGGCTTTATTTCTATGATCTACTTAATATACTTTTCACTGACGCAAATTTAGGGAAGGTGCGTATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTTATTTCATTTCTTCATGTTTGAGCCGTTTTTTTCTCCCGTAAATGCCTTGCATCAGCATATTTAGACCGTTTCTCCGCCCTTTATGTCTTTATcccc                                                                                 >  1:139866/1‑249 (MQ=18)
                                                                                                                                                            acaggacagATCCTGTCCGGGTTCGTTCTAACCTTGTCAAGGTTTTAATTCGAAGAGCTACATAATAGAAAGTTCGCTGACGCAAATTTAGGGATGGTTCGAATAAGCGGGGTAATGCTTCTCGGCTGACTGAGTGATATCATGTCTTCATGGTTGAGCCGATTGTTGCTCCCGTAAATGCCTTGAATCAGCCTAGTTAGACCGTTTCTGCTCCGTTTAAGGCGGTAACCCCAgcattgagttggatctct                                                              <  2:128352/251‑18 (MQ=255)
                                                                                                                                                                                         aaCCTTGTCAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCCGTTTTTAGTGAGATCTCTCCCACTGACGTATCATTTGGTCCGCCCg                                  >  2:211494/1‑250 (MQ=255)
                                                                                                                                                                                         aaCCGTGTCAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTGAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATTTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCTTATTTAGGCCGTTTCTTCGCCATTTAAGGCGTTAGCCCCAGTTTTTAGTGAGATCTCTCCCACTGACGTATCATTTGGTCCGCCCg                                  <  2:143587/250‑1 (MQ=255)
                                                                                                                                                                                                 cAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCTAATAAGCGGTGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCACTGACGTATCATTTGGTCCGCCCGAATCAGGtt                         <  2:111191/251‑1 (MQ=255)
                                                                                                                                                                                                                aaGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCACTGACGTATCATTTGGTCCGCCCGAAACAGGTTGGCCAGCGTGAATa           >  1:73902/1‑250 (MQ=255)
                                                                                                                                                                                                                         ttAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCACTGACGTATCATTTGGTCCGCCCGAAACAGGTTGGCCAGCGTGAATAACATCGcc   >  1:55147/1‑249 (MQ=40)
                                                                                                                                                                                                                          tAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCACTGACGTATCATTTGGTCCGCCCGAAACAGGTTGGCCAGCGTGAATAACATCGCCa  <  1:13327/249‑1 (MQ=40)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                           
ACTAATTTTATTGCTGAACAACGAGCTGACTTACCAACTAATATAAAAGATATTCAAAGCCTGTATCAAAAAATGACGAAGTTATATATTGAGCACAGTGAAAATAAAAACAGGATGAAAGTCTTTGCCGGAACCAACTTTATTGATTTCAATATGACAGGACAGAACCTGTCCGGGTTCGTTCTAACCTTGTCAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCACTGACGTATCATTTGGTCCGCCCGAAACAGGTTGGCCAGCGTGAATAACATCGCCA  >  NC_000913/3651789‑3652255

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