Predicted mutation
evidence seq id position mutation annotation gene description
JC JC NC_000913 751,781 IS5 (–) +4 bp coding (386‑389/2448 nt) ybgQ ← putative fimbrial usher protein YbgQ

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 687851 =NA (NA)46 (1.290) 34/278 0.1 100% noncoding (1195/1195 nt) IS5 repeat region
?NC_000913 = 751784 0 (0.000)coding (386/2448 nt) ybgQ putative fimbrial usher protein YbgQ
* ? NC_000913 = 689045NA (NA)37 (1.040) 32/278 0.1 100% noncoding (1/1195 nt) IS5 repeat region
?NC_000913 751781 = 0 (0.000)coding (389/2448 nt) ybgQ putative fimbrial usher protein YbgQ

TCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTAAAGATATTGCTGTTTTATTGTTTGTTTTTGCGAGATGCGCCGCACCATTCCGAAGCAAAATTCTTAAAATGCACTCTTTTAGTGCTACCGCTGGATTACTGTGGTGCAA  >  NC_000913/688918‑689183
                                                                                                                               |                                                                                                                                          
tCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACATGTTCGCACCTTCCTTAGGTAACAtt                                                                                                                                <  1:495692/140‑1 (MQ=21)
 cTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACAttt                                                                                                                               >  1:332994/1‑140 (MQ=34)
    gTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGt                                                                                                                            <  1:214402/140‑1 (MQ=38)
       tGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTg                                                                                                                         >  1:90119/1‑140 (MQ=40)
         aCGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGc                                                                                                                       >  1:439482/1‑140 (MQ=255)
         aCGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGc                                                                                                                       >  1:341113/1‑140 (MQ=255)
           ggCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTa                                                                                                                     <  2:212915/140‑1 (MQ=255)
           ggCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTa                                                                                                                     <  1:362848/140‑1 (MQ=255)
             cgcTTACTGCTGAATTCACTGTCGTCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTaaa                                                                                                                   <  2:599894/140‑1 (MQ=255)
                   ctgctgAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTaaa                                                                                                             >  1:121953/1‑140 (MQ=255)
                          aTTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTAAAGATATTg                                                                                                      <  1:62961/140‑1 (MQ=255)
                              aCTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTAAAGATATTGCTGt                                                                                                  <  1:220914/140‑1 (MQ=255)
                                   cGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTAAAGATATTGCTGTTTTAt                                                                                             >  2:236644/1‑140 (MQ=255)
                                    ggCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTAAAGATATTGCTGTTTTAtt                                                                                            >  1:29233/1‑140 (MQ=255)
                                                   atgaCTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTAAAGATATTGCTGTTTTATTGTTTGTTTTTGCgag                                                                             <  2:592055/140‑1 (MQ=255)
                                                     gaCTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTAAAGATATTGCTGTTTTATTGTTTGTTTTTGCGAGAt                                                                           <  2:434093/140‑1 (MQ=255)
                                                     gaCTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTAAAGATATTGCTGTTTTATTGTTTGTTTTTGCGAGAt                                                                           >  1:417774/1‑140 (MQ=255)
                                                     gaCTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTAAAGATATTGCTGTTTTATTGTTTGTTTTTGCGAGAt                                                                           >  2:311474/1‑140 (MQ=255)
                                                      aCTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTAAAGATATTGCTGTTTTATTGTTTGTTTTTGCGAGATg                                                                          <  2:327093/140‑1 (MQ=255)
                                                      aCTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTAAAGATATTGCTGTTTTATTGTTTGTTTTTGCGAGATg                                                                          <  1:144313/140‑1 (MQ=255)
                                                                             ggCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTAAAGATATTGCTGTTTTATTGTTTGTTTTTGCGAGATGCGCCGCACCATTCCGAAGCaaaa                                                   <  2:332994/140‑1 (MQ=255)
                                                                               cTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTAAAGATATTGCTGTTTTATTGTTTGTTTTTGCGAGATGCGCCGCACCATTCCGAAGCAAAAtt                                                 >  1:563364/1‑140 (MQ=255)
                                                                                 ccAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTAAAGATATTGCTGTTTTATTGTTTGTTTTTGCGAGATGCGCCGCACCATTCCGAAGCAAAATTCt                                               >  2:455451/1‑140 (MQ=255)
                                                                                 ccAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTAAAGATATTGCTGTTTTATTGTTTGTTTTTGCGAGATGCGCCGCACCATTCCGAAGCAAATTCtt                                              >  1:124018/1‑140 (MQ=255)
                                                                                    gATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTAAAGATATTGCTGTTTTATTGTTTGTTTTTGCGAGATGCGCCGCACCATTCCGAAGCAAAATTCTTaa                                            >  1:145852/1‑140 (MQ=255)
                                                                                    gATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTAAAGATATTGCTGTTTTATTGTTTGTTTTTGCGAGATGCGCCGCACCATTCCGAAGCAAAATTCTTaa                                            >  2:539749/1‑140 (MQ=255)
                                                                                         cAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTAAAGATATTGCTGTTTTATTGTTTGTTTTTGCGAGATGCGCCGCACCATTCCGAAGCAAAATTCTTAAAATGc                                       >  2:106947/1‑140 (MQ=255)
                                                                                                 aTCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTAAAGATATTGCTGTTTTATTGTTTGTTTTTGCGAGATGCGCCGCACCATTCCGAAGCAAAATTCTTAAAATGCACTCtttt                               >  2:289866/1‑140 (MQ=255)
                                                                                                       tatCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTAAAGATATTGCTGTTTTATTGTTTGTTTTTGCGAGATGCGCCGCACCATTCCGAAGCAAAATTCTTAAAATGCACTCTTTTAGTGCt                         >  1:178615/1‑140 (MQ=255)
                                                                                                           aGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTAAAGATATTGCTGTTTTATTGTTTGTTTTTGCGAGATGCGCCGCACCATTCCGAAGCAAAATTCTTAAAATGCACTCTTTTAGTGCTACCg                     <  2:121953/140‑1 (MQ=255)
                                                                                                             ggACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTAAAGATATTGCTGTTTTATTGTTTGTTTTTGCGAGATGCGCCGCACCATTCCGAAGCAAAATTCTTAAAATGCACTCTTTTAGTGCTACCGCt                   >  2:121945/1‑140 (MQ=255)
                                                                                                             ggACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTAAAGATATTGCTGTTTTATTGTTTGTTTTTGCGAGATGCGCCGCACCATTCCGAAGCAAAATTCTTAAAATGCACTCTTTTAGTGCTACCGCt                   >  1:556421/1‑140 (MQ=255)
                                                                                                              gACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTAAAGATATTGCTGTTTTATTGTTTGTTTTTGCGAGATGCGCCGCACCATTCCGAAGCAAAATTCTTAAAATGCACTCTTTTAGTGCTACCGCTg                  <  1:568759/140‑1 (MQ=255)
                                                                                                                cTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTAAAGATATTGCTGTTTTATTGTTTGTTTTTGCGAGATGCGCCGCACCATTCCGAAGCAAAATTCTTAAAATGCACTCTTTTAGTGCTACCGCTGGa                <  2:326208/140‑1 (MQ=255)
                                                                                                                cTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTAAAGATATTGCTGTTTTATTGTTTGTTTTTGCGAGATGCGCCGCACCATTCCGAAGCAAAATTCTTAAAATGCACTCTTTTAGTGCTACCGCTGGa                >  2:373195/1‑140 (MQ=255)
                                                                                                                cTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTAAAGATATTGCTGTTTTATTGTTTGTTTTTGCGAGATGCGCCGCACCATTCCGAAGCAAAATTCTTAAAATGCACTCTTTTAGTGCTACCGCTGGa                >  2:429626/1‑140 (MQ=255)
                                                                                                                    ttCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTAAAGATATTGCTGTTTTATTGTTTGTTTTTGCGAGATGCGCCGCACCATTCCGAAGCAAAATTCTTAAAATGCACTCTTTTAGTGCTACCGCTGGATTAc            >  1:41215/1‑140 (MQ=255)
                                                                                                                         aCCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTAAAGATATTGCTGTTTTATTGTTTGTTTTTGCGAGATGCGCCGCACCATTCCGAAGCAAAATTCTTAAAATGCACTCTTTTAGTGCTACCGCTGGATTACTgtgg       <  2:90119/140‑1 (MQ=255)
                                                                                                                         aCCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTAAAGATATTGCTGTTTTATTGTTTGTTTTTGCGAGATGCGCCGCACCATTCCGAAGCAAAATTCTTAAAATGCACTCTTTTAGTGCTACCGCTGGATTACTgtgg       <  2:145852/140‑1 (MQ=255)
                                                                                                                           cttccttAGGTAACATTTAGTTTGGCTAAATGTAAAGATATTGCTGTTTTATTGTTTGTTTTTGCGAGATGCGCCGCACCATTCCGAAGCAAAATTCTTAAAATGCACTCTTTTAGTGCTACCGCTGGATTACTgtggtg     <  1:289866/140‑1 (MQ=255)
                                                                                                                             tccttAGGTAACATTTAGTTTGGCTAAATGTAAAGATATTGCTGTTTTATTGTTTGTTTTTGCGAGATGCGCCGCACCATTCCGAAGCAAAATTCTTAAAATGCACTCTTTTAGTGCTACCGCTGGATTACTgtggtgta   >  2:193256/1‑138 (MQ=255)
                                                                                                                              ccttAGGTAACATTTAGTTTGGCTAAATGTAAAGATATTGCTGTTTTATTGTTTGTTTTTGCGAGATGCGCCGCACCATTCCGAAGCAAAATTCTTAAAATGCACTCTTTTAGTGCTACCGCTGGATTACTGTGGTGCaa  >  1:363994/1‑140 (MQ=255)
                                                                                                                               |                                                                                                                                          
TCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTAAAGATATTGCTGTTTTATTGTTTGTTTTTGCGAGATGCGCCGCACCATTCCGAAGCAAAATTCTTAAAATGCACTCTTTTAGTGCTACCGCTGGATTACTGTGGTGCAA  >  NC_000913/688918‑689183

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