New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 = 382440426 (0.890)3 (0.120) 3/148 2.0 9.7% coding (588/1710 nt) yicH AsmA family protein
?NZ_CP009273 = 3824452 33 (1.300)coding (636/1710 nt) yicH AsmA family protein
Rejected: Frequency below/above cutoff threshold.

TACTGATTGAAGGCAGTATTGATAACGATCGCGTTACGCTGACTAACCTGGGTGCCGACATCGCCCGCGGGACATTAACCGGAAACGCGCAGCGTAACGCCGACGGCAGCTGGCAAGTGGAAAACCTGCGCATGGCGGATATACGTCTACAAAGCGAAAAATCGCTAACCGACTTCTT  >  NZ_CP009273/3824364‑3824541
                                                                                        |                                                                                         
tACTGATTGAAGGCAGTATTGATAACGATCGCGTTACGCTGACTAACCTGGGTGCCGACATCGCCCGCGGGACATTAACCGGAAACGCGc                                                                                          >  1:542145‑M2/1‑89 (MQ=255)
    gATTGAAGGCAGTATTGATAACGATCGCGTTACGCTGACTAACCTGGGTGCCGACATCGCCCGCGGGACATTAACCGGAAACGCGcagcg                                                                                      <  1:469620‑M2/90‑6 (MQ=255)
    gATTGAAGGCAGTATTGATAACGATCGCGTTACGCTGACTAACCTGGGTGCCGACATCGCCCGCGGGACATTAACCGGAAACGCGcagcg                                                                                      <  1:804499‑M2/90‑6 (MQ=255)
    gATTGAAGGCAGTATTGATAACGATCGCGTTACGCTGACTAACCTGGGTGCCGACATCGCCCGCGGGACATTAACCGGAAACGCGcagcg                                                                                      <  2:422923‑M2/90‑6 (MQ=255)
      ttGAAGGCAGTATTGATAACGATCGCGTTACGCTGACTAACCTGGGTGCCGACATCGCCCGCGGGACATTAACCGGAAACGCGcagcgta                                                                                    <  2:753704‑M2/90‑8 (MQ=255)
      ttGAAGGCAGTATTGATAACGATCGCGTTACGCTGACTAACCTGGGTGCCGACATCGCCCGCGGGACATTAACCGGAAACGCGcagcgta                                                                                    >  1:807033‑M2/1‑83 (MQ=255)
         aaGGCAGTATTGATAACGATCGCGTTACGCTGACTAACCTGGGTGCCGACATCGCCCGCGGGACATTAACCGGAAACGCGcagcgtaacg                                                                                 <  2:216097‑M2/90‑11 (MQ=255)
         aaGGCAGTATTGATAACGATCGCGTTACGCTGACTAACCTGGGTGCCGACATCGCCCGCGGGACATTAACCGGAAACGCGcagcgtaacg                                                                                 <  2:611961‑M2/90‑11 (MQ=255)
              aGTATTGATAACGATCGCGTTACGCTGACTAACCTGGGTGCCGACATCGCCCGCGGGACATTAACCGGAAACGCGCAGCGTAACGCCGAc                                                                            >  2:10264/1‑90 (MQ=255)
                    gATAACGATCGCGTTACGCTGACTAACCTGGGTGCCGACATCGCCCGCGGGACATTAACCGGAAACGCGCAGCGTAACGCCGACGGCAGc                                                                      >  2:144396/1‑90 (MQ=255)
                      tAACGATCGCGTTACGCTGACTAACCTGGGTGCCGACATCGCCCGCGGGACATTAACCGGAAACGCGCAGCGTAACGCCGACGGCAGCTg                                                                    >  2:720290/1‑90 (MQ=255)
                       aaCGATCGCGTTACGCTGACTAACCTGGGTGCCGACATCGCCCGCGGGACATTAACCGGAAACGCGCAGCGTAACGCCGACGGCAGCTgg                                                                   >  1:71320/1‑90 (MQ=255)
                         cGATCGCGTTACGCTGACTAACCTGGGTGCCGACATCGCCCGCGGGACATTAACCGGAAACGCGCAGCGTAACGCCGACGGCAGCTGGCa                                                                 >  1:28374/1‑90 (MQ=255)
                         cGATCGCGTTACGCTGACTAACCTGGGTGCCGACATCGCCCGCGGGACATTAACCGGAAACGCGCAGCGTAACGCCGACGGCAGCTGGCa                                                                 >  2:718866/1‑90 (MQ=255)
                           aTCGCGTTACGCTGACTAACCTGGGTGCCGACATCGCCCGCGGGACATTAACCGGAAACGCGCAGCGTAACGCCGACGGCAGCTGGCAAg                                                               >  2:626047/1‑90 (MQ=255)
                                    cGCTGACTAACCTGGGTGCCGACATCGCCCGCGGGACATTAACCGGAAACGCGCAGCGTAACGCCGACGGCAGCTGGCAAGTGGAAAAcc                                                      <  2:721269/90‑1 (MQ=255)
                                        gACTAACCTGGGTGCCGACATCGCCCGCGGGACATTAACCGGAAACGCGCAGCGTAACGCCGACGGCAGCTGGCaa                                                                >  1:190547/1‑76 (MQ=255)
                                        gACTAACCTGGGTGCCGACATCGCCCGCGGGACATTAACCGGAAACGCGCAGCGTAACGCCGACGGCAGCTGGCaa                                                                <  2:190547/76‑1 (MQ=255)
                                        gACTAACCTGGGTGCCGACATCGCCCGCGGGACATTAACCGGAAACGCGCAGCGTAACGCCGACGGCAGCTGGCAAGTGGAAAACCTgcg                                                  >  2:260758/1‑90 (MQ=255)
                                           tAACCTGGGTGCCGACATCGCCCGCGGGACATTAACCGGAAACGCGCAGCGTAACGCCGACGGCAGCTGGCAAGTGGAAAACCTGCGCAt                                               <  1:753704/90‑1 (MQ=255)
                                             aCCTGGGTGCCGACATCGCCCGCGGGACATTAACCGGAAACGCGCAGCGTAACGCCGACGGCAGCTGGCaa                                                                >  1:815863/1‑71 (MQ=255)
                                             aCCTGGGTGCCGACATCGCCCGCGGGACATTAACCGGAAACGCGCAGCGTAACGCCGACGGCAGCTGGCaa                                                                <  2:815863/71‑1 (MQ=255)
                                             aCCTGGGTGCCGACATCGCCCGCGGGACATTAACCGGAAACGCGCAGCGTAACGCCGACGGCAGCTGGCAAGTGGAAAACCTGCGCATgg                                             <  1:720290/90‑1 (MQ=255)
                                              ccTGGGTGCCGACATCGCCCGCGGGACATTAACCGGAAACGCGCAGCGTAACGCCGACGGCAGCTGGCAAGTGGAAAACCTGCGCATGGc                                            <  1:51769/90‑1 (MQ=255)
                                                tGGGTGCCGACATCGCCCGCGGGACATTAACCGGAAACGCGCAGCGTAACGCCGACGGCAGCTGGCAAGTGGAAAACCTGCGCATGGCgg                                          >  1:568614/1‑90 (MQ=255)
                                                  ggTGCCGACATCGCCCGCGGGACATTAACCGGAAACGCGcagcgtaacgcgatcgttatcaatactgccttcaatcagtacattggtggc                                        <  2:564497‑M2/90‑52 (MQ=255)
                                                  ggTGCCGACATCGCCCGCGGGACATTAACCGGAAACGCGCAGCGTAACGCCGACGGCAGCTGGCAAGTGGAAAACCTGCGCATGGCGGat                                        >  1:783463/1‑90 (MQ=255)
                                                  ggTGCCGACATCGCCCGCGGGACATTAACCGGAAACGCGCAGCGTAACGCCGACGGCAGCTGGCAAGTGGAAAACCTGCGCATGGCGGat                                        >  2:439981/1‑90 (MQ=255)
                                                  ggTGCCGACATCGCCCGCGGGACATTAACCGGAAACGCGCAGCGTAACGCCGACGGCAGCTGGCAAGTGGAAAACCTGCGCATGGCGGat                                        >  1:46164/1‑90 (MQ=255)
                                                   gTGCCGACATCGCCCGCGGGACATTAACCGGAAACGCGcagcgtaacgcgatcgttatcaatactgccttcaatcagtacattggtggca                                       >  1:292515‑M2/1‑38 (MQ=255)
                                                    tGCCGACATCGCCCGCGGGACATTAACCGGAAACGCGCAGCGTAACGCCGACGGCAGCTGGCAAGTGGAAAACCTGCGCATGGCGGatat                                      <  2:56832/90‑1 (MQ=255)
                                                    tGCCGACATCGCCCGCGGGACATTAACCGGAAACGCGCAGCGTAACGCCGACGGCAGCTGGCAAGTGGAAAACCTGCGCATGGCGGatat                                      >  1:15531/1‑90 (MQ=255)
                                                     gCCGACATCGCCCGCGGGACATTAACCGGAAACGCGCAGCGTAACGCCGACGGCAGCTGGCAAGTGGAAAACCTGCGCATGGCGGAtata                                     <  2:196642/90‑1 (MQ=255)
                                                      ccGACATCGCCCGCGGGACATTAACCGGAAACGCGc                                                                                          <  1:833988‑M2/36‑2 (MQ=255)
                                                      ccGACATCGCCCGCGGGACATTAACCGGAAACGCGc                                                                                          >  2:833988‑M2/1‑35 (MQ=255)
                                                      ccGACATCGCCCGCGGGACATTAACCGGAAACGCGCAGCGTAACGCCGACGGCAGCTGGCAAGTGGAAAACCTGCGCATGGCGGATATAc                                    >  1:237564/1‑90 (MQ=255)
                                                      ccGACATCGCCCGCGGGACATTAACCGGAAACGCGCAGCGTAACGCCGACGGCAGCTGGCAAGTGGAAAACCTGCGCATGGCGGATATAc                                    >  1:234352/1‑90 (MQ=255)
                                                        gACATCGCCCGCGGGACATTAACCGGAAACGCGCAGCGTAACGCCGACGGCAGCTGGCAAGTGGAAAACCTGCGCATGGCGGATATACGt                                  >  1:690049/1‑90 (MQ=255)
                                                             cGCCCGCGGGACATTAACCGGAAACGCGCAGCGTAACGCCGACGGCAGCTGGCAAGTGGAAAACCTGCGCAt                                               >  2:274300/1‑72 (MQ=255)
                                                             cGCCCGCGGGACATTAACCGGAAACGCGCAGCGTAACGCCGACGGCAGCTGGCAAGTGGAAAACCTGCGCAt                                               <  1:274300/72‑1 (MQ=255)
                                                                      gACATTAACCGGAAACGCGCAGCGTAACGCCGACGGCAGCTGGCAAGTGGAAAACCTGCGCATGGCGGATATACGTCTACAAAGCGaaaa                    <  1:144396/90‑1 (MQ=255)
                                                                         aTTAACCGGAAACGCGCAGCGTAACGCCGACGGCAGCTGGCAAGTGGAAAACCTGCGCATGGCGGATATACGTCTACAAAGCGAAAAATc                 >  1:370371/1‑90 (MQ=255)
                                                                               cGGAAACGCGCAGCGTAACGCCGACGGCAGCTGGCAAGTGGAAAACCTGCGCATGGCGGATATACGTCTACAAAGCGAAAAATCGCTAAc           >  1:388322/1‑90 (MQ=255)
                                                                                   aaCGCGcagcgtaacgcgatcgttatcaatactgccttcaatcagtacattggtggcaggaacatcattgagcgaaagcgatccggcact       >  1:564497‑M2/1‑6 (MQ=255)
                                                                                   aaCGCGCAGCGTAACGCCGACGGCAGCTGGCAAGTGGAAAACCTGCGCATGGCGGATATACGTCTACAAAGCGAAAAATCGCTAACCGAc       <  1:636423/90‑1 (MQ=255)
                                                                                        gcAGCGTAACGCCGACGGCAGCTGGCAAGTGGAAAACCTGCGCATGGCGGATATACGTCTACAAAGCGAAAAATCGCTAACCGActtctt  <  2:237564/90‑1 (MQ=255)
                                                                                        |                                                                                         
TACTGATTGAAGGCAGTATTGATAACGATCGCGTTACGCTGACTAACCTGGGTGCCGACATCGCCCGCGGGACATTAACCGGAAACGCGCAGCGTAACGCCGACGGCAGCTGGCAAGTGGAAAACCTGCGCATGGCGGATATACGTCTACAAAGCGAAAAATCGCTAACCGACTTCTT  >  NZ_CP009273/3824364‑3824541

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