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

ACGCATAGCTTCACATAATTCTGGTTTATGACTTACCCTTATCGCACTACAATGGCACTCAACACCTATCAGTACAGGGAAACTACAATGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGGAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGTCCTGTTACGGACTCGTGAAAAACTGGCG  >  CP009273/3178112‑3178399
                                                                                                                                              |                                                                                                                                                 
aCGCATAGCTTCACATAATTCTGGTTTATGACTTACCCTTATCGCACTACAATGGCACTCAACACCTATCAGTACAGGGAAACTACAATGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTaa                                                                                                                                                >  2:72182/1‑146 (MQ=255)
          gccttaaatggcgaagaaacggtctaaataggctgattcaaggcatttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAAAGGAATCAgg                                                                                                                                      <  2:7529‑M2/14‑1 (MQ=255)
                 atggcgaagaaacggtctaaataggctgattcaaggcatttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAAAGGAATCAGGGAGTTCGGGGa                                                                                                                           >  1:15850‑M2/126‑150 (MQ=255)
                       gTTTATGACTTACCCTTATCGCACTACAATGGCACTCAACACCTATCAGTACAGGGAAACTACAATGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAg                                                                                                                         >  2:232811/1‑146 (MQ=255)
                           aacggtctaaataggctgattcaaagcatgtacggaagaaaaaatcggctcaaagatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAAAGGAATCAGGGAGTTCGGGGAAGATGTGgag                                                                                                                 >  1:145736‑M2/116‑150 (MQ=255)
                             cggtctaaataggcagattcaaggcatttacgggagcatcaatagggtcaaccatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAAAGGAATCAGGGAGTTCGGGGAAGATGTGGAGaa                                                                                                               <  1:115206‑M2/37‑1 (MQ=255)
                                    aataggctgattcaaggcatttacgggagaaaaaatcgactcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAAAGGAATCAGGGAGTTCGGGGAAGATGTGGAGaaaaa                                                                                                            <  2:80597‑M2/40‑1 (MQ=255)
                                      ttCTCTCACTACAATGGCACTCAACACCTCTCAGTACAGGGAAACTACAATGATTGACCCGAAAAAAATTGAGCAAATCGCCCGCCAGGTTCCCGAATCAATGCCTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtca                                                                                                      >  1:5057‑M2/1‑108 (MQ=255)
                                         gctgattcaaggcattgacgggagaaaaaatcggcccaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAAAGGAATCAGGGAGTTCGGGGAAGATGTGGAGAAAAAAATcc                                                                                                       <  2:89178‑M2/45‑1 (MQ=255)
                                         gccgaatcaaggcatttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAAAGGAATCAGGGAGTTCGGGGAAGATGTGGAGAAAAAAATCCGCCa                                                                                                   <  1:78954‑M2/49‑1 (MQ=255)
                                                caaggcatttacgggagaaaaaatcagctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAAAGGAATCAGGGAGTTCGGGGAAGATGTGGAGATAAAAATCCGCCAaacccgt                                                                                            >  1:195446‑M2/95‑148 (MQ=255)
                                                         aacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAAAGGAATCCGGGAGTTCGGGGGAGATGTGGAGACAAAAATCCGCCAAACCCTACAAGCGCAg                                                                                   >  1:184651‑M2/86‑150 (MQ=255)
                                                          tCAACACCTATTAGTACAGGGAAACTACAATGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGc                                                                                  <  1:223154/150‑1 (MQ=255)
                                                             acagctgacagtaaaggGAAGCGACAATGAATGACCCGAAAAAAATTGAGCAAAGCGCTCGCCAGGGTCACGAATCAATGCCTAAggaaggtgcgaacaagtccctgatatgagatcctgtttgtcatctggagccatagaacagg                                                                                   >  2:26846‑M2/9‑85 (MQ=255)
                                                                 cTATCAGTACAGGGAAACTACAATGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCAAAGGTCAGGAATCAATGCCTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatca                                                                           >  1:192376‑M2/1‑81 (MQ=255)
                                                                      atcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAAAGGAATCAGGGAGTTCGGGGAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCg                                                                      <  1:84019‑M2/78‑1 (MQ=255)
                                                                         aCAGGGAAACTACAATGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGAcc                                                                   <  1:82285/150‑1 (MQ=255)
                                                                          cAGGGAAACTACAATGATTGACCCGTAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCAAAACCCTACAAGCGCAGCTGACGCGCCTCGACCt                                                                  <  1:167719/150‑1 (MQ=255)
                                                                              gAAACTACAATGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCt                                                                  <  2:67231/146‑1 (MQ=255)
                                                                                          aTTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAACAATCCGCCAACCCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAGGCCGTGa                                                      >  2:69820/1‑146 (MQ=255)
                                                                                             gACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAggaaggtgcgaacaagtgcctgatatgagatcatgtttgtcaccaggagccatagaacagggtgcatcatgagtcatccacttaccttcgccg                                                   <  2:127190‑M2/146‑94 (MQ=255)
                                                                                                    aaaaaaTTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGAAATCCGGAAGTTCGGGAAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTcc                                        >  1:117313/1‑150 (MQ=255)
                                                                                                                    cgctcgcCAGGTTCACGAATCAATGCCTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtcaacagcttaccttccccgaccgtgaattcagcagtaagcgccgtc                        <  1:245863‑M2/150‑121 (MQ=255)
                                                                                                                          ccAGGTTCTCGAAGCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAACCCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGTCCTGTTAc                      <  2:174557/146‑1 (MQ=255)
                                                                                                                            cgcttattcgccccttccCTAAAGGAATGAGGGAGTTCGGGGAAGATGAGCAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGTCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGTCCTGTTACgg                    <  2:127249‑M2/128‑1 (MQ=255)
                                                                                                                               ttCAAGAATCAATGCCTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtcatcaacttaccttcgccgacagtgaattcagcagtaagcgccgtcagaccagaaaa             >  1:236021‑M2/1‑19 (MQ=255)
                                                                                                                                  ttcgcaccttccCTAAAGGAATCAGGGAGTTCGGGGAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGTCCTGTTACGGACTCGTGaaa          >  1:208442‑M2/13‑150 (MQ=255)
                                                                                                                                         aaTGCCTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtcatcaacttaccttcgccgacagtgaattcagcagtaagcgccgtcagaccagaaaagagattttct   >  1:241289‑M2/1‑9 (MQ=255)
                                                                                                                                              cTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCACCTCGACCTTGTAACCCGTGAAGAGTTCGACGTCCAAACGCCAGTCCTGTTACGGACCCGGGAACACCTGgcc  >  2:98758/1‑145 (MQ=255)
                                                                                                                                              |                                                                                                                                                 
ACGCATAGCTTCACATAATTCTGGTTTATGACTTACCCTTATCGCACTACAATGGCACTCAACACCTATCAGTACAGGGAAACTACAATGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGGAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGTCCTGTTACGGACTCGTGAAAAACTGGCG  >  CP009273/3178112‑3178399

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.