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)60 (1.030) 37/294 0.6 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)76 (1.300) 52/296 0.1 100% noncoding (1/1195 nt) IS5 repeat region
?CP009273 = 3178257 0 (0.000)coding (59/291 nt) yqiC BMFP family putative fusogenic protein

AAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCG  >  CP009273/2282251‑2282530
                                                                                                                                                   |                                                                                                                                    
aaaTTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATCTaagg                                                                                                                                     >  2:418137/1‑149 (MQ=255)
aaaTTCTTTTATATGCTGTTCTCCGACAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATacgg                                                                                                                                     >  1:533931/1‑149 (MQ=255)
  aTTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGGGAATATaaggaa                                                                                                                                   >  1:106826/1‑149 (MQ=255)
  aTTCTTTCATATGCGGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACCGTTGATGCCAATGCTATGGGTGAATATaaggaa                                                                                                                                   <  1:467099/149‑1 (MQ=255)
  aTTCTGTTACATGCTGTTCGCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATaaggaa                                                                                                                                   <  2:236495/149‑1 (MQ=255)
     cTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGt                                                                                                                                <  2:301656/149‑1 (MQ=255)
        ttATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCg                                                                                                                             <  1:930222/149‑1 (MQ=255)
          atatGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGGGAATAAAAGGAAGGTGCGac                                                                                                                           >  2:792648/1‑148 (MQ=255)
                tggtctcCGCCAGGAGTGATTTTCGTATTATTGGCGGTCCCTAAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCg                                                                                                                     <  2:196120/146‑1 (MQ=255)
                tGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGACTAAGCg                                                                                                                     >  1:820541/1‑149 (MQ=255)
                 gTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCCTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCgg                                                                                                                    <  1:352909/149‑1 (MQ=255)
                 gTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTAAATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCgg                                                                                                                    <  1:220802/149‑1 (MQ=255)
                  ttCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACAACATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCggg                                                                                                                   >  2:296526/1‑149 (MQ=255)
                         ccAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAAttc                                                                                                            >  2:350970/1‑149 (MQ=255)
                           aGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTCGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAAttctt                                                                                                          <  2:718954/149‑1 (MQ=255)
                               ttgatcTTCGTTTTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCg                                                                                                      <  1:424880/149‑1 (MQ=255)
                                   ttttCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTg                                                                                                  >  1:785128/1‑149 (MQ=255)
                                    tttCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCCGCAGTTGATGCCAATGCTAGAGGTGCATGTAAGGCAGGTGCGAATAAGCGGAGAAAATCTTCTCGGCGGa                                                                                                 >  1:409333/1‑149 (MQ=255)
                                    tttCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGa                                                                                                 <  1:26277/149‑1 (MQ=255)
                                         tattatTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAACCCTTCCCGGCTGACtgcg                                                                                            >  2:867505/1‑146 (MQ=255)
                                                  cGGTCCCATATACCAGCTGGATGTCTTTTTCCCGGACAACACCCTCAACTGTCTCACCATACTCAGCAGTTGACGACGTTCACATGGTTGAGTAGATGCAAGACGCAAATAAGCGGGGAAATTCTTCTCGGCTGACGGAGtcagttcat                                                                                   <  2:676026/149‑1 (MQ=255)
                                                       ccAAAAACCAGCTGGATATCTTTTTCCAGGACACCACCATCAACCGTTTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATttctt                                                                              <  1:672884/149‑1 (MQ=255)
                                                         aTAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCa                                                                            >  2:59268/1‑149 (MQ=255)
                                                          ttactcAGCTGGTTATCTTGTTCCAGGACAACACCATCAACTGTCTCACCATCCTCAACGGTTGATGCCAATGCCATGGGTGAATATAAGGAAGGTGCGAAGAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCAt                                                                           <  2:416105/144‑1 (MQ=25)
                                                               cAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTg                                                                      <  2:639288/149‑1 (MQ=255)
                                                               cAGCTGGATATCTTTTTCCAGGACAACACCAGCAAATGTGTCACCAGACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTg                                                                      <  1:414457/149‑1 (MQ=255)
                                                                 gCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGATCAAGTCATTTCATTTCTTAATGTTTggg                                                                    >  2:419646/1‑147 (MQ=255)
                                                                  cTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGc                                                                   >  2:913540/1‑149 (MQ=255)
                                                                           tttttCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGAtttttt                                                          <  1:59268/149‑1 (MQ=255)
                                                                            ttttCCAGGACAACACCATCAACTGTCTCAACATACTCAACGGTTGATGCCAATGCTATGGGCGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTtc                                                         <  1:648203/149‑1 (MQ=255)
                                                                                ccAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTccc                                                     >  1:511947/1‑149 (MQ=255)
                                                                                  aGGACAACACCATCAACTGTCGCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGt                                                   <  2:783500/149‑1 (MQ=255)
                                                                                      caacaCCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGTGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAAt                                               <  1:893615/149‑1 (MQ=255)
                                                                                         cacCATCAACTGTCTCACCATACTCAACAGTTGATGTCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGcc                                            >  1:369448/1‑149 (MQ=255)
                                                                                          acCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCt                                           <  1:767441/149‑1 (MQ=255)
                                                                                            cATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTAAGTAATTTCATTTCTTCATGGTTGAGCCGATTTTTTCTCCCGTAAATGCCTTg                                         >  2:229954/1‑149 (MQ=255)
                                                                                                   tGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGc                                  <  1:878777/149‑1 (MQ=255)
                                                                                                                      gTTGATGCCAATGCTNTGGNTNCNTNTAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCTTAAATGCCTTGAATCAGCCTATTTAGACCGTttcttc               <  1:608381/149‑1 (MQ=255)
                                                                                                                               aaTGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTaa      >  2:375537/1‑149 (MQ=255)
                                                                                                                               aaTGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTaa      >  2:682121/1‑149 (MQ=255)
                                                                                                                                aTGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTAGTTAGACCGTTTCTTCGCCATTTAAg     <  1:476631/149‑1 (MQ=255)
                                                                                                                                 tGCCATGGGTGAATATAAGGAAGGTGCGGATAAGCGGGGAAATGCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAgg    <  2:318104/149‑1 (MQ=255)
                                                                                                                                   cTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCg  >  2:232342/1‑149 (MQ=40)
                                                                                                                                                   |                                                                                                                                    
AAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCG  >  CP009273/2282251‑2282530

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.