Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
JC	NC_000913	3,998,555	Δ3 bp	9.5%	coding (573‑575/2163 nt)	uvrD →	DNA‑dependent ATPase I and helicase II

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	3998558 =	42 (0.960)	4 (0.100)	4/240	NT	9.5%	coding (576/2163 nt)	uvrD	DNA‑dependent ATPase I and helicase II
	?	NC_000913	= 3998554	37 (0.900)					coding (572/2163 nt)	uvrD	DNA‑dependent ATPase I and helicase II

TGGTACATCAACAGCCAGAAAGATGAAGGCCTGCGTCCGCATCATATTCAAAGCTACGGTAATCCGGTGGAGCAGACCTGGCAGAAGGTGTATCAGGCGTATCAGGAAGCGTGTGACCGCGCGGGCCTGGTGGACTTCGCCGAGCTGCTGCTGCGCGCTCACGAGTTGTGGCTTAACAAGCCGCATATCCTGCAACACTACCGCGAACGTTTTACCAATATCCTGGTGGACGAATTCCAGGATACCAACAACATTCAGTACGCGTGGATCCGCCTGCTGGCGGGCGA  >  NC_000913/3998412‑3998698                                                                                                                                               |                                                                                                                                                 tgGTACATCAACAGCCAGAAAGATGAAGGCCTGCGTCCGCATCATATTCAAAGCTACGGTAATCCGGTGGAGCAGACCTGGCAGAAGGTGTATCAGGCGTATCAGGAAGCGTGCGACCGCGCGGGCCTGGTGGACTTCGCCGAgctgct                                                                                                                                            <  2:77517‑M2/149‑7 (MQ=255)   ggacaTCAACAGCCAGAAAGATGAAGGCCTGCGTCCGCATCATATTCAAAGCTACGGTAATCCGGTGGAGCAGACCTGGCAGAAGGTGTATCAGGCGTATCAGGAAGCGTGTGACCGCGCGGGCCTGGTGGACTTCGCCGAgctgctgc                                                                                                                                          <  1:359663‑M2/147‑9 (MQ=255)     gcATCAACAGCCAGAAAGATGAAGGCCTGCGTCCGCATCATATTCAAAGCTACGGTAATCCGGTGGAGCAGACCTGGCAGAAGGTGTATCAGGCGTATCAGGAAGCGTGTGACCGCGCGGGCCTGGTGGACTTCGCCGAGctgctgctg                                                                                                                                        >  1:23085/2‑149 (MQ=255)        tCAACAGCCAGAAAGATGAAGGCCTGCGTCCGCATCACATTCAAAGCTACGGTAATCCGGTGGAGCAGACCTGGCAGAAGGTGTATCAGGCGTATCAGGAAGCGTGTGACCGCGCGGGCCTGGTGGACTTCGCCGAGCTGCTGCTgcg                                                                                                                                      >  3:104935/1‑148 (MQ=255)           aCAGCCAGAAAGATGAAGGCCTGCGTCCGCATCATATTCAAAGCTACGGTAATCCGGTGGAGCAGACCTGGCAGAAGGTGTATCAGGCGTATCAGGAAGCGTGTGACCGCGCGGGCCTGGTGGACTTCGCCGAGCTGCTGCTGCGCGCt                                                                                                                                  >  1:45531/1‑149 (MQ=255)                 aGAAAGATGAAGGCCTGCGTCCGCATCATATTCAAAGCTACGGTAATCCGGTGGAGCAGACCTGGCAGAAGGTGTATCAGGCGTATCAGGAAGCGTGTGACCGCGCGGGCCTGGTGGACTTCGCCGAGCTGCTGCTGCGCGCTCACGAg                                                                                                                            >  1:201038/1‑149 (MQ=255)                  gAAAGATGAAGGCCTGCGTCCGCATCATATTCAAAGCTACGGTAATCCGGTGGAGCAGACCTGGCAGAAGGTGTATCAGGCGTATCAGGAAGCGTGTGACCGCGCGGGCCTGGTGGACTTCGCCGAGCTGCTGCTGCGCGCTCACGAGt                                                                                                                           >  2:287155/1‑149 (MQ=255)                     aGATGAAGGCCTGCGTCCGCATCATATTCAAAGCTACGGTAATCCGGTGGAGCAGACCTGGCAGAAGGTGTATCAGGCGTATCAGGAAGCGTGTGACCGCGCGGGCCTGGTGGACTTCGCCGAGCTGCTGCTGCGCGCTCACGAGTtgt                                                                                                                        <  2:51687/149‑1 (MQ=255)                      gATGAAGGCCTGCGTCCGCATCATATTCAAAGCTACGGTAATCCGGTGGAGCAGACCTGGCAGAAGGTGTATCAGGCGTATCAGGAAGCGTGTGACCGCGCGGGCCTGGTGGACTTCGCCGAGCTGCTGCTGCGCGCTCACGAGTtgt                                                                                                                        <  4:237914/148‑1 (MQ=255)                         gAAGGCCTGCGTCCGCATCATATTCAAAGCTACGGTAATCCGGTGGAGCAGACCTGGCAGAAGGTGTATCAGGCGTATCAGGAAGCGTGTGACCGCGCGGGCCTGGTGGACTTCGCCGAGCTGCTGCTGCGCGCTCACGAGTTGTGGCt                                                                                                                    >  2:275171/1‑149 (MQ=255)                          aaGGCCTGCGTCCGCATCATATTCAAAGCTACGGTAATCCGGTGGATCAGACCTGGCAGAAGGTGTATCAGGCGTATCAGGAAGCGTGTGACAGCGCGGGCATGGTGGACTTCGCCGAGCTGCTGCTGCGCGCTCACGAGTTGTGGCt                                                                                                                    >  4:241672/1‑148 (MQ=255)                                        catcatATTCAAAGCTACGGTAATCCGGTGGAGCAGACCTGGCAGAAGGTGTATCAGGCGTATCAGGAAGCGTGTGACCGCGCGGGCCTGGTGGACTTCGCCGAGCTGCTGCTGCGCGCTCACGAGTTGTGGCTTAACAAGCCGCata                                                                                                      >  3:110164/1‑148 (MQ=255)                                         atcatATTCAAAGCTACGGTAATCCGGTGGAGCAGACCTGGCAGAAGGTGTATCAGGCGTATCAGGAAGCGTGTGACCGCGCGGGCCTGGTGGACTTCGCCGAGCTGCTGCTGCGCGCTCACGAGTTGTGGCTTAACAAGCCGCatat                                                                                                     <  4:183680/148‑1 (MQ=255)                                         atcatATTCAAAGCTACGGTAATCCGGTGGAGCAGACCTGGCAGAAGGTGTATCAGGCGTATCAGGAAGCGTGTGACCGCGCGGGCCTGGTGGACTTCGCCGAGCTGCTGCTGCGCGCTCACGAGTTGTGGCTTAACAAGCCGCatat                                                                                                     >  3:400212/1‑148 (MQ=255)                                         atcatATTCAAAGCTACGGTAATCCGGTGGAGCAGACCTGGCAGAAGGTGTATCAGGCGTATCAGGAAGCGTGTGACCGCGCGGGCCTGGTGGACTTCGCCGAGCTGCTGCTGCGCGCTCACGAGTTGTGGCTTAACAAGCCGCatat                                                                                                     >  4:323234/1‑148 (MQ=255)                                           catATTCAAAGCTACGGTAATCCGGTGGAGCAGACCTGGCAGAAGGTGTATCAGGCGTATCAGGAAGCGTGTGACCGCGCGGGCCTGGTGGACTTCGCCGAgctg                                                                                                                                              <  4:241132‑M2/105‑5 (MQ=255)                                           catATTCAAAGCTACGGTAATCCGGTGGAGCAGACCTGGCAGAAGGTGTATCAGGCGTATCAGGAAGCGTGTGACCGCGCGGGCCTGGTGGACTTCGCCGAgctg                                                                                                                                              >  3:241132‑M2/1‑101 (MQ=255)                                                tCAAAGCTACGGTAATCCGGTGGAGCAGACCTGGCAGAAGGTGTATCAGGCGTATCAGGAAGCGTGTGACCGCGCGGGCCTGGTGGACTTCGCCGAgctgctgcgcgctcacgagttgtggcttaacaagccgcatatcctgcaacac                                                                                              <  4:135874‑M2/148‑53 (MQ=255)                                                  aaaGCTACGGTAATCCGGTGGAGCAGACCTGGCAGAAGGTGTATCAGGCGTATCAGGAAGCGTGTGACCGCGCGGGCCTGGTGGACTTCGCCGAGCTGCTGCTGCGCGCTCACGAGTTGTGGCTTAACAAGCCGCATATCCTGCAacac                                                                                           >  1:256168/1‑149 (MQ=255)                                                     gCTACGGTAATCCGGTGGAGCAGACCTGGCAGAAGGTGTATCAGGCGTATCAGGAAGCGTGTGACCGCGCGGGCCTGGTGGACTTCGCCGAGCTGCTGCTGCGCGCTCACGAGTTGTGGCTTAACAAGCCGCATATCCTGCAACACTa                                                                                         >  4:126966/1‑148 (MQ=255)                                                               tCCGGTGGAGCAGACCTGGCAGAAGGTGTATCAGGCGTATCAGGAAGCGTGTGACCGCGCGGGCCTGGTGGACTTCGCCGAgctgctgcgcgctcacgagttgtggcttaacaagccgcatatcctgcaacactaccgc                                                                                        >  1:79034‑M2/1‑81 (MQ=255)                                                               tCCGGTGGAGCAGACCTGGCAGAAGGTGTATCAGGCGTATCAGGAAGCGTGTGACCGCGCGGGCCTGGTGGACTTCGCCGAgctgctgcgcgctcacgagttgtggcttaacaagccgcatatcctgcaacactaccgc                                                                                        <  2:79034‑M2/139‑59 (MQ=255)                                                                ccGGTGGAGCAGACCTGGCAGAAGGTGTATCAGGCGTATCAGGAAGCGTGTGACCGCGCGGGCCTGGTGGACTTCGCCGAGCTGCTGCTGCGCGCTCACGAGTTGTGGCTTAACAAGCCGCATATCCTGCAACACTACCGCGAACGttt                                                                             <  1:198444/149‑1 (MQ=255)                                                                ccGGTGGAGCAGACCTGGCAGAAGGTGTATCAGGCGTATCAGGAAGCGTGTGACCGCGCGGGCCTGGTGGACTTCGCCGAGCTGCTGCTGCGCGCTCACGAGTTGTGGCTTAACAAGCCGCATATCCTGCAACACTACCGCGAACGtt                                                                              <  3:126966/148‑1 (MQ=255)                                                                          aGACCTGGCAGAAGGTGTATCAGGCGTATCAGGAAGCGTGTGACCGCGCGGGCCTGGTGGACTTCGCCGAGCTGCTGCTGCGCGCTCACGAGTTGTGGCTTAACAAGCCGCATATCCTGCAACACTACCGCGAACGTTTTACCAATATc                                                                   >  2:101317/1‑149 (MQ=255)                                                                          aGACCTGGCAGAAGGTGTATCAGGCGTATCAGGAAGCGTGCGACCGCGCGGGCCTGGTGGACTTCGCCGAGCTGCTGCTGCGCGCTCACGAGTTGTGGCTTAACAAGCCGCATATCCTGCAACACTACCGCGAACGTTTTACCAATATc                                                                   <  1:77517/149‑1 (MQ=255)                                                                             ccTGGCAGAAGGTGTATCAGGCGTATCAGGAAGCGTGTGACCGCGCGGGCCTGGTGGACTTCGCCGAGCTGCTGCTGCGCGCTCACGAGTTGTGGCTTAACAAGCCGCATATCCTGCAACACTACCGCGAACGTTTTACCAATATCCtg                                                                <  2:201038/149‑1 (MQ=255)                                                                                   aGAAGGTGTATCAGGCGTATCAGGAAGCGTGTGACCGCGCGGGCCTGGTGGACTTCGCCGAGCTGCTGCTGCGCGCTCACGAGTTGTGGCTTAACAAGCCGCATATCCTGCAACACTACCGCGAACGTTTTACCAATATCCTGGTGGAc                                                          >  1:68944/1‑149 (MQ=255)                                                                                    gAAGGTGTATCAGGCGTATCAGGAAGCGTGTGACCGCGCGGGCCTGGTGGACTTCGCCGAgctgctgcgcgctcacgagttgtggcttaacaagccgcatatcctgcaacactaccgcgaacgttttaccaatatcctggtggacgaa                                                          >  3:135874‑M2/1‑60 (MQ=255)                                                                                           tATCAGGCGTATCAGGAAGCGTGTGACCGCGCGGGCCTGGTGGACTTCGCCGAGCTGCTGCTGCGCGCTCACGAGTTGTGGCTTAACAAGCCGCATATCCTGCAACACTACCGCGAACGTTTTACCAATATCCTGGTGGACGAATTcc                                                   <  3:5543/148‑1 (MQ=255)                                                                                           tATCAGGCGTATCAGGAAGCGTGTGACCGCGCGGGCCTGGTGGACTTCGCCGAGCTGCTGCTGCGCGCTCACGAGTTGTGGCTTAACAAGCCGCATATCCTGCAACACTACCGCGAACGTTTTACCAATATCCTGGTGGACGAATTCCa                                                  >  1:360003/1‑149 (MQ=255)                                                                                                   gTATCAGGAAGCGTGTGACCGCGCGGGCCTGGTGGACTTCGCCGAGCTGCTGCTGCGCGCTCACGAGTTGTGGCTTAACAAGCCGCATATCCTGCAACACTACCGCGAACGTTTTACCAATATCCTGGTGGACGAATTCCAGGATAcc                                           <  4:110164/148‑1 (MQ=255)                                                                                                   gTATCAGGAAGCGTGTGACCGCGCGGGCCTGGTGGACTTCGCCGAGCTGCTGCTGCGCGCTCACGAGTTGTGGCTTAACAAGCCGCATATCCTGCAACACTACCGCGAACGTTTTACCAATATCCTGGTGGACGAATTCCAGGATACca                                          <  1:101317/149‑1 (MQ=255)                                                                                                    tATCAGGAAGCGTGTGACCGCGCGGGCCTGGTGGACTTCGCCGAGCTGCTGCTGCGCGCTCACGAGTTGTGGCTTAACAAGCCGCATATCCTGCAACACTACCGCGAACGTTTTACCAATATCCTGGTGGACGAATTCCAGGATACcaa                                         <  1:320718/149‑1 (MQ=255)                                                                                                    tATCAGGAAGCGTGTGACCGCGCGGGCCTGGTGGACTTCGCCGAGCTGCTGCTGCGCGCTCACGAGTTGTGGCTTAACAAGCCGCATATCCTGCAACACTACCGCGAACGTTTTACCAATATCCTGGTGGACGAATTCCAGGATACca                                          <  4:249060/148‑1 (MQ=255)                                                                                                       cAGGAAGCGTGTGACCGCGCGGGCCTGGTGGACTTCGCCGAGCTGCTGCTGCGCGCTCACGAGTTGTGGCTTAACAAGCCGCATATCCTGCAACACTACCGCGAACGTTTTACCAATATCCTGGTGGACGAATTCCAGGATACcaaca                                       <  3:341334/148‑1 (MQ=255)                                                                                                          gAAGCGTGTGACCGCGCGGGCCTGGTGGACTTCGCCGAGCTGCTGCTGCGCGCTCACGAGTTGTGGCTTAACAAGCCGCATATCCTGCAACACTACCGCGAACGTTTTACCAATATCCTGGTGGACGAATTCCAGGATACCAACAACAt                                   >  2:5333/1‑149 (MQ=255)                                                                                                            aGCGTGTGACCGCGCGGGCCTGGTGGACTTCGCCGAGCTGCTGCTGCGCGCTCACGAGTTGTGGCTTAACAAGCCGCATATCCTGCAACACTACCGCGAACGTTTTACCAATATCCTGGTGGACGAATTCCAGGATACCAACAACATTc                                 <  2:256168/149‑1 (MQ=255)                                                                                                                tgtgACCGCGCGGGCCTGGTGGACTTCGCCGAGCTGCTGCTGCGCGCTCACGAGTTGTGGCTTAACAAGCCGCATATCCTGCAACACTACCGCGAACGTTTTACCAATATCCTGGTGGACGAATTCCAGGATACCAACAACATTCAGt                              >  4:60760/1‑148 (MQ=255)                                                                                                                   gACCGCGCGGGCCTGGTGGACTTCGCCGAGCTGCTGCTGCGCGCTCACGAGTTGTGGCTTAACAAGCCGCATATCCTGCAACACTACCGCGAACGTTTTACCAATATCCTGGTGGACGAATTCCAGGATACCAACAACATTCAGTAcg                           <  3:60760/148‑1 (MQ=255)                                                                                                                     ccGCGCGGGCCTGGTGGACTTCGCCGAGCTGCTGCTGCGCGCTCACGAGTTGTGGCTTAACAAGCCGCATATCCTGCAACACTACCGCGAACGTTTTACCAATATCCTGGTGGACGAATTCCAGGATACCAACAACATTCAGTACGCGt                        <  1:269105/149‑1 (MQ=255)                                                                                                                             gCCTGGTGGACTTCGCCGAGCTGCTGCTGCGCGCTCACGAGTTGTGGCTTAACAAGCCGCATATCCTGCAACACTACCGCGAACGTTTTACCAATATCCTGGTGGACGAATTCCAGGATACCAACAACATTCAGTACGCGTGGATCCGc                <  2:359663/149‑1 (MQ=255)                                                                                                                                tggtggACTTCGCCGAGCTGCTGCTGCGCGCTCACGAGTTGTGGCTTAACAAGCCGCATATCCTGCAACACTACCGCGAACGTTTTACCAATATCCTGGTGGACGAATTCCAGGATACCAACAACATTCAGTACGCGTGGATCCGCctg             >  2:254819/1‑149 (MQ=255)                                                                                                                                     gACTTCGCCGAGCTGCTGCTGCGCGCTCACGAGTTGTGGCTTAACAAGCCGCATATCCTGCAACACTACCGCGAACGTTTTACCAATATCCTGGTGGACGAATTCCAGGATACCAACAACATTCAGTACGCGTGGATCCGCCTGCTgg         >  4:396174/1‑148 (MQ=255)                                                                                                                                        ttCGCCGAGCTGCTGCTGCGCGCTCACGAGTTGTGGCTTAACAAGCCGCATATCCTGCAACACTACCGCGAACGTTTTACCAATATCCTGGTGGACGAATTCCAGGATACCAACAACATTCAGTACGCGTGGATCCGCCTGCTggcggg     >  2:147554/1‑149 (MQ=255)                                                                                                                                        ttCGCCGAGCTGCTGCTGCGCGCTCACGAGTTGTGGCTTAACAAGCCGCATATCCTGCAACACTACCGCGAACGTTTTACCAATATCCTGGTGGACGAATTCCAGGATACCAACAACATTCAGTACGCGTGGATCCGCCTGCTggcggg     >  1:43103/1‑149 (MQ=255)                                                                                                                                         tCGCCGAGCTGCTGCTGCGCGCTCACGAGTTGTGGCTTAACAAGCCGCATATCCTGCAACACTACCGCGAACGTTTTACCAATATCCTGGTGGACGAATTCCAGGATACCAACAACATTCAGTACGCGTGGATCCGCCTGCTggcgggc    <  1:275171/149‑1 (MQ=255)                                                                                                                                           gCCGAGCTGCTGCTGCGCGCTCACGAGTTGTGGCTTAACAAGCCGCATATCCTGCAACACTACCGCGAACGTTTTACCAATATCCTGGTGGACGaa                                                       >  3:111094/1‑96 (MQ=255)                                                                                                                                           gCCGAGCTGCTGCTGCGCGCTCACGAGTTGTGGCTTAACAAGCCGCATATCCTGCAACACTACCGCGAACGTTTTACCAATATCCTGGTGGACGaa                                                       <  4:111094/96‑1 (MQ=255)                                                                                                                                           gCCGAGCTGCTGCTGCGCGCTCACGAGTTGTGGCTTAACAAGCCGCATATCCTGCAACACTACCGCGAACGTTTTACCAATATCCTGGTGGACGAATTCCAGGATACCAACAACATTCAGTACGCGTGGATCCGCCTGCTGGCGGGCGa  <  1:402773/149‑1 (MQ=255)                                                                                                                                               |                                                                                                                                                 TGGTACATCAACAGCCAGAAAGATGAAGGCCTGCGTCCGCATCATATTCAAAGCTACGGTAATCCGGTGGAGCAGACCTGGCAGAAGGTGTATCAGGCGTATCAGGAAGCGTGTGACCGCGCGGGCCTGGTGGACTTCGCCGAGCTGCTGCTGCGCGCTCACGAGTTGTGGCTTAACAAGCCGCATATCCTGCAACACTACCGCGAACGTTTTACCAATATCCTGGTGGACGAATTCCAGGATACCAACAACATTCAGTACGCGTGGATCCGCCTGCTGGCGGGCGA  >  NC_000913/3998412‑3998698

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