Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,802,309	G→C	P88A (CCG→GCG)	thrS ←	threonine‑‑tRNA ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,802,309	0	G	C	100.0%	77.9 / NA	22	P88A (CCG→GCG)	thrS	threonine‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base C (14/8);  total (14/8)

CGTAGTTTTTCTCAGCAAGCTCATGCATCCGCTTCTCGAGTGCTTCGACATCTTCCTGGGTTAACGTGCGGTCAAGATCAACGTCGTAATAAAAACCGTTGTCAATAACCGGGCCGATTGCCATTTTGGTATGCGGCCAAAGTTGTTTAATCGCGTGCCCTAACAGGTGCGCACAGGAGTGACGAATGATCTCCAGACCTTCTTCGTCTTTGGCGGTAATGATCGACAGTTGTGCGTCGTTTTCAATCAGATCGC  >  NC_000913/1802174‑1802428                                                                                                                                        |                                                                                                                        cgtAGTTTTTCTCAGCAAGCTCATGCATCCGCTTCTCGAGTGCTTCGACATCTTCCTGGGTTAACGTGCGGTCAAGATCAACGTCGTAATAAAAACCGTTGTCAATAACCGGGCCGATTGCCATTTTGGTATGCGCCCaa                                                                                                                     >  1:1239590/1‑140 (MQ=255)              aGCAAGCTCATGCATCCGCTTCTCGAGTGCTTCGACATCTTCCTGGGTTAACGTGCGGTCAAGAACAACGTCGTAATAAAAACCGTTGTCAATAACCGGGCCGATTGCCATTTTGGTATGCGCCCAAAGTTGTTTAATcg                                                                                                        <  2:1434806/140‑1 (MQ=255)                                  tctcGAGTGCTTCGACATCTTCCTGGGTTAACGTGCGGTCAAGATCAACGTCGTAATAAAAACCGTTGTCAATAACCGGGCCGATTGCCATTTTGGTATGCGCCCAAAGTTGTTTAATCGCGTGCCCTAACAGGTGCGca                                                                                    <  2:1399134/140‑1 (MQ=255)                                    tcGAGTGCTTCGACATCTTCCTGGGTTAACGTGCGGTCAAGATCAACGTCGTAATAAAAACCGTTGTCAATAACCGGGCCGATTGCCATTTTGGTATGCGCCCAAAGTTGTTTAATCGCGTGCCCTAACAGGTGCGcaca                                                                                  <  2:1239590/140‑1 (MQ=255)                                         tGCTTCGACATCTTCCTGGGTTAACGTGCGGTCAAGATCAACGTCGTAATAAAAACCGTTGTCAATAACCGGGCCGATTGCCATTTTGGTATGCGCCCAAAGTTGTTTAATCGCGTGCCCTAACAGGTGCGCACAGGAGt                                                                             <  1:1441169/140‑1 (MQ=255)                                           cTTCGACATCTTCCTGGGTTAACGTGCGGTCAAGATCAACGTCGTAATAAAAACCGTTGTCAATAACCGGGCCGATTGCCATTTTGGTATGCGCCCAAAGTTGTTTAATCGCGTGCCCTAACAGGTGCGCACAGGAGTGa                                                                           >  1:1402335/1‑140 (MQ=255)                                                      tCCTGGGTTAACGTGCGGTCAAGATCAACGTCGTAATAAAAACCGTTGTCAATAACCGGGCCGATTGCCATTTTGGTATGCGCCCAAAGTTGTTTAATCGCGTGCCCTAACAGGTGCGCACAGGAGTGACGAATGAtctc                                                                >  2:1525794/1‑140 (MQ=255)                                                                     cGGTCAAGATCAACGTCGTAATAAAAACCGTTGTCAATAACCGGGCCGATTGCCATTTTGGTATGCGCCCAAAGTTGTTTAATCGCGTGCCCTAACAGGTGCGCACAGGAGTGACGAATGATCTCCAGACCTTCTTCGTc                                                 >  2:1282361/1‑140 (MQ=255)                                                                             aTCAACGTCGTAATAAAAACCGTTGTCAATAACCGGGCCGATTGCCATTTTGGTATGCGCCCAAAGTTGTTTAATCGCGTGCCCTAACAGGTGCGCACAGGAGTGACGAATGATCTCCAGACCTTCTTCGTCTTTGGCgg                                         >  2:1535123/1‑140 (MQ=255)                                                                                  cgtcgtAATAAAAACCGTTGTCAATAACCGGGCCGATTGCCATTTTGGTATGCGCCCAAAGTTGTTTAATCGCGTGCCCTAACAGGTGCGCACAGGAGTGACGAATGATCTCCAGACCTTCTTCGTCTTTGGCGGTAATg                                    >  1:1449200/1‑140 (MQ=255)                                                                                  cgtcgtAATAAAAACCGTTGTCAATAACCGGGCCGATTGCCATTTTGGTATGCGCCCAAAGTTGTTTAATCGCGTGCCCTAACAGGTGCGCACAGGAGTGACGAATGATCTCCAGACCTTCTTCGTCTTTGGCGGTAATg                                    >  2:1649591/1‑140 (MQ=255)                                                                                  cgtcgtAATAAAAACCGTTGTCAATAACCGGGCCGATTGCCATTTTGGTATGCGCCCAAAGTTGTTTAATCGCGTGCCCTAACAGGTGCGCACAGGAGTGACGAATGATCTCCAGACCTTCTTCGTCTTTGGCGGTAATg                                    >  2:1448787/1‑140 (MQ=255)                                                                                      gtAATAAAAACCGTTGTCAATAACCGGGCCGATTGCCATTTTGGTATGCGCCCAAAGTTGTTTAATCGCGTGCCCTAACAGGTGCGCACAGGAGTGACGAATGATCTCCAGACCTTCTTCGTCTTTGGCGGTAATGATCg                                >  1:1639883/1‑140 (MQ=255)                                                                                       taataaAAACCGTTGTCAATAACCGGGCCGATTGCCATTTTGGTATGCGCCCAAAGTTGTTTAATCGCGTGCCCTAACAGGTGCGCACAGGAGTGACGAATGATCTCCAGACCTTCTTCGTCTTTGGCGGTAATGATCGa                               >  1:1363492/1‑140 (MQ=255)                                                                                                  gTTGTCAATAACCGGGCCGATTGCCATTTTGGTATGCGCCCAAAGTTGTTTAATCGCGTGCCCTAACAGGTGCGCACAGGAGTGACGAATGATCTCCAGACCTTCTTCGTCTTTGGCGGTAATGATCGACAGTTGTGcgt                    <  1:1282361/140‑1 (MQ=255)                                                                                                  gTTGTCAATAACCGGGCCGATTGCCATTTTGGTATGCGCCCAAAGTTGTTTAATCGCGTGCCCTAACAGGTGCGCACAGGAGTGACGAATGATCTCCAGACCTTCTTCGTCTTTGGCGGTAATGATCGACAGTTGTGcgt                    <  2:1402335/140‑1 (MQ=255)                                                                                                     gTCAATAACCGGGCCGATTGCCATTTTGGTATGCGCCCAAAGTTGTTTAATCGCGTGCCCTAACAGGTGCGCACAGGAGTGACGAATGATCTCCAGACCTTCTTCGTCTTTGGCGGTAATGATCGACAGTTGTGcgtcgt                 >  1:1771869/1‑140 (MQ=255)                                                                                                             ccGGGCCGATTGCCATTTTGGTATGCGCCCAAAGTTGTTTAATCGCGTGCCCTAACAGGTGCGCACAGGAGTGACGAATGATCTCCAGACCTTCTTCGTCTTTGGCGGTAATGATCGACAGTTGTGCGTCGTTTTCAATc         <  1:1594791/140‑1 (MQ=255)                                                                                                                ggCCGATTGCCATTTTGGTATGCGCCCAAAGTTGTTTAATCGCGTGCCCTAACAGGTGCGCACAGGAGTGACGAATGATCTCCAGACCTTCTTCGTCTTTGGCGGTAATGATCGACAGTTGTGCGTCGTTTTCAATCAGa      >  2:1556272/1‑140 (MQ=255)                                                                                                                ggCCGATTGCCATTTTGGTATGCGCCCAAAGTTGTTTAATCGCGTGCCCTAACAGGTGCGCACAGGAGTGACGAATGATCTCCAGACCTTCTTCGTCTTTGGCGGTAATGATCGACAGTTGTGCGTCGTTTTCAATCAGa      >  2:1608052/1‑140 (MQ=255)                                                                                                                ggCCGATTGCCATTTTGGTATGCGCCCAAAGTTGTTTAATCGCGTGCCCTAACAGGTGCGCACAGGAGTGACGAATGATCTCCAGACCTTCTTCGTCTTTGGCGGTAATGATCGACAGTTGTGCGTCGTTTTCAATCAGa      >  2:1718287/1‑140 (MQ=255)                                                                                                                    gATTGCCATTTTGGTATGCGCCCAAAGTTGTTTAATCGCGTGCCCTAACAGGTGCGCACAGGAGTGACGAATGATCTCCAGACCTTCTTCGTCTTTGGCGGTAATGATCGACAGTTGTGCGTCGTTTTCAATCAGATCGc  <  1:1525794/140‑1 (MQ=255)                                                                                                                                        |                                                                                                                        CGTAGTTTTTCTCAGCAAGCTCATGCATCCGCTTCTCGAGTGCTTCGACATCTTCCTGGGTTAACGTGCGGTCAAGATCAACGTCGTAATAAAAACCGTTGTCAATAACCGGGCCGATTGCCATTTTGGTATGCGGCCAAAGTTGTTTAATCGCGTGCCCTAACAGGTGCGCACAGGAGTGACGAATGATCTCCAGACCTTCTTCGTCTTTGGCGGTAATGATCGACAGTTGTGCGTCGTTTTCAATCAGATCGC  >  NC_000913/1802174‑1802428

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