New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 4511910 =83 (0.950)4 (0.050) 4/254 NT 4.9% coding (505/957 nt) fecD ferric citrate ABC transporter permease
?NC_000913 4511933 = 77 (0.920)coding (482/957 nt) fecD ferric citrate ABC transporter permease
Rejected: Frequency below/above cutoff threshold.

CGAGCGTGGTGGCGCGCGCATCGCCGAGTGCAAGGAGGTCGAGATCGCGGCAAAAACTCAGGCTCAGCGGCAGAAATAAAATCATCAGCGGGATGGCAATCTTCACAAAGCTCCAGTCACGGCCCCATAAGCTGCCGGTCAGCCACAGCAGGGCGTTGTTCACATCCTGCGGGCGCGAGAGCATCAGATAATCCGTCAGGCTGGCCCAGCATGCAGAAAGCGCCACGCCGGTGAGCGCCAGCTTCATCGGCTGGTGGGTCTTTGCCAGCATCTTCAGCAATATCAACCCCGCCAT  >  NC_000913/4511763‑4512057
                                                                                                                                                   |                                                                                                                                                   
gctgctggcctttgcgggcggcatggcggggttgatattgctgaagatgctggcaaagacccaccagccgatgaagctggcgctcaccggcgtggcgctttctgcatgctgggccagcctgacggattatctgatgctctcgcgcccGc                                                                                                                                                    >  2:36060‑M1/148‑149 (MQ=255)
cGAGCGTGGTGGCGCGCGCATCGCCGAGTGCAAGGAGGTCGAGATCGCGGCAAAAACTCAGGCTCAGCGGCAGAAATAAAATCATCAGCGGGATGGCAATCTTCACAAAGCTCCAGTCACGGCCCCATAAGCTGCCGGTCAGCCAcag                                                                                                                                                     >  7:19194/1‑148 (MQ=255)
 ctgcttgcctttgcgggcggcatggcggggttgatattgctgaagatgctggcaaagacccaccagccgatgaagctggcgctcaccggcgtggcgctttctgcatgctgggccagcctgacggattatctgatgctctcgcgcccGCa                                                                                                                                                   >  1:339515‑M1/147‑149 (MQ=255)
  tgctggcctttgcgggcggcatggcggggttgatattgctgaagatgctggcaaagacccaccagccgatgaagctggcgctcaccggcgtggcgctttctgcatgctgggccagcctgacggattatctgatgctctcgcgcccGCAg                                                                                                                                                  >  8:145062‑M1/146‑149 (MQ=255)
    ctggcctttgcgggcggcatggcggggttgatattgctgaagatgctggcaaagacccaccagccgatgaagctggcgctcaccggcgtggcgctttctgcatgctgggccagcctgacggattatctgatgctctcgcgcccGCAg                                                                                                                                                  >  8:131985‑M1/144‑147 (MQ=255)
    cGTGGTGGCGCGCGCATCGCCGAGTGCAAGGAGGTCGAGATCGCGGCAAAAACTCAGGCTCAGCGGCAGAAATAAAATCATCAGCGGGATGGCAATCTTCACAAAGCTCCAGTCACGGCCCCATAAGCTGCCGGTCAGCCACAGCAggg                                                                                                                                                <  7:96329/149‑1 (MQ=255)
    cGTGGTGGCGCGCGCATCGCCGAGTGCAAGGAGGTCGAGATCGCGGCAAAAACTCAGGCTCAGCGGCAGAAATAAAATCATCAGCGGGATGGCAATCTTCACAAAGCTCCAGTCACGGCCCCATAAGCTGCCGGTCAGCCACAGCAgg                                                                                                                                                 <  2:160203/148‑1 (MQ=255)
      tggtggCGCGCGCATCGCCGAGTGCAAGGAGGTCGAGATCGCGGCAAAAACTCAGGCTCAGCGGCAGAAATAAAATCATCAGCGGGATGGCAATCTTCACAAAGCTCCAGTCACGGCCCCATAAGCTGCCGGTCAGCCACAGCAGGGCg                                                                                                                                              >  2:357624/1‑149 (MQ=255)
      tggtggCGCGCGCATCGCCGAGTGCAAGGAGGTCGAGATCGCGGCAAAAACTCAGGCTCAGCGGCAGAAATAAAATCATCAGCGGGATGGCAATCTTCACAAAGCTCCAGTCACGGCCCCATAAGCTGCCGGTCAGCCACAGCAGGGCg                                                                                                                                              >  3:308752/1‑149 (MQ=255)
        gtggCGCGCGCATCGCCGAGTGCAAGGAGGTCGAGATCGCGGCAAAAACTCAGGCTCAGCGGCAGAAATAAAATCATCAGCGGGATGGCAATCTTCACAAAGCTCCAGTCACGGCCCCATAAGCTGCCGGTCAGCCACAGCAGGGCgtt                                                                                                                                            <  2:11356/149‑1 (MQ=255)
          ggCGCGCGCATCGCCGAGTGCAAGGAGGTCGAGATCGCGGCAAAAACTCAGGCTCAGCGGCAGAAATAAAATCATCAGCGGGATGGCAATCTTCACAAAGCTCCAGTCACGGCCCCATAAGCTGCCGGTCAGCCACAGCAGGGCgttgt                                                                                                                                          <  4:117212/149‑1 (MQ=255)
                   aTCGCCGAGTGCAAGGAGGTCGAGATCGCGGCAAAAACTCAGGCTCAGCGGCAGAAATAAAATCATCAGCGGGATGGCAATCTTCACAAAGCTCCAGTCACGGCCCCATAAGCTGCCGGTCAGCCACAGCAGGGCGTTGTTCACATCCt                                                                                                                                 >  4:225711/1‑149 (MQ=255)
                      gCCGAGTGCAAGGAGGTCGAGATCGCGGCAAAAACTCAGGCTCAGCGGCAGAAATAAAATCATCAGCGGGATGGCAATCTTCACAAAGCTCCAGTCACGGCCCCATAAGCTGCCGGTCAGCCACAGCAGGGCGTTGTTCACATCCTGc                                                                                                                               <  7:2356/148‑1 (MQ=255)
                       ccGAGTTCAAGGAGGTCGAGATCTCTGCAAAAACTCAGGCTCAGCGGCAGAAATAAAATCATCATCTGGATGGCAATCTTCACAAAGCTCCAGTCACGGCCCCATAAGCTGCCGGTCAGCCACAGCAGGTCGTTGTTCACATCCTGCgg                                                                                                                             <  7:222271/149‑1 (MQ=255)
                       ccGAGTGCAAGGAGGTCGAGATCGCGGCAAAAACTCAGGCTCAGCGGCAGAAATAAAATCATCAGCGGGATGGCAATCTTCACAAAGCTCCAGTCACGGCCCCATAAGCTGCCGGTCAGCCACAGCAGGGCGTTGTTCACATCCTGCgg                                                                                                                             >  4:174607/1‑149 (MQ=255)
                       ccGAGTGCAAGGAGGTCGAGATCGCGGCAAAAACTCAGGCTCAGCGGCAGAAATAAAATCATCAGCGGGATGGCAATCTTCACAAAGCTCCAGTCACGGCCCCATAAGCTGCCGGTCAGCCACAGCAGGGCGTTGTTCACATCCTGCgg                                                                                                                             <  5:341555/149‑1 (MQ=255)
                         gcggggttgatattgctgacgatgctggccaagacccaccagccgatgaagctggcgctcaccggcgtggcgctttctgcatgctggaccagcctgacggattatctgatgcgctcgcgcccGCAGGGCGTTGTTCACATCCTGCGGgc                                                                                                                           >  4:7861‑M1/123‑149 (MQ=255)
                         gcggggttgatattgctgaagatgctggcaaagacccagcagccgatgaagctggcgctcaccggcgtggcgctttctgcatgcttggccagcctgacggattatctgatgctctcgcgcccGCAGGGCGTTGTTCACATCCTGCggg                                                                                                                            >  2:243558‑M1/123‑148 (MQ=255)
                              cAAGGAGGTCGAGATCGCGGCAAAAACTCAGGCTCAGCGGCAGAAATAAAATCATCAGCGGGATGGCAATCTTCACAAAGCTCCAGTCACGGCCCCATAAGCTGCCGGTCAGCCACAGCAGGGCGTTGTTCACATCCTGCGGGCGCga                                                                                                                       <  2:133809/148‑1 (MQ=255)
                               aaGGAGGTCGAGATCGCGGCAAAAACTCAGGCTCAGCGGCAGAAATAAAATCATCAGCGGGATGGCAATCTTCACAAAGCTCCAGTCACGGCCCCATAAGCTGCCGGTCGGCCACAGCAGGGCGTTGTTCACATCCTGCGGGCGCga                                                                                                                       >  5:209962/1‑147 (MQ=255)
                               aaGGAGGTCGAGATCGCGGCAAAAACTCAGGCTCAGCGGCAGAAATAAAATCATCAGCGGGATGGCAATCTTCACAAAGCTCCAGTCACGGCCCCATAAGCTGCCGGTCAGCCACAGCAGGGCGTTGTTCACATCCTGCGGGCGCga                                                                                                                       >  1:8658/1‑147 (MQ=255)
                                    ggTCGAGATCGCGGCAAAAACTCAGGCTCAGCGGCAGAAATAAAATCATCAGCGGGATGGCAATCTTCACAAAGCTCCAGTCACGGCCCCATAAGCTGCCGGTCAGCCACAGCAGGGCGTTGTTCACATCCTGCGGGCGCGAGAGCATc                                                                                                                <  8:376295/149‑1 (MQ=255)
                                     gTCGAGATCGCTGCAAAAACTTAGGCTCAGCGGCAGAAATAAAATCATCAGCGGGATGGCAATCTTCACAAAGCTCCAGTCACGGCCCCATAAGCTGCCGGTCAGCCACAGCAGGGCGTTGTTCACATCCTGCGGGCGCGAGAGCAt                                                                                                                 >  2:34100/1‑147 (MQ=255)
                                     gTCGAGATCGCGGCAAAAACTCAGGCTCAGCGGCAGAAATAAAATCATCAGCGGGATGGCAAGCTTCACAAAGCTCCAGTCACGGCCCCATAAGCTGCCGGTCAGCCACAGCAGGGCGTTGTTCACATCCTGCGGGCGCGAGAGCATCa                                                                                                               <  2:368007/149‑1 (MQ=255)
                                      tCGAGATCGCGGCAAAAACTCAGGCTCAGCGGCAGAAATAAAATCATCAGCGGGATGGCAATCTTCACAAAGCTCCAGTCACGGCCCCATAAGCTGCCGGTCAGCCACAGCAGGGCGTAGTTCACATCCAGCGGGCGCGAGAGCAt                                                                                                                 >  4:66441/1‑146 (MQ=255)
                                        gagaTCGCGGCAAAAACTCAGGCTCAGCGGCAGAAATAAAATCATCAGCGGGATGGCAATCTTCACAAAGCTCCAGTCACGGCCCCATAAGCTGCCGGTCAGCCACAGCAGGGCGTTGTTCACATCCTGCGGGCGCGAGAGCATCAGa                                                                                                             <  3:225711/148‑1 (MQ=255)
                                              gcgGCAAAAACTCAGGCTCAGCGGCAGAAATAAAATCATCAGCGGGATGGCAATCTTCACAAAGCTCCAGTCACGGCCCCATAAGCTGCCGGTCAGCCACAGCAGGGCGTTGTTCACATCCTGCGGGCGCGAGAGCATCAGATAATcc                                                                                                       <  1:245358/148‑1 (MQ=255)
                                              gcgGCAAAAACTCAGGCTCAGCGGCAGAAATAAAATCATCAGCGGGATGGCAATCTTCACAAAGCTCCAGTCACGGCCCCATAAGCTGCCGGTCAGCCACAGCAGGGCGTTGTTCACATCCTGCGGGCGCGAGAGCATCAGATAATCCg                                                                                                      >  5:365233/1‑149 (MQ=255)
                                                ggCAAAAACTCAGGCTCAGCGGCAGAAATAAAATCATCAGCGGGATGGCAATCTTCACAAAGCTCCAGTCACGGCCCCATAAGCTGCCGGTCAGCCACAGCAGGGCGTTGTTCACATCCTGCGGGCGCGAGAGCATCAGATAATCCGt                                                                                                     >  7:306614/1‑148 (MQ=255)
                                                   aaaaaCTCAGGCTCAGCGGCAGAAATAAAATCATCAGCGGGATGGCAATCTTCACAAAGCTCCAGTCACGGCCCCATAAGCTGCCGGTCAGCCACAGCAGGGCGTTGTTCACATCCTGCGGGCGCGAGAGCATCAGATAATCCGTCAgg                                                                                                 >  2:364747/1‑149 (MQ=255)
                                                     aaaCTCAGGCTCAGCGGCAGTAATAAAAGCATGAGCGGGATGGCAATCTTCACAAAGCTCCAGTCACGGCCCCATAAGCTGCCGGTCAGCCACAGCAGGGCGTTGTTCACATCCTGCGGGCGCGAGAGCATCAGATAATcc                                                                                                       <  2:306829/141‑1 (MQ=255)
                                                     aaaCTCAGGCTCAGCGGCAGAAATAAAATCATCAGCGGGATGGCAATCTTCACAAAGCTCCAGTCACGGCCCCATAAGCTGCCGGTCAGCCACAGCAGGGCGTTGTTCACATCCTGCGGGCGCGAGAGCATCAGATAATCCGACatcg                                                                                                >  1:306829/1‑145 (MQ=255)
                                                         tCAGGCTCAGCGGCAGAAATAAAATCATCAGCGGGATGGCAATCTTCACAAAGCTCCAGTCACGGCCCCATAAGCTGCCGGTCAGCCACAGCAGGGCGTTGTTCACATCCTGCGGGCGCGAGAGCATCAGATAATcc                                                                                                       <  3:1018/137‑1 (MQ=255)
                                                         tCAGGCTCAGCGGCAGAAATAAAATCATCAGCGGGATGGCAATCTTCACAAAGCTCCAGTCACGGCCCCATAAGCTGCCGGTCAGCCACAGCAGGGCGTTGTTCACATCCTGCGGGCGCGAGAGCATCAGATAATcc                                                                                                       >  4:1018/1‑137 (MQ=255)
                                                         tCAGGCTCAGCGGCAGAAATAAAATCATCAGCGGGATGGCAATCTTCACAAAGCTCCAGTCACGGCCCCATAAGCTGCCGGTCAGCCACAGCAGGGCGTTGTTCACATCCTGCGGGCGCGAGAGCATCAGATAATCCGTCAGGCTGGcc                                                                                           >  4:129406/1‑149 (MQ=255)
                                                          cAGGCTCAGCGGCAGAAATAAAATCATCAGCGGGATGGCAATCTTCACAAAGCTCCAGTCACGGCCCCATAAGCTGCCGGTCAGCCACAGCAGGGCGTTGTTCACATCCTGCGGGCGCGAGAGCATCAGATAATCCGTCAGGCTGGcc                                                                                           >  4:151056/1‑148 (MQ=255)
                                                           aGGCTCAGCGGCAGAAATAAAATCATCAGCGGGATGGCAATCTTCACAAAGCTCCAGTCACGGCCCCATAAGCTGCCGGTCAGCCACAGCAGGGCGTTGTTCACATCCTGCGGGCGCGAGAGCATCAGATAATCCGTCAGGCTGGcc                                                                                           <  1:364747/147‑1 (MQ=255)
                                                                cAGCGGCAGAAATAAAATCATCCGCGGGATGGCAATCTTCCCAAAGCTCCAGTCACGGCCCCATAAGCTGCCGGTCAGCCACAGCAGGGCGTTGTTCACATCCTGCGGGCGCGAGAGCATCAGATAATCCGTCAGGCTGGCCCAGCATg                                                                                    <  4:223631/149‑1 (MQ=255)
                                                                 aGCGGCAGAAATAAAATCATCAGCGGGATGGCAATCTTCACAAAGCTCCAGTCACGGCCCCATAAGCTGCCGGTCAGCCACAGCAGGGCGTTGTTCACATCCTGCGGGCGCGAGAGCATCAGATAATCCGTCAGGCTGGCCCAGCAt                                                                                     <  2:8658/147‑1 (MQ=255)
                                                                 aGCGGCAGAAATAAAATCATCAGCGGGATGGCAATCTTCACAAAGCTCCAGTCACGGCCCCATAAGCTGCCGGTCAGCCACAGCAGGGCGTTGTTCACATCCTGCGGGCGCGAGAGCATCAGATAATCCGTCAGGCTGGCCCAGCATGc                                                                                   <  3:129406/149‑1 (MQ=255)
                                                                    ggCAGAAATAAAATCATCAGCGGGATGGCAATCTTCACAAAGCTCCAGTCACGGCCCCATAAGCTGCCGGTCAGCCACAGCAGGGCGTTGTTCACATCCTGCGGGCGCGAGAGCATCAGATAATCCGTCAGGCTGGCCCAGCATGCa                                                                                  <  5:302044/147‑1 (MQ=255)
                                                                    ggCAGAAATAAAATCATCAGCGGGATGGCAATCTTCACAAAGCTCCAGTCACGGCCCCATAAGCTGCCGGTCAGCCACAGCAGGGCGTTGTTCACATCCTGCGGGCGCGAGAGCATCAGATAATCCGTCAGGCTGGCCCAGCATGCa                                                                                  >  6:302044/1‑147 (MQ=255)
                                                                     gCAGAAATAAAATCATCAGCGGGATGGCAATCTTCACAAAGCTCCAGTCACGGCCCCATAAGCTGCCGGTCAGCCACAGCAGGGCGTTGTTCACATCCTGCGGGCGCGAGAGCATCAGATAATCCGTCAGGCTGGCCCAGCATGCAGaa                                                                               <  3:174607/149‑1 (MQ=255)
                                                                     gCAGAAATAAAATCATCAGCGGGATGGCAATCTTCACAAAGCTCCAGTCACGGCCCCATAAGCTGCCGGTCAGCCACAGCAGGGCGTTGTTCACATCCTGCGGGCGCGAGAGCATCAGATAATCCGTCAGGCTGGCCCAGCATGCAGaa                                                                               <  4:308752/149‑1 (MQ=255)
                                                                       aGAAATAAAATCATCAGCGGGATGGCAATCTTCACAAAGCTCCAGTCACGGCCCCATAAGCTGCCGGTCAGCCACAGCAGGGCGTTGTTCACATCCTGCGGGCGCGAGAGCATCa                                                                                                               <  1:24392/115‑1 (MQ=255)
                                                                          aaTAAAATCATCAGCGGGATGGCAATCTTCACAAAGCTCCAGTCACGGCCCCATAAGCTGCCGGTCAGCCACAGCAGGGCGTTGTTCACATCCTGCGGGCGCGAGAGCATCAGATAATCCGTCAGGCTGGCCCAGCATGCAGAAAgcgc                                                                          <  1:357624/149‑1 (MQ=255)
                                                                          aaTAAAATCATCAGCGGGATGGCAATCTTCACAAAGCTCCAGTCACGGCCCCATAAGCTGCCGGTCAGCCACAGCAGGGCGTTGTTCACATCCTGCGGGCGCGAGAGCATCAGATAATCCGTCAGGCTGGCCCAGCATGCAGAAAgcgc                                                                          >  4:89431/1‑149 (MQ=255)
                                                                          aaTAAAATCATCAGCGGGATGGCAATCTTCACAAAGCTCCAGTCACGGCCCCATAAGCTGCCGGTCAGCCACAGCAGGGCGTTGTTCACATCCTGCGGGCGCGAGAGCATCAGATAATCCGTCAGGCTGGCCCAGCATGCAGAAAgcgc                                                                          <  7:98212/149‑1 (MQ=255)
                                                                            tAAAATCATCAGCGGGATGGCAATCTTAACAAAGCTCCAGTCACGGCCCCATAAGCTGCCGGTCAGCCACAGCAGGGCGTTGTTCACATCCTGCGGGCGCGCGAGCAGCAGATAATCCGTCAGGCTGGCCCAGCATGCAGAAAGCGCCa                                                                        >  6:347211/1‑149 (MQ=255)
                                                                              aaaTCATCAGCGGGATGGCAATCTTCACAAAGCTCCAGTCACGGCCCCATAAGCTGCCGGTCCGCCACAGCAGGGCGTTGTTCACATCCTGCGGGCGCGAGAGCATCAGATAATCCGTCAGGCTGGCCCAGCATGCAGAAAGCGCCAc                                                                       <  6:209962/148‑1 (MQ=255)
                                                                               aaTCATCAGCGGGATGGCAATCTTCACAAAGCTCCAGTCACGGCCCCATAAGCTGCCGGTCAGCCACAGCAGGGCGTTGTTCACATCCTGCGGGCGCGAGAGCATCAGATAATCCGTCAGGCTGGCCCAGCATGCAGAAAGCGCCACGc                                                                     >  3:11710/1‑149 (MQ=255)
                                                                                 tcatcaGCGGGATGGCAATCTTCACAAAGCTCCAGTCACGGCCCCATAAGCTGCCGGTCAGCCACAGCAGGGCGTTGTTCACATCCTGCGGGCGCGAGAGCATCAGATAATCCGTCAGGCTGGCCCAGCATGCAGAAAGCGCCACGCCg                                                                   <  1:34100/149‑1 (MQ=255)
                                                                                   atcaGCGGGATGGCAATCTTCACAAAGCTCCAGTCACGGCCCCATAAGCTGCCGGTCAGCCACAGCAGGGCGTTGTTCACATCCTGCGGGCGCGAGAGCATCAGATAATCCGTCAGGCTGGCCCAGCATGCAGAAAGCGCCACGCCGGt                                                                 >  2:97686/1‑149 (MQ=255)
                                                                                   atcaGCGGGATGGCAATCTTCACAAAGCTCCAGTCACGGCCCCATAAGCTGCCGGTCAGCCACAGCAGGGCGTTGTTCACATCCTGCGGGCGCGAGAGCATCAGATAATCCGTCAGGCTGGCCCAGCATGCAGAAAGCGCCACGCCGGt                                                                 >  5:261651/1‑149 (MQ=255)
                                                                                    tcaGCGGGATGGCAATCTTCACAAAGCTCCAGTCACGGCCCCATAAGCTGCCGGTCAGCCACAGCAGGGCGTTGTTCACATCCTGCGGGCGCGAGAGCATCAGATAATCCGTCAGGCTGGCCCAGCATGCAGAAAGCGCCACGCCGGt                                                                 >  3:93620/1‑148 (MQ=255)
                                                                                    caccggcgtggcgctttctgcatgctgggccagcctgacggattatctgatgctctcgcgcccGCAGGGCGTTGTTCACATCCTGCGGgcgc                                                                                                                         <  1:396800‑M1/29‑1 (MQ=255)
                                                                                    caccggcgtggcgctttctgcatgctgggccagcctgacggattatctgatgctctcgcgcccGCAGGGCGTTGTTCACATCCTGCGGgcgc                                                                                                                         >  2:396800‑M1/64‑92 (MQ=255)
                                                                                      aGCGGGATGGCAATCTTCACAAAGCTCCAGTCACGGCCCCATAAGCTGCCGGTCAGCCACAGCAGGGCGTTGTTCACATCCTGCGGGCGCGAGAGCATCAGATAATCCGTCAGGCTGGCCCAGCATGCAGAAAGCGCCACGCCGGTGa                                                               <  6:197292/148‑1 (MQ=255)
                                                                                       gCGGGATGGCAATCTTCACAAAGCTCCAGTCACGGCCCCATAAGCTGCCGGTCAGCCACAGCAGGGCGTTGTTCACATCCTGCGGGCGCGAGAGCATCAGATAATCCGTCAGGCTGGCCCAGCATGCAGAAAGCGCCACGCCGGTGAgc                                                             >  1:101757/1‑149 (MQ=255)
                                                                                           gATGGCAATCTTCACAAAGCTCCAGTCACGGCCCCATAAGCTGCCGGTCAGCCACAGCAGGGCGTTGTTCACATCCTGCGGGCGCGAGAGCATCAGATAATCCGTCAGGCTGGCCCAGCATGCAGAACGCGCCACGCCGGTGAgcgc                                                           >  7:152637/1‑147 (MQ=255)
                                                                                             tGGCAATCTTCACAAAGCTCCAGTCACGGCCCCATAAGCTGCCGGTCAGCCACAGCAGGGCGTTGTTCACATCCTGCGGGCGCGAGAGCATCAGATAATCCGTCAGGCTGGCCCAGCATGCAGAAAGCGCCACGCCGGTGAGCGCCAGc                                                       <  3:89431/149‑1 (MQ=255)
                                                                                               gCAATCTTCACAAAGCTCCAGTCACGGCCCCATAAGCTGCCGGTCAGCCACAGCAGGGCGTTGTTCACATCCTGCGGGCGCGAGAGCATCAGATAATCCGTCAGGCTGGCCCAGCATGCAGAAAGCGCCACGCCGGTGAGCGCCAGCt                                                      <  4:93620/148‑1 (MQ=255)
                                                                                               gCAATCTTCACAAAGCTCCAGTCACGGCCCCATAACCTGCCGGTCAGCCACAGCAGGGCGTAGTTCACATCCTGCGGGCGCGAGAGCATCAGATAATCCGACAGGCTGCCCCAGCATGCAGAAAGCGCCACGCCGGTGAGCGCCCGCtt                                                     >  1:331255/1‑149 (MQ=255)
                                                                                                  aTCTTCACAACGCTCCAGTCACGGCCCCCTCCGCTGCCGGGCAGCCACCGCAGGGCGTTGTTCACATCCTGCGGGCGCGAGAGCAGCAGATAATCCGTCAGGCTGGCCCAGCATGCCGAAAGCGCCACGCCGGTGAGCGCCAGCTTCAt                                                  <  6:155840/149‑1 (MQ=255)
                                                                                                      tCACAAAGCTCCAGTCACGGCCCCATAAGCTGCCGGTCAGCCACAGCAGGGCGTTGTTCACATCCTGCGGGCGCGAGAGCATCAGATAATCCGTCAGGCTGGCCCAGCATGCAGAAAGCGCCACGCCGGTGAGCGCCAGCTTCATCGGc                                              >  2:103282/1‑149 (MQ=255)
                                                                                                         caAAGCTCCAGTCACGGCCCCATAAGCTGCCGGTCAGCCACAGCAGGGCGTTGTTCACATCCTGCGGGCGCGAGAGCATCAGATAATCCGTCAGGCTGGCCCAGCATGCAGAAAGCGCCACGCCGGTGa                                                               >  5:25595/1‑129 (MQ=255)
                                                                                                         caAAGCTCCAGTCACGGCCCCATAAGCTGCCGGTCAGCCACAGCAGGGCGTTGTTCACATCCTGCGGGCGCGAGAGCATCAGATAATCCGTCAGGCTGGCCCAGCATGCAGAAAGCGCCACGCCGGTGa                                                               <  6:25595/129‑1 (MQ=255)
                                                                                                          aaaGCTCCAGTCACGGCCCCATAAGCTGCCGGTCAGCCACAGCAGGGCGTTGTTCACATCCTGCGGGCGCGAGAGCATCATATAATCCGTCAGGCTGGCCCAGCATGCAGAAAGCGCCACGCCGGTGAGCGCCAGCTTCATCGGCtggt                                          <  6:261651/149‑1 (MQ=255)
                                                                                                          aaaGCTCCAGTCACGGCCCCATAAGCTGCCGGTCAGCCACAGCAGGGCGTTGTTCACATCCTGCGGGCGCGAGAGCATCAGATAATCCGTCAGGCTGGCCCAGCAt                                                                                     >  4:370436/1‑106 (MQ=255)
                                                                                                          aaaGCTCCAGTCACGGCCCCATAAGCTGCCGGTCAGCCACAGCAGGGCGTTGTTCACATCCTGCGGGCGCGAGAGCATCAGATAATCCGTCAGGCTGGCCCAGCAt                                                                                     <  3:370436/106‑1 (MQ=255)
                                                                                                          aaaGCTCCAGTCACGGCCCCATAAGCTGCCGGTCAGCCACAGCAGGGCGTTGTTCACATCCTGCGGGCGCGAGAGCATCAGATAATCCGTCAGGCTGGCCCAGCATGCAGAAAGCGCCACGCCGGTGAGCGCCCGCTTCATCGGCtggt                                          >  7:109416/1‑149 (MQ=255)
                                                                                                          aaaGCTCCAGTCACGGCCCCATAAGCTGCCGGTCAGCCACAGCAGGGCGTTGTTCACATCCTGCGGGCGCGAGAGCATCAGATAATCCGTCAGGCTGGCCCAGCATGCAGAAAGCGCCACGCCGGTGAGCGCCAGCTTCATCGGCtggt                                          <  1:103282/149‑1 (MQ=255)
                                                                                                           aaGCTCCAGTCACGGCCCCATAAGCTGCCGGTCAGCCACAGCAGGGCGTTGTTCACATCCTGCGGGCGCGAGAGCATCAGATAATCCGTCAGGCTGGCCCAGCATGCAGAAAGCGCCACGCCGGTGAGCGCCAGCTTCATCGGCtggtg                                         <  1:97686/149‑1 (MQ=255)
                                                                                                             gCTCCAGTCACGGCCCCATAAGCTGCCGGTCAGCCACAGCAGGGCGTTGTTCACATCCTGCGGGCGCGAGAGCATCAGATAATCCGTCAGGCTGGCCCAGCATGCAGAAAGCGCCACGCCGGTGAGCGCCAGCTTCATCGGCTGGTggg                                       <  8:19194/149‑1 (MQ=255)
                                                                                                                ccAGTCACGGCCCCATAAGCTGCCGGTCAGCCACAGCAGGGCGTTGTTCACATCCTGCGGGCGCGAGAGCATCAGATAATCCGTCAGGCTGGCCCAGCATGCAGAAAGCGCCACGCCGGTGAGCGCCAGCTTCATCGGCTGGTGGGTCt                                    >  3:25246/1‑149 (MQ=255)
                                                                                                                   gTCACGGCCCCATAAGCTGCCGGTCAGCCACAGCAGGGCGGTGTTCACATCCTGCGGGCGCGAGAGCATCAGTTAATCCGTCAGGCTGGCCCAGCATTCAGAACGCGCCACGCCGGTGAGCGCCAGCTTCATCGGCTGGTGGGTCTTTg                                 <  5:383121/149‑1 (MQ=255)
                                                                                                                    tCACGGCCCCATAAGCTGCCGGTCAGCCACAGCAGGGCGTTGTTCACATCCTGCGGGCGCGAGAGCATCAGATAATCCGTCAGGCTGGCCCAGCATGCAGAAAGCGCCACGCCGGTGAGCGCCAGCTTCATCGGCTGGTGGGTCTTTGc                                >  2:18548/1‑149 (MQ=255)
                                                                                                                    tCACGGCCCCATAAGCTGCCGGTCAGCCACAGCAGGGCGTTGTTCACATCCTGCGGGCGCGAGAGCATCAGATAATCCGTCAGGCTGGCCCAGCATGCAGAAAGCGCCACGCCGGTGAGCGCCAGCTTCATCGGCTGGTGGGTCTTTGc                                <  6:365233/149‑1 (MQ=255)
                                                                                                                    tCACGGCCCCATAAGCTGCCGGTCAGCCACAGCAGGGCGTTGTTCACATCCTGCGGGCGCGAGAGCATCAGATAATCCGTCAGGCTGGCCCAGCATGCAGAAAGCGCCACGCCGGTGAGCGCCAGCTTCATCGGCTGGTGGGTCTTTGc                                >  6:366624/1‑149 (MQ=255)
                                                                                                                      aCGGCCCCATAAGCTGCCGGTCAGCCACAGCAGGGCGTTGTTCACATCCTGCGGGCGCGAGAGCATCAGATAATCCGTCAGGCTGGCCCAGCATGCAGAAAGCGCCACGCCGGTGAGCGCCAGCTTCATCGGCTGGTGGGTCTTTGCCa                              <  8:232909/149‑1 (MQ=255)
                                                                                                                       cGGCCCCATAAGCTGCCGGTCAGCCACAGCAGGGCGTTGTTCACATCCTGCGGGCGCGAGAGCATCAGATAATCCGTCAGGCTGGCCCAGCATGCAGAAAGCGCCACGCCGGTGAGCGCCAGCTTCATCGGCTGGTGGGTCTTTGCCAg                             >  3:63095/1‑149 (MQ=255)
                                                                                                                        ggCCCCATAAGCTGCCGGTCAGCCACAGCAGGGCGTTGTTCACATCCTGCGGGCGCGAGAGCATCAGATAATCCGTCTGGCTGGCCCAGCATGCAGAAAGCGCCACGCCGGTGAGCGCCAGCTTCATCGGCTGGTGGGTCTTTGCCAGc                            <  7:374337/149‑1 (MQ=255)
                                                                                                                            ccATAAGCTGCCGGTCAGCCACAGCAGGGCGTTGTTCACATCCTGCGGGCGCGAGAGCATCAGATAATCCGTGAGGCTGGCCCAGCATGCAGAAAGCGCCACGCCGGTGAGCGCCAGCTTCATCGGCTGGTGGGTCTTTGCCAGCATCt                        >  3:344260/1‑149 (MQ=255)
                                                                                                                             cATAAGCTGCCGGTCAGCCACAGCAGGGCGTTGTTCACATCCTGCGGGCGCGAGAGCATCAg                                                                                                              >  4:30268/1‑62 (MQ=255)
                                                                                                                             cATAAGCTGCCGGTCAGCCACAGCAGGGCGTTGTTCACATCCTGCGGGCGCGAGAGCATCAg                                                                                                              <  3:30268/62‑1 (MQ=255)
                                                                                                                              ctAAGCTGCCGGTCAGCCACAGCAGGGCGTTGTTCACATCCTGCGGGCGCGAGAGCATCAGATAATCCGTCAGGCTGGCCCAGCATGCAGAAAGCGCCACGCCGGTGAGCGCCAGCTTCATCGGCTGGTGGGTCTTTGCCAGCATCtt                       >  5:261828/2‑148 (MQ=255)
                                                                                                                                 aGCTGCCGGTCAGCCACAGCAGGGCGTTGTTCACATCCTGCGGGCGCGAGAGCATCAGATAATCCGTCAGGCTGGCCCAGCATGCAGAAAGCGCCACGCCGGTGAGCGCCAGCTTCATCGGCTGGTGGGTCTTTGCCAGCATCTTCAg                    <  7:132752/148‑1 (MQ=255)
                                                                                                                                     gCCGGTCAGCCACAGCAGGGCGTTGTTCACATCCTGCGGGCGCGAGAGCATCAGATAATCCGTCAGGCTGGCCCAGCATGCAGAAAGCGCCACGCCGGTGAGCGCCAGCTTCATCGGCTGGTGGGTCTTTGCCAGCATCTTCAGCAata               >  2:52010/1‑149 (MQ=255)
                                                                                                                                         gTCAGCCACAGCAGGGCGTTGTTCACATCCTGCGGGCGCGAGAGCATCAGATAATCCGTCAGGCTGGCCCAGCATGCAGAAAGCGCCACGCCGGTGAGCGCCAGCTTCATCGGCTGGTGGGTCTTTGCCAGCATCTTCAGCAATATCaa           <  5:366624/149‑1 (MQ=255)
                                                                                                                                          tCAGCCACAGCAGGGCGTTGTTCACATCCTCCGGGCGCGAGAGCATCAGATAATCCGTCAGGCTGGCCCAGCATGCAGAAAGCGCCACGCCGGTGAGCGCCAGCTTCATCGGCTGGTGGGTCTTTGCCAGCATCTTCAGCAATATCAAc          >  2:40544/1‑149 (MQ=255)
                                                                                                                                             gCCACAGCAGGGCGTTGTTCACATCCTGCGGGCGCGAGAGCATCAGATAATCCGTCAGGCTGGCCCAGCATGCAGAAAGCGCCACGCCGGTGAGCGCCAGCTTCATCGGCTGGTGGGTCTTTGCCAGCATCTTCAGCAata               >  3:371613/1‑141 (MQ=255)
                                                                                                                                             gCCACAGCAGGGCGTTGTTCACATCCTGCGGGCGCGAGAGCATCAGATAATCCGTCAGGCTGGCCCAGCATGCAGAAAGCGCCACGCCGGTGAGCGCCAGCTTCATCGGCTGGTGGGTCTTTGCCAGCATCTTCAGCAata               <  4:371613/141‑1 (MQ=255)
                                                                                                                                             gCCACAGCAGGGCGTTGTTCACATCCTGCGGGCGCGAGAGCATCAGATAATCCGTCAGGCTGGCCCAGCATGCAGAAAGCGCCACGCCGGTGAGCGCCAGCTTCATCGGCTGGTGGGTCTTTGCCAGCATCTTCAGCAATATCAAcccc       >  4:362148/1‑149 (MQ=255)
                                                                                                                                              ccACAGCAGGGCGTTGTTCACATCCTGCGGGCGCGAGAGCATCAGATAATCCGTCAGGCTGGCCCAGCATGCAGAAAGCGCCACGCCGGTGAGCGCCAGCTTCATCGGCTGGTGGGTCTTTGCCAGCATCTTCAGCAATATCAACCCCg      >  6:71901/1‑149 (MQ=255)
                                                                                                                                                acaGCAGGGCGTTGTTCACATCCTGCGGGCGCGAGAGCATCAGATAATCCGTCAGGCTGGCCCAGCATGCAGAAAGCGCCACGCCGGTGAGCGCCAGCTTCATCGGCTGGTGGGTCTTTGCCAGCATCTTCAGCAATATCAACCCCGcc    <  4:11710/149‑1 (MQ=255)
                                                                                                                                                acaGCAGGGCGTTGTTCACATCCTGCGGGCGCGAGAGCATCAGATAATCCGTCAGGCTGGCCCAGCATGCAGAAAGCGCCACGCCGGTGAGCGCCAGCTTCATCGGCTGGTGGGTCTTTGCCAGCATCTTCAGCAATATCAACCCCGc     >  8:219422/1‑148 (MQ=255)
                                                                                                                                                  agcagGGCGTTGTTCACATCCTGCGGGCGCGAGAGCATCAGATAATCCGTCAGGCTGGCCCAGCATGCAGAAAGCGCCACGCCGGTGAGCGCCAGCTTCATCGGCTGGTGGGTCTTTGCCAGCATCTTCAGCAATATCAACCCCGCCAt  <  1:18548/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CGAGCGTGGTGGCGCGCGCATCGCCGAGTGCAAGGAGGTCGAGATCGCGGCAAAAACTCAGGCTCAGCGGCAGAAATAAAATCATCAGCGGGATGGCAATCTTCACAAAGCTCCAGTCACGGCCCCATAAGCTGCCGGTCAGCCACAGCAGGGCGTTGTTCACATCCTGCGGGCGCGAGAGCATCAGATAATCCGTCAGGCTGGCCCAGCATGCAGAAAGCGCCACGCCGGTGAGCGCCAGCTTCATCGGCTGGTGGGTCTTTGCCAGCATCTTCAGCAATATCAACCCCGCCAT  >  NC_000913/4511763‑4512057

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