New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	4873089 =	45 (1.520)	4 (0.140)	4/236	NT	9.0%	coding (576/693 nt)	aqpZ	aquaporin Z
	?	NC_002947	4873108 =	38 (1.340)					coding (557/693 nt)	aqpZ	aquaporin Z

GTGCACCAGCGTTAAGGTTCTTTCGCCAGGCCCCGATACAGCGCGCCGCCGATAGCAGCACCGATCAGCGGCGCCACCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCGGCCCCGTGC‑TACGCGCGGGGTTGACAGAGGTATTGGTCACCGGGATCGAGATCAGGTGGATCAGGGTCAGGGCCAGGCCGATGGCGATGGGTGCAAAACCTGCCGGTGCGCGGGCGTCGGTCGCCCCCATGATGACCACCAGGAACATCGCCGTCATCACCACTTCACTG  >  NC_002947/4872961‑4873257                                                                                                                                                     |                                                                                                                                                      gTGCACCAGCGTTAAGGGTCGTGCTCCAGGCCCCGATACATCGCGCCGCCGTGAGCAGCACCGATCAGCGGCGCCACCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCGGCCCCGTGC‑TACGCGCGGGGTTg                                                                                                                                                     >  2:628580/1‑150 (MQ=255)      ccAGCGTTAAGGTTCTTTCGCCAGGCCCCGATACAGCGCGCCGCCGATAGCAGCACCGATCAGCGGCGCCACCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCGGCCCCGTGC‑TACGCGCGGGGTTGACaga                                                                                                                                                <  2:518753/150‑1 (MQ=255)      ccAGCGTTAAGGTTCTTTCGCCAGGCCCCGATACAGCGCGCCGCCGATAGCAGCACCGATCAGCGGCGCCACCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCGGCCCCGTGC‑TACGCGCGGGGTTGACaga                                                                                                                                                >  2:181003/1‑150 (MQ=255)       cAGCGTTAAGGTTCTTTCGCCAGGCCCCGATACAGCGCGCCGCCGATAGCAGCACCGATCAGCGGCGCCACCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCGGCCCCGTGC‑TACGCGCGGGGTTGACagag                                                                                                                                               <  1:381150/150‑1 (MQ=255)            ttAATGTTCTTGCGCCAGGCCCCGATACAGCGCGCCGCCGATAGCAGCACCGATCAGCGGCGCCACCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCGGCCCCGTGC‑TACGCGCGGGGTTGACAGAGGTAtt                                                                                                                                          <  1:39628/150‑1 (MQ=255)                 gTTCTTTCGCCAGGCCCCGATACAGCGCGCCGCCGATAGCAGCACCGATCAGCGGCGCCACCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCGGCCCCGTGC‑TACGCGCGGGGTTGACAGCGGTATTGGTCa                                                                                                                                     >  2:697833/1‑150 (MQ=255)                     tttCGCCAGGCCCCGATACAGCGCGCCGCCGATAGCAGCACCGATCAGCGGCGCCACCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCGGCCCCGTGC‑TACGCGCGGGGTTGACAGAGGTATTgg                                                                                                                                        <  2:227538/143‑1 (MQ=255)                     tttCGCCAGGCCCCGATACAGCGCGCCGCCGATAGCAGCACCGATCAGCGGCGCCACCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCGGCCCCGTGC‑TACGCGCGGGGTTGACAGAGGTATTgg                                                                                                                                        >  1:227538/1‑143 (MQ=255)                          ccAGGCCCCGATACAGCGCGCCGCCGATAGCAGCACCGATCAGCGGCGCCACCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCGGCCCCGTGC‑TACGCGCGGGGTTGACAGAGGTATTGGTCACCGGGATCg                                                                                                                            >  1:340694/1‑150 (MQ=255)                           cAGGCCCCGATACAGCGCGCCGCCGATAGCAGCACCGATCAGCGGCGCCACCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCGGCCCCGTGC‑TACGCGCGGGGTTGACAGAGGTATTGGTCACCGGGATCga                                                                                                                           <  2:510914/150‑1 (MQ=255)                                   gATACAGCGCGCCGCCGATAGCAGCACCGATCAGCGGCGCCACCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCGGCCCCGTGC‑TACGCGCGGGGTTGACAGAGGTATTGGTCACCGGGATCGAGATCAGGt                                                                                                                   >  2:773431/1‑150 (MQ=255)                                   gATACAGAGCGCCGCCGATAGCAGCACCGATCAGCGGCGCCACCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCGGCCCCGTGC‑TACGCGCGGGGTTGACAGAGGTATTGGTCACCGGGATCGAGATCAGGt                                                                                                                   >  2:131927/1‑150 (MQ=255)                                      aCAGCGCGCCGCCGATAGCAGCACCGATCAGCGGCGCCACCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCGGCCCCGTGC‑TACGCGCGGGGTTGACAGAGGTATTGGTCACCGGGATCGAGATCAGGTGGa                                                                                                                <  1:773431/150‑1 (MQ=255)                                        aGCGCGCCGCCGATAGCAGCACCGATCAGCGGCGCCACCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCGGCCCCGTGC‑TACGCGCGGGGTTGACAGAGGTATTGGTCACCGGGATCGAGATCAGGTGGAtc                                                                                                              <  2:340694/150‑1 (MQ=255)                                          cgcgcCGCCGATAGCAGCACCGATCAGCGGCGCCACCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCGGCCCCGTGC‑TACGCGCGGGGTTGACAGAGGTATTGGTCACCGGGATCGAGATCAGGTGGAtcag                                                                                                            <  1:375003/150‑1 (MQ=255)                                                gccgATAGCAGCACCGATCAGCGGCGCCACCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCGGCCCCGTGC‑TACGCGCGGGGTTGACAGAGGTATTGGTCACCGGGATCGAGATCAGGTGGAtcagggtcag                                                                                                      >  2:783110/1‑150 (MQ=255)                                                 ccgATAGCAGCACCGATCAGCGGCGCCACCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCGGCCCCGTGC‑TACGCGCGGGGTTGACAGAGGTATTGGTCACCGGGATCGAGATCAGGTGGAtcagggtcagg                                                                                                     <  1:521861/150‑1 (MQ=255)                                                      agccgcACCGATCCGCGGCGCCACACAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCGGCCCCGTGC‑TACGCGCGGGGTTGACAGAGGTATTGGTCACCGGGATCGAGATCAGGTGGATCAGGGTCAGGGCCAg                                                                                                <  2:770066/150‑1 (MQ=255)                                                        cagcaCCGATCAGCGGCGCCACCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCGGCCCCGTGC‑TACGCGCGGGGTTGACAGAGGTATTGGTCACCGGGATCGAGATCAGGTGGATCAGGGTCAGGGCCAGGc                                                                                              <  1:783110/150‑1 (MQ=255)                                                             ccGATCAGCGGCGCCACCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCAACCGACGAACAGGGCCGGCCCCGTGC‑TACGAGCGGGGTTGACAGCGGTATTGGTCACCGGGATCGAGATCAGGTGGATCAGGGTCAGGGCCAGGCcgatg                                                                                         <  1:181003/150‑1 (MQ=255)                                                                 gccatcggcctggccctgaccctgatccacctgatctcgatcccggtgaccaatacctctgtcaaccccgtgctacgcgcggggTTGACAGAGGTATTgg                                                                                                                                        <  2:140362‑M2/16‑1 (MQ=255)                                                                 gccatcggcctggccctgaccctgatccacctgatctcgatcccggtgaccaatacctctgtcaaccccgtgctacgcgcggggTTGACAGAGGTATTgg                                                                                                                                        >  1:140362‑M2/85‑100 (MQ=255)                                                                      ggcgCCACCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCGGCCCCGTGC‑TACGCGCGGGGTTGACAGAGGTATTGGTCACCGGGATCGAGATCAGGTGGATCAGGGTCAGGGCCAGGCCGATGGCGATGGGt                                                                                <  1:131927/150‑1 (MQ=255)                                                                      ggcgCCACCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCGGCCCCGTGC‑TACGCGCGGGGTTGACAGAGGTATTGGTCACCGGGATCGAGATCAGGTGGATCAGGGTCAGGGCCAGGCCGATGGCGATGGGt                                                                                <  1:628580/150‑1 (MQ=255)                                                                            aCCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCGGCCCCGTGC‑TACGCGCGGGGTTGACAGAGGTATTGGTCACCGGGATCgag                                                                                                                          <  2:199063/102‑1 (MQ=255)                                                                            aCCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCGGCCCCGTGC‑TACGCGCGGGGTTGACAGAGGTATTGGTCACCGGGATCgag                                                                                                                          >  1:199063/1‑102 (MQ=255)                                                                            aCCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCGGCCCCGTGC‑TACGCGCGGGGTTGACAGAGGTATTGGTCACCGGGATCGAGATCAGGTGGATCAGGGTCAGGGCCAGGCCGATGGCGATGGGTGCaaaa                                                                          >  1:583568/1‑150 (MQ=255)                                                                              cctgaccctgatccacctgatctcgatcccggtgaccaatacctctgtcaaCCCCGTGC‑TACGCGCGGGGTTGa                                                                                                                                                    >  1:658291‑M1/52‑74 (MQ=255)                                                                              cctgaccctgatccacctgatctcgatcccggtgaccaatacctctgtcaaCCCCGTGC‑TACGCGCGGGGTTGa                                                                                                                                                    <  2:658291‑M1/23‑1 (MQ=255)                                                                                       ccACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCGGCCCCGTGC‑TACGCGCGGGGTTGACAGAGGTATTGGTCACCGGGATCGAGATCAGGTGGATCAGGGTCAGGGCCAGGCcgatggcgatgg                                                                                  >  1:28192/1‑131 (MQ=255)                                                                                       ccACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCGGCCCCGTGC‑TACGCGCGGGGTTGACAGAGGTATTGGTCACCGGGATCGAGATCAGGTGGATCAGGGTCAGGGCCAGGCcgatggcgatgg                                                                                  <  2:28192/131‑1 (MQ=255)                                                                                       ccACAGCTGCTGCAGCGCCCAGCCACCGACGAACAGGCCCGGCCCCGTGC‑TACGCGCGGGGTTGACAGAGGGATTGGTCACCGGGATCGAGGTCAGGTGGCTCAGGGTCAGGGCCAGGCCGATGGCGATGGGTGCAAAACCTGCCGGtcc                                                               >  2:377358/1‑148 (MQ=255)                                                                                          caGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCGGCCCCGTGC‑TACGCGCGGGGTTGACAGAGGTATTGGTCACCGGGATCGAGATCAGGTGGATCAGGGTCAGGGCCAGGCCGATGGCGATGGGTGCAAAACCTGCCGGTGcgcg                                                            <  2:583568/150‑1 (MQ=255)                                                                                                   cAGGGCCCAGCCACCGACGAACAGGGCCGGCCCCGTGC‑TACGCGCGGGGTTGACAGAGGTATTGGTCACCGGGATCGAGATCAGGTGGATCAGGGTCAGGGCCAGGCcgatgg                                                                                        >  2:792336/1‑113 (MQ=255)                                                                                                   cAGGGCCCAGCCACCGACGAACAGGGCCGGCCCCGTGC‑TACGCGCGGGGTTGACAGAGGTATTGGTCACCGGGATCGAGATCAGGTGGATCAGGGTCAGGGCCAGGCcgatgg                                                                                        <  1:792336/113‑1 (MQ=255)                                                                                                     gggCCCAGCCACCGACGAACAGGGCCGGCCCCGTGC‑TACGCGCGGGGTTGACAGAGGTATTGGTCACCGGGATCGAGAt                                                                                                                        <  2:41810/79‑1 (MQ=255)                                                                                                     gggCCCAGCCACCGACGAACAGGGCCGGCCCCGTGC‑TACGCGCGGGGTTGACAGAGGTATTGGTCACCGGGATCGAGAt                                                                                                                        >  1:41810/1‑79 (MQ=255)                                                                                                               aCCGACGAACAGGGCCGGCCCCGTGC‑TACGCGCGGGGTTGACAGAGGTATTGGTCACCGGGATCg                                                                                                                            >  2:521626/1‑65 (MQ=255)                                                                                                               aCCGACGAACAGGGCCGGCCCCGTGC‑TACGCGCGGGGTTGACAGAGGTATTGGTCACCGGGATCg                                                                                                                            <  1:521626/65‑1 (MQ=255)                                                                                                               aCCGACGAACAGGGCCGGCCCCGTGC‑TACGCGCGGGGTTGACAGAGGTATTGGTCACCGGGATCGAGATCAGGTGGATCAGGGTCAGGGCCAGGCCGATGGCGATGGGTGCAAAACCTGCCGGTGCGCGGGCGTCGGTCGCCCCCatgat                                       >  1:487911/1‑150 (MQ=255)                                                                                                                 ccaatacctctgtcaaCCCCGCGCGTA‑GCACGGGGTTGACAGAGGTATTGGTCACCGGGATCGAGATCAGGTGGATCAGGGTCAGGGCCAGGCcgatggc                                                                                       <  2:140362‑M1/84‑1 (MQ=255)                                                                                                                 ccaatacctctgtcaaCCCCGCGCGTA‑GCACGGGGTTGACAGAGGTATTGGTCACCGGGATCGAGATCAGGTGGATCAGGGTCAGGGCCAGGCcgatggc                                                                                       >  1:140362‑M1/17‑100 (MQ=255)                                                                                                                  gacgaACAGGGCCGGCCCCGTGC‑TACGCGCGGGGTTGACAGAGGTATTGGTCACCGGGATCGAGATCAGGTGGATCAGGGTCAGGGCCAGGCCGAGGGCGATGGGTGCAAAACCTGCCGGTGCGCGGGCGTCGGTCGCCCCCATGATGac                                    <  1:221004/150‑1 (MQ=255)                                                                                                                         aGGGCCGGCCCCGTGC‑TACGCGCGGGGTTGACAGAGGTATTGGTCACCGGGATCGAGATCAGGTGGAtcagggtcagg                                                                                                     >  2:155311/1‑78 (MQ=255)                                                                                                                         aGGGCCGGCCCCGTGC‑TACGCGCGGGGTTGACAGAGGTATTGGTCACCGGGATCGAGATCAGGTGGAtcagggtcagg                                                                                                     <  1:155311/78‑1 (MQ=255)                                                                                                                              tcaaccccgcgcgtagcacggggTTGACAGAGGTATTGGTCACCGGGATCGAGATCAGGTGGAtcagggtcagg                                                                                                     >  1:658291‑M2/24‑74 (MQ=255)                                                                                                                              tcaaccccgcgcgtagcacggggTTGACAGAGGTATTGGTCACCGGGATCGAGATCAGGTGGAtcagggtcagg                                                                                                     <  2:658291‑M2/51‑1 (MQ=255)                                                                                                                                  CCCGTGC‑TACGCGCGGGGTTGACAGAGGTATTGGTCACCGGGATCGAGATCAGGTGGATCAGGGTCAGGGCCAGGCCGATGGCGATGGGTGCAAAACCTGCCGGTGCGCGGGCGTCGGTCGCCCCCATGATGACCACCAGGAACATCGcc                    <  2:23258‑M1/150‑1 (MQ=255)                                                                                                                                                     ttGACAGAGGTATCGGTCACCGGGATCGAGATCAGGTGGATCAGGGTCAGGGCCAGGCCGATGGCGATGGGTGCAAAACCTGCCGGTGCGCGGGCGTCGGTCGCCCCCATGATGACCACCAGGAACATCGCCGTCATCACCACTTCACTg  >  2:713920/1‑150 (MQ=255)                                                                                                                                                     |                                                                                                                                                      GTGCACCAGCGTTAAGGTTCTTTCGCCAGGCCCCGATACAGCGCGCCGCCGATAGCAGCACCGATCAGCGGCGCCACCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCGGCCCCGTGC‑TACGCGCGGGGTTGACAGAGGTATTGGTCACCGGGATCGAGATCAGGTGGATCAGGGTCAGGGCCAGGCCGATGGCGATGGGTGCAAAACCTGCCGGTGCGCGGGCGTCGGTCGCCCCCATGATGACCACCAGGAACATCGCCGTCATCACCACTTCACTG  >  NC_002947/4872961‑4873257

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