New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? W3110S.gb = 230464016 (0.700)3 (0.150) 3/92 NT 15.2% coding (1448/2487 nt) napA nitrate reductase, periplasmic, large subunit
?W3110S.gb = 2304649 19 (0.920)coding (1439/2487 nt) napA nitrate reductase, periplasmic, large subunit

CGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTACACATGGTCCAGTAAACATTCAGTTTGCCGTCTTTCAGCGCACGGTCTTGTGCCACCGCGTGCAGACCGA  >  W3110S.gb/2304592‑2304718
                                                         |                                                                     
cgcCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTACACATGGTCCAGTa                                                          <  1:1704372/71‑1 (MQ=255)
   cAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTACACATGGTCCAGTaaa                                                        >  1:2073764/1‑70 (MQ=255)
      cccGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTACACATGGTCCAGTAAACAt                                                     >  1:246729/1‑70 (MQ=255)
         ggCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTACACATGGTCCAGTAAACATTCAg                                                 >  1:176655/1‑71 (MQ=255)
          gCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTACACATGGTCCAGTAAACATTCAGt                                                >  1:1189747/1‑71 (MQ=255)
            aTACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTACACATGGTCCAGTAAACATTCAGttt                                              <  1:1156851/71‑1 (MQ=255)
              aCGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTACACATGGTCCAGTAAACATTCAGTTTGc                                            >  1:1575234/1‑71 (MQ=255)
               cGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTACACATGGTCCAGTaa                                                         <  1:1764347/57‑1 (MQ=255)
                      cATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTACACATGGTCCAGTAAACATTCAGTTTGCCGTCt                                       <  1:619015/68‑1 (MQ=255)
                       aTTGATGTTCGGCCCGGCCTGCATGTTGTTGGTACACATGGTCCAGTAAACATTCAGTTTGCCGTCttt                                     >  1:1783679/1‑69 (MQ=255)
                         tGATGTTCGGCCCGGCCTGCATGTTGTTGGTACACATGGTCCAGTAAACATTCAGTTTGCCGTCTTTCAgc                                 >  1:1326454/1‑71 (MQ=255)
                                   cccggccTGCATGTTGTTGGTACACATGGTCCAGTAAACATTCAGTTTGCCGTCTTTCAGCGCACGGTCtt                       <  1:1219208/71‑1 (MQ=255)
                                   cccggccTGCATGTTGTTGGTACACATGGTCCAGTAAACATTCAGTTTGCCGTCTTTCAGCGCACGGTCtt                       <  1:412728/71‑1 (MQ=255)
                                    ccggccTGCATGTTGTTGGTACACATGGTCCAGTAAACATTCAGTTTGCCGTCTTTCAGCGCACGGTCTtg                      >  1:1620251/1‑71 (MQ=255)
                                          tGCATGTTGTTGGTACacatgcaggccgggccgaacatcaatgaagagcgtatgccgggctggcgcgatcc                <  1:1262788‑M2/71‑56 (MQ=255)
                                          tGCATGTTGTTGGTACACATGGTCCAGTAAACATTCAGTTTGCCGTCTTTCAGCGCACGGTCTTGTGCCAc                >  1:68704/1‑71 (MQ=255)
                                            cATGTTGTTGGTACACATGGTCCAGTAAACATTCAGTTTGCCGTCTTTCAGCGCACGGTCTTGTGCCAc                <  1:1332479/69‑1 (MQ=255)
                                               gttgttGGTACacatgcaggccgggccgaacatcaatg                                            >  1:1496368‑M2/1‑11 (MQ=255)
                                                ttgttgGTACACATGGTCCAGTAAACATTCAGTTTGCCGTCttt                                     >  1:1590109/1‑44 (MQ=255)
                                                ttgttgGTACACATGGTCCAGTAAACATTCAGTTTGCCGTCTTTCAGCGCACGGTCTTGTGCCACCGCGTg          <  1:1251176/71‑1 (MQ=255)
                                                  gttgGTACacatgcaggccgggccgaacatcaatgaagagcgtatgccgggctggcgctatccgc              <  1:1268810‑M2/65‑58 (MQ=255)
                                                  gttgGTACACATGGTCCAGTAAACATTCAGTTTGCCGTCTTTCAGCGCACGGTCTTGTGCCACCGCGTg          >  1:363296/1‑69 (MQ=255)
                                                       tACACATGGTCCAGTAAACATTCAGTTTGCCGTCTTTCAGCGCACGGTCTTGTGCCAcc               <  1:1338256/59‑1 (MQ=255)
                                                        acacATGGTCCAGTAAACATTCAGTTTGCCGTCTTTCAGCGCACGGTCTTGTGCCACCGCGTGCAGACCg   <  1:640045/70‑1 (MQ=255)
                                                         cacaTGGTCCAGTAAACATTCAGTTTGCCGTCTTTCAGCGCACGGTCTTGTGCCACCGCGTGCAGACCGa  <  1:1639553/70‑1 (MQ=255)
                                                         cacaTGGTCCAGTAAACATTCAGTTTGCCGTCTTTCAGCGCACGGTCTTGTGCCACCGCGTGCAGACCGa  <  1:1577464/70‑1 (MQ=255)
                                                         cacaTGGTCCAGTAAACATTCAGTTTGCCGTCTTTCAGCGCACGGTCTTGTGCCACCGCGTGCAGACCGa  <  1:407785/70‑1 (MQ=255)
                                                         |                                                                     
CGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTACACATGGTCCAGTAAACATTCAGTTTGCCGTCTTTCAGCGCACGGTCTTGTGCCACCGCGTGCAGACCGA  >  W3110S.gb/2304592‑2304718

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