Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,328,811	A→G	Q559Q (CAA→CAG)	crfC →	clamp‑binding sister replication fork colocalization protein, dynamin‑related

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,328,811	0	A	G	100.0%	109.8 / NA	30	Q559Q (CAA→CAG)	crfC	clamp‑binding sister replication fork colocalization protein, dynamin‑related
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (13/17);  total (13/17)

ATTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTGAGCAACAAGTCAAAATGGAATTGAGCGAGTCAGGGTTTCGTCCTGGGTTTCATTTTCCTGTTTTTCACGGCGTAGTTCCCCACTTCAACACTCGCCAGCTGTTCAGTGAAGTCATTTCGCGCCAGGAAGCAACGG  >  NC_000913/4328680‑4328941                                                                                                                                    |                                                                                                                                   aTTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTGAGCAGCAAGTCa                                                                                                                             >  2:1020556/1‑139 (MQ=255)              gtgAGACACTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTGAGCAGCAAGTCAAAATGGAATTGAg                                                                                                                <  1:121185/139‑1 (MQ=255)              gtgAGACACTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTGAGCAGCAAGTCAAAATGGAATTGAg                                                                                                                <  2:371301/139‑1 (MQ=255)              gtgAGACACTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTGAGCAGCAAGTCAAAATGGAATTGAg                                                                                                                <  1:238145/139‑1 (MQ=255)              gtgAGACACTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTGAGCAGCAAGTCAAAATGGAATTGAg                                                                                                                <  1:593886/139‑1 (MQ=255)                 agaCACTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTGAGCAGCAAGt                                                                                                                               >  2:1008702/1‑121 (MQ=255)                 agaCACTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTGAGCAGCAAGt                                                                                                                               <  1:1008702/121‑1 (MQ=255)                       tGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTGAGCAGCAAGTCAAAATGGAATTGAGCGAGTCAgg                                                                                                       >  1:29406/1‑139 (MQ=255)                                 ggCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTGAGCAGCAAGTCAAAATGGAATTGAGCGAGTCAGGGTTTCGTCCt                                                                                             >  2:498949/1‑139 (MQ=255)                                   cgcAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTGAGCAGCAAGTCAAAATGGAATTGAGCGAGTCAGGGTTTCGTCCTgg                                                                                           >  1:513129/1‑139 (MQ=255)                                                cAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTGAGCAGCAAGTCAAAATGGAATTGAGCGAGTCAGGGTTTCGTCCTGGGTTTCATTTTCCt                                                                              <  2:859175/139‑1 (MQ=255)                                                      cgACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTGAGCAGCAAGTCAAAATGGAATTGAGCGAGTCAGGGTTTCGTCCTGGGTTTCATTTTCCTGttttt                                                                        >  1:652358/1‑139 (MQ=255)                                                      cgACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTGAGCAGCAAGTCAAAATGGAATTGAGCGAGTCAGGGTTTCGTCCTGGGTTTCATTTTCCTGttttt                                                                        >  2:157225/1‑139 (MQ=255)                                                        aCTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTGAGCAGCAAGTCAAAATGGAATTGAGCGAGTCAGGGTTTCGTCCTGGGTTTCATTTTCCTGTTTTTCa                                                                      <  2:1288245/139‑1 (MQ=255)                                                          tttACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTGAGCAGCAAGTCAAAATGGAATTGAGCGAGTCAGGGTTTCGTCCTGGGTTTCATTTTCCTGTTTTTCACg                                                                    <  1:498949/139‑1 (MQ=255)                                                             aCGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTGAGCAGCAAGTCAAAATGGAATTGAGCGAGTCAGGGTTTCGTCCTGGGTTTCATTTTCCTGTTTTTCACGGCg                                                                 <  2:623676/139‑1 (MQ=255)                                                              cGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTGAGCAGCAAGTCAAAATGGAATTGAGCGAGTCAGGGTTTCGTCCTGGGTTTCATTTTCCTGTTTTTCACGGCGt                                                                <  1:523197/139‑1 (MQ=255)                                                                    gTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTGAGCCGCAAGTCAAAATGGAATTGAGCGAGTCAGGGTTTCGTCCTGGGTTTCATTTTCCTGTTTTTCACGGCGTAGTTcc                                                          >  2:1222246/1‑139 (MQ=255)                                                                    gTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCAATTGAGCAGCAAGTCAAAATGGAATTGAGCGAGTCAGGGTTTCGTCCTGGGTTTCATTTTCCTGTTTTTCACGGCGTAGTTcc                                                          >  2:1266665/1‑139 (MQ=255)                                                                     ttttGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTGAGCAGCAAGTCAAAATGGAATTGAGCGAGTCAGGGTTTCGTCCTGGGTTTCATTTTCCTGTTTTTCACGGCGTAGTTccc                                                         <  1:157225/139‑1 (MQ=255)                                                                                     tCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTGAGCAGCAAGTCAAAATGGAATTGAGCGAGTCAGGGTTTCGTCCTGGGTTTCATTTTCCTGTTTTTCACGGCGTAGTTCCCCACTTCAACACTCGcc                                         <  1:259923/139‑1 (MQ=255)                                                                                       cACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTGAGCAGCAAGTCAAAATGGAATTGAGCGAGTCAGGGTTTCGTCCTGGGTTTCATTTTCCTGTTTTTCACGGCGTAGTTCCCCACTTCAACACTCGCCAg                                       <  2:145229/139‑1 (MQ=255)                                                                                                    gtgtTAACCGATGTTATTCGCCCCATTGAGCAGCAAGTCAAAATGGAATTGAGCGAGTCAGGGTTTCGTCCTGGGTTTCATTTTCCTGTTTTTCACGGCGTAGTTCCCCACTTCAACACTc                                            <  1:933844/121‑1 (MQ=255)                                                                                                    gtgtTAACCGATGTTATTCGCCCCATTGAGCAGCAAGTCAAAATGGAATTGAGCGAGTCAGGGTTTCGTCCTGGGTTTCATTTTCCTGTTTTTCACGGCGTAGTTCCCCACTTCAACACTc                                            >  2:933844/1‑121 (MQ=255)                                                                                                                    ttCGCCCCATTGAGCAGCAAGTCAAAATGGAATTGAGCGAGTCAGGGTTTCGTCCTGGGTTTCATTTTCCTGTTTTTCACGGCGTAGTTCCCCACTTCAACACTCGCCAGCTGTTCAGTGAAGTCATTTCGCGCCAGGa          >  1:204133/1‑139 (MQ=255)                                                                                                                       gCCCCATTGAGCAGCAAGTCAAAATGGAATTGAGCGAGTCAGGGTTTCGTCCTGGGTTTCATTTTCCTGTTTTTCACGGCGTAGTTCCCCACTTCAACACTCGCCAGCTGTTCAGTGAAGTCATTTCGCGCCAGGAAGc       <  2:513129/139‑1 (MQ=255)                                                                                                                       cgcccATTGAGCAGCAAGTCAAAATGGAATTGAGCGAGTCAGGGTTTCGTCCTGGGTTTCATTTTCCTGTTTTTCACGGCGTAGTTCCCCACTTCAACACTCGCCAGCTGTTCAGTGAAGCCATTTCGCGCCAGGAAGc       <  1:732228/137‑1 (MQ=255)                                                                                                                         cccATTGAGCAGCAAGTCAAAATGGAATTGAGCGAGTCAGGGTTTCGTCCTGGGTTTCATTTTCCTGTTTTTCACGGCGTAGTTCCCCACTTCAACACTCGCCAGCTGt                                   <  2:980219/109‑1 (MQ=255)                                                                                                                         cccATTGAGCAGCAAGTCAAAATGGAATTGAGCGAGTCAGGGTTTCGTCCTGGGTTTCATTTTCCTGTTTTTCACGGCGTAGTTCCCCACTTCAACACTCGCCAGCTGt                                   >  1:980219/1‑109 (MQ=255)                                                                                                                            aTTGAGCAGCAAGTCAAAATGGAATTGAGCGAGTCAGGGTTTCGTCCTGGGTTTCATTTTCCTGTTTTTCACGGCGTAGTTCCCCACTTCAACACTCGCCAGCTGTTCAGTGAAGTCATTTCGCGCCAGGAAGCAACgg  >  1:781877/1‑139 (MQ=255)                                                                                                                                    |                                                                                                                                   ATTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTGAGCAACAAGTCAAAATGGAATTGAGCGAGTCAGGGTTTCGTCCTGGGTTTCATTTTCCTGTTTTTCACGGCGTAGTTCCCCACTTCAACACTCGCCAGCTGTTCAGTGAAGTCATTTCGCGCCAGGAAGCAACGG  >  NC_000913/4328680‑4328941

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