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

TGCCTCAGATTCTCAGTATGTTAGGGTAGAAAAAAGTGACTATTTCCATTGGGTAATATATCGACATAGACAAATAAAGGAATCTTTCTATTGCATGGCAATTAAATTAGAAATTAAAAATCTTTATAAAATATTTGGCGAGCATCCACAGCGAGCGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTAAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGTCATCATGGGATTATCCGGCTCGGGTA  >  NC_000913/2804721‑2805013                                                                                                                                               |                                                                                                                                                       tGCCTCAGATTCTCAGTATGTTAGGGTAGAAAAAAGTGACTATTTCCATTGGGTAATATATCGACATAGACAAATAAAGGAATCTTTCTATTGCATGGCAATTAAATTAGAAATTAAAAATCTTTATAAAATATTTGGCGAGCgttcaaat                                                                                                                                                >  1:744902‑M1/1‑143 (MQ=255) tGCCTCAGATTCTCAGTATGTTAGGGTAGAAAAAAGTGACTATTTCCATTGGGTAATATATCGACATAGACAAATAAAGGAATCTTTCTATTGCATGGCAATTAAATTAGAAATTAAAAATCTTTATAAAATATTTGGCGAGCgttcaaa                                                                                                                                                 <  1:176672‑M1/150‑8 (MQ=255)      cAGATTCTCAGTATGTTAGGGTAGAAAAAAGTGACTATTTCCATTGGGTAATATATCGACATAGACAAATAAAGGAATCTTTCTATTGCATGGCAATTAAATTAGAAATTAAAAATCTTTATAAAATATTTGGCGAGCgttcaaatatatc                                                                                                                                           >  1:663949‑M1/1‑138 (MQ=255)             tcAGTATGTTAGGGTAGAAAAAAGTGACTATTTCCATTGGGTAATATATCGACATAGACAAATAAAGGAATCTTTCTATTGCATGGCAATTAAATTAGAAATTAAAAATCTTTATAAAATATTTGGCGAGCgttcaaatatatcgaaca                                                                                                                                      >  1:89594‑M1/1‑131 (MQ=255)              cAGTATGTTAGGGTAGAAAAAAGTGACTATTTCCATTGGGTAATATATCGACATAGACAAATAAAGGAATCTTTCTATTGCATGGCAATTAAATTAGAAATTAAAAATCTTTATAAAATATTTGGCGAGCgttcaaatatatcgaacaagg                                                                                                                                   <  2:489132‑M1/151‑22 (MQ=255)              cAGTATGTTAGGGTAGAAAAAAGTGACTATTTCCATTGGGTAATATATCGACATAGACAAATAAAGGAATCTTTCTATTGCATGGCAATTAAATTAGAAATTAAAAATCTTTATAAAATATTTGGCGAGCgttcaaatatatcgaaca                                                                                                                                      <  2:169654‑M1/148‑19 (MQ=255)              cAGTATGTTAGGGTAGAAAAAAGTGACTATTTCCATTGGGTAATATATCGACATAGACAAATAAAGGAATCTTTCTATTGCATGGCAATTAAATTAGAAATTAAAAATCTTTATAAAATATTTGGCGAGCgttcaaatatatcgaaca                                                                                                                                      <  2:900385‑M1/148‑19 (MQ=255)                  aTGTTAGGGTAGAAAAAAGTGACTATTTCCATTGGGTAATATATCGACATAGACAAATAAAGGAATCTTTCTATTGCATGGCAATTAAATTAGAAATTAAAAATCTTTATAAAATATTTGGCGAGCgttcaaatatatcgaacaaggactt                                                                                                                               >  1:297192‑M1/1‑126 (MQ=255)                     ttAGGGTAGAAAAAAGTGACTATTTCCATTGGGTAATATATCGACATAGACAAATAAAGGAATCTTTCTATTGCATGGCAATTAAATTAGAAATTAAAAATCTTTATAAAATATTTGGCGAGCgttcaaatatatcgaacaaggactttc                                                                                                                             >  1:280445‑M1/1‑123 (MQ=255)                        gggTAGAAAAAAGTGACTATTTCCATTGGGTAATATATCGACATAGACAAATAAAGGAATCTTTCTATTGCATGGCAATTAAATTAGAAATTAAAAATCTTTATAAAATATTTGGCGAGCgttcaaatatatcgaacaaggactttcaaaa                                                                                                                         <  2:22163‑M1/151‑32 (MQ=255)                            aGAAAAAAGTGACTATTTCCATTGGGTAATATATCGACATAGACAAATAAAGGAATCTTTCTATTGCATGGCAATTAAATTAGAAATTAAAAATCTTTATAAAATATTTGGCGAGCgttcaaatatatcgaacaaggactttcaaaagaac                                                                                                                     <  1:727405‑M1/151‑36 (MQ=255)                             gAAAAAAGTGACTATTTCCATTGGGTAATATATCGACATAGACAAATAAAGGAATCTTTCTATTGCATGGCAATTAAATTAGAAATTAAAAATCTTTATAAAATATTTGGCGAGCgttcaaatatatcgaacaaggactttcaaaagaac                                                                                                                     >  1:540163‑M1/1‑115 (MQ=255)                              aaaaaaGTGACTATTTCCATTGGGTAATATATCGACATAGACAAATAAAGGAATCTTTCTATTGCATGGCAATTAAATTAGAAATTAAAAATCTTTATAAAATATTTGGCGAGCgttcaaatatatcgaacaaggactttcaaaagaacaa                                                                                                                   <  1:544622‑M1/151‑38 (MQ=255)                              aaaaaaGTGACTATTTCCATTGGGTAATATATCGACATAGACAAATAAAGGAATCTTTCTATTGCATGGCAATTAAATTAGAAATTAAAAATCTTTATAAAATATTTGGCGAGCgttcaaatatatcgaacaaggactttcaaaagaacaa                                                                                                                   <  2:427684‑M1/151‑38 (MQ=255)                                aaaaGTGACTATTTCCATTGGGTAATATATCGACATAGACAAATAAAGGAATCTTTCTATTGCATGGCAATTAAATTAGAAATTAAAAATCTTTATAAAATATTTGGCGAGCgttcaaatatatcgaacaaggactttcaaaagaacaaat                                                                                                                 >  2:85696‑M1/1‑112 (MQ=255)                                 aaaGTGACTATTTCCATTGGGTAATATATCGACATAGACAAATAAAGGAATCTTTCTATTGCATGGCAATTAAATTAGAAATTAAAAATCTTTATAAAATATTTGGCGAGCgttcaaatatatcgaacaaggactttcaaaagaacaaatt                                                                                                                >  1:362734‑M1/1‑111 (MQ=255)                                  aaGTGACTATTTCCATTGGGTAATATATCGACATAGACAAATAAAGGAATCTTTCTATTGCATGGCAATTAAATTAGAAATTAAAAATCTTTATAAAATATTTGGCGAGCgttcaaatatatcgaacaaggactttcaaaagaacaaattc                                                                                                               <  1:370912‑M1/151‑42 (MQ=255)                                  aaGTGACTATTTCCATTGGGTAATATATCGACATAGACAAATAAAGGAATCTTTCTATTGCATGGCAATTAAATTAGAAATTAAAAATCTTTATAAAATATTTGGCGAGCgttcaaatatatcgaacaaggactttcaaaagaacaaattc                                                                                                               <  1:74259‑M1/151‑42 (MQ=255)                                      gACTATTTCCATTGGGTAATATATCGACATAGACAAATAAAGGAATCTTTCTATTGCATGGCAATTAAATTAGAAATTAAAAATCTTTATAAAATATTTGGCGAGCgttcaaatatatcgaacaaggactttcaaaagaacaaattctgga                                                                                                           <  2:887122‑M1/151‑46 (MQ=255)                                          aTTTCCATTGGGTAATATATCGACATAGACAAGCAAAGGAATCTTTCTATTGCATGGCAATTAAATTAGAAATTAAAAATCTTTATAAAATATTTGGCGAGCgttcaaatatatcgaacaaggactttcaaaagaacaaattctggaaa                                                                                                         <  1:454393‑M1/149‑48 (MQ=255)                                          aTTTCCATTGGGTAATATATCGACATAGACAAATAAAGGAATCTTTCTATTGCATGGCAATTAAATTAGAAATTAAAAATCTTTATAAAATATTTGGCGAGCgttcaaatatatcgaacaaggactttcaaaagaacaaattctggaaaa                                                                                                        >  2:924904‑M1/1‑102 (MQ=255)                                            ttNCATTGGGNNNTNNATNNNCATAGACAAATAAAGGAATCTTTNTATTGCATNGCAATTAAATTAGAAATTAAAAATCTTTATANAANATTTGGCGAGCgttcnaatatatcgaacaaggactttcaaaagaacaaattctggaana                                                                                                        <  1:192561‑M1/148‑49 (MQ=255)                                             tCCATTGGGTAATATATCGACATAGACAAATAAAGGAATCTTTCTATTGCATGGCAATTAAATTAGAAATTAAAAATCTTTATAAAATATTTGGCGAGCgttcaaatatatcgaacaaggactttcaaaagaacaaattctggaaaaaact                                                                                                    >  2:722570‑M1/1‑99 (MQ=255)                                              ccATTGGGTAATATATCGACATAGACAAATAAAGGAATCTTTCTATTGCATGGCAATTAAATTAGAAATTAAAAATCTTTATAAAATATTTGGCGAGCgttcaaatatatcgaacaaggactttcaaaagaacaaattctggacaaaact                                                                                                    >  2:783468‑M1/1‑98 (MQ=255)                                              ccATTGGGTAATATATCGACATAGACAAATAAAGGAATCTTTCTATTGCATGGCAATTAAATTAGAAATTAAAAATCTTTATAAAATATTTGGCGAGCgttcaaatatatcgaacaaggactttcaaaagaacaaattctggaaaaaact                                                                                                    <  1:596203‑M1/150‑53 (MQ=255)                                               cATTGGGTAATATATCGACATAGACAAATAAAGGAATCTTTCTATTGCATGGCAATTAAATTAGAAATTAAAAATCTTTATAAAATATTTGGCGAGCgttcaaatatatcgaacaaggactttcaaaagaacaaattctggaaaaaactgg                                                                                                  >  1:161881‑M1/1‑97 (MQ=255)                                                aTTGGGTAATATATCGACATAGACAAATAAAGGAATCTTTCTATTGCATGGCAATTAAATTAGAAATTAAAAATCTTTATAAAATATTTGGCGAGCgttcaaatatatcgaacaaggactttcaaaagaacaaattctggaaaaaactggg                                                                                                 <  1:1006635‑M1/151‑56 (MQ=255)                                                aTTGGGTAATATATCGACATAGACAAATAAAGGAATCTTTCTATTGCATGGCAATTAAATTAGAAATTAAAAATCTTAATAAAATATTTGGCGAGCgttcaaatatatcgaacaaggactttcaaaagaacaaattctggaaaaaactggg                                                                                                 <  2:574432‑M1/151‑56 (MQ=255)                                                       aaTATATCGACATAGACAAATAAAGGAATCTTTCTATTGCATGGCAATTAAATTAGAAATTAAAAATCTTTATAAAATATTTGGCGAGCgttcaaatatatcgaacaaggactttcaaaagaacaaattctggaaaaaactgggctatcg                                                                                           <  2:980820‑M1/150‑62 (MQ=255)                                                          atatCGACATAGACAAATAAAGGAATCTTTCTATTGCATGGCAATTAAATTAGAAATTAAAAATCTTTATAAAATATTTGGCGAGCgttcaaatatatcgaacaaggactttcaaaagaacaaattctggaaaaaactgggctatcgcttg                                                                                       >  1:989354‑M1/1‑86 (MQ=255)                                                            atCGACATAGACAAATAAAGGAATCTTTCTATTGCATGGCAATTAAATTAGAAATTAAAAATCTTTATAAAATATTTGGCGAGCgttcaaatatatcgaacaaggactttcaaaagaacaaattctggaaaaaactgggctatcgcttggc                                                                                     <  2:516877‑M1/151‑68 (MQ=255)                                                            atCGACATAGACAAATAAAGGAATCTTTCTATTGCATGGCAATTAAATTAGAAATTAAAAATCTTTATAAAATATTTGGCGAGCgttcaaatatatcgaacaaggactttcaaaagaacaaattctggaaaaaactgggctatcgcttggc                                                                                     <  2:398264‑M1/151‑68 (MQ=255)                                                               gACATAGACAAATAAAGGAATCTTTCTATTGCATGGCAATTAAATTAGAAATTAAAAATCTTTATAAAATATTTGGCGAGCgttcaaatatatcgaacaaggactttcaaaagaacaaattctggaaaaaactgggctatcgcttggcgta                                                                                  <  1:1024621‑M1/151‑71 (MQ=255)                                                                 cATAGACAAATAAAGGAATCTTTCTATTGCATGGCAATTAAATTAGAAATTAAAAATCTTTATAAAATATTTGGCGAGCgttcaaatatatcgaacaaggactttcaaaagaacaaattctggaaaaaactgggctatcgcttggcgtaaa                                                                                <  2:529141‑M1/151‑73 (MQ=255)                                                                      aCAAATAAAGGAATCTTTCTATTGCATGGCAATTAAATTAGAAATTAAAAATCTTTATAAAATATTTGGCGAGCgttcaaatatatcgaacaaggactttcaaaagaacaaattctggaaaaaactgggctatcgcttggcgtaaaagacg                                                                           >  2:477301‑M1/1‑74 (MQ=255)                                                                       cAAATAAAGGAATCTTTCTATTGCATGGCAATTAAATTAGAAATTAAAAATCTTTATAAAATATTTGGCGAGCgttcaaatatatcgaacaaggactttcaaaagaacaaattctggaaaaaactgggctatcgcttggcgtaaaagacgc                                                                          >  1:659202‑M1/1‑73 (MQ=255)                                                                        aaaTAAAGGAATCTTTCTATTGCATGGCAATTAAATTAGAAATTAAAAATCTTTATAAAATATTTGGCGAGCgttcaaatatatcgaacaaggactttcaaaagaacaaattctggaaaaaactgggctatcgcttggcgtaaaagacgcc                                                                         >  2:132429‑M1/1‑72 (MQ=255)                                                                                      ttCTATTGCATGGCAATTAAATTAGAAATTAAAAATCTTTATAAAATATTTGGCGAGCgttcaaatatatcgaacaaggactttcaaaagaacaaattctggaaaaaactgggctatcgcttggcgtaaaagacgccagtctggccat                                                              <  2:228123‑M1/148‑91 (MQ=255)                                                                                             gCATGGCAATTAAATTAGAAATTAAAAATCTTTATAAAATATTTGGCGAGCgttcaaatatatcgaacaaggactttcaaaagaacaaattctggaaaaaactgggctatcgcttggcgtaaaagacgccagtctggccattgaagaagg                                                     <  1:294066‑M1/150‑100 (MQ=255)                                                                                              cATGGCAATTAAATTAGAAATTAAAAATCTTTATAAAATATTTGGCGAGCgttcaaatatatcgaacaaggactttcaaaagaacaaattctggaaaaaactgggctatcgcttggcgtaaaagacgccagtctggccattgaagaaggcg                                                   >  1:683339‑M1/1‑50 (MQ=255)                                                                                                  gCAATTAAATTAGAAATTAAAAATCTTTATAAAATATTTGGCGAGCgttcaaatatatcgaacaaggactttcaaaagaacaaattctggaaaaaactgggctatcgcttggcgtaaaagacgccagtctggccattgaagaaggcg                                                   >  2:741909‑M1/1‑46 (MQ=255)                                                                                                      ttaaattaGAAATTAAAAATCTTTATAAAATATTTGGCGAGCgttcaaatatatcgaacaaggactttcaaaagaacaaattctggaaaaaactgggctatcgcttggcgtaaaagacgccagtctggccattgaagaaggcgagatattt                                           >  2:636328‑M1/1‑42 (MQ=255)                                                                                                       taaattaGAAATTAAAAATCTTTATAAAATATTTGGCGAGCgttcaaatatatcgaacaaggactttcaaaagaacaaattctggaaaaaactgggctatcgcttggcgtaaaagacgccagtctggccattgaagaaggcgagatattt                                           <  1:308150‑M1/150‑110 (MQ=255)                                                                                                       taaattaGAAATTAAAAATCTTTATAAAATATTTGGCGAGCgttcaaatatatcgaacaaggactttcaaaagaacaaattctggaaaaaactgggctatcgcttggcgtaaaagacgccagtctggccattgaagaaggcgagatattt                                           <  1:161483‑M1/150‑110 (MQ=255)                                                                                                          attaGAAATTAAAAATCTTTATAAAATATTTGGCGAGCgttcaaatatatcgaacaaggactttcaaaagaacaaattctggaaaaaactgggctatcgcttggcgtaaaagacgccagtctggccattgaagaaggcgagatatttgtca                                       <  1:353219‑M1/151‑114 (MQ=255)                                                                                                          attaGAAATTAAAAATCTTTATAAAATATTTGGCGAGCgttcaaatatatcgaacaaggactttcaaaagaacaaattctggaaaaaactgggctatcgcttggcgtaaaagacgccagtctggccattgaagaaggcgagatatttgtc                                        >  2:98467‑M1/1‑38 (MQ=255)                                                                                                             aGAAATTAAAAATCTTTATAAAATATTTGGCGAGCgttcaaatatatcgaacaaggactttcaaaagaacaaattctggaaaaaactgggctatcgcttggcgtaaaagacgccagtctggccattgaagaaggcgagatatttgtcatc                                     >  1:697186‑M1/1‑35 (MQ=255)                                                                                                                   tAAAAATCTTTATAAAATATTTGGCGAGCgttcaaatatatcgaacaaggactttcaaaagaacaaattctggaaaaaactgggctatcgcttggcgtaaaagacgccagtctggccattgaagaaggcgagatatttgtcatcatgggat                              >  2:97502‑M1/1‑29 (MQ=255)                                                                                                                             tataAAATATTTGGCGAGCgttcaaatatatcgaacaaggactttcaaaagaacaaattctggaaaaaactgggctatcgcttggcgtaaaagacgccagtctggccattgaagaaggcgagatatttgtcatcatgggattatccggct                     <  2:290592‑M1/150‑132 (MQ=255)                                                                                                                                aaaaTATTTGGCGAGCgttcaaatatatcgaacaaggactttcaaaagaacaaattctggaaaaaactgggctatcgcttggcgtaaaagacgccagtctggccattgaagaaggcgagatatttgtcatcatgggattatccggctcgg                  <  2:64148‑M1/150‑135 (MQ=255)                                                                                                                                    tatTTGGCGAGCgttcaaatatatcgaacaaggactttcaaaagaacaaattctggaaaaaactgggctatcgcttggcgtaaaagacgccagtctggccattgaagaaggcgagatatttgtcatcatgggattatccggctcgggtaa              >  1:1071409‑M1/1‑12 (MQ=255)                                                                                                                                     atTTGGCGAGCgttcaaatatatcgaacaaggactttcaaaagaacaaattctggaaaaaactgggctatcgcttggcgtaaaagacgccagtctggccattgaagaaggcgagatatttgtcatcatgggattatccggctcgggtaaat            >  2:415092‑M1/1‑11 (MQ=255)                                                                                                                                     atTTGGCGAGCgttcaaatatatcgaacaaggactttcaaaagaacaaattctggaaaaaactgggctatcgcttggcgtaaaagacgccagtctggccattgaagaaggcgagatatttgtcatcatgggattatccggctcgggtaaat            >  1:525627‑M1/1‑11 (MQ=255)                                                                                                                                     atTTGGCGAGCgttcaaatatatcgaacaaggactttcaaaagaacaaattctggaaaaaactgggctatcgcttggcgtaaaagacgccagtctggccattgaagaaggcgagatatttgtcatcatgggattatccggctcgggtaaa             <  1:238031‑M1/150‑140 (MQ=255)                                                                                                                                       ttGGCGAGCgttcaaatatatcgaacaaggactttcaaaagaacaaattctggaaaaaactgggctatcgcttggcgtaaaagacgccagtctggccattgaagaaggcgagatatttgtcatcatgggattatccggctcgggtaaatcc          >  2:198142‑M1/1‑9 (MQ=255)                                                                                                                                          gCGAGCgttcaaatatatcgaacaaggactttcaaaagaacaaattctggaaaaaactgggctatcgcttggcgtaaaagacgccagtctggccattgaagaaggcgagatatttgtcatcatgggattatccggctcgggtaaatccac        >  2:837425‑M1/1‑6 (MQ=255)                                                                                                                                               cgttcaaatatatcgaacaaggactttcaaaagaacaaattctggaaaaaactgggctatcgcttggcgtaaaagacgccagtctggccattgaagaaggcgagatatttgtcatcatgggattatccggctcgggtaaatccacaatggt  <  1:742372‑M1/151‑151 (MQ=255)                                                                                                                                               cgttcaaatatatcgaacaaggactttcaaaagaacaaattctggaaaaaactgggctatcgcttggcgtaaaagacgccagtctggccattgaagaaggcgagatatttgtcatcatgggattatccggctcgggtaaatccacaatgg   >  1:236586‑M1/1‑1 (MQ=255)                                                                                                                                               |                                                                                                                                                       TGCCTCAGATTCTCAGTATGTTAGGGTAGAAAAAAGTGACTATTTCCATTGGGTAATATATCGACATAGACAAATAAAGGAATCTTTCTATTGCATGGCAATTAAATTAGAAATTAAAAATCTTTATAAAATATTTGGCGAGCATCCACAGCGAGCGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTAAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGTCATCATGGGATTATCCGGCTCGGGTA  >  NC_000913/2804721‑2805013

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: ‑