New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 337397332 (0.980)6 (0.200) 4/224 3.5 19.3% coding (517/792 nt) nanR sialic acid‑inducible nan operon repressor
?NC_000913 = 3374015 21 (0.710)coding (475/792 nt) nanR sialic acid‑inducible nan operon repressor
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

GTAACCGTTGGGCGTGCGGCAATAAGCCAGTCGAGCAGGGCAACGTGGATCGCCATGAAGATTGGGTTACCGGGGATCTCCGCCAGCACGCGGTGGAAATCAACGTCTGAACGAATGAATGCCGCGTTGTTATCCAGCGACTGACTGTTGATTTCCAGTGCTTTTGCCAGCAAATCGATTTGCTCATCGGTGGCATGTTCAGCCGCATAGCGCACCAGACTGGATTCAAAGAACAGACGTAATTGTTCGAAATGGGCAATCCCACCGGGATGAGAAA  >  NC_000913/3373870‑3374146
                                                                                                                                                 |                                                                                                                                   
gTAACCGTTGGGCGTGCGGCAATAAGCCAGTCGAGCAGGGCAACGTGGATCGCCATGAAGATTGGGTTACCGGGGATCTCCGCCAGCACGCGGTGGAAATCAACGTCTGAACGAATGAATGCCGCGTTGTTATCCAGCGACTGACTGtt                                                                                                                                  <  5:106394/149‑1 (MQ=255)
    ccGTTGGGCGTGCGGCAATAAGCCAGTCGAGCAGGGCAACGTGGATCGCCATGAAGATTGGGTTACCGGGGATCTCCGCCAGCACGCGGTGGAAATCAACGTCTGAACGAATGAATGCCGCGTTGTTATCCAGCGACTGACTGTTGAtt                                                                                                                              >  2:118097/1‑149 (MQ=255)
     cGTTGGGCGTGCGGCAATAAGCCAGTCGAGCAGGGCAACGTGGATCGCCATGAAGATTGGGTTACCGGGGATCTCCGCCAACACGCGGTGGAAATCAACGTCTGAACGAATGAATGCCGCGTTGTTATCCAGCGACTGACTGTTGAttt                                                                                                                             >  5:101250/1‑149 (MQ=255)
          ggCGTGCGGCAATAAGCCAGTCGAGCAGGGCAACGTGGATCGCCATGAAGATTGGGTTACCGGGGATCTCCGCCAACACGCGGTGGAAATCAACGTCTGAACGAATGAATGCCGCGTTGTTATCCAGCGACTGACTGTTGATTTCCAGt                                                                                                                        <  6:101250/149‑1 (MQ=255)
           gCGTGCGGCAATAAGCCAGTCGAGCAGGGCAACGTGGATCGCCATGAAGATTGGGTTACCGGGGATCTCCGCCAGCACGCGGTGGAAATCAACGTCTGAACGAATGAATGCCGCGTTGTTATCCAGCGACTGACTGTTGATTTc                                                                                                                            <  6:84245/144‑1 (MQ=255)
           gCGTGCGGCAATAAGCCAGTCGAGCAGGGCAACGTGGATCGCCATGAAGATTGGGTTACCGGGGATCTCCGCCAGCACGCGGTGGAAATCAACGTCTGAACGAATGAATGCCGCGTTGTTATCCAGCGACTGACTGTTGATTTc                                                                                                                            >  5:84245/1‑144 (MQ=255)
                  gCAATAAGCCAGTCGAGCAGGGCAACGTGGATCGCCATGAAGATTGGGTTACCGGGGATCTCCGCCAGCACGCGGTGGAAATCAACGTCTGAACGAATGAATGCCGCGTTGTTATCCAGCGACTGACTGTTGATTTCCAGTGCTTTTGc                                                                                                                <  1:118097/149‑1 (MQ=255)
                    aaTAAGCCAGTCGAGCAGGGCAACGTGGATCGCCATGAAGATTGGGTTACCGGGGATCTCCGCCAGCACGCGGTGGAAATCAACGTCTGAACGAATGAATGCCGCGTTGTTATCCAGCGACTGACTGTTGATTTCCAGTGCTTTTGc                                                                                                                >  3:108294/1‑147 (MQ=255)
                    aaTAAGCCAGTCGAGCAGGGCAACGTGGATCGCCATGAAGATTGGGTTACCGGGGATCTCCGCCAGCACGCGGTGGAAATCAACGTCTGAACGAATGAATGCCGCGTTGTTATCCAGCGACTGACTGTTGATTTCCAGTGCTTTTGc                                                                                                                <  4:108294/147‑1 (MQ=255)
                                          aCGTGGATCGCCATGAAGATTGGGTTACCGGGGATCTCCGCCAGCACGCGGTGGAAATCAACGTCTGAACGAATGAATGCCGCGTTGTTATCCAGCGACTGACTGTTGATTTCCAGTGCTTTTGCCAGCAAATCGATTTGCTCATCGGt                                                                                        >  5:83511/1‑149 (MQ=255)
                                           cGTGGATCGCCATGAAGATTGGGTTACCGGGGATCTCCGCCAGCACGCGGTGGAAATCAACGTCTGAACGAATGAATGCCGCGTTGTTATCCAGCGACTGACTGTTGATTTCCAGTGCTTTTGCCAGCAAATCGATTTGCTCATCGGTg                                                                                       <  2:156096/149‑1 (MQ=255)
                                           cGTGGATCGCCATGAAGATTGGGTTACCGGGGATCTCCGCCAGCACGCGGTGGAAATCAACGTCTGAACGAATGAATGCCGCGTTGTTATCCAGCGACTGACTGTTGATTTCCAGTGCTTTTGCCAGCAAATCGATTTGCTCATCGGTg                                                                                       >  4:10379/1‑149 (MQ=255)
                                                 tCGCCATGAAGATTGGGTTACCGGGGATCTCCGCCAGCACGCGGTGGAAATCAACGTCTGAACGAATGAATGCCGCGTTGTTATCCAGCGACTGACTGTTGATTTCCAGTGCTTTTGCCAGCAAATCGATTTGCTCATCGGTGGCATGt                                                                                 >  1:97359/1‑149 (MQ=255)
                                                                ggTTACCGGGGATCTCCGCCAGCACGCGGTGGAAATCAACGTCTGAACGAATGAATGCCGCGTTGTTATCCAGCGACTGACTGTTGATTTCCAGTGCTTTTGCCAGCAAATCGATTTGCTCATCGGTGGCATGTTCAGCCGCATAgcg                                                                   >  2:35650/1‑148 (MQ=255)
                                                                  ttACCGGGGATCTCCGCCAGCACGCGGTGGAAATCAACGTCTGAACGAATGAATGCCGCGTTGTTATCCAGCGACTGACTGTTGAATTCCAGTGCTTTTGCCAGCAAATCGATTTGCTCATCGGTGGCATGTTCAGCCGCATAGCGCAc                                                                >  7:38201/1‑149 (MQ=255)
                                                                         ggATCTGAGCCAGCTCGCTGTGGAAATCTACGTGTGAACGAATGAATGGCGCGTTGTTATCCAGCGACTGACTGTTGATTTCAAGTGCTTTTGCCAGGTAATCGATTTGCTCATCTGTTGCTTGTTCAGCCGCTTTGCGCACCAGACTg                                                         <  8:38201/149‑1 (MQ=255)
                                                                                cGCCAGCACGCGGTGGAAATCAACGTCTGAACGAATGAATGCCGCGTTGTTATCCAGCGACTGACTGTTGATTTCCAGTGCTTTTGCCAGCAAATCGATTTGCTCATCGGTGGCATGTTCAGCCGCATAGCGCACCAGACTGGa                                                       <  3:10379/144‑1 (MQ=255)
                                                                                cGCCAGCACGCGGTGGAAATCAACGTCTGAACGAATGAATGCCGCGTTGTTATCCAGCGACTGACTGTTGATTTCCAGTGCTTTTGCCAGCAAATCGATTTGCTCATCGGTGGCATGTTCAGCCGCATAGCGCACCAGACTGGATTc                                                    <  2:97359/147‑1 (MQ=255)
                                                                                 gCCAGCACGCGGTGGAAATCAACGTCTGAACGAATGAATGCCGCGTTGTTATCCAGCGACTGACTGTTGATTTCCAGTGCTTTTGCCAGCAAATCGATTTGCTCATCGGTGGCATGTTCAGCCGCATAGCGCACCAGACTGGATTCaaa                                                 <  4:141918/149‑1 (MQ=255)
                                                                                         gcgGTGGAAATCAACGTCTGAACGAATGAATGCCGCGTTGTTATCCAGCGACTGACTGTTGATTTCCAGTGCTTTTGCCAGCAAATCGATTTGCTCATCGGTGGCATGTTCAGCCGCATAGCGCACCAGACTGGATTCAAAGAACAGAc                                         <  6:83511/149‑1 (MQ=255)
                                                                                                   tCAACGTCTGAACGAATGAATGCCGCGTTGTTATCCAGCGACTGACTGTTGATTTCCAGTGCTTTTGCCAGCAAATCGATTTGCTCATCGGTGGCATGTTCAGCCGCATAGCGCACCAGACTGGATTCAAAGAACAGACGTAATTGTTc                               <  1:35650/149‑1 (MQ=255)
                                                                                                      aCGTCTGAACGAATGAATGCCGCGTTGTTATCCAGCGACTGACTGTTGATTTCCAGTGCTTTTGCCAGCAAATCGATTTGCTCATCGGTGGCATGTTCAGCCGCATAGCGCACCAGACTGGATTCAAAGAACAGACGTaa                                     <  5:8877/140‑1 (MQ=255)
                                                                                                      aCGTCTGAACGAATGAATGCCGCGTTGTTATCCAGCGACTGACTGTTGATTTCCAGTGCTTTTGCCAGCAAATCGATTTGCTCATCGGTGGCATGTTCAGCCGCATAGCGCACCAGACTGGATTCAAAGAACAGACGTaa                                     >  6:8877/1‑140 (MQ=255)
                                                                                                                          cgcgTTGTTATCCAGCGACTGACTGTTGATTTCCAGTGCTTTTGCCAGCAAATCAATTTGCTCATCGGTGGCATGTTCAGCCGCATAGGGCACCAGACTGGATTCAAAGAACAGACGTAATTGGTCGAACTGGGCAATCCCACCGGGAt        >  8:111478/1‑149 (MQ=255)
                                                                                                                                   aTCCAGCGACTGACTGTTGATTTCCAGTGCTTTTGCCAGCAAATCGATTTGCTCATCGGTGGCATGTTCAGCCGCATAGCGCACCAGACTGGATTCAAAGAACAGACGTAATTGTTCGAAATGGGCAATCCCACCGGGATGAGaaa  >  5:125846/1‑146 (MQ=255)
                                                                                                                                   aTCCAGCGACTGACTGTTGATTTCCAGTGCTTTTGCCAGCAAATCGATTTGCTCATCGGTGGCATGTTCAGCCGCATAGCGCACCAGACTGGATTCAAAGAACAGACGTAATTGTTCGAAATGGGCAATCCCACCGGGATGAGaaa  <  6:125846/146‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                   
GTAACCGTTGGGCGTGCGGCAATAAGCCAGTCGAGCAGGGCAACGTGGATCGCCATGAAGATTGGGTTACCGGGGATCTCCGCCAGCACGCGGTGGAAATCAACGTCTGAACGAATGAATGCCGCGTTGTTATCCAGCGACTGACTGTTGATTTCCAGTGCTTTTGCCAGCAAATCGATTTGCTCATCGGTGGCATGTTCAGCCGCATAGCGCACCAGACTGGATTCAAAGAACAGACGTAATTGTTCGAAATGGGCAATCCCACCGGGATGAGAAA  >  NC_000913/3373870‑3374146

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: 
Reads not counted as support for junction
read_name Not counted due to insufficient overlap past the breakpoint.
read_name Not counted due to not crossing MOB target site duplication.