Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,018,153	A→G	Y274C (TAC→TGC)	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,153	0	A	G	100.0%	129.3 / NA	38	Y274C (TAC→TGC)	rmuC	putative recombination limiting protein RmuC
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (18/20);  total (18/20)

GGTCAGCATCGAAAATGACGCCCGCTCGCGGATGCAGCCGGATGTCATCGTGCGCCTGCCGCAGGGAAAAGATGTGGTGATCGACGCCAAAATGACGCTGGTCGCCTATGAACGCTATTTTAACGCCGAAGACGACTACACCCGCGAAAGCGCGCTACAGGAACATATCGCGTCGGTGCGTAACCATATCCGTTTGCTGGGACGCAAAGATTATCAACAGCTGCCGGGGCTGCGAACTCTGGATTACGTGCTGATGTTTATTCCCGTTGAAC  >  NC_000913/4018016‑4018287                                                                                                                                          |                                                                                                                                       ggTCAGCATCGAAAATGACGCCCGCTCGCGGATGCAGCCGGATGTCATCGTGCGCCTGCCGTAGGGAAAAGATGTGGTGATCGACGCCAAAATGACGCTGGTCGCCTATGAACGCTATTTTAACGCCGAAGACGACTgca                                                                                                                                      >  1:532481/1‑140 (MQ=255) ggTCAGCATCGAAAATGACGCCCGCTCGCGGATGCAGCCGGATGTCATCGTGCGCCTGCCGCAGGGAAAAGATGTGGTGATCGACGCCAAAATGACGCTGGTCGCCTATGAACGCTATTTTAACGCCGAAGACGACTgca                                                                                                                                      <  2:20186/140‑1 (MQ=255)    cAGCATCGAAAATGACGCCCGCTCGCGGATGCAGCCGGATGTCATCGTGCGCCTGCCGCAGGGAAAAGATGTGGTGATCGACGCCAAAATGACGCTGGTCGCCTATGAACGCTATTTTAACGCCGAAGACGACTGCAccc                                                                                                                                   >  1:126771/1‑140 (MQ=255)        aTCGAAAATGACGCCCGCTCGCGGATGCAGCCGGATGTCATCGTGCGCCTGCCGCAGGGAAAAGATGTGGTGATCGACGCCAAAATGACGCTGGTCGCCTATGAACGCTATTTTAACACCGAAGCCGACTGCACCCGCGa                                                                                                                               >  1:260307/1‑140 (MQ=255)            aaaaTGACGCCCGCTCGCGGATGCAGCCGGATGTCATCGTGCGCCTGCCGCAGGGAAAAGATGTGGTGATCGACGCCAAAATGACGCTGGTCGCCTATGAACGCTATTTTAACGCCGAAGACGACTGCACCCGCGAAAgc                                                                                                                           >  2:405383/1‑140 (MQ=255)             aaaTGACGCCCGCTCTCGGATGTAGCCGGCTGTCATCGTGCGCCTGCCGCAGGGAAAAGATGTGGTGATCGACGCCAAAATGACGCTGGTCGCCTATGAACGCTATTTTAACGCCGAAGACGACTGCACCCGCGAAAgcg                                                                                                                          <  1:226990/140‑1 (MQ=255)             aaaTGACGCCCGCTCGCGGATGCAGCCGGATGTCATCGTGCGCCTGCCGCAGGGAAAAGATGTGGTGATCGACGCCAAAATGACGCTGGTCGCCTATGAACGCTATTTTAACGCCGAAGACGACTGCACCCGCGAAAgcg                                                                                                                          <  2:109984/140‑1 (MQ=255)                 gACGCCCGCTCGCGGATGCAGCCGGATGTCATCGTGCGCCTGCCGCAGGGAAAAGATGTGGTGATCGACGCCAAAATGACGCTGGTCGCCTATGAACGCTATTTTAACGCCGAAGACGAATGCACCCGCGAAAGCGCGCt                                                                                                                      >  2:24007/1‑140 (MQ=255)                  aCGCCCGCTCGCGGATGCAGCCGGATGTCATCGTGCGCCTGCCGCAGGGAAAAGATGTGGTGATCGACGCCAAAATGACGCTGGTCGCCTATGAACGCTATTTTAACGCCGAAGACGACTGCACCCGCGAAAGCGCGCTa                                                                                                                     <  1:168542/140‑1 (MQ=255)                    gCCCGCTCGCGGATGCAGCCGGATGTCATCGTGCGCCTGCCGCAGGGAAAAGATGTGGTGATCGACGCCAAAATGACGCTGGTCGCCTATGAACGCTATTTTAACGCCGAAGACGACTGCACCCGCGAAAGCGCGCTACa                                                                                                                   <  2:382843/140‑1 (MQ=255)                                  gCAGCCGGATGTCATCGTGCGCCTGCCGCAGGGAAAAGATGTGGTGATCGACGCCAAAATGACGCTGGTCGCCTATGAACGCTATTTTAACGCCGAAGACGACTGCACCCGCGAAAGCGCGCTACAGGAACATATCGCGt                                                                                                     <  1:148677/140‑1 (MQ=255)                                        ggATGTCATCGTGCGCCTGCCGCAGGGAAAAGATGTGGTGATCGACGCCAAAATGACGCTGGTCGCCTATGAACGCTATTTTAACGCCGAAGACGACTGCACCCGCGAAAGCGCGCTACAGGAACATATCGCGTCGGTGc                                                                                               <  1:577208/140‑1 (MQ=255)                                          aTGTCATCGTGCGCCTGCCGCAGGGAAAAGATGTGGTGATCGACGCCAAAATGACGCTGGTCGCCTATGAACGCTATTTTAACGCCGAAGACGACTGCACCCGCGAAAGCGCGCTACAGGAACATATCGCGTCGGTGCGt                                                                                             <  2:664140/140‑1 (MQ=255)                                                tCGTGCGCCTGCCGCAGGGAAAAGATGTGGTGATCGACGCCAAAATGACGCTGGTCGCCTATGAACGCTATTTTAACGCCGAAGACGACTGCACCCGCGAAAGCGCGCTACAGGAACATATCGCGTCGGTGCGTAACCat                                                                                       >  1:714901/1‑140 (MQ=255)                                                 cGTGCGCCTGCCGCAGGGAAAAGATGTGGTGATCGACGCCAAAATGACGCTGGTCGCCTATGAACGCTATTTTAACGCCGAAGACGACTGCACCCGCGAAAGCGCGCTACAGGAACATATCGCGTCGGTGCGTAACCata                                                                                      <  2:128463/140‑1 (MQ=255)                                                      gcCTGCCGCAGGGAAAAGATGTGGTGATCGACGCCAAAATGACGCTGGTCGCCTATGAACGCTATTTTAACGCCGAAGACGACTGCACCCGCGAAAGCGCGCTACAGGAACATATCGCGTCGGTGCGTAACCATATCCGt                                                                                 >  2:685783/1‑140 (MQ=255)                                                                    aaaGATGTGGTGATCCACGCCAAAACGACCCTGGTCGTCTATGAACGCTATTTTAACGCCGAAGACGACTGCACCCGCGAAAGCGCGCTACAGGAACATATCGCGTCGGTGCGTAACCATATCCGTTTGCTGGGACGCaa                                                                   <  1:206399/140‑1 (MQ=255)                                                                     aaGATGTGGTGATCGACGCCAAAATGACGCTGGTCGCCTATGAACGCTATTTTAACGCCGAAGACGACTGCACCCGCGAAAGCGCGCTACAGGAACATATCGCGTCGGTGCGTAACCATATCCGTTTGCTGGGACGCaaa                                                                  >  1:2498/1‑140 (MQ=255)                                                                           tggtgATCGACGCCAAAATGACGCTGGTCGCCTATGAACGCTATTTTAACGCCGAAGACGACTGCACCCGCGAAAGCGCGCTACAGGAACACATCGCGTCAGTGCGTAGCCATATCCATTTGGTGGGACGCAAAGATTAt                                                            >  2:206399/1‑140 (MQ=255)                                                                             gtgATCGACGCCAAAATGACGCTGGTCGCCTATGAACGCTATTTTAACGCCGAAGACGACTGCACCCGCGAAAGCGCGCTACAGGAACATATCGCGTCGGTGCGTAACCATATCCGTTTGCTGGGACGCAAAGATTATCa                                                          >  1:418422/1‑140 (MQ=255)                                                                                aTCGACGCCAAAATGACGCTGGTCGCCTATGAACGCTATTTTAACGCCGAAGACGACTGCACCCGCGAAAGCGCGCTACAGGAACATATCGCGTCGGTGCGTAACCATATCCGTTTGCTGGGACGCAAAGATTATCAACa                                                       >  2:72433/1‑140 (MQ=255)                                                                                     cGCCAAAATGACGCTGGTCGCCTATGAACGCTATTTTAACGCCGAAGACGACTGCACCCGCGAAAGCGCGCTACAGGAACATATCGCGTCGGTGCGTAACCATATCCGTTTGCTGGGACGCAAAGATTATCAACAGCTGc                                                  <  1:405383/140‑1 (MQ=255)                                                                                           aaTGACGCTGGTCGCCTATGAACGCTATTTTAACGCCGAAGACGACTGCACCCGCGAAAGCGCGCTACAGGAACATATCGCGTCGGTGCGTAACCATATCCGTTTGCTGGGACGCAAAGATTATCAACAGCTGCCGGGGc                                            >  2:602412/1‑140 (MQ=255)                                                                                              gACGCTGGTCGCCTATGAACGCTATTTTAACGCCGAAGACGACTGCACCCGCGAAAGCGCGCTACAGGAACATATCGCGTCGGTGCGTAACCATATCCGTTTGCTGGGACGCAAAGATTATCAACAGCTGCCGGGGCTGc                                         <  1:602412/140‑1 (MQ=255)                                                                                                 gCTGGTCGCCTATGAACGCTATTTTAACGCCGAAGACGACTGCACCCGCGAAAGCGCGCTACAGGAACATATCGCGTCGGTGCGTAACCATATCCGTTTGCTGGGACGCAAAGATTATCAACAGCTGCCGGGGCTGCGaa                                      <  1:685783/140‑1 (MQ=255)                                                                                                     gTCGCCTATGAACGCTACTTTATCGCCGAAGTCGACTGCGCCCGCGAAGGCGCGCTACAGGAACATATCGCGTCGGTGCGTAACCATATCCGTTTGCTGGGACGCAAAGATTATCAACAGCTGCCGGGGCTGCGAActct                                  <  2:260307/140‑1 (MQ=255)                                                                                                        gCCTATGAACGCTATTTTAACGCCGAAGACGACTGCACCCGCGAAAGCGCGCTACAGGAACATATCGCGTCGGTGCGTAACCATATCCGTTTGCTGGGACGCAAAGATTATCAACAGCTGCCGGGGCTGCGAACTCTGGa                               <  2:115540/140‑1 (MQ=255)                                                                                                         ccTATGAACGCTATTTTAACGCCGAAGACGACTGCACCCGCGAAAGCGCGCTACAGGAACATATCGCGTCGGTGCGTAACCATATCCGTTTGCTGGGACGCAAAGATTATCAACAGCTGCCGGGGCTGCGAACTCTGGAt                              >  2:726138/1‑140 (MQ=255)                                                                                                               aaCGCTATTTTAACGGCGAAGAGGACTGCTCCCGCGATAGCGCGCTACAGGAACATATCGCGTCGGTGCGTAACCATATCCGTTTGCTGGGACGCAAAGATTATCAACAGCTGCCGGGGCTGCGAACTCTGGATTACGTg                        <  1:236998/140‑1 (MQ=255)                                                                                                                 cGCTATTTTAACGCCGAAGACGACTGCACCCGCGAAAGCGCGCTACAGGAACATATCGCGTCGGTGCGTAACCATATCCGTTTGCTGGGACTCAAAGATTTTCAGCAGTTGCCGGGGGTGCGAACTCTGGATTACGTGCt                      >  2:308469/1‑140 (MQ=255)                                                                                                                 cGCTATTTTAACGCCGAAGACGACTGCACCCGCGAAAGCGCGCTACAGGAACATATCGCGTCGGTGCGTAACCATATCCGTTTGCTGGGACGCAAAGATTATCAACAGCTGCCGGGGCTGCGAACTCTGGATTACGTGCt                      <  1:726138/140‑1 (MQ=255)                                                                                                                  gCTATTTTAACGCCGAAGACGACTGCACCCGCGAAAGCGCGCTACAGGAACATATCGCGTCGGTGCGTAACCATATCCGTTTGCTGGGACGCAAAGATTATCAACAGCTGCCGGGGCTGCGAACTCTGGATTACGTGCTg                     <  1:660632/140‑1 (MQ=255)                                                                                                                  gCTATTTTAACGCCGAAGACGACTGCACCCGCGAAAGCGCGCTACAGGAACATATCGCGTCGGTGCGTAACCATATCCGTTTGCTGGGACGCAAAGATTATCAACAGCTGCCGGGGCTGCGAACTCTGGATTACGTGCTg                     <  2:532481/140‑1 (MQ=255)                                                                                                                    tATTTTATCGCCGAAGACGGCTGCACCCGCGAAAGCGCGCTACAGGAACATATCGCGTCGGTGCGTAACCATATCCGTTTGCTGGGACGCAAAGATTATCAACAGCTGCCGGGGCTGCGAACTCTGGATTACGTGCTGAt                   >  2:422891/1‑140 (MQ=255)                                                                                                                    tATTTTAACGCCGAAGACGACTGCACCCGCGAAAGCGCGCTACAGGAACATATCGCGTCGGTGCGTAACCATATCCGTTTGCTGGGACGCAAAGATTATCAACAGCTGCCGGGGCTGCGAACTCTGGATTACGTGCTGAt                   >  1:501372/1‑140 (MQ=255)                                                                                                                       tttAACGCCGAAGACGACTGCACCCGCGAAAGCGCGCTACAGGAACATATCGCGTCGGTGCGTAACCATATCCGTTTGCTGGGACGCAAAGATTATCAACAGCTGCCGGGGCTGCGAACTCTGGATTACGTGCTGATGtt                <  2:126771/140‑1 (MQ=255)                                                                                                                         tAACGCCGAAGACGACTGCACCCGCGAAAGCGCGCTACAGGAACATATCGCGTCGGTGCGTAACCATATCCGTTTGCTGGGACGCAAAGATGATCAACAGATGCCGGGGCTGCGAACGCTGGATTACGTGCTGATGTTTa              >  2:200349/1‑140 (MQ=255)                                                                                                                           aCGCCGAAGACGACTGCACCCGCGAAAGCGCGCTACAGGAACATATCGCGTCGGTGCGTAACCATATCCGTTTGCTGGGACGCAAAGATTATCAACAGCTGCCGGGGCTGCGAACTCTGGATTACGTGCTGATGTTTAtt            >  1:161102/1‑140 (MQ=255)                                                                                                                           aCGCCGAAGACGAATGCACCCGCGAAAGCGCGCTACAGGAACATATCGCGTCGGTGCGTAACCATATCCGTTTGCTGGGACGCAAAGATTATCAACAGCTGCCGGGGCTGCGAACTCTGGATTACGTGCTGATGTTTAtt            <  1:24007/140‑1 (MQ=255)                                                                                                                                     cgacTGCACCCGCGAAAGCGCGCTACAGGAACATATCGCGTCGGTGCGTAACCATATCCGTTTGCTGGGACGCAAAGATTATCAACAGCTGCCGGGGCTGCGCACTCTGGATTACGTGCTGATGTTTATTCCCGTTGAAc  >  1:258320/1‑140 (MQ=255)                                                                                                                                          |                                                                                                                                       GGTCAGCATCGAAAATGACGCCCGCTCGCGGATGCAGCCGGATGTCATCGTGCGCCTGCCGCAGGGAAAAGATGTGGTGATCGACGCCAAAATGACGCTGGTCGCCTATGAACGCTATTTTAACGCCGAAGACGACTACACCCGCGAAAGCGCGCTACAGGAACATATCGCGTCGGTGCGTAACCATATCCGTTTGCTGGGACGCAAAGATTATCAACAGCTGCCGGGGCTGCGAACTCTGGATTACGTGCTGATGTTTATTCCCGTTGAAC  >  NC_000913/4018016‑4018287

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