Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	2,804,864	Δ13 bp	coding (50‑62/1203 nt)	proV →	glycine betaine/proline ABC transporter periplasmic binding protein

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	2804864	2804876	13	45 [0]	[0] 45	proV	glycine betaine/proline ABC transporter periplasmic binding protein

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 2804863	0 (0.000)	40 (0.940)	36/286	0.4	100%	coding (49/1203 nt)	proV	glycine betaine/proline ABC transporter periplasmic binding protein
	?	NC_000913	2804877 =	0 (0.000)					coding (63/1203 nt)	proV	glycine betaine/proline ABC transporter periplasmic binding protein

AGTATGTTAGGGTAGAAAAAAGTGACTATTTCCATTGGGTAATATATCGACATAGACAAATAAAGGAATCTTTCTATTGCATGGCAATTAAATTAGAAATTAAAAATCTTTATAAAATATTTGGCGAGCATCCACAGCGAGCGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTAAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGTCATCATGGGATTATCCGGCTCGGGTAAATCCACAATGGTA  >  NC_000913/2804735‑2805027                                                                                                                                               |                                                                                                                                                       gcctcagattctcagtatgttagggtagaaaaaagtgactatttccattgggtaatatatcgacatagacaaataaaggaatctttctattgcatggcaattaaattagaaattaaaaatctttataaaatatttggcgagcGTTCAAata                                                                                                                                                <  1:614812‑M2/9‑1 (MQ=255)     cagattctcagtatgttagggtagaaaaaagtgactatttccattgggtaatatatcgacatagacaaataaaggaatctttctattgcatggcaattaaattagaaattaaaaatctttataaaatatttggcgagcGTTCAAatata                                                                                                                                              <  2:604079‑M2/11‑1 (MQ=255)      agattctcagtatgttagggtagaaaaaagtgactatttccattgggtaatatatcgacatagacaaataaaggaatctttctattgcatggcaattaaattagaaattaaaaatctttataaaatatttggcgagcGTTCAAATATATCg                                                                                                                                           >  1:117153‑M2/138‑151 (MQ=255)             cagtatgttagggtagaaaaaagtgactatttccattgtgtaatatatcgacatagacaaataaaggaatctttctattgcatggcaattaaattagaaattaaaaatctttataaaatatttggcgagcGTTCAAATATATCGAACaa                                                                                                                                      <  2:499910‑M2/19‑1 (MQ=255)                   gttagggtagaaaaaagtgactatttccattgggtaatatatcgacatagacaaataaaggaatctttctattgcatggcaattaaattagaaattaaaaatctttataaaatatttggcgagcGTTCAAATATATCGAACAAGGACTTTc                                                                                                                              <  1:60110‑M2/27‑1 (MQ=255)                       gggtagaaaaaagtgactatttccattgggtaatatatcgacatagacaaataaaggaatctttctattgcatggcaattaaattagaaattaaaaatctttataaaatatttggcgagcGTTCAAATATATCGAACAAGGACTTTCaaa                                                                                                                           >  1:647431‑M2/121‑150 (MQ=255)                           agaaaaaagtgactatttccattgggtaatatatcgacatagacaaataaaggaatctttctattgcatggcaattaaattagaaattaaaaatctttataaaatatttggcgagcGTTCAAATATATCGAACAAGGACTTTCAAAAGaa                                                                                                                       <  1:195961‑M2/34‑1 (MQ=255)                             aaaaaagtgactatttccattgggtaatatatcgacatagacaaataaaggaatctttctattgcatggcaattaaattagaaattaaaaatctttataaaatatttggcgagcGTTCAAATATATCGAACAAGGACTTTCAAAAGAACa                                                                                                                     >  2:30761‑M2/115‑150 (MQ=255)                               aaaagtgactatttccattgggtaatatatcgacatagacaaataaaggaatctttctattgcatggcaattaaattagaaattaaaaatctttataaaatatttggcgagcGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAAt                                                                                                                  >  2:612758‑M2/113‑151 (MQ=255)                               aaaagtgactatttccattgggtaatatatcgacatagacaaataaaggaatctttctattgcatggcaattaaattagaaattaaaaatctttataaaatatttggcgagcGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAAt                                                                                                                  >  2:144104‑M2/113‑151 (MQ=255)                                   gtgactatttccattgggtaatatatcgacatagacaaataaaggaatctttctattgcatggcaattaaattagaaattaaaaatctttataaaatatttggcgagcGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTg                                                                                                              <  1:444545‑M2/43‑1 (MQ=255)                                      actatttccattgggtaatatatcgacatagacaaataaaggaatctttctattgcatggcaattaaattagaaattaaaaatctttataaaatatttggcgagcGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGaa                                                                                                           <  2:29645‑M2/46‑1 (MQ=255)                                                 tggttaatatatcgacatagacaaataaaggaatctttctattgcatggcaattaaattagaaattaaaaatctttataaaataattggcgagcGTTCAAATATATCGAACAACGACTTTCAAAAGAACAAATTCTGGAAAAAACTgcgc                                                                                                 <  2:626702‑M2/56‑1 (MQ=255)                                                  gggtaatatatcgacatagacaaataaaggaatctttctattgcatggcaattaaattagaaattaaaaatctttataaaatatttggcgagcGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCt                                                                                                <  2:311123‑M2/57‑1 (MQ=255)                                                       atatatcgacatagacaaataaaggaatctttctattgcatggcaattaaattagaaattaaaaatctttataaaatatttggcgagcGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCt                                                                                          >  1:346784‑M2/89‑151 (MQ=255)                                                           atcgacatagacaaataaaggaatctttctattgcatggcaattaaattagaaattaaaaatctttataaaatatttggcgagcGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCtt                                                                                         >  2:110877‑M2/85‑148 (MQ=255)                                                           atcgacatagacaaataaaggaatctttctattgcatggcaattaaattagaaattaaaaatctttataaaatatttggcgagcGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTgg                                                                                       <  1:690671‑M2/66‑1 (MQ=255)                                                              gacatagacaaataaaggaatctttctattgcatggcaattaaattagaaattaaaaatctttataaaatatttggcgagcGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTa                                                                                   <  1:227709‑M2/70‑1 (MQ=255)                                                               acatagacaaataaaggaatctttctattgcatggcaattaaattagaaattaaaaatctttataaaatatttggcgagcGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTa                                                                                   <  2:117131‑M2/70‑1 (MQ=255)                                                                   agacaaataaaggaatctttctattgcatggcaattaaattagaaattaaaaatctttataaaatatttggcgagcGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTaaa                                                                                 >  2:551499‑M2/77‑148 (MQ=255)                                                                        aataaaggaatctttctattgcatggcaattaaattagaaattaaaaatctttataaaatatttggcgagcGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTAAAAGACGCCa                                                                         <  2:292775‑M2/80‑1 (MQ=255)                                                                           aaaggaatctttctattgcatggcaattaaattagaaattaaaaatctttataaaatatttggcgagcGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTAAAAGACGCCAGt                                                                       <  1:201356‑M2/82‑1 (MQ=255)                                                                           aaaggaatctttctattgcatggcaattaaattagaaattaaaaatctttataaaatatttggcgagcGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTAAAAGACGCCAGt                                                                       <  2:383393‑M2/82‑1 (MQ=255)                                                                             aggaatctttctattgcatggcaattaaattagaaattaaaaatctttataaaatatttggcgagcGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTAAAAGACGCCAGTCTg                                                                    <  1:34694‑M2/85‑1 (MQ=255)                                                                                    tttctattgcatggcaattaaattagaaattaaaaatctttataaaatatttggcgagcGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTAAAAGACGCCAGTCTGGCCATTg                                                             >  2:456939‑M2/60‑151 (MQ=255)                                                                                     ttctattgcatggcaattaaattagaaattaaaaatctttataaaatatttggcgagcGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTAAAAGACGCCAGTCTGGCCAtt                                                              <  1:348496‑M2/91‑1 (MQ=255)                                                                                     ttctattgcatggcaattaaattagaaattaaaaatctttataaaatatttggcgagcGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTAAAAGACGCCAGTCTGGCCATTg                                                             <  2:247937‑M2/92‑1 (MQ=255)                                                                                         attgcatggcaattaaattagaaattaaaaatctttataaaatatttggcgagcGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTAAAAGACGCCAGTCTGGCCATTgaagaa                                                        <  1:56337‑M2/97‑1 (MQ=255)                                                                                            gcatggcaattaaattagaaattaaaaatctttataaaatatttggcgagcGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTAAAAGACGCCAGTCTGGCCATTGAAGAAGGc                                                     >  1:691512‑M2/52‑151 (MQ=255)                                                                                                 gcaattaaattagaaattaaaaatctttataaaatatttggcgagcGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTAAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGat                                                <  1:649832‑M2/105‑1 (MQ=255)                                                                                                         attagaaattaaaaatctttataaaatatttggcgagcGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTAAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGtc                                         >  2:212891‑M2/39‑150 (MQ=255)                                                                                                          ttagaaattaaaaatctttataaaatatttggcgagcGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTAAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGtca                                        >  1:351192‑M2/38‑150 (MQ=255)                                                                                                              aaattaaaaatctttataaaatatttggcgagcGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTAAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGTCATCATg                                   >  2:280348‑M2/34‑151 (MQ=255)                                                                                                                attaaaaatctttataaaatatttggcgagcGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTAAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGTCATCATggg                                 >  1:585257‑M2/32‑151 (MQ=255)                                                                                                                              taaaatatttggcgagcGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTAAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGTCATCATGGGATTATCCGGCTCgg                   >  1:230717‑M2/18‑151 (MQ=255)                                                                                                                                aaatatttggcgagcGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTAAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGTCATCATGGGATTATCCGGCTCGGGt                 >  1:18573‑M2/16‑151 (MQ=255)                                                                                                                                 aatatttggcgagcGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTAAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGTCATCATGGGATTATCCGGCTCGGGTa                <  1:331046‑M2/137‑1 (MQ=255)                                                                                                                                  atgtttggcgagcGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTAAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGTCATCATGGGATTATCCGGCTCGGGTa                >  2:339754‑M2/14‑150 (MQ=255)                                                                                                                                    atttggcgagcGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTAAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGTCATCATGGGATTATCCGGCTCGGGTAAAt             >  1:30384‑M2/12‑151 (MQ=255)                                                                                                                                      ttggcgagcGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTAAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGTCATCATGGGATTATCCGGCTCGGGTAAATc            <  1:375325‑M2/141‑1 (MQ=255)                                                                                                                                          cgagcGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTAAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGTCATCATGGTATTATCCGGCTCGGGTAAATCcaca        >  2:67408‑M2/6‑150 (MQ=255)                                                                                                                                           gagcGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTAAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGTCATCATGGGATTATCCGGCTCGGGTAAATCCACAAt      <  1:236210‑M2/147‑1 (MQ=255)                                                                                                                                              cGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTAAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGTCATCATGGGATTATCCGGCTCGGGTAAATCCACAATGGt   >  1:200565‑M2/2‑151 (MQ=255)                                                                                                                                               gTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTAAAAGACGGCAGTCTGGCCATTGAAGAAGGCGAGATATTTGTCATCATGGGATTATCCGGCTCGGGTAAATCCACAATGGTa  <  2:592939/151‑1 (MQ=255)                                                                                                                                               gTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTAAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGTCATCATGGGATTATCCGGCTCGGGTAAATCCACAATGGt   <  1:144131/150‑1 (MQ=255)                                                                                                                                               |                                                                                                                                                       AGTATGTTAGGGTAGAAAAAAGTGACTATTTCCATTGGGTAATATATCGACATAGACAAATAAAGGAATCTTTCTATTGCATGGCAATTAAATTAGAAATTAAAAATCTTTATAAAATATTTGGCGAGCATCCACAGCGAGCGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTAAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGTCATCATGGGATTATCCGGCTCGGGTAAATCCACAATGGTA  >  NC_000913/2804735‑2805027

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑