Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
JC JC	NC_000913	2,967,576	IS5 (+) +4 bp	coding (856‑859/2247 nt)	ptsP ←	PEP‑protein phosphotransferase enzyme I; GAF domain containing protein

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	2967576 =	0 (0.000)	84 (1.260)	69/296	0.0	100%	coding (859/2247 nt)	ptsP	PEP‑protein phosphotransferase enzyme I; GAF domain containing protein
	?	NC_000913	= 3131340	NA (NA)					noncoding (1195/1195 nt)	IS5	repeat region
*	?	NC_000913	= 2967579	0 (0.000)	63 (0.940)	57/298	0.3	100%	coding (856/2247 nt)	ptsP	PEP‑protein phosphotransferase enzyme I; GAF domain containing protein
	?	NC_000913	3130146 =	NA (NA)					noncoding (1/1195 nt)	IS5	repeat region

CCACGCCTGCCAGATGATATGGTCCGTCCCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGATCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGAAAAGAGATTTTCTTGTC  >  NC_000913/3129997‑3130290                                                                                                                                                      |                                                                                                                                                 ccacGCCTGCCAGATGATATGGTCCGTCCCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGATCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTaagg                                                                                                                                                 >  1:1012860/1‑151 (MQ=255) ccacGCCTGCCAGATGATATGGTCCGTCCCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGATCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTaagg                                                                                                                                                 <  2:704885/151‑1 (MQ=255)    cGCCTGCCAGATGATATGGTCCGTCCCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGATCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTaaggaag                                                                                                                                              <  2:726991/151‑1 (MQ=255)    cGCCTGCCAGATGATATGGTCCGTCCCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGATCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTaaggaag                                                                                                                                              <  2:1014006/151‑1 (MQ=255)    cGCCTGCCAGATGATATGGTCCGTCCCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGATCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTaagga                                                                                                                                                >  2:1043418/1‑149 (MQ=255)      ccTGCCAGATGATATGGTCCGTCCCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGATCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGt                                                                                                                                            <  1:171599/151‑1 (MQ=255)       cTGCCAGATGATATGGTCCGTCCCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGATCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTg                                                                                                                                           >  2:614217/1‑151 (MQ=255)         gCCAGATGATATGGTCCGTCCCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGATCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGc                                                                                                                                          <  2:665777/150‑1 (MQ=255)             gatgatATGGTCCGTCCCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGATCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAAc                                                                                                                                      >  1:196415/1‑150 (MQ=255)                   atGGTCCGTCCCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGATCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAGTcc                                                                                                                                <  2:269892/150‑1 (MQ=255)                   atGGTCCCTCCCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGACTAAAGGGGGTATCCAGTTTGGATCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAGTccc                                                                                                                               >  2:647129/1‑151 (MQ=255)                     ggTCCGTCCCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGATCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAGTCCCTg                                                                                                                             >  2:541744/1‑151 (MQ=255)                      gtccgtccCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGATCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAGTCCCTGa                                                                                                                            <  2:486720/151‑1 (MQ=255)                        ccgtccCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGATCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAGTCCCTGata                                                                                                                          <  2:444815/151‑1 (MQ=255)                         cgtccCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGATCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAGTCCCTGatat                                                                                                                         <  2:716500/151‑1 (MQ=255)                           tccCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGATCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAGTCCCTGATATg                                                                                                                        >  1:60941/1‑150 (MQ=255)                               cGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGATCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAGTCCCTGATATg                                                                                                                        >  1:562038/1‑146 (MQ=255)                                ggATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGATCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAGTCCCTGATATGAGAt                                                                                                                    >  1:221326/1‑149 (MQ=255)                                                      cGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGATCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCAt                                                                                              <  1:710424/149‑1 (MQ=255)                                                       gTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGATCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCAt                                                                                              >  1:862275/1‑148 (MQ=255)                                                         tGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGATCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGaa                                                                                          <  1:967461/150‑1 (MQ=255)                                                          gCGTTCAATGAATAAAGGGGGTATCCAGTTTGGATCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAAc                                                                                         >  1:986942/1‑150 (MQ=255)                                                                aaTGAATAAAGGGGGTATCCAGTTTGGATCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTtc                                                                                  <  1:69275/151‑1 (MQ=255)                                                                  tgcaTAAAGGGGGTATCCAGTTTGGATCGTCGGGATAAAGTGCTGACGGAAAAGAGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAGTACCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTtcat                                                                                >  1:56781/4‑151 (MQ=255)                                                                   gAATAAAGGGGGTATCCAGTTTGGATCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTtcat                                                                                <  1:500843/150‑1 (MQ=255)                                                                    aaTAAAGGGGGTATCCAGTTTGGATCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTtcatca                                                                              <  1:984818/151‑1 (MQ=255)                                                                           ggggTATCCAGTTTGGATCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGtca                                                                       <  1:614110/151‑1 (MQ=255)                                                                                 tCCAGTTTGGATCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACt                                                                 >  1:143856/1‑151 (MQ=255)                                                                                  ccAGTTTGGATCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACtt                                                                >  2:70294/1‑151 (MQ=255)                                                                                     gTTTGGATCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTAcc                                                             >  2:187101/1‑151 (MQ=255)                                                                                       ttGGATCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCtt                                                           >  2:1034187/1‑151 (MQ=255)                                                                                         ggATCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTc                                                          <  2:551510/150‑1 (MQ=255)                                                                                             cgtcgGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGGGTCATCAACTTACCTTCGCCg                                                      <  1:643811/150‑1 (MQ=255)                                                                                             cgtcgGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCg                                                      >  1:434706/1‑150 (MQ=255)                                                                                                cgGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGAc                                                    >  2:471111/1‑149 (MQ=255)                                                                                                cgGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACa                                                   <  2:21380/150‑1 (MQ=255)                                                                                                 gggATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGt                                                 >  2:619089/1‑151 (MQ=255)                                                                                                 gggATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGt                                                 >  1:113474/1‑151 (MQ=255)                                                                                                   gATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTg                                                >  1:544422/1‑150 (MQ=255)                                                                                                   gATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGa                                               <  1:1037170/151‑1 (MQ=255)                                                                                                       aaGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTc                                           >  2:718825/1‑151 (MQ=255)                                                                                                         gTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTc                                           >  1:266554/1‑149 (MQ=255)                                                                                                              gaCGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGt                                     >  1:710728/1‑150 (MQ=255)                                                                                                               aCGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTa                                    >  1:113697/1‑150 (MQ=255)                                                                                                                 ggAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAgc                                 <  1:649870/151‑1 (MQ=255)                                                                                                                  gAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAgcg                                <  1:672594/151‑1 (MQ=255)                                                                                                                  gAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAgcg                                <  1:407993/151‑1 (MQ=255)                                                                                                                     aaGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCg                             <  2:886055/151‑1 (MQ=255)                                                                                                                       gCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTc                           >  2:468783/1‑151 (MQ=255)                                                                                                                        cGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTc                           >  1:805161/1‑150 (MQ=255)                                                                                                                        cGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCa                          <  2:480086/151‑1 (MQ=255)                                                                                                                        cGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCa                          <  1:618980/151‑1 (MQ=255)                                                                                                                        cGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCa                          <  1:563761/151‑1 (MQ=255)                                                                                                                          aTAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTc                           >  2:524076/1‑148 (MQ=255)                                                                                                                                gCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGa                   <  2:380695/150‑1 (MQ=255)                                                                                                                                    tAAGGTCATTTTTATTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGaaa                 >  1:150351/1‑148 (MQ=39)                                                                                                                                      aGGTCATTTTTATTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGAAAag               >  2:894701/1‑148 (MQ=38)                                                                                                                                          cATTTTTATTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGAAAAGAGAtttt        >  2:456137/1‑151 (MQ=34)                                                                                                                                               ttattaAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGAAAAGAGATTTTCTTg    <  1:934285/150‑1 (MQ=31)                                                                                                                                                tattaAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGAAAAGAGATTTTCTTGTc  >  2:1013809/1‑151 (MQ=11)                                                                                                                                                      |                                                                                                                                                 CCACGCCTGCCAGATGATATGGTCCGTCCCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGATCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGAAAAGAGATTTTCTTGTC  >  NC_000913/3129997‑3130290

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