Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
JC JC	NC_000913	3,668,878	IS2 (+) +5 bp	coding (307‑311/1398 nt)	yhjA ←	putative cytochrome C peroxidase

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	2996361 =	NA (NA)	45 (0.820)	37/294	0.9	100%	noncoding (1331/1331 nt)	IS2	repeat region
	?	NC_000913	3668878 =	0 (0.000)					coding (311/1398 nt)	yhjA	putative cytochrome C peroxidase
*	?	NC_000913	3186096 =	NA (NA)	51 (0.920)	45/298	0.4	100%	noncoding (1/1331 nt)	IS2	repeat region
	?	NC_000913	= 3668882	0 (0.000)					coding (307/1398 nt)	yhjA	putative cytochrome C peroxidase

CAAGATTGCTTACGAGTAATCGCATTTGGCCGGGTTCGCCCGGCCACATCATTAACGGATTAATGATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGAT  >  NC_000913/3185947‑3186242                                                                                                                                                      |                                                                                                                                                   cAAGATTGCTTACGAGTAATCGCATTTGGCCGGGTTCGCCCGGCCACATCATTAACGGATTAATGATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGttatta                                                                                                                                                   <  1:521569/151‑1 (MQ=255)  aaGATTGCTTACGAGTAATCGCATTTGGCCGGGTTCGCCCGGCCACATCATTAACGGATTAATGATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGttatt                                                                                                                                                    >  2:492223/1‑149 (MQ=255)              gAGTAATCGCATTTGGCCGGGTTCGCCCGGCCACATCATTAACGGATTAATGATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTg                                                                                                                                      >  2:432011/1‑151 (MQ=255)                   aTCGCATTTTGCCGGGTTCGCCCTGCCACATCATTAACGGATTAATGATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGACTTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAAtct                                                                                                                                 >  1:352162/1‑151 (MQ=255)                       cATTTGGCCGGGTTCGCCCGGCCACATCATTAACGGATTAATGATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAg                                                                                                                             <  2:71982/151‑1 (MQ=255)                        atgtGGCCGAGTTCGCCCGACCACATCATAAACGGATTAATGATAAGTGGATCAGATGTATAAAAAATGAAAATTAACCGCAATAAGCGTATTAATTAAAAATATTTATTGTTCATTATCCTTTCTTAGACTGGCTCCCTGAATCTCCAg                                                                                                                             <  1:413693/147‑1 (MQ=255)                                gggTTCGCCCGGCCACATCATTAACGGATTAATGATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAat                                                                                                                    >  2:75268/1‑151 (MQ=255)                                   ttCGCCCGGCCACATCATTAACGGATTAATGATAAGTGGATCAGATGAATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAatat                                                                                                                  >  1:28560/1‑150 (MQ=255)                                    tCGCCCGGCCACATCATTAACGGATTAATGATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATc                                                                                                                 >  2:257546/1‑150 (MQ=255)                                       cccGGCCACATCATTAACGGATTAATGATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACtt                                                                                                             <  2:642016/151‑1 (MQ=255)                                        ccGGCCACATCATTAACGGATTAATGATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCATTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTa                                                                                                            <  1:153821/151‑1 (MQ=255)                                         cGGCCACATCATTAACGGATTAATGATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTa                                                                                                            >  2:295813/1‑150 (MQ=255)                                                atcatTAACGGATTCATGATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCCAGTTGATTAAACATATTTATCGTTCAGCATCCGGTATTAGACTGGCCCCCTGAATCTCCAGTCATCCTATATCGCTGaagtaagtt                                                                                                    <  1:169468/151‑2 (MQ=255)                                                  catTAACGGATTAATGATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGAt                                                                                                  >  1:823112/1‑151 (MQ=255)                                                    ttAACGGATTAATGATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGAt                                                                                                  >  1:666174/1‑149 (MQ=255)                                                     tAACGGATTAATGATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGt                                                                                               >  2:414488/1‑151 (MQ=255)                                                         ggATTAATGATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTa                                                                                           <  1:842123/151‑1 (MQ=255)                                                          gATTAATGATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTaa                                                                                          <  2:112150/151‑1 (MQ=255)                                                           aTTAATGATAAGTGGATCAGATGTATAAAAAATTAAAATTAAACACAATAAGAGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGACCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAAt                                                                                         >  1:372363/1‑151 (MQ=255)                                                            ttAATGATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAAt                                                                                         >  1:31079/1‑150 (MQ=255)                                                              aaTGATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACt                                                                                      <  2:729096/151‑1 (MQ=255)                                                               aTGATAAGTGGATCAGATGTATATAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAGATATTTATTGTTCATTATCCGTTATTAGACTGGTCCCCTGAATCTCAACACAACCAACATCTCGTAAATAAGTGATACTCTTAATACTa                                                                                     <  1:756929/151‑1 (MQ=255)                                                               aTGATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACt                                                                                      >  1:199071/1‑150 (MQ=255)                                                                 gATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAg                                                                                    <  1:828574/150‑1 (MQ=255)                                                                  aTAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGGCTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGtt                                                                                  >  2:440154/1‑151 (MQ=255)                                                                   tAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGt                                                                                               <  1:87010/137‑1 (MQ=255)                                                                   tAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGt                                                                                               >  2:87042/1‑137 (MQ=255)                                                                    aaGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGtttt                                                                                >  1:444843/1‑151 (MQ=255)                                                                        ggATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAg                                                                             <  1:837570/150‑1 (MQ=255)                                                                        ggATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTAGTAGACTGGCCCCCTGAATCTCCAGACATCCCATATCACTTAAATAAGTGATAGTCTTAATCCTAGTTTTTTGa                                                                            <  1:133302/151‑1 (MQ=255)                                                                           tCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACt                                                                          >  1:550576/1‑150 (MQ=255)                                                                            cAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTa                                                                         <  2:458799/150‑1 (MQ=255)                                                                            cAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAg                                                                        <  2:760067/151‑1 (MQ=255)                                                                             aGATGTATAAAAAATTAAAATTATCCACAATAATCGAATAGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCGGAATTTCCAGAGAACTAATATCACTAAAATAAGTGATAGTCTTAATACTAGTTGTTAGGCTAGt                                                                       <  1:669822/151‑1 (MQ=25)                                                                              gATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTc                                                                      <  2:25160/151‑1 (MQ=255)                                                                               aTGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTc                                                                      >  2:867614/1‑150 (MQ=255)                                                                               aTGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTc                                                                      >  2:105449/1‑150 (MQ=255)                                                                                tGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCAt                                                                    >  1:35486/1‑151 (MQ=255)                                                                                 gTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCAtt                                                                   >  1:391658/1‑151 (MQ=255)                                                                                  tataAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTg                                                                  >  1:390331/1‑151 (MQ=255)                                                                                   ataAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTgg                                                                 <  2:220851/151‑1 (MQ=255)                                                                                   ataAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTgg                                                                 >  2:515665/1‑151 (MQ=255)                                                                                   ataAAAAATTAAAAGTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTg                                                                  <  2:625617/150‑1 (MQ=255)                                                                                      aaaaaTTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTgg                                                                 >  1:236113/1‑148 (MQ=255)                                                                                      aaaaaTTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTgg                                                                 >  2:68861/1‑148 (MQ=255)                                                                                       aaaattaaaattAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGgaga                                                              >  2:284023/1‑150 (MQ=255)                                                                                        aaattaaaattaACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAAc                                                            >  1:398024/1‑151 (MQ=255)                                                                                         aattaaaattaaCCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACa                                                           <  1:105417/151‑1 (MQ=255)                                                                                           ttaaaattaaCCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAAc                                                            >  1:92282/1‑148 (MQ=255)                                                                                             aaaattaaCCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAgatg                                                       >  1:286443/1‑151 (MQ=255)                                                                                              aaattaaCCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAgatg                                                       >  1:573952/1‑150 (MQ=255)                                                                                                attaaCCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGAtgatt                                                    >  2:745486/1‑151 (MQ=255)                                                                                                attaaCCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGAtgatt                                                    <  2:343111/151‑1 (MQ=255)                                                                                                   aaCCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGATAACAGAtgatt                                                    >  1:524094/1‑148 (MQ=255)                                                                                                    aCCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATg                                                <  2:823209/151‑1 (MQ=255)                                                                                                    aCCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATg                                                >  2:16351/1‑151 (MQ=255)                                                                                                    aCCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATg                                                >  2:536547/1‑151 (MQ=255)                                                                                                     ccACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGt                                               <  2:539311/151‑1 (MQ=255)                                                                                                     ccACAATAAGAGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGt                                               >  1:91427/1‑151 (MQ=255)                                                                                                      cacaATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTc                                              >  2:308122/1‑151 (MQ=255)                                                                                                            aaGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAggg                                        <  1:560713/151‑1 (MQ=255)                                                                                                              gCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGcc                                      <  2:146121/151‑1 (MQ=255)                                                                                                                gAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCgg                                    >  1:637710/1‑151 (MQ=255)                                                                                                                  aTTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCgg                                    >  2:127626/1‑149 (MQ=255)                                                                                                                  aTTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGgag                                  <  1:79866/151‑1 (MQ=255)                                                                                                                  aTTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCACCCTGAATCTCCAGACAACCACTATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCg                                     >  1:414558/1‑148 (MQ=255)                                                                                                                    tGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGgaga                                 <  2:417953/150‑1 (MQ=255)                                                                                                                         aaaaaTATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGc                            >  1:148797/1‑150 (MQ=255)                                                                                                                         aaaaaTATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGc                            <  2:463938/150‑1 (MQ=255)                                                                                                                         aaaaaTATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGc                            >  2:523492/1‑150 (MQ=255)                                                                                                                           aaaTATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGc                            >  1:140455/1‑148 (MQ=255)                                                                                                                           aaaTATTTATTGTTCATTATCCGTTATTAGACTGACCCCCTACAACTCCAGACAACCAATATCACTAACATAAGTGATAGTCTTACTACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGCACCGc                            >  2:844788/1‑148 (MQ=17)                                                                                                                             atatTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACg                       >  2:136706/1‑151 (MQ=255)                                                                                                                             atatTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACg                       <  1:790067/151‑1 (MQ=255)                                                                                                                               atttattGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGc                     <  2:743955/151‑1 (MQ=255)                                                                                                                                  tattGTTCATTATCCGTTATTAGACTGGCCCCCTGTATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGTGAACTGATGATTGATGTCTTAGGTCCGGAGAAACGCAGACGGc                     >  1:546439/1‑148 (MQ=255)                                                                                                                                        tCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACg                       >  1:824990/1‑140 (MQ=40)                                                                                                                                        tCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACg                       <  2:825087/140‑1 (MQ=40)                                                                                                                                          aTTATCCGTTATTAGACTGGCCCACTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAAGTCCGGAGAAACGCAGACGGCGTACcac              >  1:437126/1‑147 (MQ=255)                                                                                                                                            tATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAgg         <  1:599015/150‑1 (MQ=35)                                                                                                                                                cGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGaaaa     >  1:862108/1‑150 (MQ=33)                                                                                                                                                   tattaGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGAt  <  2:273224/150‑1 (MQ=31)                                                                                                                                                      |                                                                                                                                                   CAAGATTGCTTACGAGTAATCGCATTTGGCCGGGTTCGCCCGGCCACATCATTAACGGATTAATGATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGAT  >  NC_000913/3185947‑3186242

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