Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
JC JC	NC_000913	2,957,831	IS5 (–) +4 bp	coding (1051‑1054/2889 nt)	ptrA ←	protease III

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	2957831 =	0 (0.000)	54 (1.060)	46/298	0.2	100%	coding (1054/2889 nt)	ptrA	protease III
	?	NC_000913	3130146 =	NA (NA)					noncoding (1/1195 nt)	IS5	repeat region
*	?	NC_000913	= 2957834	0 (0.000)	53 (1.040)	44/298	0.3	100%	coding (1051/2889 nt)	ptrA	protease III
	?	NC_000913	= 3131340	NA (NA)					noncoding (1195/1195 nt)	IS5	repeat region

GGGATAACGCCTTAAATGGCGAAGAAACGGTCTAAATAGGCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCT‑AAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGT  >  NC_000913/3131201‑3131487                                                                                                                                            |                                                                                                                                                     gggATAACGCCTTAAATGGCGAAGAAACGGTCTAAATAGGCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCT‑AAATCAGGt                                                                                                                                          >  1:525469/1‑151 (MQ=33)  ggATAACGCCTTAAATGGCGAAGAAACGGTCTAAATAGGCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCT‑AAATCAGGTc                                                                                                                                         <  1:88263/151‑1 (MQ=33)       aCGCCTTAAATGGCGAAGAAACGGTCTAAATAGGCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCT‑AAATCAGGTCATa                                                                                                                                      <  1:160139/149‑1 (MQ=35)           cTTAAATGGCGAAGAAACGGTCTAAATAGGCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCT‑AAATCAGGTCATACGCTTc                                                                                                                                >  1:586108/1‑151 (MQ=40)              aaaTGGCGAAGAAACGGTCTAAATAGGCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCT‑AAATCAGGTCATACGCTTCgag                                                                                                                             <  1:555113/151‑1 (MQ=255)              aaaTGGCGAAGAAACGGTCTAAATAGGCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCT‑AAATCAGGTCATACGCTTCg                                                                                                                               >  2:98878/1‑149 (MQ=255)                aTGGCGAAGAAACGGTCTAAATAGGCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCT‑AAATCAGGTCATACGCTTCGAGAt                                                                                                                           <  1:407738/151‑1 (MQ=255)                    cGAAGAAACGGTCTAAATAGGCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCT‑AAATCAGGTCATACGCTTCGAGATACtt                                                                                                                       >  2:738985/1‑151 (MQ=255)                     gaagaaACGGTCTAAATAGGCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCT‑AAATCAGGTCATACGCTTCGAGATACTTa                                                                                                                      <  1:494958/151‑1 (MQ=255)                       agaaACGGTCTAAATAGGCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCT‑AAATCAGGTCATACGCTTCGAGATACTTAAc                                                                                                                    <  2:523313/151‑1 (MQ=255)                        gaaACGGTCTAAATAGGCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCT‑AAATCAGGTCATACGCTTCGAGATACTTAACg                                                                                                                   <  2:43682/151‑1 (MQ=255)                                   aaTAGGCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCT‑AAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACAcc                                                                                                          >  1:610891/1‑149 (MQ=255)                                     tAGGCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCT‑AAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCg                                                                                                      >  2:781901/1‑151 (MQ=255)                                      aGGCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCT‑AAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCg                                                                                                      <  2:125969/150‑1 (MQ=255)                                      aGGCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCT‑AAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCg                                                                                                      >  2:252418/1‑150 (MQ=255)                                                  ggCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCT‑AAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCtt                                                                                         <  1:395460/151‑1 (MQ=255)                                                       ttACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCT‑AAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCcg                                                                                     <  2:557638/150‑1 (MQ=255)                                                        tACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCT‑AAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCcgac                                                                                   >  1:254325/1‑151 (MQ=255)                                                         aCGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCT‑AAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCcgacg                                                                                  <  2:199540/151‑1 (MQ=255)                                                                 aaaagTCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTAGCACCTTACCT‑AAATCAGGACATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGAGCaa                                                                             >  2:790145/1‑148 (MQ=255)                                                                  aaaaTCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCT‑AAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACAcc                                                                         <  2:68951/151‑1 (MQ=255)                                                                       cGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCT‑AAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACACcgcg                                                                      >  2:265779/1‑149 (MQ=255)                                                                       cGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCT‑AAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGa                                                                     <  2:80882/150‑1 (MQ=255)                                                                             aaaCATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCT‑AAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGccac                                                              <  2:650200/151‑1 (MQ=255)                                                                                cATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCT‑AAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGccacc                                                             >  1:120561/1‑149 (MQ=255)                                                                                  tGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCT‑AAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGa                                                          >  2:659192/1‑150 (MQ=255)                                                                                         atgaaatgaCTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCT‑AAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGAt                                                  <  1:11395/151‑1 (MQ=255)                                                                                            aaatgaCTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCT‑AAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGATTcc                                               >  1:61285/1‑151 (MQ=255)                                                                                            aaatgaCTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCT‑AAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGATTcc                                               <  1:667051/151‑1 (MQ=255)                                                                                                gaCTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCT‑AAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGATTccgccc                                           >  1:462741/1‑151 (MQ=255)                                                                                                gaCTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCT‑AAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGATTccgccc                                           >  2:572257/1‑151 (MQ=255)                                                                                                gaCTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCT‑AAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGATTccgccc                                           >  2:433516/1‑151 (MQ=255)                                                                                                     aGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCT‑AAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGATTCCGCCCGCCa                                       <  1:506802/150‑1 (MQ=255)                                                                                                      gTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCT‑AAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGATTCCGCCCGCCAGt                                     >  2:414361/1‑151 (MQ=255)                                                                                                        cAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCT‑AAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGATTCCGCCCGCCAGTa                                    >  1:764219/1‑150 (MQ=255)                                                                                                                   aaTTTCCCCGCTTATTCGCACCTTCCCT‑AAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGATTCCGCCCGCCAGTAAACCGGACTTc                         >  1:273161/1‑150 (MQ=255)                                                                                                                        ccccttttcctCGCCCCTCCCAT‑AAACGATGTCATACGCTACGAGATACTAAACCCAAACCACCGGCAAAAGGAGGGCCTCCCCGACAACCACCCCCGCGCCCCCCCCCAAAATCCCGCTCCCGCCCGCCAGTAAACCGGACTTcac                       <  2:727335/137‑1 (MQ=255)                                                                                                                            gCTTATTCGCACCTTCCCT‑AAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTaac                >  2:607080/1‑150 (MQ=255)                                                                                                                                   ccccccTCCCCA‑AAATAGGGTCATCCGCTCCGAAATACTAACCCCAAAACACCAGAAAAAGGGGCGGTTTGCCGCCAACCACCCCCGCGAGCCCCACAAAAATCCCGATCCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCAc          <  1:369173/145‑1 (MQ=255)                                                                                                                                   cGCCCCTCCCCA‑AAAGCAGGTCATGCGCTTGGATATACTAACCGCCAAACACCAGCAAAAGGGGCGGCTTGCGGACAATAACCCCCGCGACCCCCCCAAGAACCCCGATCCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACg         <  1:61786/150‑1 (MQ=255)                                                                                                                                   cGCACCTTCCCT‑AAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACgg        >  2:383839/1‑151 (MQ=255)                                                                                                                                    cccccTTCCCAAAAAGCAGGTCAACCGCTTCGATATACTAACCGCCAACCACCAGAAAAAGGAGCGGTTAGCCGACTAAAACCCCCGCGCGCCCCACCAAACCACCCATCCCCCCCGCCAGTAAACCGGACTTCACACCTaaca               <  1:242871/141‑1 (MQ=255)                                                                                                                                      aCCTTCCCT‑AAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCt     >  2:61847/1‑151 (MQ=255)                                                                                                                                       ccTTCCAA‑AAATCAGGCCACACGCTTCAATATACACAACCCCAAACACAAGCAAATTGGGCGCCTTGCCGACGACCACCCCCGCGACCCCCACAAGAACCCCGGTCCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCtg    <  2:317336/151‑1 (MQ=255)                                                                                                                                           cccT‑AAACCAGGTCATACGCTTAGAGATACTAAACCCCAACCCCCAGAAAAAGGGGCGGCTGGCCGACGACCACCCCCGCGAGCCCCACAAAAACACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCtggt  <  1:421913/149‑1 (MQ=255)                                                                                                                                            |                                                                                                                                                     GGGATAACGCCTTAAATGGCGAAGAAACGGTCTAAATAGGCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCT‑AAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGT  >  NC_000913/3131201‑3131487

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