New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	= 3926238	39 (0.770)	3 (0.060)	3/240	NT	7.9%	coding (253/1629 nt)	peaA	quinohemoprotein amine dehydrogenase subunit alpha
	?	NC_002947	= 3926260	34 (0.730)					coding (231/1629 nt)	peaA	quinohemoprotein amine dehydrogenase subunit alpha

AGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGCGAC‑GATCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGATGCTCAT‑‑CAGCCAGCCTTCCGGTGTCTTGCGCTGGTGGCTGATGCGGCTATAGGTGTCGTTGCCTTCGGGGATATGGCACCCCATGCACTTGTTCTGCAACAGGCTCGGGCCCTGCTCGGCCGCCATGGCCTGGGTGCTCAGCAGCGCGGCGGCGA  >  NC_002947/3926090‑3926381                                                                                                                                                        |                                                                                                                                                aGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGCGAC‑GATCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGATGCTCAT‑‑CAGCCA                                                                                                                                                 <  8:278728‑M1/149‑1 (MQ=255)   ccGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGCGAC‑GATCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGATGCTCAT‑‑cagccagc                                                                                                                                               <  4:54284/149‑1 (MQ=255)             aTGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGCGAC‑GATCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGATGCTCAT‑‑CAGCCAGCCTTCCGgtgt                                                                                                                                     >  5:119753/1‑149 (MQ=255)                gCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGCGAC‑GATCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGATGCTCAT‑‑CAGCCAGCCTTCCGGTGTCtt                                                                                                                                  <  1:304971/149‑1 (MQ=255)                  gTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGCGAC‑GATCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGATGCTCAT‑‑CAGCCAGCCTTCCGGTGTCTTgc                                                                                                                                >  8:19521/1‑149 (MQ=255)                     gcgcACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGCGAC‑GATCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGATGCTCAT‑‑CAGCCAGCCTTCCGGTGTCTTGCGCt                                                                                                                             >  1:124464/1‑149 (MQ=255)                             cATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGCGAC‑GATCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGATGCTCAT‑‑CCGCCAGCGTTCCGGTGTCTTGCGCTGGtggct                                                                                                                      <  8:284253‑M2/148‑6 (MQ=255)                               tCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGTGAC‑GATCGTCGTCGCTGATCTGAAGTCCATGCATCACCTGCATGCGGGCGATGCTCAT‑‑CAGCCAccagcgcaagacaccggaaggctggctgat                                                                                                                   <  7:65513‑M1/149‑31 (MQ=255)                                cGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGCGAAAGATCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGATGCTCAT‑‑CAGCCAGCCTTCCGGTGTCTTGCGCTGGtggctgat                                                                                                                   <  1:88349‑M2/149‑9 (MQ=255)                                   tttCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGCGAC‑GATCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGATGCTCAT‑‑CAGCCAGCCTTCCGGTGTCTTGCGCTGGTGGCTGATGCgg                                                                                                               <  8:182341/149‑1 (MQ=255)                                      cGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGCGAC‑GATCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGATGCTCAT‑‑CAGCCAGCCTTCCGGTGTCTTGCGCTGGTGGCTGATGCGGCt                                                                                                             <  7:230126/148‑1 (MQ=255)                                       gCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGCGAC‑GATCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGATGCTCAT‑‑CAGCCAGCCTTCCGGTGTCTTGCGCTGGtggctg                                                                                                                     <  5:265124‑M2/139‑7 (MQ=255)                                       gCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGCGAC‑GATCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGATGCTCAT‑‑CAGCCAGCCTTCCGGTGTCTTGCGCTGGtggctg                                                                                                                     >  6:265124‑M2/1‑133 (MQ=255)                                        cTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGCGAC‑GATCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGATGCTCAT‑‑CAGCCAGCCTTCCGGTGTCTTGCGCTGGTGGCTGATGCGGCtata                                                                                                          <  5:131052/149‑1 (MQ=255)                                         tGGGCGCCAGGCCCTGTTTCTCGGCCAGGTATTTGACCAGCGTGCGAC‑GATCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGATGCTCAT‑‑CAGCCAGCATTCCGGTGTCTTGCGCTGGTGGCTGATGCGGCTATAg                                                                                                         >  4:171244/1‑149 (MQ=255)                                            gcgcCAGGCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGCGAC‑GATCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGATGCTCAT‑‑CAGCCAGCCTTCCGGTGTCTTGCGCTGGTGGCTGATGCGGCTATAGgtg                                                                                                      <  1:222570/149‑1 (MQ=255)                                            gcgcCAGGCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGCGAC‑GATCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGATGCTCAT‑‑CAGCCAGCCTTCCGGTGTCTTGCGCTGGTGGCTGATGCGGCTATAGgtg                                                                                                      <  2:88349/149‑1 (MQ=255)                                                cAGGCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGCGAC‑GATCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGATGCTCAT‑‑CAGCCAccagcgcaagacaccggaaggctggctgatga                                                                                                                 <  2:18422‑M1/134‑33 (MQ=255)                                                cAGGCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGCGAC‑GATCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGATGCTCAT‑‑CAGCCAccagcgcaagacaccggaaggctggctgatga                                                                                                                 >  1:18422‑M1/1‑102 (MQ=255)                                                         tttGTCGGCCAGGTATTTGACCAGCGTGCGAC‑GATCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGATGCTCAT‑‑CAGCCAGCCTTCCGGTGTCTTGCGCTGGTGGCTGATGCGGCTATAGGTGTCGTTGCCTTCg                                                                                          <  5:36253/148‑1 (MQ=255)                                                                  cAGGTATTTGAACAGCGTGCGAC‑GATGGTCGTCGCTGATCTGAAGGCCATGCATCACCTGGATGCGGGCGATGCTCAT‑‑CAGCCAGCCTTCCGGTGTCTTGCGCTGGTGGCTGATGCGGCTATAGGTGTGGTTGCCTTCGGGGATATGGc                                                                                <  1:323649/149‑1 (MQ=255)                                                                     gTATTTGACCAGCGTGCGAC‑GATCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGATGCTCAT‑‑CAGCCAGCCTTCCGGTGTCTTGCGCTGGTGGCTGATGCGGCTATAGGTGTCGTTGCCTTCGGGGATATGGCa                                                                               <  2:179191/147‑1 (MQ=255)                                                                     gTATTTGACCAGCGTGCGAC‑GATCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGATGCTCAT‑‑CAGCCAGCCTTCCGGTGTCTTGCGCTGGTGGCTGATGCGGCTATAGGTGTCGTTGCCTTCGGGGATATGGCa                                                                               >  1:179191/1‑147 (MQ=255)                                                                        tttGACCAGCGTGCGAC‑GATCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGATGCTCAT‑‑CAGCCAGCCTTCCGGTGTCTTGCGCTGGTGGCTGATGCGGCTATAGGTGTCGTTGCCTTCGGGGATATGGCAcccc                                                                           <  2:124464/148‑1 (MQ=255)                                                                        tttGACCAGCGTGCGAC‑GATCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGATGCTCAT‑‑CAGCCAGCCTTCCGGTGTCTTGCGCTGGTGGCTGATGCGGCTATAGGTGTCGTTGCCTTCGGGGATATGGCACCCCa                                                                          >  2:7816/1‑149 (MQ=255)                                                                              cAGCGTGCGAC‑GATCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGATGCTCATAGCACCCCCCATACGGGTGTTTTGCGGGGTTGGATGAGGCGGATATAGGTGACTTTGCATGGGGGGATATGGCCCCCCATGCa                                                                      >  1:96135/1‑149 (MQ=255)                                                                              cAGCGTGCGAC‑GATCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGATGCTCAT‑‑CAGCCAGCCTTCCGGTGTCTTGCGCTGGTGGCTGATGCGGCTATAGGTGTCGTTGCCTTCGGGGATATGGCACCCCATGCACt                                                                    >  7:60372/1‑149 (MQ=255)                                                                                      gac‑gaTCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGATGCTCAT‑‑CAGCCAGCCTTCCGGTGTCTTGCGCTGGTGGCTGATGCGGCTATAGGTGTCGTTGCCTTCGGGGATATGGCACCCCATGCACTTGTTCTg                                                             >  1:43976/1‑148 (MQ=255)                                                                                      gac‑gaTCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGATGCTCAT‑‑CAGCCAGCCTTCCGGTGTCTTGCGCTGGTGGCTGATGCGGCTATAGGTGTCGTTGCCTTCGGGGATATGGCACCCCATGCACTTGTTCTg                                                             >  5:309135/1‑148 (MQ=255)                                                                                        c‑gaTCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGATGCTCAT‑‑CAGCCAGCCTTCCGGTGTCTTGCGCTGGTGGCTGATGCGGCTATAGGTGTCGTTGCCTTCGGGGATATGGCACCCCATGCACTTGTTCTGCaa                                                          >  1:86394/1‑149 (MQ=255)                                                                                           aTCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGATGCTCAT‑‑CAGCCAGCCTTCCGGTGTCTTGCGCTGGTGGCTGATGCGGCTATAGGTGTCGTTGCCTTCGGGGATATGGCACCCCATGCACTTGTTc                                                               <  6:243258/142‑1 (MQ=255)                                                                                           aTCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGATGCTCAT‑‑CAGCCAGCCTTCCGGTGTCTTGCGCTGGTGGCTGATGCGGCTATAGGTGTCGTTGCCTTCGGGGATATGGCACCCCATGCACTTGTTc                                                               >  5:243258/1‑142 (MQ=255)                                                                                           aTCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGATGCTCAT‑‑CAGCCAGCCTTCCGGTGTCTTGCGCTGGTGGCTGATGCGGCTATAGGTGTCGTTGCCTTCGGGGATATGGCACCCCATGCACTTGTTCTGCaa                                                          >  7:305041/1‑147 (MQ=255)                                                                                           aTCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGATGCTCAT‑‑CAGCCAGCCTTCCGGTGTCTTGCGCTGGTGGCTGATGCGGCTATAGGTGTCGTTGCCTTCGGGGATATGGCACCCCATGCACTTGTTCTGCAACa                                                        <  1:7816/149‑1 (MQ=255)                                                                                           aTCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGATGCTCAT‑‑CAGCCAGCCTTCCGGTGTCTTGCGCTGGTGGCTGATGCGGCTATAGGTGTCGTTGCCTTCGGGGATATGGCACCCCATGCACTTGTTCTGCAACa                                                        <  2:86394/149‑1 (MQ=255)                                                                                                 gtcgCTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGATGCTCAT‑‑CAGCCAGCCTTCCGGTGTCTTGCGCTGGTGGCTGATGCGGCTATAGGTGTCGTTGCCTTCGGGGATATGGCACCCCATGCACTTGTTCTGCAACAGGCTCg                                                  <  8:60372/149‑1 (MQ=255)                                                                                                  tcgCTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGATGCTCAT‑‑CAGCCAGCCTTCCGGTGTCTTGCGCTGGTGGCTGATGCGGCTATAGGTGTCGTTGCCTTCGGGGATATGGCACCCCATGCACTTGTTCTGCAACAGGCTCgg                                                 <  6:309135/149‑1 (MQ=255)                                                                                                     cTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGATGCTCAT‑‑CAGCCAGCCTTCCGGTGTCTTGCGCTGGTGGCTGATGCGGCTATAGGTGTCGTTGCCTTCGGGGatat                                                                                   <  7:184167/112‑1 (MQ=255)                                                                                                     cTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGATGCTCAT‑‑CAGCCAGCCTTCCGGTGTCTTGCGCTGGTGGCTGATGCGGCTATAGGTGTCGTTGCCTTCGGGGatat                                                                                   >  8:184167/1‑112 (MQ=255)                                                                                                     cTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGATGCTCAT‑‑CAGCCAGCCTTCCGGTGTCTTGCGCTGGTGGCTGATGCGGCTATAGGTGTCGTTGCCTTCGGGGATATGGCACCCCATGCACTTGTTCTGCAACAGGCTCGGGcc                                              <  8:305041/149‑1 (MQ=255)                                                                                                             aaGGCCATGCATCACCTGCATGCGGGCGATGCTCAT‑‑CAGCCAGCCTTCCGGTGTCTTGCGCTGGTGGCTGATGCGGCTATAGGTGTCGTTGCCTTCGGGGATATGGCACCCCATGCACTTGTTCTGCAACAGGCTCGGGCCCTGc                                          <  2:96135/145‑1 (MQ=255)                                                                                                                  cATGCATCACCTGCATGCGGGCGATGCTCAT‑‑CAGCCAGCCTTCCGGTGTCTTGCGCTGGTGGCTGATGCGGCTATAGGTGTCGTTGCCTTCGGGGATATGGCACCCCATGCACTTGTTCTGCAACAGGCTCGGGCCCTGCTCGgccgc                                  >  5:32957/1‑148 (MQ=255)                                                                                                                    tGCATCACCTGCATGCGGGCGATGCTCAT‑‑CAGCCAGCCTTCCGGTGTCTTGCGCTGGTGGCTGATGCGGCTATAGGTGTCGTTGCCTTCGGGGATATGGCACCCCATGCACTTGTTCTGCAACAGGCTCGGGCCCTGCTCGGCCGCCAt                               <  3:171244/149‑1 (MQ=255)                                                                                                                           ccTGCATGCGGGCGATGCTCAT‑‑CAGCCAGCCTTCCGGTGTCTTGCGCTGGTGGCTGATGCGGCTATAGGTGTCGTTGCCTTCGGGGATATGGCACCCCATGCACTTGTTCTGCAACAGGCTCGGGCGCTGCTCGGCCGCCATGGCCTgg                        <  6:119753/149‑1 (MQ=255)                                                                                                                           ccTGCAGGCGGGCGCCGCTCAG‑‑CAGCCCGCCTTCCGGGGTCGTGCGCTGGTGGCTGATGCGGCTATAGGTGTCGTTGCCTTCGGGGATATGGCACCCCATGCACTTGTTCTGCAACAGGCTCGGGCCCTGCTCGGCCGCCATGGCCTgg                        <  6:231222/149‑1 (MQ=255)                                                                                                                                          tGCTCAT‑‑CAGCCAGCCTTCCGGTGTCTTGCGCTGGTGGCTGATGCGGCTATAGGTGTCGTTGCCTTCGGGGATATGGCACCCCATGCACTTGTTCTGCAACAGGCTCGGGCCCTGCTCGGCCGCCATGGCCTGGGTGCTCAGCAGcgcg         <  3:147348/149‑1 (MQ=255)                                                                                                                                             tcat‑‑caGCCAGCCTTCCGGTGTCTTGCGCTGGtggctgatgagcatcgcccgcatgcaggtgatgcatggccttcagatcagcgacgacgatcgtcgcacgctggtcaaatacctggccgacaaacagggcctg                     <  2:18422‑M2/134‑103 (MQ=255)                                                                                                                                             tcat‑‑caGCCAGCCTTCCGGTGTCTTGCGCTGGtggctgatgagcatcgcccgcatgcaggtgatgcatggccttcagatcagcgacgacgatcgtcgcacgctggtcaaatacctggccgacaaacagggcctg                     >  1:18422‑M2/1‑32 (MQ=255)                                                                                                                                               at‑‑caGCCAccagcgcaagacaccggaaggctggctgatgagcatcgcccgcatgcaggtgatgcatggccttcagatcagcgacgacgatctttcgcacgctggtcaaatacctggccgacaaacagggcctggcgcccagcgaaaccg    <  1:88349‑M1/149‑142 (MQ=255)                                                                                                                                               at‑‑caGCCAGCCTTCCGGTGTCTTGCGCTGGtggctgatgagcatcgcccgcatgcaggtgatgcatggacttcagatcagcgacgacgatcgtcacacgctggtcaaatacctggccgacaaacagggcctggcgcccagcgaaaccga    <  7:65513‑M2/149‑120 (MQ=255)                                                                                                                                                t‑‑caGCCAGCCTTCCGGTGTCTTGCGCTGGTGGCTGATGCGGCTATAGGTGTCGTTGCCTTCGGTGATATGGCACCCCATGCACTTGTTCTGCAACAGGCGCGGGCCCTGCTCGGCCGCCATGGCCTGGGTGCTCAGCAGCgcggcggcg   <  7:83977/149‑1 (MQ=255)                                                                                                                                                t‑‑caGCCAGCCTTCCGGGTTCTTGCTCTGGTGGCTGATGCGGCTATAGGTGTCGTTGCCTTCGGGGATATGGCACCCCATGCACTTGTTCTGCAACAGGCTCGGGCCCTGCTCGGCCGCCATGGCCTGGGTGCTCAGCAGCgcggcggcg   <  7:88380/149‑1 (MQ=255)                                                                                                                                                   caGCCAccagcgcaagacaccggaaggctggctgatgagcatcgcccgcatgcaggtgatgcatggccttcagatcagcgacgacgatcgtcgcacgctggtcaaatacctggccgacaaacagggcctggcgcccagc            <  5:265124‑M1/139‑134 (MQ=255)                                                                                                                                                   caGCCAccagcgcaagacaccggaaggctggctgatgagcatcgcccgcatgcaggtgatgcatggccttcagatcagcgacgacgatcgtcgcacgctggtcaaatacctggccgacaaacagggcctggcgcccagc            >  6:265124‑M1/1‑6 (MQ=255)                                                                                                                                                    aGCCAccagcgcaagacaccggaacgctggcggatgagcatcgcccgcatgcaggtgatgcatggccttcagatcagcgacgacgatcgtcgcacgctggtcaaatacctggccgacaaacagggcctggcgcccagcgaaaccgatg  <  8:284253‑M1/148‑144 (MQ=255)                                                                                                                                                        |                                                                                                                                                AGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGCGAC‑GATCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGATGCTCAT‑‑CAGCCAGCCTTCCGGTGTCTTGCGCTGGTGGCTGATGCGGCTATAGGTGTCGTTGCCTTCGGGGATATGGCACCCCATGCACTTGTTCTGCAACAGGCTCGGGCCCTGCTCGGCCGCCATGGCCTGGGTGCTCAGCAGCGCGGCGGCGA  >  NC_002947/3926090‑3926381

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