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

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

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913_3_hsa_tpiA = 19785020 (0.000)38 (0.850) 37/262 0.3 100% intergenic (‑305/+16) flhD/insB1 flagellar class II regulon transcriptional activator, with FlhC/IS1 transposase B
?NC_000913_3_hsa_tpiA 1979279 = 0 (0.000)intergenic (‑64/‑474) insA/uspC IS1 repressor TnpA/universal stress protein

AAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCA  >  NC_000913_3_hsa_tpiA/1978357‑1978647
                                                                                                                                                 |                                                                                                                                                 
aaaTGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtat                                                                                                                                                <  1:341015‑M1/149‑4 (MQ=255)
   gtggtcTTTACGCTTATCTAAGATTTTTCATAAATCGACGCAACGTGACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATGACACAAAACACTCAATTACTTAACATAAATGtataag                                                                                                                                             <  2:748727‑M1/147‑7 (MQ=255)
              acgtatcTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACGACACGCACGTACAACGGAGGGGGGTTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAACTGtataagtcatacttttg                                                                                                                                  <  2:711000‑M1/146‑18 (MQ=255)
                ttATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGgataagtcatactttggtt                                                                                                                                >  1:177029‑M1/1‑130 (MQ=255)
                      aaGATTTTTCCTAAAGCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggt                                                                                                                          <  1:487793‑M1/149‑26 (MQ=255)
                       aGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttcggtg                                                                                                                         >  2:784670‑M1/1‑123 (MQ=255)
                         aTTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgg                                                                                                                            >  2:763124‑M1/1‑121 (MQ=255)
                         aTTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgg                                                                                                                            <  1:763124‑M1/144‑24 (MQ=255)
                                        aCGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaa                                                                                                        >  1:117171‑M1/1‑106 (MQ=255)
                                          gCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACACAATTACTTGACATAAATTcataagtcatacttttgttttgggtgtatttcaatctgttaaaaa                                                                                                      <  2:562911‑M1/149‑46 (MQ=255)
                                              cTGTACTCGTCCCTACACGCACATACAACGGAGGGGGGCGGCGATTTCCAATAATGCGTGATGCATATCACACAAAACACTCAATTACTTAACATAAATGtataattcatacttttgttttgggtgtatttcaatctgttaaaaagttt                                                                                                  <  1:292062‑M1/149‑50 (MQ=255)
                                                 tACTCTTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttc                                                                                               <  1:465919‑M1/149‑53 (MQ=255)
                                                    tCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcactctgttaaaaagtttttcgct                                                                                            >  1:705548‑M1/1‑94 (MQ=255)
                                                    tCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgct                                                                                            >  1:482107‑M1/1‑94 (MQ=255)
                                                      gTCACTCCACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatct                                                                                                              >  2:399979‑M1/1‑92 (MQ=255)
                                                      gTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatct                                                                                                              <  1:399979‑M1/129‑38 (MQ=255)
                                                          cTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgcta                                                                                      <  2:117171‑M1/149‑62 (MQ=255)
                                                                  acaTACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaa                                                                              >  1:632727‑M1/1‑80 (MQ=255)
                                                                          cGGAGGGAGGCTGCTATTTTCAATAATTCGGGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaa                                                                      <  2:836772‑M1/149‑78 (MQ=255)
                                                                                 gggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaat                                                               >  2:222016‑M1/1‑65 (MQ=255)
                                                                                           tttCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtca                                                                                                                                          <  1:14300‑M1/64‑10 (MQ=255)
                                                                                           tttCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtca                                                                                                                                          >  2:14300‑M1/1‑55 (MQ=255)
                                                                                              cAATAATGCGTGATGTATATCACAAAAAAAACTCAATTAATTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatggga                                                  <  1:137851‑M1/149‑98 (MQ=255)
                                                                                               aataatGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggat                                                 >  1:392338‑M1/1‑51 (MQ=255)
                                                                                                   atGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatggggtgttg                                             >  1:599557‑M1/1‑47 (MQ=255)
                                                                                                       gTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataattcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgca                                         <  1:700287‑M1/149‑107 (MQ=255)
                                                                                                         gATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcaca                                       <  1:222016‑M1/149‑109 (MQ=255)
                                                                                                                aTCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataa                                                             >  1:820134‑M1/1‑34 (MQ=255)
                                                                                                                aTCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataa                                                             <  2:820134‑M1/120‑87 (MQ=255)
                                                                                                                   acacaAAATACTTCATTACGTAACATAAATGtataagtcatacgtgtgttttgggtgtatttcaatctgttaaaaagttttgcgctacgctagcaagcaaaaatgaaacaggaataagcgaaatgggatgttgcgcacagtcaaaataa                             <  2:170844‑M1/149‑119 (MQ=255)
                                                                                                                     acaAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaaatttttcgatacgctaataagcaaaaatcatacaggaataatcgaaatgggaggttgcgcacagtcaaagtaact                           >  2:428804‑M1/1‑29 (MQ=255)
                                                                                                                        aaaCACGCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcac                        <  1:113803‑M1/149‑124 (MQ=255)
                                                                                                                         aaCACTTAATTATTTAACATAAATGtattagtcatacttttgttttgggtgtatttcaatctgttgacaagtttttcgctacgctagcaagcaaaaatgatacaggaatagtcgaaatgggatgttgggcacattcaaaataactcacc                       <  2:168358‑M1/149‑125 (MQ=255)
                                                                                                                             cTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgtcgcgcacagtcaaaataactcaccgtaa                   <  2:632727‑M1/149‑129 (MQ=255)
                                                                                                                              tCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaa                  >  2:223260‑M1/1‑20 (MQ=255)
                                                                                                                                aaTTACTTAACATAAATGtataagtcatacttttgtgttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactacccgtaaata                >  1:386742‑M1/1‑18 (MQ=255)
                                                                                                                                       tAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatct         >  1:222424‑M1/1‑11 (MQ=255)
                                                                                                                                        aaCATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctg        >  2:781947‑M1/1‑10 (MQ=255)
                                                                                                                                              aaTGtagaagtcatacttttgttttgggtgtagttttagctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgctatac  <  2:177029‑M1/149‑146 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
AAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCA  >  NC_000913_3_hsa_tpiA/1978357‑1978647

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: