Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,018,030	Δ1 bp	coding (698/1428 nt)	rmuC →	putative recombination limiting protein RmuC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,018,027	0	A	.	100.0%	186.2 / NA	44	coding (695/1428 nt)	rmuC	putative recombination limiting protein RmuC
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base . (19/25);  total (19/25)

CGAGGTAGTATTGACGCGGGTGCTGGAGGCTTCCGGTCTGCGTGAAGGGTATGAATATGAAACCCAGGTCAGCATCGAAAATGACGCCCGCTCGCGGATGCAGCCGGATGTCATCGTGCGCCTGCCGCAGGGAAAAGATGTGGTGATCGACG  >  NC_000913/4017950‑4018101                                                                              |                                                                           cGAGGTAGTATTGACGCGGGTGCTGGAGGCTTCCGGTCTGCGTGAAGGGTATGAATATGAAACCCAGGTCAGCATCGaaat                                                                         <  1:806944/81‑2 (MQ=255)   aGGTAGTATTGACGCGGGTGCTGGAGGCTTCCGGTCTGCGTGAAGGGTATGAATATGAAACCCAGGTCAGCATCGaaatga                                                                       <  1:806945/81‑4 (MQ=255)        gtaTTGACGCGGGTGCTGGAGGCTTCCGGTCTGCGTGAAGGGTATGAATATGAAACCCAGGTCAGCATCG‑AAATGACGccc                                                                 <  2:806947/81‑1 (MQ=255)        gtaTTGACGCGGGTGCTGGAGGCTTCCGGTCTGCGTGAAGGGTATGAATATGAAACCCAGGTCAGCATCG‑AAATGACGcc                                                                  <  2:806946/80‑1 (MQ=255)         taTTGACGCGGGTGCTGGAGGCTTCCGGTCTGCGTGAAGGGTATGAATATGAAACCCAGGTCAGCATCG‑AAATGACGCCCg                                                                <  1:806948/81‑1 (MQ=255)          aTTGACGCGGGTGCTGGAGGCTTCCGGTCTGCGTGAAGGGTATGAATATGAAACCCAGGTCAGCATCG‑AAATGACGCCCGc                                                               >  2:807031/1‑81 (MQ=255)          aTTGACGCGGGTGCTGGAGGCTTCCGGTCTGCGTGAAGGGTATGAATATGAAACCCAGGTCAGCATCG‑AAATGACGCCCGc                                                               >  2:807005/1‑81 (MQ=255)          aTTGACGCGGGTGCTGGAGGCTTCCGGTCTGCGTGAAGGGTATGAATATGAAACCCAGGTCAGCATCG‑AAATGACGCCCGc                                                               >  1:807023/1‑81 (MQ=255)           ttGACGCGGGTGCTGGAGGCTTCCGGTCTGCGTGAAGGGTATGAATATGAAACCCAGGTCAGCATCG‑AAATGACGCCCGCt                                                              <  2:806950/81‑1 (MQ=255)           ttGACGCGGGTGCTGGAGGCTTCCGGTCTGCGTGAAGGGTATGAATATGAAACCCAGGTCAGCATCG‑AAATGACGCCCGCt                                                              <  2:806951/81‑1 (MQ=255)           ttGACGCGGGTGCTGGAGGCTTCCGGTCTGCGTAAAGGGTATGAATATGAAACCCAGGTCAGCATCG‑AAATGACGCCCGc                                                               <  1:806949/80‑1 (MQ=255)             gACGCGGGTGCTGGAGGCTTCCGGTCTGCGTGAAGGGTATGAATATGAAACCCAGGTCAGCATCG‑AAATGACGCCCGCTc                                                             <  1:806952/80‑1 (MQ=255)             gACGCGGGTGCTGGAGGCTTCCGGTCTGCGTGAAGGGTATGAATATGAAACCCAGGTCAGCATCG‑AAATGACGCCCGCTc                                                             <  2:806953/80‑1 (MQ=255)              aCGCGGGTGCTGGAGGCTTCCGGTCTGCGTGAAGGGTATGAATATGAAACCCAGGTCAGCATCG‑AAATGACGCCCGCTcgc                                                           <  1:806955/81‑1 (MQ=255)              aCGCGGGTGCTGGAGGCTTCCGGTCTGCGTGAAGGGTATGAATATGAAACCCAGGTCAGCATCG‑AAATGACGCCCGCTcg                                                            <  1:806954/80‑1 (MQ=255)               cgcgGGTGCTGGAGGCTTCCGGTCTGCGTGAAGGGTATGAATATGAAACCCAGGTCAGCATCG‑AAATGACGCCCGCTcgcg                                                          <  1:806956/81‑1 (MQ=255)                   ggTGCTGGAGGCTTCCGGTCTGCGTGAAGGGTATGAATATGAAACCCAGGTCAGCATCG‑AAATGACGCCCGCTCGCGGATg                                                      >  2:807010/1‑81 (MQ=255)                    gTGCTGGAGGCTTCCGGTCTGCGTGAAGGGTATGAATATGAAACCCAGGTCAGCATCG‑AAATGACGCCCGCTCGCGGATGc                                                     <  2:806957/81‑1 (MQ=255)                       cTGGAGGCTTCCGGTCTGCGTGAAGGGTATGAATATGAAACCCAGGTCAGCATCG‑AAATGACGCCCGCTCGCGGATGCAg                                                   <  2:806958/80‑1 (MQ=255)                       cTGGAGGCTTCCGGTCTGCGTGAAGGGTATGAATATGAAACCCAGGTCAGCATCG‑AAATGACGCCCGCTCGCGGATGCAGc                                                  >  2:807007/1‑81 (MQ=255)                         ggAGGCTTCCGGTCTGCGTGAAGGGTATGAATATGAAACCCAGGTCAGCATCG‑AAATGACGCCCGCTCGCGGATGCAGCCg                                                >  2:807008/1‑81 (MQ=255)                         ggAGGCTTCCGGTCTGCGTGAAGGGTATGAATATGAAACCCAGGTCAGCATCG‑AAATGACGCCCGCTCGCGGATGCAGCCg                                                >  2:807015/1‑81 (MQ=255)                                  cGGTCTGCGTGAAGGGTATGAATATGAAACCCAGGTCAGCATCG‑AAATGACGCCCGCTCGCGGATGCAGCCGGATGTCAt                                        <  1:806959/80‑1 (MQ=255)                                  cGGTCTGCGTGAAGGGTATGAATATGAAACCCAGGTCAGCATCG‑AAATGACGCCCGCTCGCGGATGCAGCCGGATGTCATc                                       <  2:806960/81‑1 (MQ=255)                                   ggTCTGCGTGAAGGGTATGAATATGAAACCCAGGTCAGCATCG‑AAATGACGCCCGCTCGCGGATGCAGCCGGATGTCATCg                                      >  1:807027/1‑81 (MQ=255)                                         cGTGAAGGGTATGAATATGAAACCCAGGTCAGCATCG‑AAATGACGCCCGCTCGCGGATGCAGCCGGATGTc                                          <  2:806961/71‑1 (MQ=255)                                         cGTGAAGGGTATGAATATGAAACCCAGGTCAGCATCG‑AAATGACGCCCGCTCGCGGATGCAGCCGGATGTCATCGTgc                                   <  2:806963/78‑1 (MQ=255)                                         cGTGAAGGGTATGAATATGAAACCCAGGTCAGCATCG‑AAATGACGCCCGCTCGCGGATGCAGCCGGATGTCATCGTGCGcc                                >  1:807013/1‑81 (MQ=255)                                         cGTGAAGGGTATGAATATGAAACCCAGGACAGCATCG‑AAATGACGCCCGCTCGCGGATGCAGCCGGATGTCATCGTgc                                   <  2:806962/78‑1 (MQ=255)                                          gTGAAGGGTATGAATATGAAACCCAGGTCAGCATCG‑AAATGACGCCCGCTCGCGGATGCAGCCGGATGTCATCGTGCGCCt                               >  1:807037/1‑81 (MQ=255)                                                ggTATGAATATGAAACCCAGGTCAGCATCG‑AAATGACGCCCGCTCGCGGATGCAGCCGGATGTCATCGTGCGCCTGCCg                           >  1:807021/1‑79 (MQ=255)                                                    tgaatatgaaACCCAGGTCAGCATCG‑AAATGACGCCCGCTCGCGGATGCAGCCGGATGTCATCGTGCGCCTGCCGCAggg                      <  1:806964/80‑1 (MQ=255)                                                      aatatgaaACCCAGGTCAGCATCG‑AAATGACGCCCGCTCGCGGATGCAGCCGGATGTCATCGTGCGCCTGCCGCAGGGaa                    >  1:807020/1‑80 (MQ=255)                                                      aatatgaaACCCAGGTCAGCATCG‑AAATGACGCCCGCTCGCGGATGCAGCCGGATGTCATCGTGCGCCTGCCGCAGGGaa                    <  1:806965/80‑1 (MQ=255)                                                       atatgaaACCCAGGTCAGCATCG‑AAATGACGCCCGCTCGCGGATGCAGCCGGATGTCATCGTGCGCCTGCCGCAGGGaaaa                  <  1:806966/81‑1 (MQ=255)                                                       atatgaaACCCAGGTCAGCATCG‑AAATGACGCCCGCTCGCGGATGCAGCCGGATGTCATCGTGCGCCTGCCGCAGGGaa                    >  1:807016/1‑79 (MQ=255)                                                        tatgaaACCCAGGTCAGCATCG‑AAATGACGCCCGCTCGCGGATGCAGCCGGATGTCATCGTGCGCCTGCCGCAGGGaa                    >  2:807030/1‑78 (MQ=255)                                                        tatgaaACCCAGGTCAGCATCG‑AAATGACGCCCGCTCGCGGATGCAGCCGGATGTCATCGTGCGCCTGCCGCAGGGaa                    >  2:807032/1‑78 (MQ=255)                                                          tgaaACCCAGGTCAGCATCG‑AAATGACGCCCGCTCGCGGATGCAGCCGGATGTCATCGTGCGCCTGCCGCAGGGAAAAGAt               <  1:806967/81‑1 (MQ=255)                                                            aaaCCCAGGTCAGCATCG‑AAATGACGCCCGCTCGCGGATGCAGCCGGATGTCATCGTGCGCCTGCCGCAGGGAAAAGAtgt             <  2:806968/81‑1 (MQ=255)                                                             aaCCCAGGTCAGCATCG‑AAATGACGCCCGCTCGCGGATGCAGCCGGATGTCATCGTGCGCCTGCCGCAGGGAAAAGAtgtg            >  2:807039/1‑81 (MQ=255)                                                                   ggTCAGCATCG‑AAATGACGCCCGCTCGCGGATGCAGCCGGATGTCATCGTGCGCCTGCCGCAGGGAAAAGATGTGGTGATc      <  2:806969/81‑1 (MQ=255)                                                                    gTCAGCATCG‑AAATGACGCCCGCTCGCGGATGCAGCCGGATGTCATCGTGCGCCTGCCGCAGGGAAAAGATGTGGTGATCg     <  2:806970/81‑1 (MQ=255)                                                                      cAGCATCG‑AAATGACGCCCGCTCGCGGATGCAGCCGGATGTCATCGTGCGCCTGCCGCAGGGAAAAGATGTGGTGATCGAc   >  1:807041/1‑81 (MQ=255)                                                                      cAGCATCG‑AAATGACGCCCGCTCGCGGATGCAGCCGGATGTCATCGTGCGCCTGCCGCAGGGAAAAGATGTGGTGATCGAc   >  2:807018/1‑81 (MQ=255)                                                                       aGCATCG‑AAATGACGCCCGCTCGCGGATGCAGCCGGATGTCATCGTGCGCCTGCCGCAGGGAAAAGATGTGGTGATCGACg  >  2:807019/1‑81 (MQ=255)                                                                              |                                                                           CGAGGTAGTATTGACGCGGGTGCTGGAGGCTTCCGGTCTGCGTGAAGGGTATGAATATGAAACCCAGGTCAGCATCGAAAATGACGCCCGCTCGCGGATGCAGCCGGATGTCATCGTGCGCCTGCCGCAGGGAAAAGATGTGGTGATCGACG  >  NC_000913/4017950‑4018101

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