Predicted mutation
evidence seq id position mutation freq annotation gene description
JC JC CP009273 3,178,254 IS5 (+) +4 bp 15.8% 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)5 (0.150) 4/268 NT 16.8% noncoding (1195/1195 nt) IS5 repeat region
?CP009273 3178254 = 24 (0.720)coding (56/291 nt) yqiC BMFP family putative fusogenic protein
* ? CP009273 = 2283592NA (NA)4 (0.120) 4/270 NT 13.9% noncoding (1/1195 nt) IS5 repeat region
?CP009273 = 3178257 24 (0.720)coding (59/291 nt) yqiC BMFP family putative fusogenic protein

ATAGCTTCACATAATTCTGGTTTATGACTTACCCTTATCGCACTACAATGGCACTCAACACCTATCAGTACAGGGAAACTACAATGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGGAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGTCCTGTTACGGACTCGTGAAAAACTGGCGCTGC  >  CP009273/3178116‑3178403
                                                                                                                                             |                                                                                                                                                  
aagcgatcatctacttCAGGTTTATGCCGTACCCGTATCGCACCACAATGGCACTCAACACCTAGCATTACAGGGCAACTACAACGATTGCTCCGAAAAAAATTGAGCGGATCGCTCGCCAGGTTCACGAATCAATGCCTaaa                                                                                                                                                   <  1:2343/129‑1 (MQ=255)
  aGCTTCACATAATTCTGGTTTATGACTTACCCTTATCGCACTACAATGGCACTCAACACCTATCAGTACTGGGAAACTACAATGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCCGGTTCACGAATCAATGCCTAAgga                                                                                                                                                 <  1:18012‑M2/143‑4 (MQ=255)
     ttCACATAATTCTGGTTTATGACTTACCCTTATCGCACTACAATGGCACTCAACACCTATCAGTACAGGGAAACTACAATGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAAt                                                                                                                                              <  1:255710/143‑1 (MQ=255)
         caTAATCCTGGTTTATGACTTACCCTTATCGCACTACAATGGCACTCAACACCTATCAGTACAGGGAAACTACCATGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAgg                                                                                                                                          <  1:154291/143‑1 (MQ=255)
               tCTGGTTTATGACTTACCCTTATCGCACTACAATGGCACTCAACACCTATCAGTACAGGGAAACTACAATGATTGACCCGGAAAAAATTGAGCAAATCGCTCGCCAGGTTAACGAAACAAGGCCTAAAGGAACCAGGGAggcc                                                                                                                                    >  1:149597/1‑140 (MQ=255)
                   gTTTATGACTTACCCTTATCGCACTACAATGGCACTCAACACCTATCAGTACAGGGAAACTACAATGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGGTCGGGGAAGAtg                                                                                                                          >  2:128450/1‑149 (MQ=255)
                    tttATGACTTACCCTTATCGCACTACAATGGCACTCAACACCTATCAGTACAGGGAAACTACAATGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGATCGGGGa                                                                                                                               >  1:81917/1‑143 (MQ=255)
                      aaacggtctaaataggctgattcaaggcatttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAAAGGAATCAGGGAGTTCGGGGAAGATGTgg                                                                                                                       <  2:117443‑M2/33‑1 (MQ=255)
                           cTTACCCTTATCGCACTACAATGGCACTCAACACCTATCAGTACAGGGAAACTACAATGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGGAAGAtgtg                                                                                                                        <  1:381868/143‑1 (MQ=255)
                                  ttATCGCACTACAATGGCACTCAACACCTATCAGTACAGGGAAACTACAATGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAggaaggtgcgaacaagtccctgatatgcgatcatg                                                                                                                 <  1:184104‑M2/143‑36 (MQ=255)
                                    aTCGCACTACAATGGCACTCAACACCTATCAGTACAGGGAAACTACAATGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGGTCGGGGAAGATGTGGAGaagaga                                                                                                               >  1:288355/1‑139 (MQ=255)
                                      cGCACTACAATGGCACTCAACACCTATCAGTACAGGGAAACTACAATGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGGAAGATGTGGAGCAAAAAATCCGCCa                                                                                                       >  2:494094/1‑149 (MQ=255)
                                         attctagccatttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAAAGGAATCAGGGAGTTCGGGGAAGATGTGGAGAAAAAAATCCg                                                                                                          <  1:110133‑M2/46‑1 (MQ=255)
                                          cgcaactgtaaCTCAACACCTATCAGTACAGGGAAACTACAGTGAATGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGGAAGATGTGGAGAAAAAAATCCGc                                                                                                         <  1:87590/133‑1 (MQ=255)
                                            aaaagcgcactcaACATCTAGCAGTACAGGGAAACTACAATGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatct                                                                                                       >  1:112029‑M2/7‑98 (MQ=255)
                                                      tCAACACCTATCAGTACAGGGAAACTACAATGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGGAAGATGTGGAGAAAAAAATCCGCCAAACCCTACaa                                                                                             >  1:159365/1‑143 (MQ=255)
                                                               aTCAGTACAGGGAAACTACAATGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGGAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTg                                                                                    >  1:576536/1‑143 (MQ=255)
                                                               aTCAGTACAGGGAAACTACAATGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGGAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTg                                                                                    >  1:363906/1‑143 (MQ=255)
                                                                     aCAGGGAAACTACAATGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGGAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGAc                                                                        <  2:322389/149‑1 (MQ=255)
                                                                      cAGGGAACCTACAATGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGACTCAATGCCTAAAGGAATCAGGGAGTTCGGGGAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGAcc                                                                       <  2:599164/149‑1 (MQ=255)
                                                                      cAGGGAAACTACAATGATTGACCCGAAAAAAACTGAGCAAATCGCTCGCCAGGTTCACGAAGCAATGCCTAAggaaggtgcgaacaagtcccggatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtcat                                                                       >  2:20662‑M2/1‑72 (MQ=255)
                                                                           aaaCTACAATGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGGAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGAc                                                                        <  1:408597/143‑1 (MQ=255)
                                                                               tACAATGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGACTAAAGGAATCAGGGAGTTCGGGGAAGATGTGGAGAAAAAAATCCGCCAAACCCTCCAAGCGCAGCTGACGCGACTCGCCGTTGTAAGcc                                                              >  2:582000/1‑149 (MQ=255)
                                                                                        gaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAAAGGAATCAGGGAGTTCGGGGAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTg                                                           >  1:545350‑M2/51‑143 (MQ=255)
                                                                                         gACCCGAAAAAAACTGAGCAAATCGCTCGCCAGGTTCACGAAACAATGCCTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctgtagccatagaacagggttcatcatgagtcatcaacttaccttcgccgaca                                                    >  2:206893‑M2/1‑53 (MQ=255)
                                                                                          aatgactgagtcagccgcgaagaatttcccccgttatcctcaccgtccCTAAAGGAATCAGGGAGTTCGGGGAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTc                                                   <  2:118174‑M2/101‑1 (MQ=255)
                                                                                          aatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAACGGAATCAGGGAGTTCGGGGAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTc                                                   <  2:52122‑M2/101‑1 (MQ=255)
                                                                                                aaaaaaTTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGGAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTc                                             >  2:546673/1‑149 (MQ=255)
                                                                                                 aaaaaTTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGGAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGACGAGTTCg                                                  >  1:426884/1‑143 (MQ=255)
                                                                                                               tcgatggcCGGGTGCACGCATCAAGGCCCAAAGGAATCGGGGAGTTCGGGGAAGATGTGGAGAAAAAGATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCaa                                    <  1:77192/137‑1 (MQ=255)
                                                                                                                      ccAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGGAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGTCCTGTTACGGa                       <  2:91682/149‑1 (MQ=255)
                                                                                                                       cAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGGAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGTCCTGtt                            >  1:187808/1‑143 (MQ=255)
                                                                                                                          gTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGGAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGTCCTGTTaag                         >  1:576809/1‑141 (MQ=255)
                                                                                                                              aCGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGGAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGTCCTGTTACGGACt                     >  1:130136/1‑143 (MQ=255)
                                                                                                                              aCGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGGAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGTCCTGTTACGGACt                     >  1:130247/1‑143 (MQ=255)
                                                                                                                                         ccTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtcatcaacttaccttcgccgacagtgaattcagcagtaagcgccgtcagaccggaaaagagagtttcttgt    <  2:216742‑M2/149‑145 (MQ=255)
                                                                                                                                           tAAAGGAATCAGGGAGTTCGGGGAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGTCCTGTTACGGACTCGTGAAAAACTGGCgctgc  >  2:291287/1‑149 (MQ=255)
                                                                                                                                             |                                                                                                                                                  
ATAGCTTCACATAATTCTGGTTTATGACTTACCCTTATCGCACTACAATGGCACTCAACACCTATCAGTACAGGGAAACTACAATGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGGAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGTCCTGTTACGGACTCGTGAAAAACTGGCGCTGC  >  CP009273/3178116‑3178403

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