Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,328,006	C→T	P291L (CCG→CTG)	crfC →	clamp‑binding sister replication fork colocalization protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,328,006	0	C	T	100.0%	96.4 / NA	29	P291L (CCG→CTG)	crfC	clamp‑binding sister replication fork colocalization protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (11/18);  total (11/18)

TGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGCGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGAAGAGGTCCGTGAGGCGATTTTGGCGGTGGGGCAATCGGTGCCGCTGTACGTGCTGGTCAATAAGTTCGATCAACAGGATCGTAACAGTGACGACGCCGACCAGGTGCGGGCACTGATTTCCGGGACGCTGATGAAAGGCTGTATTACGCCACAGCAGATATTTCCGGTGTCGTCGAT  >  NC_000913/4327877‑4328141                                                                                                                                  |                                                                                                                                        tGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGCGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGAAGAGGTCCGTGAGGCGATTTTGGCGGTGGGGCAATCGGTGCTGCTGTACGTg                                                                                                                               <  1:425164/140‑1 (MQ=255)          tGCTTAACCAGCAGCTGGCACGCGCCTCGGCGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGAAGAGGTCCGTGAGGCGATTTTGGCGGTGGGGCAATCGGTGCTGCTGTACGTGCTGGTCAAt                                                                                                                      <  1:588486/140‑1 (MQ=255)                 ccAGCAGCTGGCACGCGCCTCGGCGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGAAGAGGTCCGTGAGGCGATTTTGGCGGTGGGGCAATCGGTGCTGCTGTACGTGCTGGTCAATAAGTTCg                                                                                                               <  2:783721/140‑1 (MQ=255)                       gcTGGCACGCGCCTCGGCGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGAAGAGGTCCGTGAGGCGATTTTGGCGGTGGGGCAATCGGTGCTGCTGTACGTGCTGGTCAATAAGTTCGATCAAc                                                                                                         <  1:67038/140‑1 (MQ=255)                               gcgcCTCGGCGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGAAGAGGTCCGTGAGGCGATTTTGGCGGTGGGGCAATCGGTGCTGCTGTACGTGCTGGTCAATAAGTTCGATCAACAGGATCGt                                                                                                 <  2:235610/140‑1 (MQ=255)                                            aCTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGAAGAGGTCCGTGAGGCGATTTTGGCGGTGGTGCAATCGGTGCTGCTGTACGTGCTGGTCAATAAGTTCGATCAACAGGATCGTAACAGTGacgacg                                                                                    <  2:761540/140‑1 (MQ=255)                                                gCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGAAGAGGTCCGTGAGGCGATTTTGGCGGTGGGGCAATCGGTGCTGCTGTACGTGCTGGTCAATAAGTTCGATCAACAGGATCGTAACAGTGACGACGCCGa                                                                                <  1:333028/140‑1 (MQ=255)                                                       tGGGTTATACGCAACTGAAAACGATCTCCGATGAAGAGGGCCGTGAGGCGGTTTTGGCGGGGGGGGAATCGGGGCTGCTGTACGTGGTGGGCAATAAATTTGATCAACAGGATCGTAACAGTGACGACGCCGACCCgggg                                                                         >  2:147221/1‑138 (MQ=255)                                                           ttATACGCAACTGAAATCGATCTCCGATGAAGAGGTCCGTGAGGCGATTTTGGCGGTGGGGCAATCGGTGCTGCTGTACGTGCTGGTCAATAAGTTCGATCAACAGGATCGTAACAGTGACGACGCCGACCAGGTGCggg                                                                     <  1:459951/140‑1 (MQ=255)                                                             ataCGCAACTGAAATCGATCTCCGATGAAGAGGTCCGTGAGGCGATTTTGGCGGTGGGGCAATCGGTGCTGCTGTACGTGCTGGTCAATAAGTTCGATCAACAGGATCGTAACAGTGACGACGCCGACCAGGTGCGGGCa                                                                   <  1:147221/140‑1 (MQ=255)                                                              taCGCAACTGAAATCGATCTCCGATGAAGAGGTCCGTGAGGCGATTTTGGCGGTGGGGCAATCGGTGCTGCTGTACGTGCTGGTCAATAAGTTCGATCAACAGGATCGTAACAGTGACGACGCCGACCAGGTGCGGGCAc                                                                  >  2:515122/1‑140 (MQ=255)                                                                   aaCTGAAATCGATCTCCGATGAAGAGGTCCGTGAGGCGATTTTGGCGGTGGGGCAATCGGTGCTGCTGTACGTGCTGGTCAATAAGTTCGATCAACAGGATCGTAACAGTGACGACGCCGACCAGGTGCGGGCACTGAtt                                                             >  2:660962/1‑140 (MQ=255)                                                                   aaCTGAAATCGATCTCCGATGAAGAGGTCCGTGAGGCGATTTTGGCGGTGGGGCAATCGGTGCTGCTGTACGTGCTGGTCAATAAGTTCGATCAACAGGATCGTAACAGTGACGACGCCGACCAGGTGCGGGCACTGAtt                                                             >  2:656139/1‑140 (MQ=255)                                                                   aaCTGAAATCGATCTCCGATGAAGAGGTCCGTGAGGCGATTTTGGCGGTGGGGCAATCGGTGCTGCTGTACGTGCTGGTCAATAAGTTCGATCAACAGGATCGTAACAGTGACGACGCCGACCAGGTGCGGGCACTGAtt                                                             >  1:35716/1‑140 (MQ=255)                                                                   aaCTGAAATCGATCTCCGATGAAGAGGTCCGTGAGGCGATTTTGGCGGTGGGGCAATCGGTGCTGCTGTACGTGCTGGTCAATAAGTTCGATCAACAGGATCGTAACAGTGACGACGCCGACCAGGTGCGGGCACTGAtt                                                             >  1:860273/1‑140 (MQ=255)                                                                        aaaTCGATCTCCGATGAAGAGGTCCGTGAGGCGATTTTGGCGGTGGGGCAATCGGTGCTGCTGTACGTGCTGGTCAATAAGTTCGATCAACAGGATCGTAACAGTGACGACGCCGACCAGGTGCGGGCACTGATTTCCgg                                                        >  2:500799/1‑140 (MQ=255)                                                                           tCGATCTCCGATGAAGAGGTCCGTGAGGCGATTTTGGCGGTGGGGCAATCGGTGCTGCTGTACGTGCTGGGCAATAAGTTCGATCAACAGGATCGTAACAGTGACGACGCCGACCAGGTGCGGGCACTGATTTCCGGGAc                                                     <  1:697132/140‑1 (MQ=255)                                                                               tctcCGATGAAGAGGTCCGTGAGGCGATTTTGGCGGTGGGGCAATCGGTGCTGCTGTACGTGCTGGTCAATAAGTTCGATCAACAGGATCGTAACAGTGACGACGCCGACCAGGTGCGGGCACTGATTTCCGGGACGCtg                                                 >  2:235174/1‑140 (MQ=255)                                                                                     aTGAAGAGGTCCGTGAGGCGATTTTGGCGGTGGGGCAATCGGTGCTGCTGTACGTGCTGGTCAATAAGTTCGATCAACAGGATCGTAACAGTGACGACGCCGACCAGGTGCGGGCACTGATTTCCGGGACGCTGATGaaa                                           <  1:656139/140‑1 (MQ=255)                                                                                       gAAGAGGTCCGTGAGGCGATTTTGGCGGTGGGGCAATCGGTGCTGCTGTACGTGCTGGTCAATAAGTTCGATCAACAGGATCGTAACAGTGACGACGCCGACCAGGTGCGGGCACTGATTTCCGGGACGCTGATGAAAgg                                         <  1:761414/140‑1 (MQ=255)                                                                                             gTCCGTGAGGCGATTTTGGCGGTGGGGCAATCGGTGCTGCTGTACGTGCTGGTCAATAAGTTCGATCAACAGGATCGTAACAGTGACGACGCCGACCAGGTGCGGGCACTGATTTCCGGGACGCTGATGAAAGGCTGTAt                                   >  2:266039/1‑140 (MQ=255)                                                                                                    aGGCGATTTTGGCGGTGGGGCAATCGGTGCTGCTGTACGTGCTGGTCAATAAGTTCGATCAACAGGATCGTAACAGTGACGACGCCGACCAGGTGCGGGCACTGATTTCCGGGACGCTGATGAAAGGCTGTATTACGCca                            <  1:500799/140‑1 (MQ=255)                                                                                                       cGATTTTGGCGGTGGGGCAATCGGTGCTGCTGTACGTGCTGGTCAATAAGTTCGATCAACAGGATCGTAACAGTGACGACGCCGACCAGGTGCGGGCACTGATTTCCGGGACGCTGATGAAAGGCTGTATTACGCCAcag                         >  1:22391/1‑140 (MQ=255)                                                                                                       cGATTTTGGCGGTGGGGCAATCGGTGCTGCTGTACGTGCTGGTCAATAAGTTCGATCAACAGGATCGTAACAGTGACGACGCCGACCAGGTGCGGGCACTGATTTCCGGGACGCTGATGAAAGGCTGTATTACGCCAcag                         >  1:131382/1‑140 (MQ=255)                                                                                                         aTTTTGGCGGTGGGGCAATCGGTGCTGCTGTACGTGCTGGTCAATAAGTTCGATCAACAGGATCGTAACAGTGACGACGCCGACCAGGTGCGGGCACTGATTTCCGGGACGCTGATGAAAGGCTGTATTACGCCAcagca                       <  1:515122/140‑1 (MQ=255)                                                                                                         aTTTTGGCGGTGGGGCAATCGGTGCTGCTGTACGTGCTGGTCAATAAGTTCGATCAACAGGATCGTAACAGTGACGACGCCGACCAGGTGCGGGCACTGATTTCCGGGACGCTGATGAAAGGCTGTATTACGCCAcagca                       <  2:713820/140‑1 (MQ=255)                                                                                                                   tGGGGCAATCGGTGCTGCTGTACGTGCTGGTCAATAAGTTCGATCAACAGGATCGTAACAGTGACGACGCCGACCAGGTGCGGGCACTGATTTCCGGGACGCTGATGAAAGGCTGTATTACGCCACAGCAGATATTTCCg             >  1:264962/1‑140 (MQ=255)                                                                                                                          aTCGGTGCTGCTGTACGTGCTGGTCAATAAGTTCGATCAACAGGATCGTAACAGTGACGACGCCGACCAGGTGCGGGCACTGATTTCCGGGACGCTGATGAAAGGCTGTATTACGCCACAGCAGATATTTCCGGTgtcgt      <  2:22391/140‑1 (MQ=255)                                                                                                                             ggTGCTGCTGTACGTGCTGGTCAATAAGTTCGATCAACAGGATCGTAACAGTGACGACGCCGACCAGGTGCGGGCACTGATTTCCGGGACGCTGATGAAAGGCTGTATTACGCCACAGCAGATATTTCCGGTGTCGTCGa   <  1:235174/140‑1 (MQ=255)                                                                                                                              gTGCTGCTGTACGTGCTGGTCAATAAGTTCGATCAACAGGATCGTAACAGTGACGACGCCGACCAGGTGCGGGCACTGATTTCCGGGACGCTGATGAAAGGCTGTATTACGCCACAGCAGATATTTCCGGTGTCGTCGAt  <  1:660962/140‑1 (MQ=255)                                                                                                                                  |                                                                                                                                        TGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGCGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGAAGAGGTCCGTGAGGCGATTTTGGCGGTGGGGCAATCGGTGCCGCTGTACGTGCTGGTCAATAAGTTCGATCAACAGGATCGTAACAGTGACGACGCCGACCAGGTGCGGGCACTGATTTCCGGGACGCTGATGAAAGGCTGTATTACGCCACAGCAGATATTTCCGGTGTCGTCGAT  >  NC_000913/4327877‑4328141

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