Predicted mutation
evidence seq id position mutation freq annotation gene description
JC JC CP009273 3,178,254 IS5 (+) +4 bp 40.0% 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)9 (0.340) 9/274 NT 42.1% noncoding (1195/1195 nt) IS5 repeat region
?CP009273 3178254 = 12 (0.470)coding (56/291 nt) yqiC BMFP family putative fusogenic protein
* ? CP009273 = 2283592NA (NA)7 (0.260) 7/276 NT 36.2% noncoding (1/1195 nt) IS5 repeat region
?CP009273 = 3178257 12 (0.460)coding (59/291 nt) yqiC BMFP family putative fusogenic protein

AAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTGTTAAAAGCATTGGGTACAGAAAATACCCATAGCTCCATACCCGGAG  >  CP009273/2283451‑2283735
                                                                                                                                             |                                                                                                                                               
taaaaaaTCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTGCCATAACGc                                                                                                                                          >  2:48037/4‑149 (MQ=21)
   acaaTCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATCCACTGTCGGCGAAGGTAAGTTGATCACTCATGATCAACCCTGTTCTATGGCTCCAGATGACAAACACGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGt                                                                                                                                       <  2:330212/147‑1 (MQ=255)
    aaaTCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGTAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTa                                                                                                                                      <  1:452121/149‑1 (MQ=38)
              tCTGGTCCGACGGCGCTTACTGCTGAATTCACTGCCGGCGAAGGCAAGTCGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGaa                                                                                                                            <  1:153857/149‑1 (MQ=255)
                tGGTCTGACGGCGCTTACTGCTGAATTCCCTGTCGTCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACa                                                                                                                          >  1:345366/1‑149 (MQ=255)
                  ggcggACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGa                                                                                                                        <  2:447736/145‑1 (MQ=255)
                       aCGGCGCTTACTGCTGAATTCACTGTCGGCGAAAGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTATAGCCAGCAGAACAGATATTg                                                                                                                   >  2:121759/1‑149 (MQ=255)
                                 ctgctgAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACa                                                                                                         >  1:292374/1‑149 (MQ=255)
                                  tgctgAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGTACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACaa                                                                                                        <  1:248077/149‑1 (MQ=255)
                                  tgctgAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACaa                                                                                                        >  1:467400/1‑149 (MQ=255)
                                         ttCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGaaa                                                                                                 <  1:43728/149‑1 (MQ=255)
                                                             gTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCAtt                                                                             <  2:226352/149‑1 (MQ=255)
                                                                       cATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCAt                                                                   >  2:234371/1‑149 (MQ=255)
                                                                             gaACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATg                                                             >  1:298316/1‑149 (MQ=255)
                                                                                     ttCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAaattaat                                                     >  1:33941/1‑149 (MQ=255)
                                                                                       cTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTg                                                   >  1:336130/1‑149 (MQ=255)
                                                                                            gCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTGtt                                              >  2:394943/1‑149 (MQ=255)
                                                                                              tCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATGGAGCTTTCATCCTATGAAATTAATTGCTGTTaa                                            >  1:397835/1‑149 (MQ=255)
                                                                                                cAGATGACAAACATGAGCCCATATCAGGGACTTGTTCGCACCTTCCAGAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTGTTaaaa                                          <  1:22266/149‑1 (MQ=255)
                                                                                                        aaacatttcggcATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTGTTAAAAGCATTggg                                  >  2:122560/12‑149 (MQ=255)
                                                                                                           cATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTGTTAAAAGCATTGGGTAc                               <  1:455638/149‑1 (MQ=255)
                                                                                                                  tcATATCAGGGACTTGTGCGCACCTTCCATAACGCTGTAGCCACCAGTACAGATAGTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGATTTCATCCTATCAACACTGCGGCTGTTAAAAGCGGTGGGTACAGAAAAt                        <  1:345374/149‑1 (MQ=255)
                                                                                                                    atatCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTGTTAAAAGCATTGGGTACAGAAAATAc                      <  1:235532/149‑1 (MQ=255)
                                                                                                                    atatCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTGTTAAAAGCATTGGGTACAGAAAATAc                      >  1:213802/1‑149 (MQ=255)
                                                                                                                         aGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTGTTAAAAGCATTGGGTACAGAAAATACCCATa                 >  2:84261/1‑149 (MQ=255)
                                                                                                                          cgcaCTTGTTCCCACCTTCCATAACGCTGGAGCCACCAGAACAGATATTGCGGAACGACACAGAGAAACAGAACCAGATTGATGCATTATGCTTTCAACCTATGAAACTAACCGCTCTTAACCGCAATCAACACACAAAATACCCATAg                <  2:345158/146‑1 (MQ=255)
                                                                                                                                 ttcccaaCCTTCCATAACCCGGTAGCCACAAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATACTATGAAATTAATTGCTGTTAAAAGCATTGGGTACAGAAAATACCCATAGCTCCATa         <  1:121759/143‑1 (MQ=255)
                                                                                                                                       aCCTTCCATAACGCTGTAGCCACGAGAACAGATATTGCTGAACGACCAAGAGAAACAGAACCAGAATGATGCATTGAGCTTTCAGCCTATGCAATTAATTGCTGTTCAAAGAATTGGGTACACGAAATACGCATAGATCCAAActctta   >  2:350301/1‑144 (MQ=255)
                                                                                                                                        ccTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTGTTAAAAGCATTGGGTACAGAAAATACCCATAGCTCCATACCCGGAg  <  1:122560/149‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                               
AAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTGTTAAAAGCATTGGGTACAGAAAATACCCATAGCTCCATACCCGGAG  >  CP009273/2283451‑2283735

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.