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

AAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCC  >  CP009273/2282251‑2282537
                                                                                                                                                   |                                                                                                                                           
aaaTTCTTTTATATGCTGTTCTCCGCCAGGATTGGTTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACGGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAAGATaagg                                                                                                                                            >  2:5959/1‑149 (MQ=255)
aaaTTCTTTTATATGCTGTTCTCCGACAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATAAaagg                                                                                                                                            >  2:5923/1‑149 (MQ=255)
       tttATATGCTGTGCGCCGCCAGGATTGAGTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCAACGACGGTCTCACCATCCTCAACCGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGc                                                                                                                                     <  1:460448/149‑1 (MQ=255)
                     tcCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGaa                                                                                                                       >  1:587773/1‑149 (MQ=255)
                        gCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAAtt                                                                                                                    >  1:110868/1‑149 (MQ=255)
                        gCCAGGATTGATTTTAGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTCATGCCAATGCTCTGAGGGAATATAAGGAAGGTGCGGATAAGCGGGGAAAtt                                                                                                                    >  2:437707/1‑149 (MQ=255)
                           aGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAAttctt                                                                                                                 <  2:98972/149‑1 (MQ=255)
                           aGGATTGATTTTCGTAGTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAAttctt                                                                                                                 <  2:99011/149‑1 (MQ=255)
                              gttgatcTTCGTATTATTGCCCGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCGTCAACTGTCTCACCCTACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTtctc                                                                                                              <  2:589116/148‑1 (MQ=39)
                              attgattTTCGTATTATTGGCGGGCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTtctc                                                                                                              <  2:253927/149‑1 (MQ=255)
                              attgatcTTCGTATTATTGGCGGTCCCCTAAACACCCTGGATATCTTCTTCCAGGACACCACCATCAACGGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATACGCGGGGAAATTCTtctc                                                                                                              <  1:444475/149‑1 (MQ=255)
                                  attTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCt                                                                                                          <  1:293182/149‑1 (MQ=255)
                                        gTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCCATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCa                                                                                                    <  2:571116/149‑1 (MQ=255)
                                          attattGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGAGAACTCAGt                                                                                                  >  1:297351/1‑149 (MQ=255)
                                            tattGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCGCACCATACTCCACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGCAATTCTTCTCGGCTGACTCAGtca                                                                                                >  1:137075/1‑149 (MQ=255)
                                                        cATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATttcttc                                                                                    >  1:321353/1‑149 (MQ=255)
                                                            acccAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGt                                                                                <  2:374141/147‑1 (MQ=255)
                                                              ccAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCACCAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGttt                                                                              <  1:331484/149‑1 (MQ=255)
                                                              ccAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCAGACTCAACAGTTGATGCCAATGCTATGGGTGAATAATATGAACGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGttt                                                                              <  1:132534/149‑1 (MQ=37)
                                                                aGCTGGATATCTTTTTCCAGGACAACGCCAGTAACTGTCTCACCATTCTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGa                                                                            <  2:40594/149‑1 (MQ=255)
                                                                aGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGACTCAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTgg                                                                            >  2:381408/1‑148 (MQ=255)
                                                                aGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCGCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGa                                                                            <  2:19419/149‑1 (MQ=255)
                                                                           tttttCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGAtttttt                                                                 <  1:554421/149‑1 (MQ=255)
                                                                              ttCCAGGACAACACCATCAACTGTCTCCCCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTtctc                                                              <  1:517145/149‑1 (MQ=255)
                                                                                 cAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCg                                                           >  2:529026/1‑149 (MQ=255)
                                                                                 cAGGACAACACCATCAACTGTCTCACCATACTCAAAAGTTGATGCCAATGCTATGGGTCAATATAACGACGGAACGAATAAGCTGGGAAATTTTTCTCGGCTGACTCAGTCAATTCATTTCTTCATGTTTGGGCCGACTTTTTCTCCCg                                                           >  2:420777/1‑149 (MQ=255)
                                                                                        acacCATTAACTGGCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGc                                                    <  2:684582/149‑1 (MQ=255)
                                                                                            cATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTCTGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTg                                                >  2:80188/1‑149 (MQ=255)
                                                                                              tCAACTGTCTCACCATAATCTACAGGTGATTCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGaa                                              <  1:156940/149‑1 (MQ=255)
                                                                                                    gTCTCACCATACTCAACAGTTGATGACAATACTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGcc                                        >  1:358701/1‑149 (MQ=255)
                                                                                                       tcACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTAt                                     >  2:247324/1‑149 (MQ=255)
                                                                                                       tcACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTAt                                     >  1:524810/1‑149 (MQ=255)
                                                                                                         aCCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGACTACGCGGGGAAATTCTTCTCGGCTGACCCCGTCATTTCATTTCTTCATGTTTGAGCCGCTTTTTTCTCCCGTAAATGCCTTGATTCAGCCTCttt                                   >  1:8723/1‑149 (MQ=18)
                                                                                                         aCCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTGttt                                   >  2:177066/1‑149 (MQ=255)
                                                                                                          ccATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTa                                  <  1:110204/149‑1 (MQ=255)
                                                                                                            aTACTCAACAGTTGATTCCAATGCTATGGGTGAATATAAGGAAGGTGCGACTAAGCGGGGAAATTCTTCTCGGCTGCCTCAGTTATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTAGTTAGa                                <  1:310281/149‑1 (MQ=17)
                                                                                                             tACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGTAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGAc                               <  1:492810/149‑1 (MQ=255)
                                                                                                              aCTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGAcc                              >  1:210680/1‑149 (MQ=255)
                                                                                                                 cAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACGCAGTCATTTCATTTCTTCACGTCTGACACGCTACTTCCTTCCGTCAGTGCCTTGAATCACCATATTTAGACCGtg                           >  2:72575/1‑148 (MQ=17)
                                                                                                                   acggTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTttc                         <  2:107177/146‑1 (MQ=255)
                                                                                                                   aCAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTCGACCGTttc                         >  2:586184/1‑149 (MQ=255)
                                                                                                                       ttGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTTCCGTTAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCg                     >  2:214859/1‑149 (MQ=255)
                                                                                                                        tGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGc                    <  2:17396/149‑1 (MQ=255)
                                                                                                                          aTGCCAAGGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCa                  >  2:463700/1‑149 (MQ=255)
                                                                                                                            gCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCAtt                >  1:199024/1‑149 (MQ=255)
                                                                                                                             ccAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCAttt               <  1:581174/149‑1 (MQ=255)
                                                                                                                                  gCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGc          <  1:214997/149‑1 (MQ=40)
                                                                                                                                  gCTATGGGTGAATATAAGGAAGGTGAGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGc          >  1:347658/1‑149 (MQ=255)
                                                                                                                                    tATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCAGGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTTAATCAGCCTTTTTAGACCGTTTCTTCGCCATTTAAGGCgg        >  2:74424/1‑148 (MQ=255)
                                                                                                                                      tGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTa      >  2:147443/1‑149 (MQ=38)
                                                                                                                                       gggTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTAt     <  2:682099/149‑1 (MQ=37)
                                                                                                                                       gggTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTAt     <  2:682121/149‑1 (MQ=37)
                                                                                                                                          tGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTCAAGGCGTTATccc  <  2:358701/149‑1 (MQ=21)
                                                                                                                                          ggAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATccc  <  2:52531/148‑1 (MQ=18)
                                                                                                                                                   |                                                                                                                                           
AAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCC  >  CP009273/2282251‑2282537

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.