Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
JC JC	NC_000913	564,107	IS2 (+) +5 bp	100%	coding (370‑374/633 nt)	sfmZ ←	response regulator family protein

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	381260 =	NA (NA)	60 (0.890)	50/260	NT	100%	noncoding (1/1331 nt)	IS2	repeat region
	?	NC_000913	= 564111	0 (0.000)					coding (370/633 nt)	sfmZ	response regulator family protein
*	?	NC_000913	= 382590	NA (NA)	54 (0.800)	50/262	NT	100%	noncoding (1331/1331 nt)	IS2	repeat region
	?	NC_000913	564107 =	0 (0.000)					coding (374/633 nt)	sfmZ	response regulator family protein

CCTGGTTTGCAACTACCGAGTCGGCAATAAAACATTGCGGACCAATCGTTGCTTCCGCTTCAATAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAA  >  NC_000913/381113‑381401                                                                                                                                                    |                                                                                                                                              ccTGGTTTGCAACTACCGAGTCGGCAATAAAACATTGCGGACCAATCGTTGCTTCCGCTTCAATAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTa                                                                                                                                              <  5:346524/149‑1 (MQ=255)   tGGTTTGCAACTACCGAGTCGGCAATAAAACATTGCGGACCAATCGTTGCTTCCGCTTCAATAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGa                                                                                                                                            >  6:207214/1‑149 (MQ=255)     gTTTGCAACTACCGAGTCGGCAATAAAACATTGCGGACCAATCGTTGCTTCCGCTTCAATAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGa                                                                                                                                            <  4:206451/147‑1 (MQ=255)     gTTTGCAACTACCGAGTCGGCAATAAAACATTGCGGACCAATCGTTGCTTCCGCTTCAATAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGAc                                                                                                                                           <  2:272829/148‑1 (MQ=255)         gCAACTACCGAGTCGGCAATAAAACATTGCGGACCAATCGTTGCTTCCGCTTCAATAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGcc                                                                                                                                      >  1:31316/1‑149 (MQ=255)         gCAACTACCGAGTCGGCAATAAAACATTGCGGACCAATCGTTGCTTCCGCTTCAATAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGc                                                                                                                                       <  2:392051/148‑1 (MQ=255)           aaCTACCGAGTCGGCAATAAAACATTGCGGACCAATCGTTGCTTCCGCTTCAATAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTa                                                                                                                                              <  1:316358/139‑1 (MQ=255)           aaCTACCGAGTCGGCAATAAAACATTGCGGACCAATCGTTGCTTCCGCTTCAATAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTa                                                                                                                                              >  2:316358/1‑139 (MQ=255)              tACCGAGTCGGCAATAAAACATTGCGGACCAATCGTTGCTTCCGCTTCAATAACCGCATTTTTAATCTCGGTGTCAAAACCAAATTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAAGTACTAGACTGGCCCCCt                                                                                                                                  >  4:265408/1‑148 (MQ=255)               aCCGAGTCGGCAATAAAACATTGCGGACCAATCGTTGCTTCCGCTTCAATAACCGCATTTTTAATCTCGGTGGCAAAACCAAGTTTTACGCCATTGCTGATTATTGTGCCAGGCCGAATAATCGCATAATTACTAGACTGGCCCCCTg                                                                                                                                 >  2:330068/1‑148 (MQ=255)                ccGAGTCGGCAATAAAACATTGCGGACCAATCATTTCTTCCGCTTCAATAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCt                                                                                                                                  >  4:55523/1‑146 (MQ=255)                ccGAGTCGGCAATAAAACATTGCGGACCAATCATTGCTTCCGCTTCAATAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCt                                                                                                                                  <  3:55523/146‑1 (MQ=255)                        gCAATAAAACATTGCGGACCAATCGTTGCTTCCGCTTCAATAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAAATACTAGACTGGCCCCCTGAATCTCCAGa                                                                                                                       >  2:138975/1‑149 (MQ=255)                         cAATAAAACATTGCGGACCAATCGTTGCTTCCGCTTCAATAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGAc                                                                                                                      <  4:13118/149‑1 (MQ=255)                          aaTAAAACATTGCGGACCAATCGTTGCTTCCGCTTCAATAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACa                                                                                                                     >  8:181105/1‑149 (MQ=255)                           aTAAAACATTGCGGACCAATCGTTGCTTCCGCTTCAATAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGTCaa                                                                                                                    <  2:317233/149‑1 (MQ=255)                             aaaaCATTGCGGACCAATCGTTGCTTCCGCTTCAATAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGa                                                                                                                                            >  4:116492/1‑123 (MQ=255)                             aaaaCATTGCGGACCAATCGTTGCTTCCGCTTCAATAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGa                                                                                                                                            <  3:116492/123‑1 (MQ=255)                                aCATTGCGGACCAATCGTTGCTTCCGCTTCAATAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAat                                                                                                               <  8:195661/149‑1 (MQ=255)                                  aTTGCGGACCAATCGTTGCTTCCGCTTCAATAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAatat                                                                                                             <  5:167508/149‑1 (MQ=255)                                  aTTGCGGACCAATCGTTGCTTCCGCTTCAATAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAatat                                                                                                             <  5:240521/149‑1 (MQ=255)                                      cGGACCAATCGTTGCTTCCGCTTCAATAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACt                                                                                                         <  7:348137/149‑1 (MQ=255)                                       ggACCAATCGTTGCTTCCGCTTCAATAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACtt                                                                                                        <  4:143349/149‑1 (MQ=255)                                           cAATCGTTGCTTCCGCTTCAATAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAAt                                                                                                    <  6:326160/149‑1 (MQ=255)                                             aTCGTTGCTTCCGCTTCAATAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAAt                                                                                                    >  1:268284/1‑147 (MQ=255)                                                 ttGCTTCCGCTTCAATAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGa                                                                                              >  4:127058/1‑149 (MQ=255)                                                    cTTCCGCTTCAATAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATa                                                                                            <  5:326258/148‑1 (MQ=255)                                                        cGCTTCAATAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCtt                                                                                       <  8:130137/149‑1 (MQ=255)                                                         gCTTCAATAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCtt                                                                                       >  6:43892/1‑148 (MQ=255)                                                          cTTCAATAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTaa                                                                                     <  5:207214/149‑1 (MQ=255)                                                          cTTCAATAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTaa                                                                                     <  7:244766/149‑1 (MQ=255)                                                             cAATAACCTCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATa                                                                                   <  8:202058/148‑1 (MQ=255)                                                             cAATAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTGAATAc                                                                                  >  5:358395/1‑149 (MQ=255)                                                               aTAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTa                                                                                <  3:37279/149‑1 (MQ=255)                                                                    cGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGtttt                                                                           >  8:272310/1‑149 (MQ=255)                                                                     gCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGttttt                                                                          <  6:383427/149‑1 (MQ=255)                                                                     gCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGtttt                                                                           >  5:383427/1‑148 (MQ=255)                                                                     gCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGtttt                                                                           >  4:187991/1‑148 (MQ=255)                                                                     gCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGtttt                                                                           >  7:6930/1‑148 (MQ=255)                                                                       aTTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCtt                                                                                       >  2:35164/1‑134 (MQ=255)                                                                       aTTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCtt                                                                                       <  1:35164/134‑1 (MQ=255)                                                                       aTTTTTAATCTAGGTGTCAAAACCAATTTTTACTCCATTGCTGATTATTGTGCAAGGACTAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTa                                                                         <  1:313074/148‑1 (MQ=25)                                                                        tttttAATCTCGGTTGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTAGTAGAGTGGCCCCCTGAATCTCCCGACAACCATTATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGa                                                                       <  1:138975/149‑1 (MQ=37)                                                                        tttttAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATa                                                                                   >  8:222327/1‑137 (MQ=255)                                                                        tttttAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATa                                                                                   <  7:222327/137‑1 (MQ=255)                                                                          tttAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGAc                                                                      <  7:287261/148‑1 (MQ=255)                                                                                   ggTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCAt                                                               <  2:268284/146‑1 (MQ=255)                                                                                       gccaaaCCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAAAAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAAAATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGaa                                                        <  2:59904/146‑1 (MQ=255)                                                                                         aaaaCCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGaa                                                        >  2:315427/1‑147 (MQ=255)                                                                                         aaaaCCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAAATACTAGAATTTCCCCCTGAATCTCCAGAAAACCAATATCACTTAAATAAGTGATAGTCTTTATACTACTTTTTAGACTAGACATTGGAGaa                                                        >  3:51809/1‑147 (MQ=21)                                                                                          aaaCCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTg                                                                                               <  4:77696/107‑1 (MQ=255)                                                                                          aaaCCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTg                                                                                               >  3:77696/1‑107 (MQ=255)                                                                                          aaaCCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTAATAGACTGGCACCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGaa                                                        >  1:133764/1‑146 (MQ=255)                                                                                              cAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCGCCAGCCAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAgatga                                                 >  2:186001/1‑149 (MQ=255)                                                                                               aaTTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTTAATCTGCAGAAAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAgatgat                                                <  3:265408/149‑1 (MQ=255)                                                                                                aTTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGAtgatt                                               <  8:151167/149‑1 (MQ=255)                                                                                                         ccATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTa                                      <  3:187991/149‑1 (MQ=255)                                                                                                          cATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAg                                     <  8:6930/149‑1 (MQ=255)                                                                                                          cATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAg                                     <  5:233608/149‑1 (MQ=255)                                                                                                            ttGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTa                                                                                      <  4:314453/98‑1 (MQ=255)                                                                                                            ttGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTa                                                                                      >  3:314453/1‑98 (MQ=255)                                                                                                            ttGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAggg                                   >  2:266329/1‑149 (MQ=255)                                                                                                                          gCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACg                        >  1:172693/1‑146 (MQ=255)                                                                                                                               gACGAATAAACGCATACTTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGc                >  7:347834/1‑149 (MQ=255)                                                                                                                                   aaTAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTAc            <  5:43892/149‑1 (MQ=255)                                                                                                                                    aTAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTAc            >  3:254241/1‑148 (MQ=40)                                                                                                                                       aaCGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACcaca        >  2:308871/1‑149 (MQ=38)                                                                                                                                        aCGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGTCCGGAGAAACGCAGACGGCGTACcaca        >  7:303951/1‑148 (MQ=38)                                                                                                                                           cATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAg     >  7:318494/1‑148 (MQ=35)                                                                                                                                              aaTTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGaa  >  1:363635/1‑148 (MQ=32)                                                                                                                                                    |                                                                                                                                              CCTGGTTTGCAACTACCGAGTCGGCAATAAAACATTGCGGACCAATCGTTGCTTCCGCTTCAATAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAA  >  NC_000913/381113‑381401

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.