New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 228859	NA (NA)	9 (0.310) +AG	5/94	2.5	27.3%	noncoding (104/120 nt)	rrfH	5S ribosomal RNA
	?	NC_000913	= 1363199	25 (0.820)					coding (457/1488 nt)	puuC	gamma‑glutamyl‑gamma‑aminobutyraldehyde dehydrogenase

AAAGTGTATGGCGAAGTGGCGACCACCAGTAGCCATGAGCTGGCGATGATCGTGCGTGAACCGGTCGGCGTGATTGCCGCCATCGTGCCGTGGAACT  >  NC_000913/1363151‑1363247                                                 |                                                 aaaGTGTATGGCGAAGTGGCGACCACCAGTAGCCATGAGCTGGCgatgat                                                 >  1:1063783/1‑50 (MQ=255)   aGTGTATGGCGAAGTGGCGACCACCAGTAGCCATGAGCTGGCGATGATcg                                               >  1:895144/1‑50 (MQ=255)      gtATGGCGAAGTGGCGACCACCAGTAGCCATGAGCTGGCGATGATcgtgc                                            >  1:2246197/1‑50 (MQ=255)        aTGGCGAAGTGGCGACCACCAGTAGCCATGAGCTGGCGATGATcgtgcgt                                          >  1:1274469/1‑50 (MQ=255)             gAAGTGGCGACCACCAGTAGCCATGAGCTGGCGATGATCGTGCGTGAAcc                                     <  1:508014/50‑1 (MQ=255)                 tGGCGACCACCAGTAGCCATGAGCTGGCGATGATCGTGCGTGAACCGGTc                                 >  1:1153991/1‑50 (MQ=255)                     gACCACCAGTAGCCATGAGCTGGCGATGATCGTGCGTGAACCGGTCGGCg                             >  1:1293110/1‑50 (MQ=255)                        caccaGTAGCCATGAGCTGGCGATGATCGTGCGTGAACCGGTCGGCGTGa                          >  1:1825369/1‑50 (MQ=255)                        caccaGTAGCCATGAGCTGGCGATGATCGTGCGTGAACCGGTCGGCGTGa                          >  1:173760/1‑50 (MQ=255)                         accaGTAGCCATGAGCTGGCGATGATCGTGCGTGAACCGGTCGGCGTGAt                         >  1:984395/1‑50 (MQ=255)                         accaGTAGCCATGAGCTGGCGATGATCGTGCGTGAACCGGTCGGCGTGAt                         >  1:1420194/1‑50 (MQ=255)                              tAGCCATGAGCTGGCGATGATCGTGCGTGAACCGGTCGGCGTGATTgccg                    <  1:1681212/50‑1 (MQ=255)                               aGCCATGAGCTGGCGATGATCGTGCGTGAACCGGTCGGCGTGATTgccgc                   >  1:1910505/1‑50 (MQ=255)                                gCCATGAGCTGGCGATGATCGTGCGTGAACCGGTCGGCGTGATTgccgc                   >  1:145592/1‑49 (MQ=255)                                gCCATGAGCTGGCGATGATCGTGCGTGAACCGGTCGGCGTGATTgccgc                   >  1:233805/1‑49 (MQ=255)                                gCCATGAGCTGGCGATGATCGTGCGTGAACCGGTCGGCGTGATTgccgc                   >  1:2648114/1‑49 (MQ=255)                                   aTGAGCTGGCGATGActtactctcacatgggtaaccccacactaccatca               >  1:951916‑M2/1‑15 (MQ=255)                                   aTGAGCTGGCGATGActtactctcacatgggtaaccccacactaccatc                >  1:2443704‑M2/1‑15 (MQ=255)                                   aTGAGCTGGCGATGActtactctcacatgggtaaccccacactaccat                 <  1:536444‑M2/48‑34 (MQ=255)                                   aTGAGCTGGCGATGATCGTGCGTGAACCGGTCGGCGTGATTGCCGCCATc               >  1:1780962/1‑50 (MQ=255)                                    tGAGCTGGCGATGActtactctcacatgggtaaccccacactaccatcag              >  1:636206‑M2/1‑14 (MQ=255)                                    tGAGCTGGCGATGActtactctcacatgggtaaccccacactaccatcag              >  1:692370‑M2/1‑14 (MQ=255)                                     gAGCTGGCGATGATCGTGCGTGAACCGGTCGGCGTGATTGCCGCCATCGt             >  1:760267/1‑50 (MQ=255)                                      aGCTGGCGATGATCGTGCGTGAACCGGTCGGCGTGATTGCCGCCATCGTg            <  1:2699432/50‑1 (MQ=255)                                      aGCTGGCGATGATCGTGCGTGAACCGGTCGGCGTGATTGCCGCCATCGTg            <  1:604261/50‑1 (MQ=255)                                        cTGGCGATGActtactctcacatgggtaaccccacactaccatcagcgct          >  1:226475‑M2/1‑10 (MQ=255)                                        cTGGCGATGActtactctcacatgggtaaccccacactaccatcagcgct          >  1:98314‑M2/1‑10 (MQ=255)                                         tGGCGATGActtactctcacatgggtaaccccacactaccatcagcgct          <  1:66229‑M2/49‑41 (MQ=255)                                         tGGCGATGATCGTGCGTGAACCGGTCGGCGTGATTGCCGCCATCGTGCCg         >  1:1410781/1‑50 (MQ=255)                                         tGGCGATGATCGTGCGTGAACCGGTCGGCGTGATTGCCGCCATCGTGCCg         >  1:1260758/1‑50 (MQ=255)                                           gCGATGActtactctcacatgggtaaccccacactaccatcagcgct          <  1:2640942‑M2/47‑41 (MQ=255)                                            cGATGATCGTGCGTGAACCGGTCGGCGTGATTGCCGCCATCGTGCCGTg       >  1:2614208/1‑49 (MQ=255)                                              atgatCGTGCGTGAACCGGTCGGCGTGATTGCCGCCATCGTGCCGTGGaa    >  1:2251635/1‑50 (MQ=255)                                                gatCGTGCGTGAACCGGTCGGCGTGATTGCCGCCATCGTGCCGTGGAACt  >  1:907956/1‑50 (MQ=255)                                                 |                                                 AAAGTGTATGGCGAAGTGGCGACCACCAGTAGCCATGAGCTGGCGATGATCGTGCGTGAACCGGTCGGCGTGATTGCCGCCATCGTGCCGTGGAACT  >  NC_000913/1363151‑1363247

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