Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	643,960	C→A	R3I (AGA→ATA)	rnk ←	nucleoside diphosphate kinase regulator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	643,960	0	C	A	100.0%	80.7 / NA	25	R3I (AGA→ATA)	rnk	nucleoside diphosphate kinase regulator
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (11/14);  total (11/14)

TCTTCTGGCGAACACATTTGGGCGCGATCCAACTCTGCGTTTAACGCGTCGGCGATTGGCAAACCAGCATAGGCGGGTTGCTCCAGCAGAATATCGATGCGTTCGGCATCCAGGTCGTTAATGATGATAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTT  >  NC_000913/643825‑644096                                                                                                                                        |                                                                                                                                         tcttctGGCGAACACATTTGGGCGCGATCCAACTCTGCGTTTAACGCGTCGGCGATTGGCAAACCAGCATAGGCGGGTTGCTCCAGCAGAATATCGATGCGTTCGGCATCCAGGTCGTTAATGATGATAGTTGGTATGGa                                                                                                                                      <  1:206548/140‑1 (MQ=255) tcttctGGCGAACACATTTGGGCGCGATCCAACTCTGCGTTTAACGCGTCGGCGATTGGCAAACCAGCATAGGCGGGTTGCTCCAGCAGAATATCGATGCGTTCGGCATCCAGGTCGTTAATGATGATAGTTGGTATGGa                                                                                                                                      <  1:278139/140‑1 (MQ=255) tcttctGGCGAACACATTTGGGCGCGATCCAACTCTGCGTTTAACGCGTCGGCGATTGGCAAACCAGCATAGGCGGGTTGCTCCAGCAGAATATCGATGCGTTCGGCATCCAGGTCGTTAATGATGATAGTTGGTATGGa                                                                                                                                      <  2:316335/140‑1 (MQ=255)   ttctGGCGAACACATTTGGGCGCGATCCAACTCTGCGTTTAACGCGTCGGCGATTGGCAAACCAGCATAGGCGGGTTGCTCCAGCAGAATATCGATGCGTTCGGCATCCAGGTCGTTAATGATGATAGTTGGTATGGACa                                                                                                                                    >  2:351720/1‑140 (MQ=255)          gAACACATTTGGGCGCGATCCAACTCTGCGTTTAACGCGTCGGCGATTGGCAAACCAGCATAGGCGGGTTGCTCCAGCAGAATATCGATGCGTTCGGCATCCAGGTCGTTAATGATGATAGTTGGTATGGACATTTTTAc                                                                                                                             <  1:268065/140‑1 (MQ=255)                 tttGGGCGCGATCCAACTCTGCGTTTAACGCGTCGGCGATTGGCAAACCAGCATAGGCGGGTTGCTCCAGCAGAATATCGATGCGTTCGGCATCCAGGTCGTTAATGATGATAGTTGGTATGGACATTTTTACTCCATgt                                                                                                                      <  2:128921/140‑1 (MQ=255)                       cgcgATCCAACTCTGCGTTTAACGCGTCGGCGATTGGCAAACCAGCATAGGCGGGTTGCTCCAGCAGAATATCGATGCGTTCGGCATCCAGGTCGTTAATGATGATAGTTGGTATGGACATTTTTACTCCATGTCGTCgg                                                                                                                <  1:368153/140‑1 (MQ=255)                        gcgATCCAACTCTGCGTTTAACGCGTCGGCGATTGGCAAACCAGCATAGGCGGGTTGCTCCAGCAGAATATCGATGCGTTCGGCATCCAGGTCGTTAATGATGATAGTTGGTATGGACATTTTTACTCCATGTCGTCGGt                                                                                                               >  1:168691/1‑140 (MQ=255)                           aTCCAACTCTGCGTTTAACGCGTCGGCGATTGGCAAACCAGCATAGGCGGGTTGCTCCAGCAGAATATCGATGCGTTCGGCATCCAGGTCGTTAATGATGATAGTTGGTATGGACATTTTTACTCCATGTCGTCGGtgct                                                                                                            <  1:209916/140‑1 (MQ=255)                                   ctGCGTTTAACGCGTCGGCGATTGGCAAACCAGCATAGGCGGGTTGCTCCAGCAGAATATCGATGCGTTCGGCATCCAGGTCGTTAATGATGATAGTTGGTATGGACATCTTTACTCCATGTCGTCGGTGCTGCGAgtgt                                                                                                    <  1:11929/140‑1 (MQ=255)                                     gCGTTTAACGCGTCGGCGATTGGCAAACCAGCATAGGCGGGTTGCTCCAGCAGAATATCGATGCGTTCGGCATCCAGGTCGTTAATGATGATAGTTGGTATGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCg                                                                                                  <  2:83348/140‑1 (MQ=255)                                      cGTTTAACGCGTCGGCGATTGGCAAACCAGCATAGGCGGGTTGCTCCAGCAGAATATCGATGCGTTCGGCATCCAGGTCGTTAATGATGATAGTTGGTATGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGc                                                                                                 >  1:140209/1‑140 (MQ=255)                                      cGTTTAACGCGTCGGCGATTGGCAAACCAGCATAGGCGGGTTGCTCCAGCAGAATATCGATGCGTTCGGCATCCAGGTCGTTAATGATGATAGTTGGTATGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGc                                                                                                 >  2:62424/1‑140 (MQ=255)                                      cGTTTAACGCGTCGGCGATTGGCAAACCAGCATAGGCGGGTTGCTCCAGCAGAATATCGATGCGTTCGGCATCCAGGTCGTTAATGATGATAGTTGGTATGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGc                                                                                                 >  2:220847/1‑140 (MQ=255)                                      cGTTTAACGCGTCGGCGATTGGCAAACCAGCATAGGCGGGTTGCTCCAGCAGAATATCGATGCGTTCGGCATCCAGGTCGTTAATGATGATAGTTGGTATGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGc                                                                                                 >  2:153250/1‑140 (MQ=255)                                                    gCGATTGGCAAACCAGCATAGGCGGGTTGCTCCAGCAGAATATCGATGCGTTCGGCATCCAGGTCGTTAATGATGATAGTTGGTATGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATAccc                                                                                   <  1:220847/140‑1 (MQ=255)                                                        ttGGCAAACCAGCATAGGCGGGTTGCTCCAGCAGAATATCGATGCGTTCGGCATCCAGGTCGTTAATGATGATAGTTGGTATGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCaaaa                                                                               >  1:409788/1‑140 (MQ=255)                                                                   gCATAGGCGGGTTGCTCCAGCAGAATATCGATGCGTTCGGCATCCAGGTCGTTAATGATGATAGTTGGTATGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCa                                                                    <  2:41484/140‑1 (MQ=255)                                                                       aGGCGGGTTGCTCCAGCAGAATATCGATGCGTTCGGCATCCAGGTCGTTAATGATGATAGTTGGTATGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGt                                                                <  1:372921/140‑1 (MQ=255)                                                                         gCGGGTTGCTCCAGCAGAATATCGATGCGTTCGGCATCCAGGTCGTTAATGATGATAGTTGGTATGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCa                                                              >  2:304584/1‑140 (MQ=255)                                                                                gCTCCAGCAGAATATCGATGCGTTCGGCATCCAGGTCGTTAATGATGATAGTTGGTATGGACATTTTTACTCCTTGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGgcggcg                                                       <  1:351720/140‑1 (MQ=255)                                                                                 cTCCAGCAGAATATCGATGCGTTCGGCATCCAGGTCGTTAATGATGATAGTTGGTATGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGa                                                      >  1:168444/1‑140 (MQ=255)                                                                                                                     cataccgatgatAGTTGGTATGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCa                  >  2:152884/7‑140 (MQ=255)                                                                                                                               atAGTTGGTATGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAg        <  1:62424/140‑1 (MQ=255)                                                                                                                                tAGTTGGTATGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGt       >  1:149802/1‑140 (MQ=255)                                                                                                                                     ggTATGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTACCCGTCAGGCGGCGAGGGCTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGAttt  >  2:331128/1‑140 (MQ=255)                                                                                                                                        |                                                                                                                                         TCTTCTGGCGAACACATTTGGGCGCGATCCAACTCTGCGTTTAACGCGTCGGCGATTGGCAAACCAGCATAGGCGGGTTGCTCCAGCAGAATATCGATGCGTTCGGCATCCAGGTCGTTAATGATGATAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTT  >  NC_000913/643825‑644096

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