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	58 [0]	[0] 58	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)	55 (0.830)	48/286	0.7	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

CAGTATGTTAGGGTAGAAAAAAGTGACTATTTCCATTGGGTAATATATCGACATAGACAAATAAAGGAATCTTTCTATTGCATGGCAATTAAATTAGAAATTAAAAATCTTTATAAAATATTTGGCGAGCATCCACAGCGAGCGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTAAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGTCATCATGGGATTATCCGGCTCGGGTAAATCCACAATGGT  >  NC_000913/2804734‑2805026                                                                                                                                                |                                                                                                                                                      tgcctcagattctcagtatgttagggtagaaaaaagtgactatttccattgggtaatatatcgacatagacaaataaaggaatctttctattgcatggcaattaaattagaaattaaaaatctttataaaatatttggcgagcGTTCaaa                                                                                                                                                 <  1:176672‑M2/7‑1 (MQ=255) tgcctcagattctcagtatgttagggtagaaaaaagtgactatttccattgggtaatatatcgacatagacaaataaaggaatctttctattgcatggcaattaaattagaaattaaaaatctttataaaatatttggcgagcGTTCAAat                                                                                                                                                >  1:744902‑M2/144‑151 (MQ=255)      cagattctcagtatgttagggtagaaaaaagtgactatttccattgggtaatatatcgacatagacaaataaaggaatctttctattgcatggcaattaaattagaaattaaaaatctttataaaatatttggcgagcGTTCAAATATATc                                                                                                                                           >  1:663949‑M2/139‑151 (MQ=255)             tcagtatgttagggtagaaaaaagtgactatttccattgggtaatatatcgacatagacaaataaaggaatctttctattgcatggcaattaaattagaaattaaaaatctttataaaatatttggcgagcGTTCAAATATATCGAACa                                                                                                                                      >  1:89594‑M2/132‑149 (MQ=255)              cagtatgttagggtagaaaaaagtgactatttccattgggtaatatatcgacatagacaaataaaggaatctttctattgcatggcaattaaattagaaattaaaaatctttataaaatatttggcgagcGTTCAAATATATCGAACa                                                                                                                                      <  2:169654‑M2/18‑1 (MQ=255)              cagtatgttagggtagaaaaaagtgactatttccattgggtaatatatcgacatagacaaataaaggaatctttctattgcatggcaattaaattagaaattaaaaatctttataaaatatttggcgagcGTTCAAATATATCGAACa                                                                                                                                      <  2:900385‑M2/18‑1 (MQ=255)              cagtatgttagggtagaaaaaagtgactatttccattgggtaatatatcgacatagacaaataaaggaatctttctattgcatggcaattaaattagaaattaaaaatctttataaaatatttggcgagcGTTCAAATATATCGAACAAgg                                                                                                                                   <  2:489132‑M2/21‑1 (MQ=255)                  atgttagggtagaaaaaagtgactatttccattgggtaatatatcgacatagacaaataaaggaatctttctattgcatggcaattaaattagaaattaaaaatctttataaaatatttggcgagcGTTCAAATATATCGAACAAGGACtt                                                                                                                               >  1:297192‑M2/127‑151 (MQ=255)                     ttagggtagaaaaaagtgactatttccattgggtaatatatcgacatagacaaataaaggaatctttctattgcatggcaattaaattagaaattaaaaatctttataaaatatttggcgagcGTTCAAATATATCGAACAAGGACTTTc                                                                                                                             >  1:280445‑M2/124‑150 (MQ=255)                        gggtagaaaaaagtgactatttccattgggtaatatatcgacatagacaaataaaggaatctttctattgcatggcaattaaattagaaattaaaaatctttataaaatatttggcgagcGTTCAAATATATCGAACAAGGACTTTCaaaa                                                                                                                         <  2:22163‑M2/31‑1 (MQ=255)                            agaaaaaagtgactatttccattgggtaatatatcgacatagacaaataaaggaatctttctattgcatggcaattaaattagaaattaaaaatctttataaaatatttggcgagcGTTCAAATATATCGAACAAGGACTTTCAAAAGAAc                                                                                                                     <  1:727405‑M2/35‑1 (MQ=255)                             gaaaaaagtgactatttccattgggtaatatatcgacatagacaaataaaggaatctttctattgcatggcaattaaattagaaattaaaaatctttataaaatatttggcgagcGTTCAAATATATCGAACAAGGACTTTCAAAAGAAc                                                                                                                     >  1:540163‑M2/116‑150 (MQ=255)                              aaaaaagtgactatttccattgggtaatatatcgacatagacaaataaaggaatctttctattgcatggcaattaaattagaaattaaaaatctttataaaatatttggcgagcGTTCAAATATATCGAACAAGGACTTTCAAAAGAACaa                                                                                                                   <  1:544622‑M2/37‑1 (MQ=255)                              aaaaaagtgactatttccattgggtaatatatcgacatagacaaataaaggaatctttctattgcatggcaattaaattagaaattaaaaatctttataaaatatttggcgagcGTTCAAATATATCGAACAAGGACTTTCAAAAGAACaa                                                                                                                   <  2:427684‑M2/37‑1 (MQ=255)                                aaaagtgactatttccattgggtaatatatcgacatagacaaataaaggaatctttctattgcatggcaattaaattagaaattaaaaatctttataaaatatttggcgagcGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAAt                                                                                                                 >  2:85696‑M2/113‑151 (MQ=255)                                 aaagtgactatttccattgggtaatatatcgacatagacaaataaaggaatctttctattgcatggcaattaaattagaaattaaaaatctttataaaatatttggcgagcGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAAtt                                                                                                                >  1:362734‑M2/112‑151 (MQ=255)                                  aagtgactatttccattgggtaatatatcgacatagacaaataaaggaatctttctattgcatggcaattaaattagaaattaaaaatctttataaaatatttggcgagcGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTc                                                                                                               <  1:370912‑M2/41‑1 (MQ=255)                                  aagtgactatttccattgggtaatatatcgacatagacaaataaaggaatctttctattgcatggcaattaaattagaaattaaaaatctttataaaatatttggcgagcGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTc                                                                                                               <  1:74259‑M2/41‑1 (MQ=255)                                      gactatttccattgggtaatatatcgacatagacaaataaaggaatctttctattgcatggcaattaaattagaaattaaaaatctttataaaatatttggcgagcGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGa                                                                                                           <  2:887122‑M2/45‑1 (MQ=255)                                          atttccattgggtaatatatcgacatagacaagcaaaggaatctttctattgcatggcaattaaattagaaattaaaaatctttataaaatatttggcgagcGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGaaa                                                                                                         <  1:454393‑M2/47‑1 (MQ=255)                                          atttccattgggtaatatatcgacatagacaaataaaggaatctttctattgcatggcaattaaattagaaattaaaaatctttataaaatatttggcgagcGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGaaaa                                                                                                        >  2:924904‑M2/103‑150 (MQ=255)                                            ttncattgggnnntnnatnnncatagacaaataaaggaatctttntattgcatngcaattaaattagaaattaaaaatctttatanaanatttggcgagcGTTCNAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGaana                                                                                                        <  1:192561‑M2/48‑1 (MQ=255)                                             tccattgggtaatatatcgacatagacaaataaaggaatctttctattgcatggcaattaaattagaaattaaaaatctttataaaatatttggcgagcGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACt                                                                                                    >  2:722570‑M2/100‑151 (MQ=255)                                              ccattgggtaatatatcgacatagacaaataaaggaatctttctattgcatggcaattaaattagaaattaaaaatctttataaaatatttggcgagcGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGACAAAACt                                                                                                    >  2:783468‑M2/99‑150 (MQ=255)                                              ccattgggtaatatatcgacatagacaaataaaggaatctttctattgcatggcaattaaattagaaattaaaaatctttataaaatatttggcgagcGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACt                                                                                                    <  1:596203‑M2/52‑1 (MQ=255)                                               cattgggtaatatatcgacatagacaaataaaggaatctttctattgcatggcaattaaattagaaattaaaaatctttataaaatatttggcgagcGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTgg                                                                                                  >  1:161881‑M2/98‑151 (MQ=255)                                                attgggtaatatatcgacatagacaaataaaggaatctttctattgcatggcaattaaattagaaattaaaaatctttataaaatatttggcgagcGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTggg                                                                                                 <  1:1006635‑M2/55‑1 (MQ=255)                                                attgggtaatatatcgacatagacaaataaaggaatctttctattgcatggcaattaaattagaaattaaaaatcttaataaaatatttggcgagcGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTggg                                                                                                 <  2:574432‑M2/55‑1 (MQ=255)                                                       aatatatcgacatagacaaataaaggaatctttctattgcatggcaattaaattagaaattaaaaatctttataaaatatttggcgagcGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCg                                                                                           <  2:980820‑M2/61‑1 (MQ=255)                                                          atatcgacatagacaaataaaggaatctttctattgcatggcaattaaattagaaattaaaaatctttataaaatatttggcgagcGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTg                                                                                       >  1:989354‑M2/87‑151 (MQ=255)                                                            atcgacatagacaaataaaggaatctttctattgcatggcaattaaattagaaattaaaaatctttataaaatatttggcgagcGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGc                                                                                     <  2:516877‑M2/67‑1 (MQ=255)                                                            atcgacatagacaaataaaggaatctttctattgcatggcaattaaattagaaattaaaaatctttataaaatatttggcgagcGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGc                                                                                     <  2:398264‑M2/67‑1 (MQ=255)                                                               gacatagacaaataaaggaatctttctattgcatggcaattaaattagaaattaaaaatctttataaaatatttggcgagcGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTa                                                                                  <  1:1024621‑M2/70‑1 (MQ=255)                                                                 catagacaaataaaggaatctttctattgcatggcaattaaattagaaattaaaaatctttataaaatatttggcgagcGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTaaa                                                                                <  2:529141‑M2/72‑1 (MQ=255)                                                                      acaaataaaggaatctttctattgcatggcaattaaattagaaattaaaaatctttataaaatatttggcgagcGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTAAAAGACg                                                                           >  2:477301‑M2/75‑151 (MQ=255)                                                                       caaataaaggaatctttctattgcatggcaattaaattagaaattaaaaatctttataaaatatttggcgagcGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTAAAAGACGc                                                                          >  1:659202‑M2/74‑151 (MQ=255)                                                                        aaataaaggaatctttctattgcatggcaattaaattagaaattaaaaatctttataaaatatttggcgagcGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTAAAAGACGcc                                                                         >  2:132429‑M2/73‑151 (MQ=255)                                                                                      ttctattgcatggcaattaaattagaaattaaaaatctttataaaatatttggcgagcGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTAAAAGACGCCAGTCTGGCCAt                                                              <  2:228123‑M2/90‑1 (MQ=255)                                                                                             gcatggcaattaaattagaaattaaaaatctttataaaatatttggcgagcGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTAAAAGACGCCAGTCTGGCCATTGAAGAAgg                                                     <  1:294066‑M2/99‑1 (MQ=255)                                                                                              catggcaattaaattagaaattaaaaatctttataaaatatttggcgagcGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTAAAAGACGCCAGTCTGGCCATTGAAGAAGGCg                                                   >  1:683339‑M2/51‑151 (MQ=255)                                                                                                  gcaattaaattagaaattaaaaatctttataaaatatttggcgagcGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTAAAAGACGCCAGTCTGGCCATTGAAGAAGGCg                                                   >  2:741909‑M2/47‑147 (MQ=255)                                                                                                      ttaaattagaaattaaaaatctttataaaatatttggcgagcGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTAAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATAttt                                           >  2:636328‑M2/43‑151 (MQ=255)                                                                                                       taaattagaaattaaaaatctttataaaatatttggcgagcGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTAAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATAttt                                           <  1:308150‑M2/109‑1 (MQ=255)                                                                                                       taaattagaaattaaaaatctttataaaatatttggcgagcGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTAAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATAttt                                           <  1:161483‑M2/109‑1 (MQ=255)                                                                                                          attagaaattaaaaatctttataaaatatttggcgagcGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTAAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGtca                                       <  1:353219‑M2/113‑1 (MQ=255)                                                                                                          attagaaattaaaaatctttataaaatatttggcgagcGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTAAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGtc                                        >  2:98467‑M2/39‑150 (MQ=255)                                                                                                             agaaattaaaaatctttataaaatatttggcgagcGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTAAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGtcatc                                     >  1:697186‑M2/36‑150 (MQ=255)                                                                                                                   taaaaatctttataaaatatttggcgagcGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTAAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGTCATCATGGGAt                              >  2:97502‑M2/30‑151 (MQ=255)                                                                                                                             tataaaatatttggcgagcGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTAAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGTCATCATGGGATTATCCGGCt                     <  2:290592‑M2/131‑1 (MQ=255)                                                                                                                                aaaatatttggcgagcGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTAAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGTCATCATGGGATTATCCGGCTCgg                  <  2:64148‑M2/134‑1 (MQ=255)                                                                                                                                    tatttggcgagcGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTAAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGTCATCATGGGATTATCCGGCTCGGGTaa              >  1:1071409‑M2/13‑150 (MQ=255)                                                                                                                                     atttggcgagcGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTAAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGTCATCATGGGATTATCCGGCTCGGGTaaa             <  1:238031‑M2/139‑1 (MQ=255)                                                                                                                                     atttggcgagcGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTAAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGTCATCATGGGATTATCCGGCTCGGGTAAAt            >  2:415092‑M2/12‑151 (MQ=255)                                                                                                                                     atttggcgagcGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTAAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGTCATCATGGGATTATCCGGCTCGGGTAAAt            >  1:525627‑M2/12‑151 (MQ=255)                                                                                                                                       ttggcgagcGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTAAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGTCATCATGGGATTATCCGGCTCGGGTAAATcc          >  2:198142‑M2/10‑151 (MQ=255)                                                                                                                                          gcgagcGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTAAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGTCATCATGGGATTATCCGGCTCGGGTAAATCcac        >  2:837425‑M2/7‑150 (MQ=255)                                                                                                                                               cGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTAAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGTCATCATGGGATTATCCGGCTCGGGTAAATCCACAATgg   >  1:236586‑M2/2‑150 (MQ=255)                                                                                                                                               cGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTAAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGTCATCATGGGATTATCCGGCTCGGGTAAATCCACAATGGt  <  1:742372‑M2/150‑1 (MQ=255)                                                                                                                                                |                                                                                                                                                      CAGTATGTTAGGGTAGAAAAAAGTGACTATTTCCATTGGGTAATATATCGACATAGACAAATAAAGGAATCTTTCTATTGCATGGCAATTAAATTAGAAATTAAAAATCTTTATAAAATATTTGGCGAGCATCCACAGCGAGCGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTAAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGTCATCATGGGATTATCCGGCTCGGGTAAATCCACAATGGT  >  NC_000913/2804734‑2805026

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