Predicted mutation
evidence seq id position mutation freq annotation gene description
JC JC CP009273 3,178,254 IS5 (+) +4 bp 54.3% 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)12 (0.840) 11/272 NT 59.3% noncoding (1195/1195 nt) IS5 repeat region
?CP009273 3178254 = 8 (0.580)coding (56/291 nt) yqiC BMFP family putative fusogenic protein
* ? CP009273 = 2283592NA (NA)7 (0.490) 7/274 NT 45.9% noncoding (1/1195 nt) IS5 repeat region
?CP009273 = 3178257 8 (0.570)coding (59/291 nt) yqiC BMFP family putative fusogenic protein

CCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCC  >  CP009273/2282249‑2282536
                                                                                                                                                     |                                                                                                                                          
ccAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCCACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATaag                                                                                                                                            >  1:180262/1‑150 (MQ=255)
              atGCTGTTCTCCGCCAGGATTGATTATCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATaa                                                                                                                              >  1:195805/1‑150 (MQ=255)
                   gctctcCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAATGTGCGAATAAGCggg                                                                                                                         <  1:239090/148‑1 (MQ=255)
                                                  ggCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGtcat                                                                                              <  2:200144/146‑1 (MQ=255)
                                                     ggTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGGCTCGGtcatttc                                                                                           >  2:62047/1‑146 (MQ=255)
                                                     ggTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAAAACCATCAACTGTCTCACCATCCTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGACTAAGCGGGGAAATTTTTCTCGGATAGCTCAGtcacttt                                                                                           >  2:96459/1‑145 (MQ=255)
                                                            tAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCTAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCGGACTCAGTCATTTCATTTCTTCATg                                                                                >  1:162019/1‑150 (MQ=255)
                                                                           tCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGAttttt                                                                 >  1:121221/1‑150 (MQ=255)
                                                                                                            ccATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTAACTCCGGCAGTTCAGTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAg                                <  1:127825/150‑1 (MQ=18)
                                                                                                                             tGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTGTCTTCGc                   >  2:163177/1‑146 (MQ=255)
                                                                                                                                          ggTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATcn  <  1:242010/150‑2 (MQ=25)
                                                                                                                                                     |                                                                                                                                          
CCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCC  >  CP009273/2282249‑2282536

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.