Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
JC	NC_000913	3,176,344	Δ3 bp	17.8%	coding (488‑490/828 nt)	cpdA ←	3',5' cAMP phosphodiesterase

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	3176347 =	69 (0.750)	14 (0.160)	14/240	NT	17.8%	coding (487/828 nt)	cpdA	3',5' cAMP phosphodiesterase
	?	NC_000913	= 3176343	65 (0.760)					coding (491/828 nt)	cpdA	3',5' cAMP phosphodiesterase

CCAGTCGAGATCCAGCTCCTGATGAATATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATCCAGCAACAGGATTTGCCATTGCTCACCAATAAACACGCG  >  NC_000913/3176198‑3176491                                                                                                                                                  |                                                                                                                                                     ccAGTCGAGATCCAGCTCCTGATGAATATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATgc                                                                                                                                                    >  3:150040‑M2/1‑146 (MQ=255) ccAGTCGAGATACAGCTCCTGATGAATATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGAGGATGATGATg                                                                                                                                                     <  5:322233‑M2/147‑2 (MQ=255)     tCGAGATCCAGCTCCTGATGAATATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATgcagcag                                                                                                                                               <  7:274194‑M2/149‑8 (MQ=255)     tCGAGATCCAGCTCCTGATGAATATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATgcagcag                                                                                                                                               <  8:335151‑M2/149‑8 (MQ=255)       gagaTCCAGCTCCTGATGAATATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATgcagcagca                                                                                                                                             >  2:272036‑M2/1‑140 (MQ=255)            ccAGCTCCTGATGAATATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAAc                                                                                                                                        >  6:334130/1‑149 (MQ=255)              aGCTCCTGATGAATATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACg                                                                                                                                       >  2:86405/1‑148 (MQ=255)              aGCTCCTGATGAATATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGt                                                                                                                                      <  2:310652/149‑1 (MQ=255)               gCTCCTGATGAATATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTAGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTa                                                                                                                                     >  4:307397/1‑149 (MQ=255)                   cTGATGAATATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATg                                                                                                                                   <  1:271503/147‑1 (MQ=255)                   cTGATGAATATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGc                                                                                                                                 >  1:144554/1‑149 (MQ=255)                    tgatgaATATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCg                                                                                                                                >  3:238251/1‑149 (MQ=255)                      atgaATATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATgcagcagcaacgtatggcgttc                                                                                                                                >  3:255340‑M2/1‑125 (MQ=255)                           tatGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGgc                                                                                                                         <  7:335151/149‑1 (MQ=255)                           tatGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGgc                                                                                                                         <  5:74706/149‑1 (MQ=255)                               aCCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTTCAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATCGCGTTCTGGCGCAt                                                                                                                     >  7:362297/1‑149 (MQ=255)                               aCCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCAt                                                                                                                     >  7:326626/1‑149 (MQ=255)                                ccGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATc                                                                                                                    <  5:158101/149‑1 (MQ=255)                                  gCACAGCAAGTATTTGACGTTCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATgcagcagcaacgtatggcgttctggcgcatcggc                                                                                                                    <  7:312522‑M2/147‑35 (MQ=255)                                      aGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAg                                                                                                              <  8:158424/149‑1 (MQ=255)                                       gCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGAGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGt                                                                                                             >  4:92578/1‑149 (MQ=255)                                        cAAGTATTTGACGTGCGGAAACTTCGCCAGCTCGGTATCCAGTTCGCCCGCGTTACGCAGCCTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGAGGATGATGCAGCAGCCGCAACGTATGGCGTTCTGGCGCCGCGGCCAGtt                                                                                                            >  1:284722/1‑149 (MQ=255)                                          aGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTa                                                                                                          <  8:274194/149‑1 (MQ=255)                                          aGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTa                                                                                                          >  4:294866/1‑149 (MQ=255)                                          aGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTa                                                                                                          <  3:294866/149‑1 (MQ=255)                                          aGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTa                                                                                                          <  3:62333/149‑1 (MQ=255)                                           gTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATgcagca                                                                                                                                                <  4:215009‑M2/110‑7 (MQ=255)                                           gTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATgcagca                                                                                                                                                >  3:215009‑M2/1‑104 (MQ=255)                                              tttGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATgcagcagcaacgtatggcgttctggcgcatcggccagtttacgttc                                                                                                        <  6:393352‑M2/147‑47 (MQ=255)                                                  aCGTGCGGAAACTTCGCCAGCTCGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCTACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCa                                                                                                   <  8:175353/148‑1 (MQ=255)                                                   cGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGc                                                                                                 <  5:138172/149‑1 (MQ=255)                                                       cGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCAc                                                                                              <  7:408519/148‑1 (MQ=255)                                                         gAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTAGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTc                                                                                            >  1:320222/1‑148 (MQ=255)                                                              ttCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCt                                                                                       >  1:425589/1‑148 (MQ=255)                                                                cGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTATCGGACGATGATgcagcagcaacgtatggcgttctggcgcatcggccagtttacgttccagccactcaagctgaaa                                                                                      >  4:225500‑M2/1‑83 (MQ=255)                                                                cGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGa                                                                                     >  3:18004/1‑148 (MQ=255)                                                                 gCCAGCACGGTATCCAGTCCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCAGCGGCCAGTTTACGTTCCAGCCACTCAAGCTGaaa                                                                                   <  2:284722/149‑1 (MQ=255)                                                                        cGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCgctc                                                                            <  4:18004/149‑1 (MQ=255)                                                                             tCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGCGCCAACTACAACCCGAAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCCGCCACTCCAGCGGAAACTcgctccgctc                                                                       >  6:343087/1‑149 (MQ=255)                                                                             tCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCgctcagct                                                                        >  5:256090/1‑148 (MQ=255)                                                                               cAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCAc                                                                     >  1:302316/1‑149 (MQ=255)                                                                                aGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCAcc                                                                    >  3:179376/1‑149 (MQ=255)                                                                                aGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCAcc                                                                    <  2:144554/149‑1 (MQ=255)                                                                                 gTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATgcagcagcaacgtatggcgttctggcgcatcggccagtttacgttccagccactcaagctgaaactcgctcagctcaccgtgc                                                                   >  1:195042‑M2/1‑66 (MQ=255)                                                                                     gCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGc                                                                                                                 >  2:350918/1‑99 (MQ=255)                                                                                     gCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGc                                                                                                                 <  1:350918/99‑1 (MQ=255)                                                                                       ccGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCgg                                                              >  1:143681/1‑148 (MQ=255)                                                                                               cGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTGCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCGCAGCTCACCGGGCGGCACGCCaaa                                                     <  1:86405/149‑1 (MQ=255)                                                                                                 cAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAaca                                                   >  6:57097/1‑149 (MQ=255)                                                                                                  aGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATgcagcagcaacgtatggcgttctggcgcatcggccagtttacgtgccagccactcaagctgaaactcgctcagctcaccgtgcggcacgccaaacacctg                                                  <  1:304989‑M2/149‑101 (MQ=255)                                                                                                  aGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCTGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGGCAAacac                                                  <  8:326626/149‑1 (MQ=255)                                                                                                   gACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATgcagcagcaacgtatggcgttctggcgcatcggccagtttacgttccagccactcaagctgaaactcgctcagctcaccgtgcggcacgccaaacacctgg                                                 >  5:393352‑M2/1‑48 (MQ=255)                                                                                                       gtgtTGATCGAGCCAACTACAACCCGCAGGTAGCGCAAGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCACACCGTGCGGCACGCCAAACCCCAGGc                                             >  6:98986/1‑149 (MQ=255)                                                                                                        tgtTGATCGCGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCCAACACCTGGc                                             >  8:315129/1‑148 (MQ=255)                                                                                                        tgtTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTgg                                              <  4:319883/147‑1 (MQ=255)                                                                                                         gtTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTa                                                                                                                                     <  5:275086/59‑1 (MQ=255)                                                                                                         gtTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTa                                                                                                                                     >  6:275086/1‑59 (MQ=255)                                                                                                            gATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATc                                         >  1:159560/1‑148 (MQ=255)                                                                                                             aTCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCCCTCAAGCTGAACCTCGCTCAGCTCACCGTGCGGCACGCCAACCACCTGGCTATc                                         >  5:126802/1‑147 (MQ=255)                                                                                                             aTCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATCCa                                       <  2:320222/149‑1 (MQ=255)                                                                                                                 aGCCAACTACAACCCGCAGGTAGCGGATGATGATgcagcagcaacgtatggcgttctggcgcatcggccagtttacgttccagccactcaagctgaaactcgctcagctcaccgtgcggcacgc                                                            <  7:65582‑M2/124‑91 (MQ=255)                                                                                                                 aGCCAACTACAACCCGCAGGTAGCGGATGATGATgcagcagcaacgtatggcgttctggcgcatcggccagtttacgttccagccactcaagctgaaactcgctcagctcaccgtgcggcacgc                                                            >  8:65582‑M2/1‑34 (MQ=255)                                                                                                                   ccAACTACAACCCGCAGGTAGCGGATGATGATgcagcagcaacgtatggcgttctggcgcatcggccagtttacgttccagccactcaagctgaaactcgctcagctcaccgtgcggcacgccaaacacctggctatccagcaacagga                                 <  6:150512‑M2/149‑118 (MQ=255)                                                                                                                   ccAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCa                                                                                           >  7:137428/1‑91 (MQ=255)                                                                                                                   ccAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCa                                                                                           <  8:137428/91‑1 (MQ=255)                                                                                                                   ccAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATCCAGCAACa                                 <  2:159560/149‑1 (MQ=255)                                                                                                                     aaCTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATCCAGCAACAgg                               <  3:307397/149‑1 (MQ=255)                                                                                                                      aCTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGc                                                                >  2:12765/1‑115 (MQ=255)                                                                                                                        tACAACCCGAATGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCTCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATCCAGCAACAGGAtt                            >  6:290672/1‑149 (MQ=255)                                                                                                                          cAACCCGCAGGTAGCGGATGATGATgcagcagcaacgtatggcgttctggcgcatcggccagtttacgttccagccactcaagctgaaactcgctcagctcaccgtgcggcacgccaaacacctggctatccagcaacaggatttgc                            >  8:312522‑M2/1‑25 (MQ=255)                                                                                                                          cAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATCCAGCAACAGGAttt                           <  6:322233/148‑1 (MQ=255)                                                                                                                           aaCCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATCCAGCAACAGGATTTGc                         >  6:190817/1‑149 (MQ=255)                                                                                                                             gccGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATCCAGCAACAGGATTTGCCa                       <  8:362297/148‑1 (MQ=255)                                                                                                                             cccGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATCCAGCAACAGGATTTGCCa                       <  5:334130/149‑1 (MQ=255)                                                                                                                              ccGCAGGTAGCGTATGATGATGCAGCAGCAGCAACGTATGGAGTTCTGGCGCATCGGCCAGTTTACGTTCAAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATCCAGCAACAGGATTTGCCAt                      <  5:98986/149‑1 (MQ=255)                                                                                                                               cGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATCCAGCAACAGGATTTGCCAtt                     <  5:57097/149‑1 (MQ=255)                                                                                                                                gCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGGTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACGCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATCCAGCAACAGGATTTGCCATTg                    <  4:238251/149‑1 (MQ=255)                                                                                                                                  aGGTAGCGGATGATGATgcagcagcaacgtatggcgttctggcgcatcggccagtttacgttccagccactcaagctgaaactcgctcagctcaccgtgcggcacg                                                             >  1:99840‑M2/1‑17 (MQ=255)                                                                                                                                  aGGTAGCGGATGATGATgcagcagcaacgtatggcgttctggcgcatcggccagtttacgttccagccactcaagctgaaactcgctcagctcaccgtgcggcacg                                                             <  2:99840‑M2/106‑90 (MQ=255)                                                                                                                                  aGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATCCAGCAACAGGATTTGCCATTGCt                  <  2:143681/149‑1 (MQ=255)                                                                                                                                    gTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATCCAGCAACAGGAttt                           <  7:56200/138‑1 (MQ=255)                                                                                                                                    gTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATCCAGCAACAGGAttt                           >  8:56200/1‑138 (MQ=255)                                                                                                                                       gCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATCCAGCAACAGGATTTGCCATTGCTCAcc              <  6:256090/148‑1 (MQ=255)                                                                                                                                          gatgatgatgcagcagcaacgtatggcgttctggcgcatcggccagtttacgttcccgccactcaagctgaaactcgctcagctcaccgtgcggcacgcccaacacctggctatccagcaacaggatttgccattgctcaccaataaa           <  6:132629‑M2/148‑140 (MQ=255)                                                                                                                                          gatgatgatGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATCCAGCAACAGGATTTGCCATTGCTc                 >  8:326115/1‑142 (MQ=255)                                                                                                                                          gatgatgatGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATCCAGCAACAGGATTTGCCATTGCTc                 <  7:326115/142‑1 (MQ=255)                                                                                                                                          gatgatgatGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACAGTGCGGCACCCCAAACACCTGGCTATCCCGCAACAGGATTTGCCAATGCTCACCAAt           >  1:106661/1‑148 (MQ=255)                                                                                                                                                gatgCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATCCAGCAACa                                 <  7:37527/120‑1 (MQ=255)                                                                                                                                                gatgCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATCCAGCAACa                                 >  8:37527/1‑120 (MQ=255)                                                                                                                                                gatgCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATCCAGCAACAGGATTTGCCATTGCTCACCAATAAACAcg    >  1:275235/1‑149 (MQ=255)                                                                                                                                                  tgcagcagcaacgtatggcgttctggcgcatcggccagtttacgttccagccactcaagctgaaactcgctcagctcaccgtgcggcacgccaaacacctggctatccagcaacaggatttgcctttgctcaccaataaacacgcgctt  >  5:150512‑M2/1‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                     CCAGTCGAGATCCAGCTCCTGATGAATATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATCCAGCAACAGGATTTGCCATTGCTCACCAATAAACACGCG  >  NC_000913/3176198‑3176491

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