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

ACCTACGCCGCCGGGCGCTGCGTGCACCAGCGTTAAGGTTCTTTCGCCAGGCCCCGATACAGCGCGCCGCCGATAGCAGCACCGATCAGCGGCGCCACCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCGGCCCCGTGC‑TACGCGCGGGGTTGACAGAGGTATTGGTCACCGGGATCGAGATCAGGTGGATCAGGGTCAGGGCCAGGCCGATGGCGATGGGTGCAAAACCTGCCGGTGCGCGGGCGTCGGTCGCCCCCATGATGAC  >  NC_002947/4872940‑4873223                                                                                                                                                      |                                                                                                                                        aCCTACGCCGCCGGGCGCTGCGTGCACCAGAGTTAAGGTTCTTTCGCCAGGCCCCGATACAGCGCGCCGCCGATAGCAGCACCGATCAGCGGCGCCACCCAGAACAGCCACAGCTGCGGCAGGGCCCAGCCCCCGACGCACAGgggcggc                                                                                                                                         >  2:375003/1‑150 (MQ=255)  ccTACGCCGCCGGGCGCTGCGTGCACCAGCGTTAAGGTTCTTTCGCCAGGCCCCGATACAGCGCGCCGCCGATAGCAGCACCGATCAGCGGCGCCACCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGggccggcc                                                                                                                                        <  2:65856/150‑1 (MQ=255)         cgccgGGCGCTGCGTGCACCAGCGTTAAGGTTCTTTCGCCAGGCCCCGATACAGCGCGCCGCCGATAGCAGCACCGATCAGCGGCGCCACCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCGGCCCCGTGC‑t                                                                                                                                >  1:510914/1‑150 (MQ=255)           ccgGGCGCTGCGTGCACCAGCGTTAAGGTTCTTTCGCCAGGCCCCGATACAGCGCGCCGCCGATAGCAGCACCGATCAGCGGCGCCACCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCGGCCCCGTGC‑TAc                                                                                                                              <  1:724953/150‑1 (MQ=255)                  cTGCGTGCACCAGCGTTAAGGTTCTTTCGCCAGGCCCCGATACAGCGCGCCGCCGATAGCAGCACCGATCAGCGGCGCCACCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCGGCCCCGTGC‑TACGCGCggg                                                                                                                       <  2:608323/150‑1 (MQ=255)                   tGCGTGCACCAGCGTTAAGGTTCTTTCGCCAGGCCCCGATACAGCGCGCCGCCGATAGCAGCACCGATCAGCGGCGCCACCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCGGCCCCGTGC‑TACGCGCgggg                                                                                                                      >  1:770066/1‑150 (MQ=255)                      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‑TAcgcgcg                                                                                                                         >  2:157059/1‑124 (MQ=255)                                          tttCGCCAGGCCCCGATACAGCGCGCCGCCGATAGCAGCACCGATCAGCGGCGCCACCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCGGCCCCGTGC‑TAcgcgcg                                                                                                                         <  1:157059/124‑1 (MQ=255)                                          tttCGCCAGGCCCCGATACAGCGCGCCGCCGATAGCAGCACCGATCAGCGGCGCCACCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCGGCCCCGTGC‑TACGCGCGGGGTTGACAGAGGTATTgg                                                                                                      >  1:227538/1‑143 (MQ=255)                                          tttCGCCAGGCCCCGATACAGCGCGCCGCCGATAGCAGCACCGATCAGCGGCGCCACCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCGGCCCCGTGC‑TACGCGCGGGGTTGACAGAGGTATTgg                                                                                                      <  2:227538/143‑1 (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)                                                                                           ggcgCCACCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCGGCCCCGTGC‑TACGCGCGGGGTTGACAGAGGTATTGGTCACCGGGATCGAGATCAGGTGGATCAGGGTCAGGGCCAGGCCGATGGCGATGGGt                                              <  1:628580/150‑1 (MQ=255)                                                                                           ggcgCCACCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCGGCCCCGTGC‑TACGCGCGGGGTTGACAGAGGTATTGGTCACCGGGATCGAGATCAGGTGGATCAGGGTCAGGGCCAGGCCGATGGCGATGGGt                                              <  1:131927/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                                                                                                                  <  2:658291‑M1/23‑1 (MQ=255)                                                                                                   cctgaccctgatccacctgatctcgatcccggtgaccaatacctctgtcaaCCCCGTGC‑TACGCGCGGGGTTGa                                                                                                                  >  1:658291‑M1/52‑74 (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)                                                                                                                                                      |                                                                                                                                        ACCTACGCCGCCGGGCGCTGCGTGCACCAGCGTTAAGGTTCTTTCGCCAGGCCCCGATACAGCGCGCCGCCGATAGCAGCACCGATCAGCGGCGCCACCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCGGCCCCGTGC‑TACGCGCGGGGTTGACAGAGGTATTGGTCACCGGGATCGAGATCAGGTGGATCAGGGTCAGGGCCAGGCCGATGGCGATGGGTGCAAAACCTGCCGGTGCGCGGGCGTCGGTCGCCCCCATGATGAC  >  NC_002947/4872940‑4873223

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.