New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	1794268 =	42 (1.030)	5 (0.130)	3/250	5.2	11.4%	coding (885/981 nt)	btuC	vitamin B12 ABC transporter permease
	?	NC_000913	1794318 =	38 (0.990)					coding (835/981 nt)	btuC	vitamin B12 ABC transporter permease
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

TGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGCTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGACCACGCCAATAGGCAGCTCTGCGGCAGCTAATGCCAGGCGCGCTACAATATCGGCCAGCAGCAATGCGCTCGCCCCTGCCAGCGCGCAGCCGGGAAGTAATACGCGATGATCGGTTAAACCACACAACCGGAGAATATGGGGGATCACCAGACCAATAAAGCCGATAG  >  NC_000913/1794125‑1794398                                                                                                                                                |                                                                                                                                   tgGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGCTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGACCACGCCAATAGGCAGCTCTGCGGCAGCTAATGCCAGgcgcgc                                                                                                                               <  1:189484/149‑1 (MQ=255)      aGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGCTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGACCACGCCAATAGGCAGCTCTGCGGCAGCTAATGCCAGGCGCGCTACaa                                                                                                                          >  1:74613/1‑149 (MQ=255)       gTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGCTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGACCACGCCAATAGGCAGCTCTGCGGCAGCTAATGCCAGGCGCGCTACAat                                                                                                                         <  3:111806/149‑1 (MQ=255)                   gCTTTTAATCATAGACCGTCTTTTTGCGGCTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGACCACGCCAATAGGCAGCTCTGCGGCAGCTAATGCCAGGCGCGCTACAATATCGGCcagcag                                                                                                             >  7:28879/1‑149 (MQ=255)                   gCTTTTAATCATAGACCGTCTTTTTGCGGCTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGACCACGCCAATAGGCAGCTCTGCGGCAGCTAATGCCAGGCGCGCTACAATATCGGCcagcag                                                                                                             <  2:74613/149‑1 (MQ=255)                   gCTTTTAATCATAGACCGTCTTTTTGCGGCTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGACCACGCCAATAGGCAGCTCTGCGGCAGCTAATGCCAGGCGCGCTACAATATCGGCcagca                                                                                                              <  2:11232/148‑1 (MQ=255)                        tAATCATAGACCGTCTTTTTGCGGCTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGACCACGCCAATAGGCAGCTCTGCGGCAGCTAATGCCAGGCGCGCTACAATATCGGCCAGCAGCAAt                                                                                                         <  5:148672/148‑1 (MQ=255)                                         tttGCGGCTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGACCACGCCAATAGGCAGCTCTGCGGCAGCTAATGCCAGGCGCGCTACAATATCGGCCAGCAGCAATGCGCTCGCCCCTGCCAgc                                                                                       >  2:44468/1‑149 (MQ=255)                                                  aaCGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGACCACGCCAATAGGCAGCTCTGCGGCAGCTAATGCCAGGCGCGCTACAATATCGGCCAGCAGCAATGCGCTCGCCCCTGCCAGCGCGCAGCCg                                                                              >  8:31869/1‑149 (MQ=255)                                                        cTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGACCACGCCAATAGGCAGCTCTGCGGCAGCTAATGCCAGGCGCGCTACAATATCGGCCAGCAGCAATGCGCTCGCCCCTGCCAGCGCGCAGCCg                                                                              <  1:59334/143‑1 (MQ=255)                                                        cTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGACCACGCCAATAGGCAGCTCTGCGGCAGCTAATGCCAGGCGCGCTACAATATCGGCCAGCAGCAATGCGCTCGCCCCTGCCAGCGCGCAGCCg                                                                              >  2:59334/1‑143 (MQ=255)                                                          gCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGACCACGCCAATAGGCAGCTCTGCGGCAGCTAATGCCAGGCGCGCTACAATATCGGCCAGCAGCAATGCGCTCGCCCCTGCCAGCGCGCAGCCGGGAAGTaa                                                                      >  2:60957/1‑149 (MQ=255)                                                            ttttAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGACCACGCCAATAGGCAGCTCTGCGGCAGCTAATGCCAGGCGCGCTa                                                                                                                             >  8:69020/1‑92 (MQ=255)                                                            ttttAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGACCACGCCAATAGGCAGCTCTGCGGCAGCTAATGCCAGGCGCGCTa                                                                                                                             <  7:69020/92‑1 (MQ=255)                                                                 aCAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGACCACGCCAATAGGCAGCTCTGCGGCAGCTAATGCCAGGCGCGCTACAATATCGGCCAGCAGCAATGCGCTCGCCCCTGCCAGCGCGCAGCCGGGAAGTAATAcgcg                                                                <  1:60957/148‑1 (MQ=255)                                                                 aCAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGACCACGCCAATAGGCAGCTCTGCGGCAGCTAATGCCAGGCGCGCTACAATATCGGCCAGCAGCAATGCGCTCGCCCCTGCCAGCGCGCAGCCGGGAAGTAATACGCga                                                               <  5:106705/149‑1 (MQ=255)                                                                 aCAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGACCACGCCAATAGGCAGCTCTGCGGCAGCTAATGCCAGGCGCGCTACAATATCGGCCAGCAGCAATGCGCTCGCCCCTGCCAGCGCGCAGCCGGGAAGTAATACGCga                                                               <  5:104239/149‑1 (MQ=255)                                                                        ccAGATAAACACCGGCGCACCTAACGTTGCGGTGACCACGCCAATAGGCAGCTCTGCGGCAGCTAATGCCAGGCGCGCTACAATATCGGCCAGCAGCAATGCGCTCGCCCCTGCCAGCGCGCAGCCGGGAAGTAATACGCGATGATCgg                                                        >  3:39417/1‑149 (MQ=255)                                                                          aGATAAACACCGGCGCACCTAACGTTGCGGTGACCACGCCAATAGGCAGCTCTGCGGCAGCTAATGCCAGGCGCGCTACAATATCGGCCAGCAGCAATGCGCTCGCCCCTGCCAGCGCGCAGCCGGGAAGTAATACGCGATGATCGGtt                                                      <  7:31869/149‑1 (MQ=255)                                                                           gATAAACACCGGCGCACCTAACGTTGCGGTGACCACGCCAATAGGCAGCTCTGCGGCAGCTAATGCCAGGCGCGCTACAATATCGGCCAGCAGCAAt                                                                                                         >  3:143135/1‑97 (MQ=255)                                                                           gATAAACACCGGCGCACCTAACGTTGCGGTGACCACGCCAATAGGCAGCTCTGCGGCAGCTAATGCCAGGCGCGCTACAATATCGGCCAGCAGCAAt                                                                                                         <  4:143135/97‑1 (MQ=255)                                                                                 cacCGGCGCACCTAACGTTGCGGTGACCACGCCAATAGGCAGCTCTGCGGCAGCTAATGCCAGGCGCGCTACAATATCGGCCAGCAGCAATGCGCTCGCCCCTGCCAGCGCGCAGCCGGGAAGTAATACGCGATGATCGGTTAAACcac                                               <  2:118461/149‑1 (MQ=255)                                                                                  acCGGCGCACCTAACGTTGCGGTGACCACGCCAATAGGCAGCTCTGCGGCAGCTAATGCCAGGCGCGCTACAATATCGGCCAGCAGCAATGCGCTCGCCCCTGCCAGCGCGCAGCCGGGAAGTAATACGCgatga                                                            <  2:88645/135‑1 (MQ=255)                                                                                  acCGGCGCACCTAACGTTGCGGTGACCACGCCAATAGGCAGCTCTGCGGCAGCTAATGCCAGGCGCGCTACAATATCGGCCAGCAGCAATGCGCTCGCCCCTGCCAGCGCGCAGCCGGGAAGTAATACGCgatga                                                            >  1:88645/1‑135 (MQ=255)                                                                                   ccGGCGCACCTAACGTTGCGGTGACCACGCCAATAGGCAGCTCTGCGGCAGCTAATGCCAGGCGCGCTACAATATCGGCCAGCAGCAATGCGCTCGCCCCTGCCAGCGCGCAGCCGGGAAGTAATACGCGATGATCGGTTAAACcaca                                              >  5:158253/1‑148 (MQ=255)                                                                                       cgcACCTAACGTTGCGGTGACCACGCCAATAGGCAGCTCTGCGGCAGCTAATGCCAGGCGCGCTACAATATCGGCCAGCAGCAATGCGCTCGCCCCTGCCAGCGCGCAGCCGGGAAGTAATACGCGATGATCGGTTAAACCACACAAcc                                         <  8:28879/149‑1 (MQ=255)                                                                                            cTAACGTTGCGGTGACCACGCCAATAGGCAGCTCTGCGGCAGCTAATGCCAGGCGCGCTACAATATCGGCCAGCAGCAATGCGCTCGCCCCTGCCAGCGCGCAGCCGGGAAGTAATACGCGATGATCGGTTAAACCACACAACCGgag                                     >  2:5720/1‑148 (MQ=255)                                                                                                       gTGACCACGCCAATAGGCAGCTCTGCGGCAGCTAATGCCAGGCGCGCTACAATATCGGCCAGCAGCAATGCGCTCGCCCCTGCCAGCGCGCAGCCGGGAAGTAATACGCGATGATCGGTTAAACCACACAACCGGAGAATATGGGGGAt                         >  3:133918/1‑149 (MQ=255)                                                                                                          aCCACGCCAATAGGCAGCTCTGCGGCAGCTAATGCCAGGCGCGCTACAATATCGGCCAGCAGCAATGCGCTCGCCCCTGCCAGCGCGCAGCCGGGAAGTAATACGCGATGATCGGTTAAACCACACAACCGGAGAATATGGGGGATCAc                      <  1:44468/149‑1 (MQ=255)                                                                                                             aCGCCAATAGGCAGCTCTGCGGCAGCTAATGCCAGGCGCGCTACAATATCGGCCAGCAGCAATGCGCTCGCCCCTGCCAGCGCGCAGCCGGGAAGTAATACGCGATGATCGGTTAAACCACACAACCGGAGAata                                 <  5:149315/135‑1 (MQ=255)                                                                                                             aCGCCAATAGGCAGCTCTGCGGCAGCTAATGCCAGGCGCGCTACAATATCGGCCAGCAGCAATGCGCTCGCCCCTGCCAGCGCGCAGCCGGGAAGTAATACGCGATGATCGGTTAAACCACACAACCGGAGAata                                 >  6:149315/1‑135 (MQ=255)                                                                                                               gCCAATAGGCAGCTCTGCGGCAGCTAATGCCAGGCGCGCTACAATATCGGCCAGCAGCAATGCGCTCGCCCCTGCCAGCGCGCAGCCGGGAAGTAATACGCGATGATCGGTTAAACCACACAACCGGAGAATATGGGGGATCACCAGAc                 <  4:39417/149‑1 (MQ=255)                                                                                                                 cAATAGGCAGCTCTGCGGCAGCTAATGCCAGGCGCGCTACAATATCGGCCAGCAGCAATGCGCTCGCCCCTGCCAGCGCGCAGCCGGGAAGTAATACGCGATGATCGGTTAAACCACACAACCGgag                                     >  8:91194/1‑127 (MQ=255)                                                                                                                 cAATAGGCAGCTCTGCGGCAGCTAATGCCAGGCGCGCTACAATATCGGCCAGCAGCAATGCGCTCGCCCCTGCCAGCGCGCAGCCGGGAAGTAATACGCGATGATCGGTTAAACCACACAACCGgag                                     <  7:91194/127‑1 (MQ=255)                                                                                                                       gCAGCTCTGCGGCAGCTAATGCCAGGCGCGCTACAATATCGGCCAGCAGCAATGCGCTCGCCCCTGCCAGCGCGCAGCCGGGAAGTAATACGCGATGATCGGTTAAACCACACAACCGGAGAATATGGGGGATCACCAGACCAATAAAg         >  6:95668/1‑149 (MQ=255)                                                                                                                            tctGCGGCAGCTAATGCCAGGCGCGCTACAATATCGGCCAGCAGCAATGCGCTCGCCCCTGCCAGCGCGCAGCCGGGAAGTAATACGCGATGATCGGTTAAACCACACAACCGGAGAATATGGGGGATCACCAGACCAATAAAGCCGa     >  5:171848/1‑148 (MQ=255)                                                                                                                            tctGCGGCAGCTAATGCCAGGCGCGCTACAATATCGGCCAGCAGCAATGCGCTCGCCCCTGCCAGCGCGCAGCCGGGAAGTAATACGCGATGATCGGTTAAACCACACAACCGGAGAATATGGGGGATCACCAGACCAATAAAGCCGAt    <  6:158253/149‑1 (MQ=255)                                                                                                                              tGCGGCAGCTAATGCCAGGCGCGCTACAATATCGGCCAGCAGCAATGCGCTCGCCCCTGCCAGCGCGCAGCCGGGAAGTAATACGCGATGATCGGTTAAACCACACAACCGGAGAATATGGGGGATCACCAGACCAATAAAGCCGATAg  >  5:174226/1‑149 (MQ=255)                                                                                                                                      cTAATGCCAGGCGCGCTACAATATCGGCCAGCAGCAATGCGCTCGCCCCTGCCAGCGCGCAGCCGGGAAGTAATACGCGATGATCGGTTAAACCACACAACCGGAGAATATGGGGGATCACCAGACCAATAAAGCCGa     >  2:109992/1‑138 (MQ=255)                                                                                                                                      cTAATGCCAGGCGCGCTACAATATCGGCCAGCAGCAATGCGCTCGCCCCTGCCAGCGCGCAGCCGGGAAGTAATACGCGATGATCGGTTAAACCACACAACCGGAGAATATGGGGGATCACCAGACCAATAAAGCCGa     <  1:109992/138‑1 (MQ=255)                                                                                                                                        aaTGCCAGGCGCGCTACAATATCGGCCAGCAGCAATGCGCTCGCCCCTGCCAGCGCGCAGCCGGGAAGTAATACGCGATGATCGGTTAAACCACACAACCGGAGAATATGGGGGATCAc                      <  2:39913/119‑1 (MQ=255)                                                                                                                                        aaTGCCAGGCGCGCTACAATATCGGCCAGCAGCAATGCGCTCGCCCCTGCCAGCGCGCAGCCGGGAAGTAATACGCGATGATCGGTTAAACCACACAACCGGAGAATATGGGGGATCAc                      >  1:39913/1‑119 (MQ=255)                                                                                                                                                |                                                                                                                                   TGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGCTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGACCACGCCAATAGGCAGCTCTGCGGCAGCTAATGCCAGGCGCGCTACAATATCGGCCAGCAGCAATGCGCTCGCCCCTGCCAGCGCGCAGCCGGGAAGTAATACGCGATGATCGGTTAAACCACACAACCGGAGAATATGGGGGATCACCAGACCAATAAAGCCGATAG  >  NC_000913/1794125‑1794398

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.