Predicted mutation
evidence seq id position mutation annotation gene description
JC JC NC_000913 4,543,608 Δ1 bp :: IS186 (–) +6 bp :: Δ1 bp coding (494‑499/549 nt) fimA → major type 1 subunit fimbrin (pilin)

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 2514274 =NA (NA)61 (1.900) 50/238 0.0 100% noncoding (2/1345 nt) IS186 repeat region
?NC_000913 4543608 = 0 (0.000)coding (494/549 nt) fimA major type 1 subunit fimbrin (pilin)
* ? NC_000913 = 2515616NA (NA)20 (0.620) 19/238 0.8 100% noncoding (1344/1345 nt) IS186 repeat region
?NC_000913 = 4543613 0 (0.000)coding (499/549 nt) fimA major type 1 subunit fimbrin (pilin)

AGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGTTGCGTCTCGACAAACTCCGCGACCACACTCAATGA  >  NC_000913/2515485‑2515762
                                                                                                                                   |                                                                                                                                                  
aGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGttttt                                                                                                                                   >  5:7355/1‑149 (MQ=39)
  gTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGtttttt                                                                                                                                  >  4:102480/1‑148 (MQ=40)
        aTTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCt                                                                                                                           <  6:7355/149‑1 (MQ=255)
               ggACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAAt                                                                                                                     >  7:84995/1‑148 (MQ=255)
                  cAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGtttttt                                                                                                                                  >  1:48026/1‑132 (MQ=40)
                  cAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGtttttt                                                                                                                                  <  2:48026/132‑1 (MQ=40)
                     tATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTa                                                                                                                                <  8:72863/131‑1 (MQ=255)
                     tATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTa                                                                                                                                >  7:72863/1‑131 (MQ=255)
                       tCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAAt                                                                                                             <  8:84995/148‑1 (MQ=255)
                        cccTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATgg                                                                                                           >  6:75866/1‑149 (MQ=255)
                           tGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTg                                                                                                        >  7:33425/1‑149 (MQ=255)
                                   tACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGAc                                                                                                >  5:78431/1‑149 (MQ=255)
                                                    gaaCGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGc                                                                               <  5:75866/149‑1 (MQ=255)
                                                              cATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCa                                                                     <  6:78431/149‑1 (MQ=255)
                                                               aTCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGtttttt                                                                                                                                  <  6:94187/87‑1 (MQ=255)
                                                               aTCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGtttttt                                                                                                                                  >  5:94187/1‑87 (MQ=255)
                                                                tCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATAccc                                                                                   <  2:21869/133‑1 (MQ=255)
                                                                tCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATAccc                                                                                   >  1:21869/1‑133 (MQ=255)
                                                                   tCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAAt                                                                                            >  8:89933/1‑121 (MQ=255)
                                                                   tCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAAt                                                                                            <  7:89933/121‑1 (MQ=255)
                                                                   tCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCAcccc                                                                 >  2:67869/1‑148 (MQ=255)
                                                                                ccGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATa                                                   <  3:102480/149‑1 (MQ=255)
                                                                                       gTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATc                                                                                           >  2:53258/1‑102 (MQ=255)
                                                                                       gTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATc                                                                                           <  1:53258/102‑1 (MQ=255)
                                                                                       gTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCctg                                            >  7:36703/1‑149 (MQ=255)
                                                                                           cATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCctgctgc                                        >  7:64461/1‑149 (MQ=255)
                                                                                                ttcAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCtt                                                                                 >  1:32181/1‑103 (MQ=255)
                                                                                                ttcAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCtt                                                                                 <  2:32181/103‑1 (MQ=255)
                                                                                                    gATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGTTGCGtc                               <  8:64461/149‑1 (MQ=255)
                                                                                                         ttCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGTTGCGTCTCGAc                          <  1:67869/149‑1 (MQ=255)
                                                                                                          tCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGtt                                    <  2:59060/138‑1 (MQ=255)
                                                                                                          tCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGtt                                    >  1:59060/1‑138 (MQ=255)
                                                                                                            aaCCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGAc                                                                                                <  2:58069/76‑1 (MQ=255)
                                                                                                            aaCCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGAc                                                                                                >  1:58069/1‑76 (MQ=255)
                                                                                                               ccTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCCGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGTTGCGTCTCGACAAACTc                    >  1:184935/1‑149 (MQ=255)
                                                                                                                    aGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGTTGCGTCTCGACAAACTCCGCGa               >  4:18067/1‑149 (MQ=255)
                                                                                                                                tGGGATAATTCCCCGGTGTTTTTACGCCTGTTAATCAGGTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGTTGCGTCTCGACAAACTCCGCGACCACACTCTAt    <  2:184935/148‑1 (MQ=255)
                                                                                                                                  ggATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGTTGCGTCTCGACAAACTCCGCGACCACACTCAATGa  <  3:18067/148‑1 (MQ=255)
                                                                                                                                   |                                                                                                                                                  
AGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGTTGCGTCTCGACAAACTCCGCGACCACACTCAATGA  >  NC_000913/2515485‑2515762

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