New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 4546905 =217 (1.370)6 (0.040) 6/252 NT 2.6% coding (1810/2637 nt) fimD fimbrial usher outer membrane porin protein; FimCD chaperone‑usher
?NC_000913 4546951 = 237 (1.580)coding (1856/2637 nt) fimD fimbrial usher outer membrane porin protein; FimCD chaperone‑usher
Rejected: Frequency below/above cutoff threshold.

Only 100 of 278 total aligned reads displayed.

TACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGATATCAACTGGACGCTCAGCTATAGCCTGACGAAAAACGCCTGGCAAAAAGGACGGGATCAGATGTTAGCGCTTAACGTCAATATTCCTTTCAGCCACTGGCTGCGTTCTGACAGTAAATCTCAGTGGCGACATGCCAGTGCCAGCTACAGCATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACGGTACGTTGCTGGAAGACAACAACCTCAGCTATAGCGTGCAAACCGGCTATGCCGGGGGAGGCGATGGAAATAGCGGAAGTACAGGCTACGCCACGCTGAATTATCGCGGTGGTTACGGCAATGCCAATATCGGTTACAGCCATAGCGATGATATTAAGCAGCTCTATTACGGAGTCAGCGGTGGGGTACTGGCTCATGCCAATGGC  >  NC_000913/4546681‑4547216
                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                         
                    cAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGATATCAACTGGACGCTCAGCTATAGCCTGACGAAAAACGCCTGGCAAAAAGGACGGGATCAGATGTTAGCGCTTAACGTCAATATTCCTTTCAGCCACTGGCTGCGTTCTGACAGTAAATCTCAGTGGCGACCTGCCAGTGCCAGCTACAGCATGTCACACGAtc                                                                                                                                                                                                                                                    >  2:384608/1‑274 (MQ=255)
                                  acaaaCAATGGATTTGAGTGGTAGCCATCAAAATTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTGCGAAGACATCAACTGGACGCGCAGCTAGAGACTGACGAAAAACGCCTGGCAAAAAGGACGGGATCAGATGTTAGCGCTTAACGTCAATATTCCGTTCAGCCACTGGCTGCGTTCTGACAGTAAATCTCAGTGGCGACATGCCAGTGCCAGCTACAGCATGTCACACGATCTCAACGGTCGGATGa                                                                                                                                                                                                                                     <  4:401057/272‑1 (MQ=255)
                                        acttgaTTTGAGTGGTAGCCATCAAACGTATTGGGGAACGAGTAATGTGGAGGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGATATCAACTGGACGCTCAGCTATAGCCTGACGAAAAACGCCTGGCAAAAAGGACGGGATCAGATGTTAGCGCTTAACGTCAATATTCCTTTCAGCCACTGGCTGCGTTCTGACAGTAAATCTCAGTGGCGACATGCCAGTGCCAGCTACAGCATGTCACACGATCTCAACGGTCGGATGACCAAt                                                                                                                                                                                                                                <  4:131922/269‑1 (MQ=255)
                                          tGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGATATCAACTGGACGCTCAGCTATAGCCTGACGAAAAACGCCTGGCAAAAAGGACGGGATCAGATGTTAGCGCTTAACGTCAATATTCCTTTCAGCCACTGGCTGCGTTCTGACAGTAAATCTCAGTGGCGACATGCCAGTGCCAGCTACAGCATGTCACACGATCCCAACGGGCGGATGCCCAATctg                                                                                                                                                                                                                             >  4:245163/1‑275 (MQ=255)
                                            tATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGCTGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGATATCAACTGGACGCTCAGCTATAGCCTGACGAAAAACGCCTGGCAAAAAGGACGGGATCAGATGTTAGCGCTTAACGTCAATATTCCTTTCAGCCACTGGCTGCGTTCTGACAGTAAATCCACGTGCCGACATTCCAGTGCCAGCTACAGCCTCTCACACTCTCTCCACGGTCGGATGCCCAAttcgg                                                                                                                                                                                                                            >  4:384691/1‑270 (MQ=255)
                                            tATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGATATCAACTGGACGCTCAGCTATAGCCTGACGAAAAACGCCTGGCAAAAAGGACGGGATCAGATGTTAGCGCTTAACGTCAATATTCCTTTCAGCCACTGGCTGCGTTCTGACAGTAAATCTCAGTGGCGACATGCCAGTGCCAGCTACAGCATGTCACACGATCTCAACGGTCGGATGACCAATctgg                                                                                                                                                                                                                            <  2:123339/274‑1 (MQ=255)
                                               ttGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGATATCAACTGGACGCTCAGCTATAGCCTGACGAAAAACGCCTGGCAAAAAGGACGGGATCAGATGTTAGCGCTTAACGTCAATATTCCTTTCAGCCACTGGCTGCGTTCTGACAGNAAATCTCAGTGGCGACATGCCAGTGCCAGCGAAAACATGTCAAACGAACCAAAAGGTCGGATAGACAAActggccgg                                                                                                                                                                                                                        >  3:133591/1‑275 (MQ=255)
                                                      gTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGATATCAACTGGACGCTCAGCTATAGCCTGACGAAAAACGCCTGGCAAAAAGGACGGGATCAGATGTTAGCGCTTAACGTCAATATTCCTTTCAGCCACTGGCTGCGTTCTGACAGTAAATCTCAGTGGCGACATGCCAGTGCCAGCTACAGCATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATAc                                                                                                                                                                                                                 >  2:220179/1‑275 (MQ=255)
                                                        aGCCATCAAACTTATTGGGGAACGTGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGATATCAACTGGACGCTCAGCTATAGCCTGACGAAAAACGCCTGGCAAAAAGGACGGGATCAGATGTTAGCGCTTAACGTCAATATTCCTTTCAGCCACTGGCTGCGTTCTGACAGTAAATCTCAGTGGCGACATGCCAGTGCCAGCTACAGCATGTCACCCGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACgg                                                                                                                                                                                                               >  4:188925/1‑275 (MQ=255)
                                                        aGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGATATCAACTGGACGCTCAGCTATAGCCTGACGAAAAACGCCTGGCAAAAAGGACGGGATCAGATGTTAGCGCTTAACGTCAATATTCCTTTCAGCCACTGGCTGCGTTCTGACAGTAAATCTCAGTGGCGACATGCAAGTGCCAGCTACAGCATGTCACACGATCTCAACGGTCGGAGTACCAATc                                                                                                                                                                                                                               >  4:492144/1‑259 (MQ=255)
                                                                                aGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGATATCAACTGGACGCTCAGCTATAGCCTGACGAAAAACGCCTGGCAAAAAGGACGGGATCAGATGTTAGCGCTTAACGTCAATATTCCTTTCAGCCACTGGCTGCGTTCTGACAGTAAATCTCAGTGGCGACATGCCAGTGCCAGCTACAGCATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACGGTACGTTGCTGGAAGACAACAACCt                                                                                                                                                                                       >  1:94631/1‑275 (MQ=255)
                                                                                aGTAAAGTCGATGAGCAATTCCAGGCTGGAGTAAATCCTGCGTTCGAAGATATCAACTGGACGCTCAGCTATAGCCTGACGAAAAACGCCTGGCAAAAAGGACGTGATCAGATGTTAGCGCTTAACGTCAATATTCCTTTCAGCCACTGGCTGCGTTCTGACAGTAAATCTCAGTGGCGACATGCCAGTGCCAGCTACAGCATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACGGTACGTTGCTGGAAGACAACAACCt                                                                                                                                                                                       <  4:50407/275‑1 (MQ=255)
                                                                                                                              aaaataacgaCTGGACGCTCAGCTATAGCCTGACGAAAAACGCCTGTCAAAAAGGACGGGATCAGATGTTAGCGCTTAACGTCAATATTCCTTTCAGCCACTGGCTGCGTTCTGACAGTAAATCTCAGTGGCGACATGCCAGTGCCAGCTACAGCATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTTGTGTATACGGTACGTTGCTGGAAGACAACAACCTCAGCTATAGCGTGCAAACCGGCTATGCCGGGGGAGGCGATGGAAAt                                                                                                                                         <  1:54582/266‑1 (MQ=255)
                                                                                                                                   atCAACTGGACGCTCAGCTATAGCCTGACGAAAAACGCCTGGCAAAAAGGACGGGATCAGATGTTAGCGCTTAACGTCAATATTCCTTTCAGCCACTGGCTGCGTTCTGACAGTAAATCTCAGTGGCGACATGCCAGTGCCAGCTACAGCATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACGGTACGTTGCTGGAAGACAACAACCTCCGCCATAGCGTGCAAACCGGCTATGCCGGGGGAGGCGGTGGAAATAGCg                                                                                                                                     >  4:73544/1‑274 (MQ=255)
                                                                                                                                   atCAACTGGACGCTCAGCTATAGCCTGACGAAAAACGCCTGGCAAAAAGGACGGGATCAGATGTTAGCGCTTAACGTCAATATTCCTTTCAGCCACTGGCTGCGTTCTGACAGTAAATCTCAGTGGCGACATGCCAGTGCCAGCTACAGCATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACGGTACGTTGCTGGAAGACAACAACCTCAGCTATAGCGTGCAAACCGGCTATGCCGGGGGAGGCGATGGAAATAGCgg                                                                                                                                    <  4:454331/275‑1 (MQ=255)
                                                                                                                                                  aGCTATAGCCTGACGAAAAACGCCTGGCAAAAAGGACGGGATCAGATGTTAGCGCTTAACGTCAATATTCCTTTCAGCCACTGGCTGCGTTCTGACAGTAAATCTCAGTGGCGACATGCCAGTGCCAGCTACAGCAt                                                                                                                                                                                                                                                               <  5:291357/137‑1 (MQ=255)
                                                                                                                                                  aGCTATAGCCTGACGAAAAACGCCTGGCAAAAAGGACGGGATCAGATGTTAGCGCTTAACGTCAATATTCCTTTCAGCCACTGGCTGCGTTCTGACAGTAAATCTCAGTGGCGACATGCCAGTGCCAGCTACAGCAt                                                                                                                                                                                                                                                               >  6:291357/1‑137 (MQ=255)
                                                                                                                                                    cTATAGCCTGACGAAAAACGCCTGGCAAAAAGGACGGGATCAGATGTTAGCGCTTAACGTCAATATTCCTTTCAGCCACTGGCTGCGTTCTGACAGTAAATCTCAGTGGCGACATGCCAGTGCCAGCTACAGCATGTCACACGAtctc                                                                                                                                                                                                                                                  >  8:659318/1‑148 (MQ=255)
                                                                                                                                                              aCGAAAAACGCCTGGCAAAAAGGACGGGATCAGATGTTAGCGCTTAACGTCAATATTCCTTTCAGCCACTGGCTGCGTTCTGACAGTAAATCTCAGTGGCGACATGCCAGTGCCAGCTACAGCATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACGGTACGTTGCTGGAAGACAACAACCTCAGCTATAGCGTGCAAACCGGCTATGCCGGGGGCGGCGATGGAAATAGCGGAAGTACAGGCTACGCCACGCTGAATTa                                                                                                         >  2:278374/1‑275 (MQ=255)
                                                                                                                                                              aCGAAAAACGCCTGGCAAAAAGGACGGGATCAGATGTTAGCGCTTAACGTCAATATTCCTTTCAGCCACTGGCTGCGTTCTGACAGTAAATCTCAGTGGCGACATGCCAGTGCCAGCTACAGCATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACGGTACGTTGCTGGAAGACAACAACCTCAGCTATAGCGTGCAAACCGGCTATGCCGGGGGAGGCGATGGAAATAGCGGAAGTACAGGCTACGCCACGCTTAATTa                                                                                                         >  2:445392/1‑275 (MQ=255)
                                                                                                                                                                 aaaaaCGCCTGGCAAAAAGGACGGGATCAGATGTTAGCGCTTAACGTCAATATTCCTTTCAGCCACTGGCTGCGTTCTGACAGTAAATCTCAGTGGCGACATGCCAGTGCCAGCTACAGCATGTCACACGATCTCaa                                                                                                                                                                                                                                                >  5:353870/1‑137 (MQ=255)
                                                                                                                                                                  aaaaCGCCTGGCAAAAAGGACGGGATCAGATGTTAGCGCTTAACGTCAATATTCCTTTCAGCCACTGGCTGCGTTCTGACAGTAAATCTCAGTGGCGACAtgccagtgcca                                                                                                                                                                                                                                                                         >  8:919599/1‑111 (MQ=255)
                                                                                                                                                                   aaaCGCCTGGCAAAAAGGACGGGATCAGATGTTAGCGCTTAACGTCAATATTCCTTTCAGCCACTGGCTGCGTTCTGACAGTAAATCTCAGTGGCGACATGCCAGTGCCAGCTACAGCATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACGGTACGTTGCTGGAAGACAACAACCTCAGCTATAGCGTGCAAACCGGCTATGCCGGGGGAGGCGATGGAAATAGCGGAAGTACAGGCTACGCCACGCTGAATTATcgcg                                                                                                    <  2:157277/275‑1 (MQ=255)
                                                                                                                                                                           ggCAAAAAGGACGGGATCAGATGTTAGCGCTTAACGTCAATATTCCTTTCAGCCACTGGCTGCGTTCTGACAGTAAATCTCAGTGGCGACATGCCAGTGCCAGCTACAGCATGTCACACGATCTCAACGGTCGGAt                                                                                                                                                                                                                                       <  7:1018384/136‑1 (MQ=255)
                                                                                                                                                                                aaaGGACGGGATCAGATGTTAGCGCTTAACGTCAATATTCCTTTCAGCCACTGGCTGCGTTCTGACAGTAAATCTCAGTGGCGACATGCCAGTGCCAGCTACGGCATGTCACACGATCTCAACGGTCGGATAACCAATCTGGCTGGTGTATACGGTACGTTGCTGGAAGACAACAACCTCAGCTATAGCGTGCAAACCGGCTATGCCGGGGGAGGCGATGGAAATAGCGGAAGTACAGGCTACGCCACGCTGAATTATCGCGGTGGTTACGGc                                                                                         <  2:346188/273‑1 (MQ=255)
                                                                                                                                                                                aaaGGACGGGATCAGATGTTAGCGCTTAACGTCAATATTCCTTTCAGCCACTGGCTGCGTTCTGACAGTAAATCTCAGTGGCGACATGCCAGTGCCAGCTACAGCATGTCACACGATCTCa                                                                                                                                                                                                                                                 <  8:371194/121‑1 (MQ=255)
                                                                                                                                                                                aaaGGACGGGATCAGATGTTAGCGCTTAACGTCAATATTCCTTTCAGCCACTGGCTGCGTTCTGACAGTAAATCTCAGTGGCGACATGCCAGTGCCAGCTACAGCATGTCACACGATCTCa                                                                                                                                                                                                                                                 >  7:371194/1‑121 (MQ=255)
                                                                                                                                                                                   ggACGGGATCAGATGTTAGCGCTTAACGTCAATATTCCTTTCAGCCACTGGCTGCGTTCTGACAGTAAATCTCAGTGGCGACATGCCAGTGCCAGCTACAGCATGTCACACGATCTCAACGGTCGGATGACCAAt                                                                                                                                                                                                                                >  6:614098/1‑135 (MQ=255)
                                                                                                                                                                                         gATCAGATGTTAGCGCTTAACGTCAATATTCCTTTCAGCCACTGGCTGCGTTCTGACAGTAAATCTCAGTGGCGACATGCCAGTGCCAGCTACAGCATGTCACACGATCTCAACGGt                                                                                                                                                                                                                                            >  6:728475/1‑117 (MQ=255)
                                                                                                                                                                                          aTCAGAGGTTAGCGCTTAACGTCAATATTCCTTTCAGCCACTGGCTGCGTTCTGACAGTAAATCTCAGTGGCGACATGCCAGTGCCAGCTACAGCATGTCACACGATCTCAACGGTCGGATGACCAATCTGTCTGGTGTATACGGTAc                                                                                                                                                                                                            >  7:1112632/1‑148 (MQ=255)
                                                                                                                                                                                            cAGATGTTAGCGCTTAACGTCAATATTCCTTTCAGCCACTGGCTGCGTTCTGACAGTAAATCTCAGTGGCGACATGCCAGTGCCAGCTACAGCATGTCACACGATCTCaa                                                                                                                                                                                                                                                <  6:112001/110‑1 (MQ=255)
                                                                                                                                                                                            cAGATGTTAGCGCTTAACGTCAATATTCCTTTCAGCCACTGGCTGCGTTCTGACAGTAAATCTCAGTGGCGACATGCCAGTGCCAGCTACAGCATGTCACACGATCTCaa                                                                                                                                                                                                                                                >  5:112001/1‑110 (MQ=255)
                                                                                                                                                                                               aTGTTAGCGCTTAACGTCAATATTCCTTTCAGCCACTGGCTGCGTTCTGACAGTAAATCTCAGTGGCGACATGCCAGTGCCAGCTACAGCATGTCACACGATCTCAACGGTCGGATGa                                                                                                                                                                                                                                     <  5:560861/118‑1 (MQ=255)
                                                                                                                                                                                               aTGTTAGCGCTTAACGTCAATATTCCTTTCAGCAACTGGTGGCGTTCTGACAGTAAATCTCAGTGGCGACATGCCAGTGCCAGCTACAGCATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACGGTACGTTGCTGGAAGACAACAACCTCAGCTATAGCGTGCAAACCGGCTATGCCGGGGGAGGCGATGGAAATAGCGGAAGTACAGGCTACGCCACGCTGAATTATCGCGGTGGTTACGGCAATGCCAATATCGGTTa                                                                        <  2:501736/275‑1 (MQ=255)
                                                                                                                                                                                                tGTTAGCGCTTAACGTCAATATTCCTTTCAGCCACTGGCTGCGTTCTGACAGTAAATCTCAGTGGCGACATGCCAGTGCCAGCTACAGCATGTCACACGATCTCAACGGTCGGATGACCAATctggct                                                                                                                                                                                                                          <  8:580208/128‑1 (MQ=255)
                                                                                                                                                                                                   tAGCGCTTAACGTCAATATTCCTTTCAGCCACTGGCTGCGTTCTGACAGTAAATCTCAGTTGCGACATGCCAGTGCCAGCTACAGCATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGGNTATACGGTACGTTGCTGTAAGACAACAACCTCAGCTATAGCGGGCCAACGGGCTATGCCGGGGGAGACGNTGGAAATAGCGGAAGTCAAGGCTACGCCCCGATGAAGTATAGCGGNGNGTTCGGCAAGTCCAataaccgctacctc                                                                    >  3:52182/1‑264 (MQ=255)
                                                                                                                                                                                                   tAGCGCTTAACGTCAATATTCCTTTCAGCCACTGGCTGCGTTCTGACAGTAAATCTCAGTGGCGACATGCCAGTGCCAGCTACAGCATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACGGTACGTTGCTGGa                                                                                                                                                                                                   <  8:1112632/148‑1 (MQ=255)
                                                                                                                                                                                                        cTTAACGTCAATATTCCTTTCAGCCACTGGCTGCGTTCTGACAGTAAATCTCAGTGGCGACATGCCAGTGCCAGCTACAGCATGTCACACGATCTCAACGGTCGGATGa                                                                                                                                                                                                                                     >  6:7601/1‑109 (MQ=255)
                                                                                                                                                                                                        cTTAACGTCAATATTCCTTTCAGCCACTGGCTGCGTTCTGACAGTAAATCTCAGTGGCGACATGCCAGTGCCAGCTACAGCATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGgt                                                                                                                                                                                                                       <  3:95923/123‑1 (MQ=255)
                                                                                                                                                                                                          tAACGTCAATATTCCTTTCAGCCACTGGCTGCGTTCTGACAGTAAATCTCAGTGGCGACATGCCAGTGCCAGCTACAGCATGTCACACGATCTCAACGGTCGGATGACCAAt                                                                                                                                                                                                                                >  6:813203/1‑112 (MQ=255)
                                                                                                                                                                                                           aaCGTCAATATTCCTTTCAGCCACTGGCTGCGTTCTGACAGTAAATCTCAGTGGCGACATGCCAGTGCCAGCTACAGCATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACGGTACGTTGCTGGAAGacaaca                                                                                                                                                                                           >  1:115534/1‑148 (MQ=255)
                                                                                                                                                                                                           aaCGTCAATATTCCTTTCAGCCACTGGCTGCGTTCTGACAGTAAATCTCAGTGGCGACATGCCAGTGCCAGCTACAGCATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACGGTACGTTGCTGGAAGACAACAACCTCAGCTATAGCGTGCAAACCGGCTATGCCGGGGGAGGCGATGGAAATAGCGGAAGTACAGGCTACGCCACGCTGAATTATCGCGGTGGTTACGGCAATGCCAATATCGGTTACAGCCATAGCga                                                            <  2:188926/275‑1 (MQ=255)
                                                                                                                                                                                                           aaCGTCAATAGTCCTTTCAGCCACTGCGTGCGTTCGGACAGTAAATCTCAGTGGCGACAGGCCAGTGCCAGCTACAGCATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACGGTACGTTGCTGGAAGACAACAACATAAGCTATAGCGTGCAAACCGGCTATGCCGGGGGAGGCGATGGAAATAGCGGAAGTACAGGCTACGCCACGCTGAATTATCGCGGTGGTTACGGCAATGCCAATATCGGTTACAGCCATAGCga                                                            <  1:124294/275‑1 (MQ=255)
                                                                                                                                                                                                             cGTCAATATTCCTTTCAGCCACTGGCTGCGTTCTGACAGTAAATCTCAGTGGCGACATGCCAGTGCCAGCTACAGCATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACGGTACGTTGCTGGAAGaca                                                                                                                                                                                              <  5:647584/143‑1 (MQ=255)
                                                                                                                                                                                                               tCAATATTCCTTTCAGCCACTGGCTGCGTTCTGACAGTAAATCTCAGTGGCGACATGCCAGTGCCAGCTACAGCATGTCACACGATCTCAACGGt                                                                                                                                                                                                                                            <  5:1037123/95‑1 (MQ=255)
                                                                                                                                                                                                               tCAATATTCCTTTCAGCCACTGGCTGCGTTCTGACAGTAAATCTCAGTGGCGACATGCCAGTGCCAGCTACAGCATGTCACACGATCTCAACGGt                                                                                                                                                                                                                                            >  6:1037123/1‑95 (MQ=255)
                                                                                                                                                                                                               tCAATATTCCTTTCAGCCACTGGCTGCGTTCTGACAGTAAATCTCAGTGGCGACATGCCAGTGCCAGCTACAGCATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACGGTACGTTGCt                                                                                                                                                                                                      <  8:1097320/133‑1 (MQ=255)
                                                                                                                                                                                                               tCAATATTCCTTTCAGCCACTGGCTGCGTTCTGACAGTAAATCTCAGTGGCGACATGCCAGTGCCAGCTACAGCATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACGGTACGTTGCTGGAAGACAACAACCt                                                                                                                                                                                       <  3:115534/148‑1 (MQ=255)
                                                                                                                                                                                                                cAATATTCCTTTCAGCCACTGGCTGCGTTCTGACAGTAAATCTCAGTGGCGACATGCCAGTGCCAGCTACAGCATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACGGTACGTTGCTGGAAGACAACAACCTc                                                                                                                                                                                      <  5:908776/148‑1 (MQ=255)
                                                                                                                                                                                                                    atTCCTTTCAGCCACTGGCTGCGTTCTGACAGTAAATCTCAGTGGCGACATGCCAGTGCCAGCTACAGCATGTCACACGATCTCa                                                                                                                                                                                                                                                 <  7:789056/85‑1 (MQ=255)
                                                                                                                                                                                                                    atTCCTTTCAGCCACTGGCTGCGTTCTGACAGTAAATCTCAGTGGCGACATGCCAGTGCCAGCTACAGCATGTCACACGATCTCa                                                                                                                                                                                                                                                 >  8:789056/1‑85 (MQ=255)
                                                                                                                                                                                                                      tCCTTTCAGCCACTGGCTGCGTTCTGACAGTAAATCTCAGTGGCGACATGCCAGTGCCAGCTACAGCATGTCACACGATCTCAACGGTCGGATGACCAATc                                                                                                                                                                                                                               >  5:156737/1‑101 (MQ=255)
                                                                                                                                                                                                                      tCCTTTCAGCCACTGGCTGCGTTCTGACAGTAAATCTCAGTGGCGACATGCCAGTGCCAGCTACAGCATGTCACACGATCTCAACGGTCGGATGACCAATc                                                                                                                                                                                                                               <  6:156737/101‑1 (MQ=255)
                                                                                                                                                                                                                      tCCTTTCAGCCACTGGCTGCGTTCTGACAGTAAATCTCAGTGGCGACATGCCAGTGCCAGCTACAGCATGTCACACGATCTCAACGGTCGGATGACCAATc                                                                                                                                                                                                                               <  7:560694/101‑1 (MQ=255)
                                                                                                                                                                                                                          ttCAGCCACTGGCTGCGTTCTGACAGTAAATCTCAGTGGCGACATGCCAGTGCCAGCTACAGCATGTCACACGATCTCAACGGTCGGa                                                                                                                                                                                                                                        >  5:469066/1‑88 (MQ=255)
                                                                                                                                                                                                                          ttCAGCCACTGGCTGCGTTCTGACAGTAAATCTCAGTGGCGACATGCCAGTGCCAGCTACAGCATGTCACACGATCTCAACGGTCGGa                                                                                                                                                                                                                                        <  6:469066/88‑1 (MQ=255)
                                                                                                                                                                                                                          ttCAGCCACTGGCTGCGTTCTGACAGTAAATCTCAGTGGCGACATGCCAGTGCCAGCTACAGCATGTCACACGATCTCAACGGTCGGATGACCAATctggctg                                                                                                                                                                                                                         >  7:693488/1‑103 (MQ=255)
                                                                                                                                                                                                                          ttCAGCCACTGGCTGCGTTCTGACAGTAAATCTCAGTGGCGACATGCCAGTGCCAGCTACAGCATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACGGTACg                                                                                                                                                                                                           <  8:230132/117‑1 (MQ=255)
                                                                                                                                                                                                                          ttCAGCCACTGGCTGCGTTCTGACAGTAAATCTCAGTGGCGACATGCCAGTGCCAGCTACAGCATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACGGTACGTTGCTGGAAGAcaacaa                                                                                                                                                                                          <  8:689723/134‑1 (MQ=255)
                                                                                                                                                                                                                          ttCAGCCACTGGCTGCGTTCTGACAGTAAATCTCAGTGGCGACATGCCAGTGCCAGCTACAGCATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACGGTACGTTGCTGGAAGACAACAACCTCAGCTATAGCGTGCAAACCGGCTATGCCGGGGGAGGCGATGGAAATAGCGGAAGTACAGGCTACGCCACGCTGAATTATCGCGGTGGTTACGGCAATGCCAATATCGGTTACAGCCATAGCGATGATATTaa                                                   >  4:313782/1‑269 (MQ=255)
                                                                                                                                                                                                                          ttCAGCCACTGGCTGCGTTCTGACAGGAAATCTCAGTGGCGACATGACAGTGCCAGCTACAGCATGTCACACGATCTCAACGGTCGGATGAACAATCTGGCTGGTGTATACGGTACGTTGCTGGAAGACAACAACCTCAGCTATAGCGTGCAAACCGGCTATGCCGGGGGAGGCGATGGAAATAGCGGAAGTACAGGCTACGCCACGCTGAATTATCGCGGTGGTTACGGCAATGCCAATATCGGTTACAGCCATAGCGATGATATTAAGCAGct                                             <  2:401674/275‑1 (MQ=255)
                                                                                                                                                                                                                           tCAGCCACTGGCTGCGTTCTGACAGTAAATCTCAGTGGCGACATGCCAGTGCCAGCTACAGCATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACGGTACGTTGCTGGAAGACAACAACCTCAGCTATAGCGTGCAAACCGGCTATGCCGGGGGAGGCGATGGAAATAGCGGAAGTACAGGCTACGCCACGCTGAATTATCGCGGTGGTTACGGCAATGCCAATATCGGTTACAGCCATAGCGAGGATATTAagcagcct                                            >  2:8634/1‑273 (MQ=255)
                                                                                                                                                                                                                                CACTGGCTGCGTTCTGACAGTAAATCTCAGTGGCGACATGCCAGTGCCAGCTACAGCATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACGGTACGTTGCTGGAAGACAACAACCTCAGCta                                                                                                                                                                                 <  5:85924‑M1/137‑1 (MQ=255)
                                                                                                                                                                                                                                 ACTGGCTGCGTTCTGACAGTAAATCTCAGTGGCGACATGCCAGTGCCAGCTACAGCATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACGGTACg                                                                                                                                                                                                           <  7:550919‑M1/110‑1 (MQ=255)
                                                                                                                                                                                                                                 ACTGGCTGCGTTCTGACAGTAAATCTCAGTGGCGACATGCCAGTGCCAGCTACAGCATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACGGTACg                                                                                                                                                                                                           >  8:550919‑M1/1‑110 (MQ=255)
                                                                                                                                                                                                                                 ACTGGCTGCGTTCTGACAGTAAATCTCAGTGGCGACATGCCAGTGCCAGCTACAGCATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACGGTACGTTGCTGGAAGacaaca                                                                                                                                                                                           >  7:210195‑M1/1‑126 (MQ=255)
                                                                                                                                                                                                                                  CTGGCTGCGTTCTGACAGTAAATCTCAGTGGCGACATGCCAGTGCCAGCTACAGCATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACGGTACGTTGCTGGAAGACAACAACCTCAGCTATAGCGTGCAAACCGGCTATGCCGGGGGAGGCGATGGAAATAGCGGAAGTACAGGc                                                                                                                          <  2:561931‑M1/190‑1 (MQ=255)
                                                                                                                                                                                                                                   tggcatgtcgccactgagatttactgtcagaacgcagccagtgCCAGCTACAGCATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACGGTACGTTGCTGGAAGACAACAACCTc                                                                                                                                                                                      >  7:319846‑M2/44‑129 (MQ=255)
                                                                                                                                                                                                                                    ggcatgtcgccactgagatttactgtcagaacgcagccagtgCCAGCTACAGCa                                                                                                                                                                                                                                                                >  5:721584‑M2/43‑54 (MQ=255)
                                                                                                                                                                                                                                      ntgagTTCTGACAGTAAATCTCAGTGGCGACATGTCATTGCCAGCTACAGCATGTCAAACGATCTCAACGGTCTGATGACCAATCTGGCTGGTGTATACGGGACGGTGCTGGAAGACAACAACCTCAGCTATAGCGTGCAAACCGGCTATGCGGGGGGAGGCGATGGAAATAGCGGAAATACAGGCTACg                                                                                                                      <  4:561930/186‑1 (MQ=255)
                                                                                                                                                                                                                                      CTGCGTTCTGACAGTAAATCTCAGTGGCGACATGCCAGTGCCAGCTACAGCATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACGGTACGTTGCTGGAAGACAACAACCTCAGCTATAGCGTGCAAACCGGCt                                                                                                                                                                >  7:645619‑M1/1‑148 (MQ=255)
                                                                                                                                                                                                                                      CTGCGTTCTGACAGTAAATCTCAGTGGCGACATGCCAGTGCCAGCTACAGCATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACGGTACGTTGCTGGAAGACAACAACCTCAGCTATAGCGTGCAAACCGGCTATGCCGGGGGAGGCGATGGAAATAGCGGAAGTa                                                                                                                               >  4:97857‑M1/1‑181 (MQ=255)
                                                                                                                                                                                                                                        gCGTTCTGACAGTAAATCTCAGTGGCGACATGCCAGTGCCAGCTACAGCATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACGGTACGTTGCTGGAAGACAACAACCTCAGCTATAGCGTGCAAACCGGCTAt                                                                                                                                                              >  8:361640/1‑148 (MQ=255)
                                                                                                                                                                                                                                          tcgccactgagatttactgtcagaacgcagccagtgGCAGCTACAGCATGTCACACGATCTCAACGGTCgg                                                                                                                                                                                                                                         >  5:231625‑M2/37‑71 (MQ=255)
                                                                                                                                                                                                                                          tcgccactgagatttactgtcagaacgcagccagtgGCAGCTACAGCATGTCACACGATCTCAACGGTCgg                                                                                                                                                                                                                                         <  6:231625‑M2/35‑1 (MQ=255)
                                                                                                                                                                                                                                            tCTGACAGTAAATCTCAGTGGCGACATGCCAGTGCCAGCTACAGCATGTCACACGATCTCAACGGTCGGATGACCAAt                                                                                                                                                                                                                                >  7:878484/1‑78 (MQ=255)
                                                                                                                                                                                                                                            tCTGACAGTAAATCTCAGTGGCGACATGCCAGTGCCAGCTACAGCATGTCACACGATCTCAACGGTCGGATGACCAAt                                                                                                                                                                                                                                <  8:878484/78‑1 (MQ=255)
                                                                                                                                                                                                                                            tCTGACAGTAAATCTCAGTGGCGACATGCCAGTGCCAGCTACAGCATGTCACACGATCTCAACGGTCGGATGACCAAt                                                                                                                                                                                                                                <  7:1086195/78‑1 (MQ=255)
                                                                                                                                                                                                                                              ggACAGTAAATCTCAGTGGCGACATGCCAGGGCCAGCTACAGCATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACGGTACGTTGCTGGAAGACAACAACCTCAGCTATAGCGTGCAAACCGGCTATGCCGGGGGAGGCGATGGAAATAGCGGAAGTACAGGCTACGCCACGCTGAATTATCGCGGTGGTTACGGCAATGCCAATATCGGTTACAGCCATAGCGATGATATTAAGCAGCTCTATTACGGAGTCAGCGGTg                         <  2:474980/274‑1 (MQ=255)
                                                                                                                                                                                                                                                aCATTAAATCTCAGTGGCGACATACCAGTGCCAGCTACAGCATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACGGTACGTTGCt                                                                                                                                                                                                      >  5:956691/1‑100 (MQ=255)
                                                                                                                                                                                                                                                aCAGTAACTCTCAGTGGCGACATGCCAGTGCCAGCTACAGCATGTCACACGATCTGAACGGTCGGATGACCCATCTGGCTGGTGTATACGGTACGTTGCTGGAAGACAACAACCTCAGCTATAg                                                                                                                                                                              <  8:160289/124‑1 (MQ=255)
                                                                                                                                                                                                                                                  aGTAAATCTCAGTGGCGACATGCCAGTGCCAGCTACAGCATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACGGTACGTTGCTGGAAGACAACAACCTCa                                                                                                                                                                                     <  5:264536/115‑1 (MQ=255)
                                                                                                                                                                                                                                                  aGTAAATCTCAGTGGCGACATGCCAGTGCCAGCTACAGCATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACGGTACGTTGCTGGAAGACAACAACCTCa                                                                                                                                                                                     >  6:264536/1‑115 (MQ=255)
                                                                                                                                                                                                                                                   gTAAATCTCAGTGGCGACATGCCAGTGCCAGCTACAGCATGTCACACGATCTCAACGGTCGGATGACCAATctg                                                                                                                                                                                                                             <  6:921290/74‑1 (MQ=255)
                                                                                                                                                                                                                                                    tAAATCTCAGTGGCGACATGCCAGTGCCAGCTACAGCATGTCACACGATCTCAACGGTCGGATGACCAATctggc                                                                                                                                                                                                                           >  7:159713/1‑75 (MQ=255)
                                                                                                                                                                                                                                                     aaaTCTCAGTGGCGACATGCCAGTGCCAGCTACAGCATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACGGt                                                                                                                                                                                                              >  6:891884/1‑87 (MQ=255)
                                                                                                                                                                                                                                                     aaaTCTCAGTGGCGACATGCCAGTGCCAGCTACAGCATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACGGTACGTTGCTGGAAGACAACAACCTCAGCTATAGCg                                                                                                                                                                            >  8:1025231/1‑121 (MQ=255)
                                                                                                                                                                                                                                                        tctcAGTGGCGACATGCCAGTGCCAGCTACAGCATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACGGTACGTTGCTGGAAGtca                                                                                                                                                                                              <  6:354954/100‑1 (MQ=255)
                                                                                                                                                                                                                                                           cAGTGGCGACATGCCAGTGCCAGCTACAGCATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACgg                                                                                                                                                                                                               <  6:386970/80‑1 (MQ=255)
                                                                                                                                                                                                                                                           cAGTGGCGACATGCCAGTGCCAGCTACAGCATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACgg                                                                                                                                                                                                               >  5:386970/1‑80 (MQ=255)
                                                                                                                                                                                                                                                           cAGTGGCGACATGCCAGTGCCAGCTACAGCATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACGGTACGTTGCTGGAAGACAACAACCTCAg                                                                                                                                                                                    <  6:695821/107‑1 (MQ=255)
                                                                                                                                                                                                                                                           cAGTGGCGACATGCCAGTGCCAGCTACAGCATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACGGTACGTTGCTGGAAGACAACAACCTCAGCTATAGCGTGCAAACCGGCTATGCCGGGGGAGGCGATGGAAATAGCGGAAGTACAGGCTACGCCACGCTGAATTATCGCGGTGGTTACGGCAATGCCAATATCGGTTACAGCCATAGCGATGATATTAAGCAGCTCTATTACGGAGTCAGCGGTGGGGTACTGGCTc             >  2:31136/1‑274 (MQ=255)
                                                                                                                                                                                                                                                              tGGCGACATGCCAGTGCCAGCTACAGCATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACGGTACGTTGCTGGAAGACAACAACCTCAGCTATAGCGTGCAAACCGGCTATGCCGGGGGAGGCGAt                                                                                                                                               <  5:750072/141‑1 (MQ=255)
                                                                                                                                                                                                                                                                   aCATGCCAGTGCCAGCTACAGCATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACGGTACGTTGCTGGAAGACAACAACCTCAGCTATAGCGTGCaa                                                                                                                                                                       <  5:803842/112‑1 (MQ=255)
                                                                                                                                                                                                                                                                   aCATGCCAGTGCCAGCTACAGCATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACGGTACGTTGCTGGAAGACAACAACCTCAGCTATAGCGTGCAAACCGGCTATGCCGGGGGAGGCGATGGAAATAGCGGa                                                                                                                                   >  6:412957/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                      tgccagtgccagCTACAGCATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACGGTACGTTGCTGGAAGACAACAACCTCAGCTATAGCGTGCAAACCGGCTATGCCGCGGGAGGCGATGGAAATAGCGGAAGTACAGGCTACGCCACGCTGAATTATCGCGGTGGTTGCGTCAATGCCAATATCGGTTACAGCCATAGCGATGATATTAAGCAGCCCTATTACGGAGTCAGCGGTGGGGTACTGGCTCATGGCAATGGc  >  1:148723/1‑274 (MQ=255)
                                                                                                                                                                                                                                                                         cagtgCCAGCTACAGCATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACGGTACGTTGCTGGAAGACAACAACCTCAGCTATAGCGTGCAAACCGGCTATGCCGGGGGAGGCGATGGAAATAGCGGAAGTACa                                                                                                                             <  7:955280‑M2/143‑1 (MQ=255)
                                                                                                                                                                                                                                                                          agtgCCAGCTACAGCATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACGGTACGTTGCTGGAAGACAACAACCTCAGCTATAGCGTGCAAACCGGCTATGCCGGGGGAGGCGATGGAAATAGCGGAAGTACAg                                                                                                                            <  5:412957‑M2/144‑1 (MQ=255)
                                                                                                                                                                                                                                                                           gtgCCAGCTACAGCATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACGGTACGTTGCTGGAAGACAACAACCTCAGCTATAGCGTGCAAACCGGCTATg                                                                                                                                                             <  8:61481‑M2/111‑1 (MQ=255)
                                                                                                                                                                                                                                                                              ccagCTACAGCATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACGGTACGTTGCTGGAAGACAACAACCTCAGCTATAGCGTGCAAACCGGCTATGCCGGGGGAGGCGATGGAAATAGCGGAAGTACAGGCTa                                                                                                                        <  7:613383/148‑1 (MQ=255)
                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                         
TACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGATATCAACTGGACGCTCAGCTATAGCCTGACGAAAAACGCCTGGCAAAAAGGACGGGATCAGATGTTAGCGCTTAACGTCAATATTCCTTTCAGCCACTGGCTGCGTTCTGACAGTAAATCTCAGTGGCGACATGCCAGTGCCAGCTACAGCATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACGGTACGTTGCTGGAAGACAACAACCTCAGCTATAGCGTGCAAACCGGCTATGCCGGGGGAGGCGATGGAAATAGCGGAAGTACAGGCTACGCCACGCTGAATTATCGCGGTGGTTACGGCAATGCCAATATCGGTTACAGCCATAGCGATGATATTAAGCAGCTCTATTACGGAGTCAGCGGTGGGGTACTGGCTCATGCCAATGGC  >  NC_000913/4546681‑4547216

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