Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
JC JC	NC_000913	3,992,705	IS5 (+) +4 bp	100%	coding (1553‑1556/2547 nt)	cyaA →	adenylate cyclase

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	2288919 =	NA (NA)	100 (1.010)	81/322	NT	100%	noncoding (1195/1195 nt)	IS5	repeat region
	?	NC_000913	3992705 =	0 (0.000)					coding (1553/2547 nt)	cyaA	adenylate cyclase
*	?	NC_000913	= 2290113	NA (NA)	137 (1.390)	96/322	NT	100%	noncoding (1/1195 nt)	IS5	repeat region
	?	NC_000913	= 3992708	0 (0.000)					coding (1556/2547 nt)	cyaA	adenylate cyclase

GTCTGGTTTGCCGCGCCACCCATTTGAACAATTTGATAACCGTCATTTAATACTGAGTATTCTGCGGCGCCATTCATTTCAATGTACTGATACCCTCCGTTAATTTCAGTATTATTACTTACACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCACTGACGTATCATTTGGTCCGCCCGAAACAGGTTGGCCAGCGTGAATAACATCGCCAGTTGGTTATCGTTTTTCAGCAACCCCTTGTATCTGGCTTTCACGAAGC  >  NC_000913/2288647‑2289186                                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                                            ggttGGTTTGCAGGGCCACCCATTTGAACAATTTGATAACCGTCATTTAATACTGAGTATTCTGCGGCGCCATGCATTTAAATGTACTGATACCCTCCGTGAATTTCAGTATTATTACTTACACCAAATGCGTTTATATGCTGTTCTCCGCCAGGATTGATTTTTGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTGTTCCAGGACAACACCATCAACTGTCTAACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATaagga                                                                                                                                                                                                                                                                           <  4:64543/272‑1 (MQ=255)                                                        aGTATTCTGCGGCGCCATTCATTTCAATGTACTGATACCCTCCGTTAATTTCAGTATTATTACTTACACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGTTTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGAAATCTTTTTCCAGTACAACACCATCAACTGTCTCCCCATACTCAACAGTTCATGCCAATGCAATTGGTGAATATACGGAAGGTGCGACTACGCGGGGAACTTCTTCCCGGCTGTCTTATTCGTTTACTttcttc                                                                                                                                                                                                                    >  4:300056/1‑275 (MQ=255)                                                        aGTATTCTGCGGCGCCATTCATTTCAATGTACTGATACCCTCCGTTAATTTCAGTATTATTACTTACACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATttattc                                                                                                                                                                                                                    >  1:72221/1‑275 (MQ=255)                                                                    cgcCATTCATTTCAATGTAATGATACCCTCCGTTAATTTCAGTATTATTACTTACACTAAATTCTTTTATATTCTGTTCTCCGCCAGGATTGACTTTCGTATTTTCGCCTGTCCCATAAACCAGCTGGCTACCTGTTTCAAGGTCAACACCATCACCTGTCTCACCATACCCAACAGTTGATGCCAATGCCATGGGTGGATTTAAGGAAGGTTAGAATAAGCGGGGAATTTCTTCGCGGCCGACTCAGTCATTGCATTTATTAATGTTTGacgcg                                                                                                                                                                                                        >  4:425885/1‑271 (MQ=255)                                                                           cattTCAATGTACTGATACCCTCCGTTAATTTCAGTATTATTACTTACACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGAtttttt                                                                                                                                                                                                 <  2:300060/275‑1 (MQ=255)                                                                                                   tttAATTTCATTATTATTACTTACCCCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCTTATTATTGCCGGTCCCCTAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTTTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGc                                                                                                                                                                         <  4:301152/274‑1 (MQ=255)                                                                                                   gTTAATTTCAGTATTATTACTTACACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGt                                                                                                                                                                                                                                  <  2:141835/218‑1 (MQ=255)                                                                                                                        tacACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAACTATTCTCGGCTGACTCAGTCATTTCGTTTCTTAATGTTTGAGCCGATTTTTTCTCCCTTAAATGCCTTGAATCAGCCTTTTTAGAcctgtctcttt                                                                                                                                                     >  4:363683/1‑265 (MQ=255)                                                                                                                        tacACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGAGGGGAAATTCTTCTCTGCTGACTCAGTCCTTTCCTTTCTTCATGTTTGTGCCGATTTTTTCTCCCGTCAATCCCTTGAATCTGCCTATTTAGACCGTTTCTTCGc                                                                                                                                                    >  1:9048/1‑275 (MQ=255)                                                                                                                            ccAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAg                                                                                                                                                                          <  2:363688/249‑1 (MQ=255)                                                                                                                                aTTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATaaggaa                                                                                                                                                                                                                                                                          <  7:249943/149‑1 (MQ=255)                                                                                                                                     tttATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCCACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATACAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCCCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGCGCCGATTTTTTCTACCGTAAATGCCTTTAATCAGCTTTTTTAGACCGTTTCTTCCCCATTTAAGGCttt                                                                                                                                       >  4:489334/1‑275 (MQ=37)                                                                                                                                              tGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGaa                                                                                                                                                                              >  4:274156/1‑227 (MQ=255)                                                                                                                                                ttCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGACTATACGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTCTTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGtt                                                                                                                             >  2:459667/1‑274 (MQ=255)                                                                                                                                                ttCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGaa                                                                                                                                                                              >  2:117245/1‑225 (MQ=255)                                                                                                                                                ttCTCCGACATTATTTATTTTCTTTTTATNNCCGGTCCCATAAACCAGCTGCTTATCTTTTNCCAGGCCAACACCACAAACTGTCCAACAAAACTCAACATTTGATGCCAATGCTTTGGGTTAATATAAGGAATGTNCGAATAAGCCGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGttt                                                                                                                            <  3:72221/275‑1 (MQ=18)                                                                                                                                                  ctcCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAg                                                                                                                                                                          <  2:274160/227‑1 (MQ=255)                                                                                                                                                    ccGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGAtttttt                                                                                                                                                                                                 <  1:46423/202‑1 (MQ=255)                                                                                                                                                    acGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAg                                                                                                                                                                          <  4:117244/224‑1 (MQ=255)                                                                                                                                                     cGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCATCTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTa                                                                                                                         >  2:117800/1‑273 (MQ=255)                                                                                                                                                     cGCCAGGATTGATTTTCGTATTATTGGCGGGCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGt                                                                                                                       <  2:489339/275‑1 (MQ=255)                                                                                                                                                      gCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAAtt                                                                                                                                                                                                                                                    >  7:252930/1‑149 (MQ=255)                                                                                                                                                         aGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACGATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAAttctt                                                                                                                                                                                                                                                 <  8:271660/149‑1 (MQ=255)                                                                                                                                                            attgattTTCGTATTATTGGCGGTCCCATAAACCAGCTTGATATCTTTCTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTtctc                                                                                                                                                                                                                                              <  8:158293/149‑1 (MQ=255)                                                                                                                                                             ttgattTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGAtctc                                                                                                               <  1:90665/275‑1 (MQ=255)                                                                                                                                                                       tattatTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAAttcttc                                                                                                                                                                                                                                                <  8:198550/136‑1 (MQ=255)                                                                                                                                                                       tattatTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAAttcttc                                                                                                                                                                                                                                                >  7:198550/1‑136 (MQ=255)                                                                                                                                                                                       aTAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCa                                                                                                                                                                                                                   >  8:307814/1‑149 (MQ=255)                                                                                                                                                                                       aTAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCa                                                                                                                                                                                                                   >  5:131359/1‑149 (MQ=255)                                                                                                                                                                                        tAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCAt                                                                                                                                                                                                                  >  8:120643/1‑149 (MQ=255)                                                                                                                                                                                          aaCCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGt                                                                                                                                                                                                                <  7:307814/149‑1 (MQ=255)                                                                                                                                                                                           aCCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGtt                                                                                                                                                                                                               >  7:44743/1‑149 (MQ=255)                                                                                                                                                                                             cAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTg                                                                                                                                                                                                             <  8:373385/149‑1 (MQ=255)                                                                                                                                                                                              aGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGa                                                                                                                                                                                                            >  8:112580/1‑149 (MQ=255)                                                                                                                                                                                               gCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCACTGACGTATCATTTGGTCCGCCCGAAACa                                                                             >  4:390044/1‑275 (MQ=255)                                                                                                                                                                                                  ggATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCg                                                                                                                                                                                                        >  6:256331/1‑149 (MQ=255)                                                                                                                                                                                                         tttttCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGAtttttt                                                                                                                                                                                                 >  6:205392/1‑149 (MQ=255)                                                                                                                                                                                                         tttttCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGAtttttt                                                                                                                                                                                                 >  6:326137/1‑149 (MQ=255)                                                                                                                                                                                                          ttttCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTtc                                                                                                                                                                                                <  5:306124/149‑1 (MQ=255)                                                                                                                                                                                                          ttttCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTtc                                                                                                                                                                                                <  5:6252/149‑1 (MQ=255)                                                                                                                                                                                                          ttttCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTtc                                                                                                                                                                                                <  7:112580/149‑1 (MQ=255)                                                                                                                                                                                                              ccAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTccc                                                                                                                                                                                            <  6:34724/149‑1 (MQ=255)                                                                                                                                                                                                               cAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCg                                                                                                                                                                                           <  5:326137/149‑1 (MQ=255)                                                                                                                                                                                                                aGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGt                                                                                                                                                                                          <  7:476324/149‑1 (MQ=255)                                                                                                                                                                                                                aGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGt                                                                                                                                                                                          >  6:129719/1‑149 (MQ=255)                                                                                                                                                                                                                       cacCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGaa                                                                                                                                                                              >  4:149677/1‑154 (MQ=255)                                                                                                                                                                                                                       cacCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCACTGACGTATCATTTGGTCCGCCCGAAACAGGTTGGCCAGCGTGAATAACATCg                                                     >  2:435131/1‑275 (MQ=255)                                                                                                                                                                                                                        acCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCt                                                                                                                                                                                  <  8:44743/149‑1 (MQ=255)                                                                                                                                                                                                                         ccATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTGCCGTAAATGCCtt                                                                                                                                                                                 <  5:205392/149‑1 (MQ=255)                                                                                                                                                                                                                         ccATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCtt                                                                                                                                                                                 <  6:401055/149‑1 (MQ=255)                                                                                                                                                                                                                           aTCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAg                                                                                                                                                                          <  2:149678/154‑1 (MQ=255)                                                                                                                                                                                                                            tCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCACTg                                                                                                       >  2:157049/1‑220 (MQ=255)                                                                                                                                                                                                                             cAACTGTCTCACCATACTCAACAGATGATGCCCATGCTTTGGGTGAATATATGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTTATTTCATTTCTTCATGTTTGAGCAGATTATTTCTCCCGTAAATGCCTTGAAt                                                                                                                                                                             >  8:232733/1‑149 (MQ=18)                                                                                                                                                                                                                              aaCTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATc                                                                                                                                                                            >  8:89126/1‑149 (MQ=255)                                                                                                                                                                                                                               aCTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCa                                                                                                                                                                           <  7:89126/149‑1 (MQ=255)                                                                                                                                                                                                                                cTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTTAATATAAGGAAGGTTCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCACTGACGt                                                                                                   <  4:157047/220‑1 (MQ=255)                                                                                                                                                                                                                                   tctcACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCAt                                                                                                                                                                                                                  >  5:458505/1‑106 (MQ=255)                                                                                                                                                                                                                                   tctcACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCAt                                                                                                                                                                                                                  <  6:458505/106‑1 (MQ=255)                                                                                                                                                                                                                                   tctcACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCt                                                                                                                                                                       >  5:276286/1‑149 (MQ=255)                                                                                                                                                                                                                                     tcACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTAt                                                                                                                                                                     >  8:189662/1‑149 (MQ=255)                                                                                                                                                                                                                                        ccATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTa                                                                                                                                                                  <  8:252930/149‑1 (MQ=255)                                                                                                                                                                                                                                        ccATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTa                                                                                                                                                                  >  7:135455/1‑149 (MQ=255)                                                                                                                                                                                                                                             ttaaccactggATTCCAATGTGATGGGTGAATATAAAGAAGGTGCGAATAAGCGTGGAAATTCTTCTCGGCTGACTCGGTCATTTCATTTCTTCATGTTTGAACCGGTTTGTTCTCCCGTGAATGCCTTGAATCAGCCTATTTAGGCCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGATTTTAGTGAGATCTCTCCCACTGAAGTATCATTTGGTCCGCCCGAAACAGGTTGGCCAGCGTGAATAACATCGCCAGTTGGTTATCGTTTTTCAg                               <  4:117799/265‑1 (MQ=11)                                                                                                                                                                                                                                              tCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGt                                                                                                                                                            >  6:418372/1‑149 (MQ=255)                                                                                                                                                                                                                                              tCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGt                                                                                                                                                            >  8:323243/1‑149 (MQ=255)                                                                                                                                                                                                                                               cAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCACTGACGTATCATTTGGTCCGCCCGAAACAGGTTGGCCAGCGTGAATAACATCGCCAGTTGGTTATCGTTTTTCAGCa                             <  2:425890/275‑1 (MQ=255)                                                                                                                                                                                                                                                aaCAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGc                                                                                                                                                                                                          >  5:261071/1‑101 (MQ=255)                                                                                                                                                                                                                                                aaCAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGc                                                                                                                                                                                                          <  6:261071/101‑1 (MQ=255)                                                                                                                                                                                                                                                aaCAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGttt                                                                                                                                                          >  8:291734/1‑149 (MQ=255)                                                                                                                                                                                                                                                aaCAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTttc                                                                                                                                                         >  1:45887/1‑150 (MQ=255)                                                                                                                                                                                                                                                    gTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGttt                                                                                                                                                                                                              <  5:81676/93‑1 (MQ=255)                                                                                                                                                                                                                                                    gTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGttt                                                                                                                                                                                                              >  6:81676/1‑93 (MQ=255)                                                                                                                                                                                                                                                    gTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCg                                                                                                                                                     <  3:45887/150‑1 (MQ=255)                                                                                                                                                                                                                                                     ttGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCa                                                                                                                                                  >  2:541350/1‑152 (MQ=255)                                                                                                                                                                                                                                                     ttGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCACTGACGTATCATTTGGTCCGCCCGAAACAGGTTGGCCAGCGTGAATAACATCGCCAGTTGGTTATCGTTTTTTAGCAACCCCt                       >  1:14218/1‑275 (MQ=255)                                                                                                                                                                                                                                                     ttGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTACCACTGACGTATCATTTGGTCCGCCCGAAACGGGTTGGCCAGCGTGAATAACATCGCCCGTTGGTTATCGTTTTTCGGAAAcccc                        >  4:304234/1‑274 (MQ=255)                                                                                                                                                                                                                                                       gATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGcc                                                                                                                                                   >  6:195683/1‑149 (MQ=255)                                                                                                                                                                                                                                                         tGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGTGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTa                                                                                                                                              <  4:541345/152‑1 (MQ=255)                                                                                                                                                                                                                                                          gCCCATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCt                                                                                                                                                                                  >  3:59109/1‑115 (MQ=255)                                                                                                                                                                                                                                                             aaTGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTaa                                                                                                                                             <  5:129719/149‑1 (MQ=255)                                                                                                                                                                                                                                                              aTGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGaa                                                                                                                                                                              <  1:59109/115‑1 (MQ=255)                                                                                                                                                                                                                                                                gCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGc                                                                                                                                          <  5:256331/149‑1 (MQ=40)                                                                                                                                                                                                                                                                  tATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCt                                                                                                                                                                                                                                          >  4:244839/1‑51 (MQ=255)                                                                                                                                                                                                                                                                  tATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGt                                                                                                                                        >  6:308630/1‑149 (MQ=39)                                                                                                                                                                                                                                                                   aTGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCACTGACGTATCATTTGGTCCGCCCGAAACAGGTTGGCCAGCGTGAATAACATCGCCAGTTGGTTATCGTTTTTCAGCAACCCCTTGTATCTGGc             >  2:379071/1‑271 (MQ=34)                                                                                                                                                                                                                                                                     gggTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCACTGACGTATCATTTGGTCCGCCCGAAACAGGtt                                                                         >  2:248205/1‑209 (MQ=34)                                                                                                                                                                                                                                                                      ggTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACt                                                                                                                                                                                                                                      <  2:244843/51‑1 (MQ=255)                                                                                                                                                                                                                                                                         gAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCACTGACGTATCATTTGGTCCGCCCGAAACAGGTTGGcc                                                                     <  4:248201/209‑1 (MQ=32)                                                                                                                                                                                                                                                                          aaTATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTATTGAGATCTCTCCCACTGACGTATCATTTGTTCCGCCCGAAACAGGTTGGCCAGCTTGAATACAATCGCAAGTTGGTTAACGTTTTTCAGCAACCCCTTTTATCTGGCTTTCACGAAGc  >  4:377672/1‑275 (MQ=255)                                                                                                                                                                                                                                                                          aaTATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCACTGACGTATCATTTGGTCCGCCCGAAACAGGTTGGCCAGCGTGAATAACATCGCCAGTTGGTTATCGTTTTTCAGCAACCCCTTGTATCTGGCTTTCACGAAg   <  2:377677/274‑1 (MQ=31)                                                                                                                                                                                                                                                                           atatAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTgaga                                                                                                                   >  2:216537/1‑161 (MQ=31)                                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                                            GTCTGGTTTGCCGCGCCACCCATTTGAACAATTTGATAACCGTCATTTAATACTGAGTATTCTGCGGCGCCATTCATTTCAATGTACTGATACCCTCCGTTAATTTCAGTATTATTACTTACACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCACTGACGTATCATTTGGTCCGCCCGAAACAGGTTGGCCAGCGTGAATAACATCGCCAGTTGGTTATCGTTTTTCAGCAACCCCTTGTATCTGGCTTTCACGAAGC  >  NC_000913/2288647‑2289186

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