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)77 (1.120) 49/288 0.5 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)100 (1.440) 59/290 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

CAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCC  >  CP009273/2282250‑2282538
                                                                                                                                                    |                                                                                                                                            
cAAATTCTTTTATATGCTGTTCTCCGCCAGGGTTGATTTTCGTATTATTGGCGGTCCCATAACCCAGCTGAATCTCTTTTTCCAGAACAACACCATCAACTGTCTCACCATACTCAACAGTTACTGCCATTGCGTTGGGTGACTATaag                                                                                                                                              >  2:773778/1‑149 (MQ=255)
cAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATaag                                                                                                                                              <  2:258169/149‑1 (MQ=255)
cAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCAACAACTTATCCGCCATACTCAACAGTTGATGGCAAGGCTATGGGTGAATATaag                                                                                                                                              >  2:1153772/1‑149 (MQ=255)
   tttCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATATACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCCAATTTTTATTCCAATGCTATGGGTGAATATaaggaa                                                                                                                                           <  2:583380/148‑1 (MQ=255)
      cTTTTATATGCTGTGCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATaag                                                                                                                                              <  1:509560/143‑1 (MQ=255)
                cTGTTCTCGGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAACCCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGa                                                                                                                                    <  1:1042109/143‑1 (MQ=255)
                 tGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGaa                                                                                                                                   >  1:773456/1‑143 (MQ=255)
                  gTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAAt                                                                                                                                  <  1:800979/143‑1 (MQ=255)
                          ccAGGATTGATTTTCGTATTATCGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCACGACAAAACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAAttc                                                                                                                    <  2:1159721/149‑1 (MQ=255)
                            aGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCCACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAAttctt                                                                                                                  <  2:1130721/149‑1 (MQ=255)
                                ttgattTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGGAAttc                                                                                                                    >  1:105710/1‑143 (MQ=255)
                                 tgattTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCgg                                                                                                             >  2:342878/1‑149 (MQ=255)
                                   attTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCt                                                                                                           >  2:1043876/1‑149 (MQ=255)
                                   attTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCt                                                                                                           <  2:339729/149‑1 (MQ=255)
                                   attTTCCTATTCTTGGCGGTCCCATAAAACAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCt                                                                                                           <  2:1094082/149‑1 (MQ=255)
                                     tttCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGa                                                                                                         <  2:1036790/149‑1 (MQ=255)
                                         gTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCa                                                                                                     <  2:474428/149‑1 (MQ=255)
                                               ttGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACGGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCa                                                                                                     <  1:959820/143‑1 (MQ=255)
                                                   cGGTTCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCAGCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGtcatttcat                                                                                           <  2:530458/149‑1 (MQ=255)
                                                         cATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGGTGCCAATGCTATGGGGGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATttcctc                                                                                     >  2:173691/1‑149 (MQ=255)
                                                           tAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTcatttcattt                                                                                         <  1:664109/143‑1 (MQ=255)
                                                               ccAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGttt                                                                               <  2:713064/149‑1 (MQ=255)
                                                                 aGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGGAATTCTTCTCGGCTGACTCGGTACTTTCATTTCTTCAt                                                                                   >  1:287513/1‑143 (MQ=255)
                                                                 aGCGGGATATCTTTTTCCAGGACAACCCCATCACCTGGCCCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGa                                                                             <  2:516074/149‑1 (MQ=25)
                                                                     ggATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGttt                                                                               <  1:1008508/143‑1 (MQ=255)
                                                                     ggATATCTTTTTCAAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTCCTCGGCTGACTCAGTCATTTCATTTCTGCATGTTTGAGCCg                                                                         <  2:848537/149‑1 (MQ=255)
                                                                            tttttCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGa                                                                        >  1:159744/1‑143 (MQ=255)
                                                                             ttttCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCAGGTTTGAGCCggt                                                                       >  1:163243/1‑141 (MQ=255)
                                                                                 ccAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTccc                                                             <  2:1135805/149‑1 (MQ=255)
                                                                                     gACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCTATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTaa                                                         <  2:524468/149‑1 (MQ=255)
                                                                                        aacaCCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATg                                                      >  2:957138/1‑149 (MQ=255)
                                                                                        aacaCCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATg                                                      >  2:225367/1‑149 (MQ=255)
                                                                                             cATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTg                                                 >  2:1014772/1‑149 (MQ=255)
                                                                                                 aaCTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCt                                                   >  1:746388/1‑143 (MQ=255)
                                                                                                  aCTGTCTCGCCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCtt                                                  <  1:719918/143‑1 (MQ=255)
                                                                                                  aCTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCtt                                                  <  1:720195/143‑1 (MQ=255)
                                                                                                  aCTGTCTCACCATACTCAACAGTTGAAGCCAATGCTATGGGTGAATGTAAGGAAGGTGCGAATAAGCGGGGCAATTCTTCTCGGCTGACTCCGTCATTGCATTTCTTCATATTTGAGCGGATTTTTTCTCCCGTAAGTGCCtc                                                  >  1:643483/1‑142 (MQ=18)
                                                                                                   cTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAg                                           <  2:919091/149‑1 (MQ=255)
                                                                                                          aCCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATTCCTTGAATCAGCCTAttt                                    >  2:643655/1‑149 (MQ=255)
                                                                                                             aaaaacAACAGTGTATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGa                                 <  2:510758/144‑1 (MQ=255)
                                                                                                             aTACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTa                                       <  1:1064543/143‑1 (MQ=255)
                                                                                                                cTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTACTCGGCTGACTAAGTAATTTCATTTCTCCATGTTTGAGCCGATTTTTTCTCCCGTAAATTCCTTGAATCAGCCTAttt                                    >  1:568741/1‑143 (MQ=18)
                                                                                                                 tCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTa                                   >  1:757197/1‑143 (MQ=255)
                                                                                                                 tCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCAGAGTTTTTCTCCCCTAAATGCCTTGAATCAGCCTATTTAGACCGt                             >  2:222370/1‑149 (MQ=255)
                                                                                                                  cAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTAAAGAAATTAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGtt                            >  2:860772/1‑149 (MQ=17)
                                                                                                                   aaCAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGttt                           >  2:541271/1‑149 (MQ=255)
                                                                                                                       gTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTttcttc                       >  2:640913/1‑149 (MQ=255)
                                                                                                                          gATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCACTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTttc                          <  1:372009/143‑1 (MQ=255)
                                                                                                                           aTGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCa                   <  2:694793/149‑1 (MQ=255)
                                                                                                                             gCCAATGCTATGGGTGAATATAAGGAAGGTGCTAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTttcttc                       <  1:62219/143‑1 (MQ=255)
                                                                                                                             gCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTttcttc                       >  1:394086/1‑143 (MQ=255)
                                                                                                                              ccATTGCTATGGGTGAATAGAAGGAAGGTTCGAATAAGCGGGGAAATTCTGCTCGGCTGACTCAGTCATTTCAGTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTGAGGCCGTTTCTTCg                      <  1:266633/143‑1 (MQ=14)
                                                                                                                              ccAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCg                      <  1:1072341/143‑1 (MQ=255)
                                                                                                                                 aTGCTATGGGTGAATATAAGGAAGGCGCGAATAAGCAGGTAAATTCTTCTCGGCTGACTCAGCCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAg             <  2:804128/149‑1 (MQ=255)
                                                                                                                                    cTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAAGTCTTCGCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCg          <  2:517417/149‑1 (MQ=255)
                                                                                                                                    cTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGAACATTCTTCGCTTCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCg          <  2:517730/149‑1 (MQ=14)
                                                                                                                                        gggTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAgg            >  1:1002598/1‑143 (MQ=37)
                                                                                                                                            gAATATAAGGAAGGTGCGAATAAACGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTCTTTGGACCGTTTCTCCGCCATTTAAGGCGTTATcccc  >  2:659903/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                            
CAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCC  >  CP009273/2282250‑2282538

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.