New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 179793 =49 (1.320)7 (0.200) 7/224 3.0 13.7% coding (557/1518 nt) dgt deoxyguanosine triphosphate triphosphohydrolase
?NC_000913 179819 = 42 (1.200)coding (583/1518 nt) dgt deoxyguanosine triphosphate triphosphohydrolase
Rejected: Frequency below/above cutoff threshold.

GCCAACGTTTGCACCCGGAAGATGCCGAAAGCCAGCCTCTGACTGACGATCGCTGCAGCGTGGCGGCACTACGTTTACGGGACGGGGAAGAACCGCTTAACGAGCTGCGGCGCAAGATTCGTCAGGACTTATGTCATTTTGAGGGGAATGCACAAGGCATTCGCCTGGTGCATACATTGATGCGGATGAATCTCACCTGGGCACAGGTTGGCGGTATTTTAAAATATACCCGTCCGGCGTGGTGGCGTGGCGAAACGCCTGAGACACATCACTATTTAATGAAAAAGCCGGGTTA  >  NC_000913/179646‑179940
                                                                                                                                                   |                                                                                                                                                   
gCCAACGTTTGCACCCGGAAGATGCCGAAAGCCAGCCTCTGACTGACGATCGCTGCAGCGTGGCGGCACTACGTTTACGGGACGGGGAAGAACCGCTTAACGAGCTGCGGCGCAAGATTCGTCAGGACTTATGTCATTTTGAGGGGAAt                                                                                                                                                    <  1:109188/149‑1 (MQ=255)
gCCAACGTTTGCACCCGGAAGATGCCGAAAGCCAGCCTCTGACTGACGATCGCTGCAGCGTGGCGGCACTACGTTTACGGGACGGGGAAGAACCGCTTAACGAGCTGCGGCGCAAGATTCGTCAGGACTTATGTCATTTTGAGGGGAAt                                                                                                                                                    <  2:91354/149‑1 (MQ=255)
gCCAACGTTTGCACCCGGAAGATGCCGAAAGCCAGCCTCTGACTGACGATCGCTGCAGCGTGGCGGCACTACGTTTACGGGACGGGGAAGAACCGCTTAACGAGCTGCGGCGCAAGATTCGTCAGGACTTATGTCATTTTGAGGGGAAt                                                                                                                                                    <  2:88393/149‑1 (MQ=255)
  tataggcttcttcagaaagataataacccggctttttcattaaatagtgatgtgtctcaggcgtttcgccacgccaccacgccggacgggtatattttaaaataccgccaacctgtgcccaggtgagattcatccgcatcaatgtATGc                                                                                                                                                  >  6:107604‑M1/146‑149 (MQ=255)
       gcttcttcagaaagataataacccggctttttcattaaatagtgatgtgtctcaggcgtttcgccacgccaccacgccggacgggtatattttaaaataccgccaacctgtgcccaggtgagattcatccgcatcaatgtATGCACAAg                                                                                                                                             <  4:72258‑M1/9‑1 (MQ=255)
            aCCCGGAAGATGCCGAAAGCCAGCCTCTGACTGACGATCGCTGCAGCGTGGCGGCACTACGTTTACGGGACGGGGAAGAACCGCTTAACGAGCTGCGGCGCAAGATTCGTCAGGACTTATGTCATTTTGAGGGGAATGCACAAGGCAtt                                                                                                                                        >  7:21037/1‑149 (MQ=255)
              ccGGAAGATGCCGAAAGCCAGCCTCTGACTGACGATCGCTGCAGCGTGGCGGCACTACGTTTACGGGACGGGGAAGAACCGCTTAACGAGCTGCGGCGCAAGATTCGTCAGGACTTATGTCATTTTGAGGGGAATGCACAAGGc                                                                                                                                           <  6:131677/144‑1 (MQ=255)
              ccGGAAGATGCCGAAAGCCAGCCTCTGACTGACGATCGCTGCAGCGTGGCGGCACTACGTTTACGGGACGGGGAAGAACCGCTTAACGAGCTGCGGCGCAAGATTCGTCAGGACTTATGTCATTTTGAGGGGAATGCACAAGGc                                                                                                                                           >  5:131677/1‑144 (MQ=255)
                  aaGATGCCGAAAGCCAGCCTCTGACTGACGATCGCTGCAGCGTGGCGGCACTACGTTTACGGGACGGGGAAGAACCGCTTAACGAGCTGCGGCGCAAGATTCGTCAGGACTTATGTCATTTTGAGGGGAATGCACAAGGCATTCGCCTg                                                                                                                                  <  1:71268/149‑1 (MQ=255)
                          gAAAGCCAGCCTCTGACTGACGATCGCTGCAGCGTGGCGGCACTACGTTTACGGGACGGGGAAGAACCGCTTAACGAGCTGCGGCGCAAGATTCGTCAGGACTTATGTCATTTTGAGGGGAATGCACAAGGCATTCGCCTGGTGCATAc                                                                                                                          <  1:37082/149‑1 (MQ=255)
                            cccggctttttcattaaatagtgatgtgtctcaggcgtttcgccacgccaccacgccggacgggtatattttaaaataccgccaacctgtgcccaggtgagattcatccgcatcaatgtATGCACAAGGCATTCGCCTGGTGCATACAt                                                                                                                        >  3:76485‑M1/120‑149 (MQ=255)
                                 agccTCTGACTGACGATCGCTGCAGCGTGGCGGCACTACGTTTACGGGACGGGGAAGAACCGCTTAACGAGCTGCGGCGCAAGATTCGTCAGGACTTATGTCATTTTGAGGGGAATGc                                                                                                                                                  <  2:34674/118‑1 (MQ=255)
                                 agccTCTGACTGACGATCGCTGCAGCGTGGCGGCACTACGTTTACGGGACGGGGAAGAACCGCTTAACGAGCTGCGGCGCAAGATTCGTCAGGACTTATGTCATTTTGAGGGGAATGc                                                                                                                                                  >  1:34674/1‑118 (MQ=255)
                                     tctGACTGACGATCGCTGCAGCGTGGCGGCACTACGTTTACGGGACGGGGAAGAACCGCTTAACGAGCTGCGGAGCAAGATTCGTCAGGACTTATGTCATTTTGAGTGGAATGCACAAGGCATTCGCCTGGTGGATACATTGATGCGGa                                                                                                               >  8:38485/1‑149 (MQ=255)
                                         ttaaatagtgatgtgtctcaggcgtttcgccacgccaccacgccggacgggtatattttaaaataccgccaacctgtgcccaggtgagattcatccgcatcaatgtATGCACAAGGCATTCGCCTGGt                                                                                                                                >  1:159640‑M1/107‑128 (MQ=255)
                                         ttaaatagtgatgtgtctcaggcgtttcgccacgccaccacgccggacgggtatattttaaaataccgccaacctgtgcccaggtgagattcatccgcatcaatgtATGCACAAGGCATTCGCCTGGt                                                                                                                                <  2:159640‑M1/22‑1 (MQ=255)
                                              tagtgatgtgtctcaggcgtttcgccacgccaccacgccggacgggtatattttaaaataccgccaacctgtgcccaggtgagattcatccgcatcaatgtATGc                                                                                                                                                  <  1:79554‑M1/4‑1 (MQ=255)
                                              tagtgatgtgtctcaggcgtttcgccacgccaccacgccggacgggtatattttaaaataccgccaacctgtgcccaggtgagattcatccgcatcaatgtATGc                                                                                                                                                  >  2:79554‑M1/102‑105 (MQ=255)
                                                   gCTGCAGCGTGGCGGCACTACGTTTACGGGACGGGGAAGAACCGCTTAACGAGCTGCGGCGCAAGATTCGTCAGGACTTATGTCATTTTGAGGGGAATGCACAAGGCATTCGCCTGGTGCATACATTGATGCg                                                                                                                 >  4:143504/1‑133 (MQ=255)
                                                   gCTGCAGCGTGGCGGCACTACGTTTACGGGACGGGGAAGAACCGCTTAACGAGCTGCGGCGCAAGATTCGTCAGGACTTATGTCATTTTGAGGGGAATGCACAAGGCATTCGCCTGGTGCATACATTGATGCg                                                                                                                 >  4:141070/1‑133 (MQ=255)
                                                   gCTGCAGCGTGGCGGCACTACGTTTACGGGACGGGGAAGAACCGCTTAACGAGCTGCGGCGCAAGATTCGTCAGGACTTATGTCATTTTGAGGGGAATGCACAAGGCATTCGCCTGGTGCATACATTGATGCg                                                                                                                 <  3:143504/133‑1 (MQ=255)
                                                   gCTGCAGCGTGGCGGCACTACGTTTACGGGACGGGGAAGAACCGCTTAACGAGCTGCGGCGCAAGATTCGTCAGGACTTATGTCATTTTGAGGGGAATGCACAAGGCATTCGCCTGGTGCATACATTGATGCg                                                                                                                 <  3:141070/133‑1 (MQ=255)
                                                     tGCAGCGTGGCGGCACTACGTTTACGGGACGGGGAAGAACCGCTTAACGAGCTGCGGCGCAAGATTCGTCAGGACTTATGTCATTTTGAGGGGAATGCACAAGGCATTCGCCTGGTGCATACATTGATGCGGATGAAtctc                                                                                                       >  8:68186/1‑141 (MQ=255)
                                                     tGCAGCGTGGCGGCACTACGTTTACGGGACGGGGAAGAACCGCTTAACGAGCTGCGGCGCAAGATTCGTCAGGACTTATGTCATTTTGAGGGGAATGCACAAGGCATTCGCCTGGTGCATACATTGATGCGGATGAAtctc                                                                                                       <  7:68186/141‑1 (MQ=255)
                                                                    cTACGGTTACGGGACGGGGAAGAACCGCTTAACGAGCTGCGGCGCAAGATTCGTCAGGACTTATGTCATTTTGAGGGGAATGCACAAGGCATTCGCCTGGt                                                                                                                                <  3:46171/101‑1 (MQ=255)
                                                                    cTACGGTTACGGGACGGGGAAGAACCGCTTAACGAGCTGCGGCGCAAGATTCGTCAGGACTTATGTCATTTTGAGGGGAATGCACAAGGCATTCGCCTGGt                                                                                                                                >  4:46171/1‑101 (MQ=255)
                                                                     tACGTTTACGGGACGGGGAAGAACCGCTTAACGAGCTGCGGCGCAAGATTCGTCAGGACTTATGTCATTTTGAGGGGAATGCACAAGGCATTCGCCTGGTGCATACAt                                                                                                                        >  1:212701/1‑108 (MQ=255)
                                                                     tACGTTTACGGGACGGGGAAGAACCGCTTAACGAGCTGCGGCGCAAGATTCGTCAGGACTTATGTCATTTTGAGGGGAATGCACAAGGCATTCGCCTGGTGCATACAt                                                                                                                        <  2:212701/108‑1 (MQ=255)
                                                                            acgggacgtgGAAGAACCGCTTAACGAGCTGCGGCGCAAGATTCGTCAGGACTTATGTCATTTTGAGGGGAATGCACAAGGCATTCGCCTGGTGCATACATTGATGCGGATGAATCTCACCTGGGCACAGGTTGGCGGTATTTTaaaa                                                                         <  8:21037/148‑1 (MQ=255)
                                                                               ggacgggGAAGAACCGCTTAACGAGCTGCGGCGCAAGATTCGTCAGGACTTATGTCATTTTGAGGGGAATGCACAAGGCATTCGCCTGGTGCATGCATTGATGCGGATGAAtctc                                                                                                       >  7:161727/1‑115 (MQ=255)
                                                                               ggacgggGAAGAACCGCTTAACGAGCTGCGGCGCAAGATTCGTCAGGACTTATGTCATTTTGAGGGGAATGCACAAGGCATTCGCCTGGTGCATGCATTGATGCGGATGAAtctc                                                                                                       <  8:161727/115‑1 (MQ=255)
                                                                                     ggAAGAACCGCTTAACGAGCTGCGGCGCAAGATTCGTCAGGACTTATGTCATTTTGAGGGGAATGCACAAGGCATTCGCCTGGTGCATACATTGATGCGGATGAATCTCACCTGGGCACAGGTTGGCGGTATTTTAAAATATACCCGTc                                                               >  3:110535/1‑149 (MQ=255)
                                                                                       aagaaCCGCTTAACGAGCTGCGGCGCAAGATTCGTCAGGACTTATGTCATTTTGAGGGGAATGCACAAGGCATTCGCCTGGTGCATACATTGATGCGGATGAATCTCACCTGGGCACAGGTTGGCGGTATTTTaaaa                                                                         >  1:108021/1‑137 (MQ=255)
                                                                                       aagaaCCGCTTAACGAGCTGCGGCGCAAGATTCGTCAGGACTTATGTCATTTTGAGGGGAATGCACAAGGCATTCGCCTGGTGCATACATTGATGCGGATGAATCTCACCTGGGCACAGGTTGGCGGTATTTTaaaa                                                                         <  2:108021/137‑1 (MQ=255)
                                                                                              gCTTAACGAGCTGCGGCGCAAGATTCGTCAGGACTTATGTCATTTTGAGGGGAATGCACAAGGCATTCGCCTGGTGCATACATTGATGCGGATGAATCTCACCTGGGCACAGGTTGGCGGTATTTTaaaa                                                                         >  3:80192/1‑130 (MQ=255)
                                                                                              gCTTAACGAGCTGCGGCGCAAGATTCGTCAGGACTTATGTCATTTTGAGGGGAATGCACAAGGCATTCGCCTGGTGCATACATTGATGCGGATGAATCTCACCTGGGCACAGGTTGGCGGTATTTTaaaa                                                                         <  4:80192/130‑1 (MQ=255)
                                                                                                       ataccgccaacctgtgcccaggggagattcatccgcatcaatgtATGCACAAGGCATTCGCCTGGTGCATACATTGATGc                                                                                                                  >  3:125073‑M1/45‑80 (MQ=255)
                                                                                                         tGCGGCGCAAGATTCGTCAGGACTTATGTCATTTTGAGGGGAATGCACAAGGCATTCGCCTGGTGCATACATTGATGCGGATGAATCTCACCTGGGCACAGGTTGGCGGTATTTTAAAATATACCCGTCCGGCGTGGTGGCGTGGCGa                                            >  1:14400/1‑148 (MQ=255)
                                                                                                               gcAAGATTCGTCAGGACTTATGTCATTTTGAGGGGAATGCACAAGGCATTCGCCTGGt                                                                                                                                >  1:117790/1‑58 (MQ=255)
                                                                                                               gcAAGATTCGTCAGGACTTATGTCATTTTGAGGGGAATGCACAAGGCATTCGCCTGGTGCATACATTGATGCGGATGAATCTCACCTGGGCACAGGTTGGCGGTATTTTAAAATATACCCGTCCGGCGTGGTGGCGTGGCGAAACGCCt                                     <  7:38485/149‑1 (MQ=255)
                                                                                                                acctgtgcccaggtgagattcatccgcatcaatgtATGc                                                                                                                                                  >  1:165446‑M1/36‑39 (MQ=255)
                                                                                                                acctgtgcccaggtgagattcatccgcatcaatgtATGc                                                                                                                                                  <  2:165446‑M1/4‑1 (MQ=255)
                                                                                                                       cGTCAGGACTTATGTCATTTTGAGGGGAATGCACAAGGCATTCGCCTGGTGCATACATTGATGCGGATGAATCTCACCTGGGCACAGGTTGGCGGTATTTTAAAATATACCCGTCCGGCGTGGTGGCGTGGCGAAACGCCTGAGacac                              <  4:110535/148‑1 (MQ=255)
                                                                                                                         tCAGGACTTATGTCATTTTGAGGGGAATGCACAAGGCATTCGCCTGGTGCATACATTGATGCGGATGAATCTCACCTGGGCACAGGTTGGCGGTATTTTAAAATATACCCGTCCGGCGTGGTGGCGTGGCGAAACGCCTGAGACACAt                            >  4:61684/1‑148 (MQ=255)
                                                                                                                          aggtgagattcatccgcatcaatgtATGCACCAGGCGAATGCCTTGTGCATACATTGATGCGGATGAATCTCACCTGGGCACAGGTTGGCGGTATTTTAAAATATACCCGTCCGGCGTGGTGGCGTGGCGAAACGCCTGAGACACATCa                          <  3:72258‑M1/124‑1 (MQ=255)
                                                                                                                             gACTTATGTCATTTTGAGGGGAATGCACAAGGCATTCGCCTGGTGCATACATTGATGCGGATGAATCTCACCTgg                                                                                                 >  8:140200/1‑75 (MQ=255)
                                                                                                                             gACTTATGTCATTTTGAGGGGAATGCACAAGGCATTCGCCTGGTGCATACATTGATGCGGATGAATCTCACCTgg                                                                                                 <  7:140200/75‑1 (MQ=255)
                                                                                                                                ttATGTCATTTTGAGGGGAATGCACAAGGCATTCGCCTGGTGCATACATTGATGCGGATGAATCTCACCTGGGCACAGGTTGGCGGTATTTTAAAATATACCCGTCCGGCGTGGTGGCGTGGCGAAACGCCTGAGACACATCACTAtt                     >  1:125064/1‑148 (MQ=255)
                                                                                                                                     tCATTTTGAGGGGAATGCACAAGGCATTCGCCTGGTGCATACATTGATGCGGATGAATCTCACCTGGGCACAGGTTGGCGGTATTTTAAAATATACCCGTCCGGCGTGGTGGCGTGGCGAAACGCCTga                                   <  1:50791/129‑1 (MQ=255)
                                                                                                                                     tCATTTTGAGGGGAATGCACAAGGCATTCGCCTGGTGCATACATTGATGCGGATGAATCTCACCTGGGCACAGGTTGGCGGTATTTTAAAATATACCCGTCCGGCGTGGTGGCGTGGCGAAACGCCTga                                   >  2:50791/1‑129 (MQ=255)
                                                                                                                                       aTTTTGAGGGGAATGCACAAGGCATTCGCCTGGTGCATACATTGATGCGGATGAATCTCACCTGGGCACAGGTTGGCGGTATTTTAAAATATACCCGTCCGGCGTGGTGGCGTGGCGAAACGCCTGAGACACATCACTATTTAATGaaa             <  2:125064/149‑1 (MQ=255)
                                                                                                                                        ttttGAGGGGAATGCACAAGGCATTCGCCTGGTGCATACATTGATGCGGATGAATCTCACCTGGGCACAGGTTGGCGGTATTTTaaa                                                                          <  5:195366/87‑1 (MQ=255)
                                                                                                                                        ttttGAGGGGAATGCACAAGGCATTCGCCTGGTGCATACATTGATGCGGATGAATCTCAACTGGGCACAGGTTGGCGGTATTTTaaa                                                                          >  6:195366/1‑87 (MQ=255)
                                                                                                                                         gcatcaatgtATGCACCAGGCGAATGCCTTGTGCATACATTGATGCGGATGAATCTCACCTGGGCACAGGTTGGCGGTAt                                                                                >  4:125073‑M1/11‑80 (MQ=255)
                                                                                                                                                 gAATGCACAAGGCATTCGCCTGGTGCATACATTGATGCGGATGAATCTCACCTGGGCACAGGTTGGCGGTATTTTAAAATATACCCGTCCGGCGTGGTGGCGTGGCGAAACGCCTGAGACACATCa                          >  6:197703/1‑126 (MQ=255)
                                                                                                                                                 gAATGCACAAGGCATTCGCCTGGTGCATACATTGATGCGGATGAATCTCACCTGGGCACAGGTTGGCGGTATTTTAAAATATACCCGTCCGGCGTGGTGGCGTGGCGAAACGCCTGAGACACATCa                          <  5:197703/126‑1 (MQ=255)
                                                                                                                                                 gAATGCACAAGGCATTCGCCTGGTGCATACATTGATGCGGATGAATCTCACCTGGGCACAGGTTGGCGGTATTTTAAAATATACCCGTCCGGCGTGGTGGCGTGGCGAAACGCCTGAGACACATCACTATTTAATGAAAAAGCCgg      <  6:7425/146‑1 (MQ=255)
                                                                                                                                                 gAATGCACAAGGCATTCGCCTGGTGCATACATTGATGCGGATGAATCTCACCTGGGCACAGGTTGGCGGTATTTTAAAATATACCCGTCCGGCGTGGTGGCGTGGCGAAACGCCTGAGACACATCACTATTTAATGAAAAAGCCgg      >  5:7425/1‑146 (MQ=255)
                                                                                                                                                  aaTGCACAAGGCATTCGCCTGGTGCATACATTGATGCGGATGAATCTCACCTGGGCACAGGTTGGCGGTATTTTAAAATATACCCGTCCGGCGTGGTGGCGTGGCGAAACGCCTGAGACACATCACTATTTAATGa               <  7:11419/136‑1 (MQ=255)
                                                                                                                                                  aaTGCACAAGGCATTCGCCTGGTGCATACATTGATGCGGATGAATCTCACCTGGGCACAGGTTGGCGGTATTTTAAAATATACCCGTCCGGCGTGGTGGCGTGGCGAAACGCCTGAGACACATCACTATTTAATGa               >  8:11419/1‑136 (MQ=255)
                                                                                                                                                  aaTGCACAAGGCATTCGCCTGGTGCATACATTGATGCGGATGAATCTCACCTGGGCACAGGTTGGCGGTATTTTAAAATATACCCGTCCGGCGTGGTGGCGTGGCGAAACGCCTGAGACACATCACTATTTAATGAAAAAGCCGGGtta  <  2:14400/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GCCAACGTTTGCACCCGGAAGATGCCGAAAGCCAGCCTCTGACTGACGATCGCTGCAGCGTGGCGGCACTACGTTTACGGGACGGGGAAGAACCGCTTAACGAGCTGCGGCGCAAGATTCGTCAGGACTTATGTCATTTTGAGGGGAATGCACAAGGCATTCGCCTGGTGCATACATTGATGCGGATGAATCTCACCTGGGCACAGGTTGGCGGTATTTTAAAATATACCCGTCCGGCGTGGTGGCGTGGCGAAACGCCTGAGACACATCACTATTTAATGAAAAAGCCGGGTTA  >  NC_000913/179646‑179940

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ 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.