Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913_3_vch_tpiA 1,978,503 Δ776 bp insB1insA insB1, insA

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913_3_vch_tpiA 1978503 1979278 776 35 [0] [0] 35 insB1–insA insB1,insA

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913_3_vch_tpiA = 19785020 (0.000)32 (0.750) 26/274 0.8 100% noncoding (1978502/4641682 nt) repeat_region repeat region
?NC_000913_3_vch_tpiA 1979279 = 0 (0.000)noncoding (1979279/4641682 nt) repeat_region repeat region

TTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAA  >  NC_000913_3_vch_tpiA/1978366‑1978648
                                                                                                                                        |                                                                                                                                                  
tttACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACCCAAAACACTCAATTACTTAACATAAATTtataagtcatac                                                                                                                                        >  1:62895‑M1/1‑137 (MQ=255)
  tACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatactt                                                                                                                                      >  2:79254‑M1/1‑135 (MQ=255)
    cacTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatactttt                                                                                                                                    <  1:194549‑M1/149‑17 (MQ=255)
           cTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgg                                                                                                                             <  2:277408‑M1/149‑24 (MQ=255)
            tAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttggg                                                                                                                            >  2:441096‑M1/1‑125 (MQ=255)
            tAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttggg                                                                                                                            >  1:716287‑M1/1‑125 (MQ=255)
            tAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACAAACAACGGAGGGAGACAGCGATTTTCAATAATGCACGATACAGATCACACAAAACACTAAATTACTCAACATAAATGtaaaaatcatactcttgcttcggg                                                                                                                            >  2:51130‑M1/1‑125 (MQ=255)
                  ttttCCTAAATCGACGCAACTGTACTCGTCACTAAAAACAAATACAACGGAGGGGGGCTGCGATTTTCAATACTGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAACTGtataagtcatacttttcttttgcgtgtatt                                                                                                                      >  1:621341‑M1/1‑119 (MQ=255)
                   tttCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtattt                                                                                                                     <  1:79254‑M1/149‑32 (MQ=255)
                            tCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGAGATTTTCAATAATGCGTGATGCGGATCACACAAAACCCTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatcggtt                                                                                                            >  2:44417‑M1/1‑109 (MQ=255)
                                    aCTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtt                                                                                                    <  2:252932‑M1/149‑49 (MQ=255)
                                      tGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttt                                                                                                  <  1:77947‑M1/149‑51 (MQ=255)
                                           tCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATTtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttggct                                                                                             >  1:508766‑M1/1‑94 (MQ=255)
                                           tCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttggttgtatttcaatctgttaaaaagtttttcgct                                                                                             >  2:697363‑M1/1‑94 (MQ=255)
                                             gTCACTACACGCACATACAACGGAAGGGGGCTGCGATTTACAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctac                                                                                           <  1:282725‑M1/149‑58 (MQ=255)
                                                  tACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACATAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgataaaaagtttttcgctacgctag                                                                                      >  1:661465‑M1/1‑87 (MQ=255)
                                                           aTACAACGGAGGGGGGCTGCGATTTTCAATAATGCATAATGCAGGTCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaa                                                                             >  1:464003‑M1/1‑78 (MQ=255)
                                                                 cGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtgtaagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaa                                                                       <  1:749312‑M1/149‑78 (MQ=255)
                                                                            tGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatc                                                            >  2:425387‑M1/1‑61 (MQ=255)
                                                                             gCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcg                                                           >  1:204288‑M1/1‑60 (MQ=255)
                                                                                        taatGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgcatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgt                                                <  1:697363‑M1/149‑101 (MQ=255)
                                                                                         aatGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgtt                                               >  2:216424‑M1/1‑48 (MQ=255)
                                                                                                          acagaAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaaggttttcgccacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataa                              <  2:148027‑M1/149‑119 (MQ=255)
                                                                                                          acacaAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataa                              <  2:204288‑M1/149‑119 (MQ=255)
                                                                                                              aaaaCACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgcatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactca                          >  1:516783‑M1/1‑27 (MQ=255)
                                                                                                                   acTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgta                     >  1:588694‑M1/1‑22 (MQ=255)
                                                                                                                       aaTTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaa                                                                       <  2:61607‑M1/95‑78 (MQ=255)
                                                                                                                       aaTTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaa                                                                       >  1:61607‑M1/1‑18 (MQ=255)
                                                                                                                         ttACTTAACATAAACGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtcaataat               <  2:247327‑M1/149‑134 (MQ=255)
                                                                                                                           aCTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatca             <  1:463234‑M1/149‑136 (MQ=255)
                                                                                                                            cTTAACATAAATGtataagtcatacttttgttttgggtgtattttagtatgtttaacagttttttgctacgctagcgagcaaaaatgaatcaggaataatcgaaatggggtgttgcgcacagtcaaaatgactcagagtaaataatcat            >  2:330088‑M1/1‑13 (MQ=255)
                                                                                                                            cTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacatgaataatcgaaatgggatgttgcgcacagtcacaataactaaccgtaaataatcat            >  2:329857‑M1/1‑13 (MQ=255)
                                                                                                                                  aTAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgcta      >  1:355483‑M1/1‑7 (MQ=255)
                                                                                                                                  aTAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagttattcgctacgctagcaagcaaaaattaaaccggaataatcgaaatggcattttgcgcacagtcaaaataactcaccgtaaataatcatccgata      >  2:150872‑M1/1‑7 (MQ=255)
                                                                                                                                      aTGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgctataaa  >  1:409785‑M1/1‑3 (MQ=255)
                                                                                                                                        |                                                                                                                                                  
TTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAA  >  NC_000913_3_vch_tpiA/1978366‑1978648

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.