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

CAGGGGTTGCCTCAGATTCTCAGTATGTTAGGGTAGAAAAAAGTGACTATTTCCATTGGGTAATATATCGACATAGACAAATAAAGGAATCTTTCTATTGCATGGCAATTAAATTAGAAATTAAAAATCTTTATAAAATATTTGGCGAGCATCCACAGCGAGCGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTAAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGTCATCATGGGATTATCCGGCTCGGGTA  >  NC_000913/2804714‑2805013                                                                                                                                                      |                                                                                                                                                       cAGGGGTTGCCTCAGATTCTCAGTATGTTAGGGTAGAAAAAAGTGACTATTTCCATTGGGTAATATATCAACATAGACAAATAAAGGAATCTTTCTATTGCATGGCAATTAAATTAGAAATTAAAAATCTTTATAAAATATTTGGCGAGC                                                                                                                                                        <  2:366584‑M1/150‑1 (MQ=255)  ctgggTTGCCTCCGCTTCTCCGTCTGTTAGGGTAGACCAACGTGCCTCTTTCCCTTGGGTCCTCTCTCGCCCTCGCCCCATCCCGGAATCTTTCTATTGCATGGCAATTAAATTAGAAATTAAAAATCTTTATAAAATATTTGGCGAGcgt                                                                                                                                                      >  1:397658/3‑149 (MQ=255)         gCCTCAGATTCTCAGTATGTTAGGGTAGAAAAAAGTGACTATTTCCATTGGGTAATATATCGACATAGACAAATAAAGGAATCTTTCTATTGCATGGCAATTAAATTAGAAATTAAAAATCTTTATAAAATATTTGGCGAGCgttcaaata                                                                                                                                               <  1:614812‑M1/151‑10 (MQ=255)             cAGATTCTCAGTATGTTAGGGTAGAAAAAAGTGACTATTTCCATTGGGTAATATATCGACATAGACAAATAAAGGAATCTTTCTATTGCATGGCAATTAAATTAGAAATTAAAAATCTTTATAAAATATTTGGCGAGCgttcaaatata                                                                                                                                             <  2:604079‑M1/149‑12 (MQ=255)              aGATTCTCAGTATGTTAGGGTAGAAAAAAGTGACTATTTCCATTGGGTAATATATCGACATAGACAAATAAAGGAATCTTTCTATTGCATGGCAATTAAATTAGAAATTAAAAATCTTTATAAAATATTTGGCGAGCgttcaaatatatcg                                                                                                                                          >  1:117153‑M1/1‑137 (MQ=255)                     cAGTATGTTAGGGTAGAAAAAAGTGACTATTTCCATTGTGTAATATATCGACATAGACAAATAAAGGAATCTTTCTATTGCATGGCAATTAAATTAGAAATTAAAAATCTTTATAAAATATTTGGCGAGCgttcaaatatatcgaacaa                                                                                                                                     <  2:499910‑M1/149‑20 (MQ=255)                           gTTAGGGTAGAAAAAAGTGACTATTTCCATTGGGTAATATATCGACATAGACAAATAAAGGAATCTTTCTATTGCATGGCAATTAAATTAGAAATTAAAAATCTTTATAAAATATTTGGCGAGCgttcaaatatatcgaacaaggactttc                                                                                                                             <  1:60110‑M1/151‑28 (MQ=255)                               gggTAGAAAAAAGTGACTATTTCCATTGGGTAATATATCGACATAGACAAATAAAGGAATCTTTCTATTGCATGGCAATTAAATTAGAAATTAAAAATCTTTATAAAATATTTGGCGAGCgttcaaatatatcgaacaaggactttcaaa                                                                                                                          >  1:647431‑M1/1‑120 (MQ=255)                                   aGAAAAAAGTGACTATTTCCATTGGGTAATATATCGACATAGACAAATAAAGGAATCTTTCTATTGCATGGCAATTAAATTAGAAATTAAAAATCTTTATAAAATATTTGGCGAGCgttcaaatatatcgaacaaggactttcaaaagaa                                                                                                                      <  1:195961‑M1/150‑35 (MQ=255)                                     aaaaaaGTGACTATTTCCATTGGGTAATATATCGACATAGACAAATAAAGGAATCTTTCTATTGCATGGCAATTAAATTAGAAATTAAAAATCTTTATAAAATATTTGGCGAGCgttcaaatatatcgaacaaggactttcaaaagaaca                                                                                                                    >  2:30761‑M1/1‑114 (MQ=255)                                       aaaaGTGACTATTTCCATTGGGTAATATATCGACATAGACAAATAAAGGAATCTTTCTATTGCATGGCAATTAAATTAGAAATTAAAAATCTTTATAAAATATTTGGCGAGCgttcaaatatatcgaacaaggactttcaaaagaacaaat                                                                                                                 >  2:144104‑M1/1‑112 (MQ=255)                                       aaaaGTGACTATTTCCATTGGGTAATATATCGACATAGACAAATAAAGGAATCTTTCTATTGCATGGCAATTAAATTAGAAATTAAAAATCTTTATAAAATATTTGGCGAGCgttcaaatatatcgaacaaggactttcaaaagaacaaat                                                                                                                 >  2:612758‑M1/1‑112 (MQ=255)                                           gTGACTATTTCCATTGGGTAATATATCGACATAGACAAATAAAGGAATCTTTCTATTGCATGGCAATTAAATTAGAAATTAAAAATCTTTATAAAATATTTGGCGAGCgttcaaatatatcgaacaaggactttcaaaagaacaaattctg                                                                                                             <  1:444545‑M1/151‑44 (MQ=255)                                              aCTATTTCCATTGGGTAATATATCGACATAGACAAATAAAGGAATCTTTCTATTGCATGGCAATTAAATTAGAAATTAAAAATCTTTATAAAATATTTGGCGAGCgttcaaatatatcgaacaaggactttcaaaagaacaaattctggaa                                                                                                          <  2:29645‑M1/151‑47 (MQ=255)                                                         tGGTTAATATATCGACATAGACAAATAAAGGAATCTTTCTATTGCATGGCAATTAAATTAGAAATTAAAAATCTTTATAAAATAATTGGCGAGCgttcaaatatatcgaacaacgactttcaaaagaacaaattctggaaaaaactgcgc                                                                                                <  2:626702‑M1/150‑57 (MQ=255)                                                          gggTAATATATCGACATAGACAAATAAAGGAATCTTTCTATTGCATGGCAATTAAATTAGAAATTAAAAATCTTTATAAAATATTTGGCGAGCgttcaaatatatcgaacaaggactttcaaaagaacaaattctggaaaaaactgggct                                                                                               <  2:311123‑M1/150‑58 (MQ=255)                                                               atatatCGACATAGACAAATAAAGGAATCTTTCTATTGCATGGCAATTAAATTAGAAATTAAAAATCTTTATAAAATATTTGGCGAGCgttcaaatatatcgaacaaggactttcaaaagaacaaattctggaaaaaactgggctatcgct                                                                                         >  1:346784‑M1/1‑88 (MQ=255)                                                                   atCGACATAGACAAATAAAGGAATCTTTCTATTGCATGGCAATTAAATTAGAAATTAAAAATCTTTATAAAATATTTGGCGAGCgttcaaatatatcgaacaaggactttcaaaagaacaaattctggaaaaaactgggctatcgcttgg                                                                                      <  1:690671‑M1/150‑67 (MQ=255)                                                                   atCGACATAGACAAATAAAGGAATCTTTCTATTGCATGGCAATTAAATTAGAAATTAAAAATCTTTATAAAATATTTGGCGAGCgttcaaatatatcgaacaaggactttcaaaagaacaaattctggaaaaaactgggctatcgctt                                                                                        >  2:110877‑M1/1‑84 (MQ=255)                                                                      gACATAGACAAATAAAGGAATCTTTCTATTGCATGGCAATTAAATTAGAAATTAAAAATCTTTATAAAATATTTGGCGAGCgttcaaatatatcgaacaaggactttcaaaagaacaaattctggaaaaaactgggctatcgcttggcgta                                                                                  <  1:227709‑M1/151‑71 (MQ=255)                                                                       aCATAGACAAATAAAGGAATCTTTCTATTGCATGGCAATTAAATTAGAAATTAAAAATCTTTATAAAATATTTGGCGAGCgttcaaatatatcgaacaaggactttcaaaagaacaaattctggaaaaaactgggctatcgcttggcgta                                                                                  <  2:117131‑M1/150‑71 (MQ=255)                                                                           aGACAAATAAAGGAATCTTTCTATTGCATGGCAATTAAATTAGAAATTAAAAATCTTTATAAAATATTTGGCGAGCgttcaaatatatcgaacaaggactttcaaaagaacaaattctggaaaaaactgggctatcgcttggcgtaaa                                                                                >  2:551499‑M1/1‑76 (MQ=255)                                                                                aaTAAAGGAATCTTTCTATTGCATGGCAATTAAATTAGAAATTAAAAATCTTTATAAAATATTTGGCGAGCgttcaaatatatcgaacaaggactttcaaaagaacaaattctggaaaaaactgggctatcgcttggcgtaaaagacgcca                                                                        <  2:292775‑M1/151‑81 (MQ=255)                                                                                   aaaGGAATCTTTCTATTGCATGGCAATTAAATTAGAAATTAAAAATCTTTATAAAATATTTGGCGAGCgttcaaatatatcgaacaaggactttcaaaagaacaaattctggaaaaaactgggctatcgcttggcgtaaaagacgccagt                                                                      <  1:201356‑M1/150‑83 (MQ=255)                                                                                   aaaGGAATCTTTCTATTGCATGGCAATTAAATTAGAAATTAAAAATCTTTATAAAATATTTGGCGAGCgttcaaatatatcgaacaaggactttcaaaagaacaaattctggaaaaaactgggctatcgcttggcgtaaaagacgccagt                                                                      <  2:383393‑M1/150‑83 (MQ=255)                                                                                     aGGAATCTTTCTATTGCATGGCAATTAAATTAGAAATTAAAAATCTTTATAAAATATTTGGCGAGCgttcaaatatatcgaacaaggactttcaaaagaacaaattctggaaaaaactgggctatcgcttggcgtaaaagacgccagtctg                                                                   <  1:34694‑M1/151‑86 (MQ=255)                                                                                            tttCTATTGCATGGCAATTAAATTAGAAATTAAAAATCTTTATAAAATATTTGGCGAGCgttcaaatatatcgaacaaggactttcaaaagaacaaattctggaaaaaactgggctatcgcttggcgtaaaagacgccagtctggccattg                                                            >  2:456939‑M1/1‑59 (MQ=255)                                                                                             ttCTATTGCATGGCAATTAAATTAGAAATTAAAAATCTTTATAAAATATTTGGCGAGCgttcaaatatatcgaacaaggactttcaaaagaacaaattctggaaaaaactgggctatcgcttggcgtaaaagacgccagtctggccattg                                                            <  2:247937‑M1/150‑93 (MQ=255)                                                                                             ttCTATTGCATGGCAATTAAATTAGAAATTAAAAATCTTTATAAAATATTTGGCGAGCgttcaaatatatcgaacaaggactttcaaaagaacaaattctggaaaaaactgggctatcgcttggcgtaaaagacgccagtctggccatt                                                             <  1:348496‑M1/149‑92 (MQ=255)                                                                                                 aTTGCATGGCAATTAAATTAGAAATTAAAAATCTTTATAAAATATTTGGCGAGCgttcaaatatatcgaacaaggactttcaaaagaacaaattctggaaaaaactgggctatcgcttggcgtaaaagacgccagtctggccattgaagaa                                                       <  1:56337‑M1/151‑98 (MQ=255)                                                                                                    gCATGGCAATTAAATTAGAAATTAAAAATCTTTATAAAATATTTGGCGAGCgttcaaatatatcgaacaaggactttcaaaagaacaaattctggaaaaaactgggctatcgcttggcgtaaaagacgccagtctggccattgaagaaggc                                                    >  1:691512‑M1/1‑51 (MQ=255)                                                                                                         gCAATTAAATTAGAAATTAAAAATCTTTATAAAATATTTGGCGAGCgttcaaatatatcgaacaaggactttcaaaagaacaaattctggaaaaaactgggctatcgcttggcgtaaaagacgccagtctggccattgaagaaggcgagat                                               <  1:649832‑M1/151‑106 (MQ=255)                                                                                                                 attaGAAATTAAAAATCTTTATAAAATATTTGGCGAGCgttcaaatatatcgaacaaggactttcaaaagaacaaattctggaaaaaactgggctatcgcttggcgtaaaagacgccagtctggccattgaagaaggcgagatatttgtc                                        >  2:212891‑M1/1‑38 (MQ=255)                                                                                                                  ttaGAAATTAAAAATCTTTATAAAATATTTGGCGAGCgttcaaatatatcgaacaaggactttcaaaagaacaaattctggaaaaaactgggctatcgcttggcgtaaaagacgccagtctggccattgaagaaggcgagatatttgtca                                       >  1:351192‑M1/1‑37 (MQ=255)                                                                                                                      aaaTTAAAAATCTTTATAAAATATTTGGCGAGCgttcaaatatatcgaacaaggactttcaaaagaacaaattctggaaaaaactgggctatcgcttggcgtaaaagacgccagtctggccattgaagaaggcgagatatttgtcatcatg                                  >  2:280348‑M1/1‑33 (MQ=255)                                                                                                                        aTTAAAAATCTTTATAAAATATTTGGCGAGCgttcaaatatatcgaacaaggactttcaaaagaacaaattctggaaaaaactgggctatcgcttggcgtaaaagacgccagtctggccattgaagaaggcgagatatttgtcatcatggg                                >  1:585257‑M1/1‑31 (MQ=255)                                                                                                                                      taAAATATTTGGCGAGCgttcaaatatatcgaacaaggactttcaaaagaacaaattctggaaaaaactgggctatcgcttggcgtaaaagacgccagtctggccattgaagaaggcgagatatttgtcatcatgggattatccggctcgg                  >  1:230717‑M1/1‑17 (MQ=255)                                                                                                                                        aaaTATTTGGCGAGCgttcaaatatatcgaacaaggactttcaaaagaacaaattctggaaaaaactgggctatcgcttggcgtaaaagacgccagtctggccattgaagaaggcgagatatttgtcatcatgggattatccggctcgggt                >  1:18573‑M1/1‑15 (MQ=255)                                                                                                                                         aaTATTTGGCGAGCgttcaaatatatcgaacaaggactttcaaaagaacaaattctggaaaaaactgggctatcgcttggcgtaaaagacgccagtctggccattgaagaaggcgagatatttgtcatcatgggattatccggctcgggta               <  1:331046‑M1/151‑138 (MQ=255)                                                                                                                                          atgtttggcGAGCgttcaaatatatcgaacaaggactttcaaaagaacaaattctggaaaaaactgggctatcgcttggcgtaaaagacgccagtctggccattgaagaaggcgagatatttgtcatcatgggattatccggctcgggta               >  2:339754‑M1/4‑13 (MQ=255)                                                                                                                                            atTTGGCGAGCgttcaaatatatcgaacaaggactttcaaaagaacaaattctggaaaaaactgggctatcgcttggcgtaaaagacgccagtctggccattgaagaaggcgagatatttgtcatcatgggattatccggctcgggtaaat            >  1:30384‑M1/1‑11 (MQ=255)                                                                                                                                              ttGGCGAGCgttcaaatatatcgaacaaggactttcaaaagaacaaattctggaaaaaactgggctatcgcttggcgtaaaagacgccagtctggccattgaagaaggcgagatatttgtcatcatgggattatccggctcgggtaaatc           <  1:375325‑M1/150‑142 (MQ=255)                                                                                                                                                  cGAGCgttcaaatatatcgaacaaggactttcaaaagaacaaattctggaaaaaactgggctatcgcttggcgtaaaagacgccagtctggccattgaagaaggcgagatatttgtcatcatggtattatccggctcgggtaaatccaca       >  2:67408‑M1/1‑5 (MQ=255)                                                                                                                                                   gAGCgttcaaatatatcgaacaaggactttcaaaagaacaaattctggaaaaaactgggctatcgcttggcgtaaaagacgccagtctggccattgaagaaggcgagatatttgtcatcatgggattatccggctcgggtaaatccacaat     <  1:236210‑M1/151‑148 (MQ=255)                                                                                                                                                      cgttcaaatatatcgaacaaggactttcaaaagaacaaattctggaaaaaactgggctatcgcttggcgtaaaagacgccagtctggccattgaagaaggcgagatatttgtcatcatgggattatccggctcgggtaaatccacaatggt  >  1:200565‑M1/1‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                       CAGGGGTTGCCTCAGATTCTCAGTATGTTAGGGTAGAAAAAAGTGACTATTTCCATTGGGTAATATATCGACATAGACAAATAAAGGAATCTTTCTATTGCATGGCAATTAAATTAGAAATTAAAAATCTTTATAAAATATTTGGCGAGCATCCACAGCGAGCGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTAAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGTCATCATGGGATTATCCGGCTCGGGTA  >  NC_000913/2804714‑2805013

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