Predicted mutation
evidence seq id position mutation annotation gene description
JC JC CP009273 3,178,254 IS5 (+) +4 bp coding (56‑59/291 nt) yqiC → BMFP family putative fusogenic protein

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? CP009273 2282398 =NA (NA)24 (0.570) 23/294 1.0 100% noncoding (1195/1195 nt) IS5 repeat region
?CP009273 3178254 = 0 (0.000)coding (56/291 nt) yqiC BMFP family putative fusogenic protein
* ? CP009273 = 2283592NA (NA)30 (0.700) 28/296 0.6 100% noncoding (1/1195 nt) IS5 repeat region
?CP009273 = 3178257 0 (0.000)coding (59/291 nt) yqiC BMFP family putative fusogenic protein

GACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTGTTAAAAGCATTGGGTACAGAAAATACCCATAGCTCCATACCCGGAGTCA  >  CP009273/2283448‑2283738
                                                                                                                                                |                                                                                                                                                  
gACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACAAAATCAGAACGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATaa                                                                                                                                                >  1:160666/1‑149 (MQ=11)
  cAAGAAAATCTCTTGTCTGGTCTTACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGAGGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACg                                                                                                                                              <  1:480358/149‑1 (MQ=255)
    aGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGCCAGTGAATGTCTCGTGCTACACCCTTTTCTATGGCTCCAGATGACAAACATGATCTCATAACAGGGACTTGGTCGCACCTTCCATAACGCt                                                                                                                                            >  1:163468/1‑149 (MQ=11)
           ctctTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGcca                                                                                                                                     <  1:130704/149‑1 (MQ=255)
           ctctTTTCTGGTCTCACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGcca                                                                                                                                     <  1:539200/149‑1 (MQ=255)
            tctTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAAATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGccac                                                                                                                                    >  1:359937/1‑149 (MQ=255)
                 ccTGGTCTTGCGGCCGTTACGGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGaa                                                                                                                               <  2:584624/148‑1 (MQ=14)
                     gTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGa                                                                                                                           >  1:387684/1‑149 (MQ=255)
                         gACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACGTGATCTCATATCAGGGAATTGTTGGCACCTTCCATAACGCTGTAGCCACCAGAACAGATAtt                                                                                                                       >  1:471655/1‑149 (MQ=255)
                                cTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACTCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAAc                                                                                                                <  2:101419/149‑1 (MQ=255)
                                       ctgAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACaagg                                                                                                         >  1:209842/1‑147 (MQ=255)
                                           aTTCACTGTCGGCGACGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCCGATGACTAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGATGTAGTCACCATAACAGATGTTGCGGAACGACAAAGAGaa                                                                                                     <  2:43046/149‑1 (MQ=17)
                                              cACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCCCCCCTTCCATAACGCTGTACCCCCCCGAACCGATACCGCGGAACGACAAAGCGAAACa                                                                                                  >  2:357252/1‑149 (MQ=17)
                                                 tGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGaa                                                                                               >  2:386935/1‑149 (MQ=255)
                                                 tGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAAAAGac                                                                                               >  1:675374/1‑148 (MQ=255)
                                                   tCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAAcc                                                                                             <  1:255868/149‑1 (MQ=255)
                                                          aGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGGAACAGAACAAGATTGa                                                                                      >  2:662475/1‑149 (MQ=255)
                                                                gTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCAtt                                                                                >  2:68827/1‑149 (MQ=255)
                                                                   gatgaCTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAg                                                                             <  1:613543/149‑1 (MQ=255)
                                                                    atgaCTCATGATGAACCCTGTTCTATGGATCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGATCGACAAAGAGAAACAGAACCAGATTGATGCATTGAGc                                                                            >  1:453015/1‑149 (MQ=255)
                                                                         gttccggaatttcTGTTCTATGGCTCCAGATGACAAACATGATCTCCTATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTct                                                                       <  1:524879/137‑2 (MQ=37)
                                                                          cATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCAt                                                                      <  1:360596/149‑1 (MQ=255)
                                                                           atgatgAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGGAAAAGAGCCAGATTGATGCATTGAGCCTTCATc                                                                     >  2:315624/1‑149 (MQ=255)
                                                                           atgatgAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCGGATTGATGCATTGAGCTTTCATc                                                                     <  1:454462/149‑1 (MQ=255)
                                                                           atgatgAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCGGTAGCCACCAGAACAGATATTGCGGAACGACAAAAAGAAACAGAACCAGAATGATGCATTGGGCTTTCATc                                                                     >  2:96515/1‑149 (MQ=255)
                                                                              atgaACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCCCCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTa                                                                  <  1:490970/149‑1 (MQ=255)
                                                                                    ccTGTTCTATGTCTCCAGATCAAAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAaat                                                            <  1:88647/149‑1 (MQ=255)
                                                                                    ccTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAaat                                                            >  2:291760/1‑149 (MQ=255)
                                                                                    ccTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTAAGaaaa                                                            >  1:125242/1‑148 (MQ=255)
                                                                                    ccTGTTCTATGGCTCCAGATGACAAACATGACCTCATATCAGGGACTTGTTCGCACCTTCCATGACGCTGTCTCGACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCCACCTATGAaat                                                            >  2:177763/1‑149 (MQ=37)
                                                                                         tCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAaattaatt                                                       <  1:593570/149‑1 (MQ=255)
                                                                                          cTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTg                                                      >  1:154193/1‑149 (MQ=255)
                                                                                            aTGGCTCCTGATGACAAACATGATCTCATATCAGCGACTTATTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGCACGACAAAGAGAAACAGAACCAGATTTATGCATTGAGCTTTCATCCTATGAAATTAATTGCt                                                    <  2:137225/149‑1 (MQ=38)
                                                                                                           aaaCATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTGTTAAAAGCATTggg                                     >  1:598614/1‑149 (MQ=255)
                                                                                                           aaaCATGAGCTCATATCAGGGACTTGTGCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGAATGATGCATTGAGCTTTCATCCGATGCAATCAATTGCTGTTAGAAGCATTggg                                     >  2:611322/1‑149 (MQ=255)
                                                                                                              cATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGAAAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAAGTAATGGCGGTTAAAAGCATTGGGTAc                                  >  1:397363/1‑149 (MQ=255)
                                                                                                                    ctcATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTGTTAAAAGCATTGGGTACAGaaaa                            >  1:213001/1‑149 (MQ=255)
                                                                                                                     tcATATCAGGGACTTGTTCGCACATTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTGTTAAAAGCATTGGGTCCAGCAAAt                           >  1:384317/1‑149 (MQ=255)
                                                                                                                      cATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTGTTAAAAGCATTGGGTACAGAAAATa                          >  2:113809/1‑149 (MQ=255)
                                                                                                                         atCAGGGACTTGTTCGCACCTTCCAAAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTGTTAAAAGCATTGGGTACAGAAAATAccc                       >  1:415683/1‑149 (MQ=255)
                                                                                                                          tCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGGGCTTTCATCCTATGAAATTAATTGCTGTTAAAAGCATTGGGTACAGAAAATACCCa                      >  2:593330/1‑149 (MQ=255)
                                                                                                                                aCTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTGTTAAAAGCATTGGGTACAGAAAATACCCATAGCTc                <  1:180430/149‑1 (MQ=255)
                                                                                                                                        gCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTGTTAAAAGCATTGGGTACAGAAAATACCCATAGCTCCATACCCg        <  1:632585/149‑1 (MQ=255)
                                                                                                                                             ttCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTGTTAAAAGCATTGGGTACAGAAAATACCCATAGCCCCATACCCGGAGTc   >  2:23945/1‑149 (MQ=255)
                                                                                                                                              tCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTGTTAAAAGCATTGGGTACAGAAAATACCCATAGCTCCATACCCGGAGTCa  >  1:330861/1‑149 (MQ=255)
                                                                                                                                              tCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTGTTAAAAGCATTGGGTACAGAAAATACCCATAGCTCCATACCCGGAGGCa  >  2:370327/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                  
GACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTGTTAAAAGCATTGGGTACAGAAAATACCCATAGCTCCATACCCGGAGTCA  >  CP009273/2283448‑2283738

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