Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,046,391	T→G	Q884P (CAG→CCG)	gcvP ←	glycine decarboxylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,046,391	0	T	G	100.0%	110.4 / NA	29	Q884P (CAG→CCG)	gcvP	glycine decarboxylase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (15/14);  total (15/14)

CCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGCGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTCAGTCGGTTCAACCATCAGCGTACCCGCCACCGGGAACGACATCGTCGGCG  >  NC_000913/3046263‑3046511                                                                                                                                 |                                                                                                                         ccACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGCGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCGGGTCAATTTCt                                                                                                               >  1:354369/1‑140 (MQ=255) ccACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGCGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCGGGTCAATTTCt                                                                                                               >  2:217442/1‑140 (MQ=255)     acCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGCGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCGGGTCAACTTCTgcgc                                                                                                           <  2:555872/140‑1 (MQ=255)            gggAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGAGGCGCGTTCACCAGCGGGTTATCTTCCAGCCGCCAGACGCCGGCTTTCACCGGGTCAATTTCTGCGCGGATAGc                                                                                                    <  2:98654/140‑1 (MQ=255)                  aCCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGCGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCGGGTCAATTTCTGCGCGGATAGCCAGCAt                                                                                              <  2:42333/140‑1 (MQ=255)                    cGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGCGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCGGGTCAATTTCTGCGCGGATAGCCAGCATcg                                                                                            <  1:319910/140‑1 (MQ=255)                         cTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGCGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCGGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCg                                                                                       >  1:340904/1‑140 (MQ=255)                               ggCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGCGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCGGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAg                                                                                 <  2:706152/140‑1 (MQ=255)                                gCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGCGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCGGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGc                                                                                <  1:714036/140‑1 (MQ=255)                                  tGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGCGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCGGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGa                                                                              <  1:630852/140‑1 (MQ=255)                                    tACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGCGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCGGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATc                                                                            >  1:295321/1‑140 (MQ=255)                                    tACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGCGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCGGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATc                                                                            >  2:708061/1‑140 (MQ=255)                                           gcgcCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGCGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCGGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTcc                                                                     >  2:658262/1‑140 (MQ=255)                                            cgcCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGCGTTCACCAGCGGGTTATCTTCCAGCGGTCAGACACCGGCTTTCACCGGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCa                                                                    >  1:398501/1‑140 (MQ=255)                                              cccACTCGGCGCCCAGTTCGCTCTGAATGTGCGGCGCGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCGGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACt                                                                  <  2:254041/140‑1 (MQ=255)                                              cccACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGCGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCGGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACt                                                                  >  1:707336/1‑140 (MQ=255)                                              cccACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGCGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCGGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACt                                                                  >  2:390896/1‑140 (MQ=255)                                                 aCTCGGCGACCAGTTCGCTCTGAATGTGCGGCGCGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCGGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTg                                                               <  1:443841/140‑1 (MQ=255)                                                         aCCAGTTCGCTCTGAATGTGCGGCGCGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCGGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGa                                                       <  1:658262/140‑1 (MQ=255)                                                           cAGTTCGCTCTGAATGTGCGGCGCGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCGGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGAtt                                                     >  2:597506/1‑140 (MQ=255)                                                                           tgCGGCGCGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCGGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTCAGTCGGTTCAACCAt                                     <  1:752245/140‑1 (MQ=255)                                                                                    ttCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCGGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTCAGTCGGTTCAACCATCAGCGTAcc                            >  1:22344/1‑140 (MQ=255)                                                                                           gCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCGGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTCAGTCGGTTCAACCATCAGCGTACCCGCCAcc                     <  2:295321/140‑1 (MQ=255)                                                                                              ggTTATCTTCCAGCGGCCAGACACCGGCTTTCACCGGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTCAGTCGGTTCAACCATCAGCGTACCCGCCACCggg                  <  2:354369/140‑1 (MQ=255)                                                                                                      tCCAGCGGCCAGACACCGGCTTTCACCGGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTCAGTCGGTTCAACCATCAGCGTACCCGCCACCGGGAACGACAt          >  2:156957/1‑140 (MQ=255)                                                                                                         aGCGGCCAGACACCGGCTTTCACCGGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTCAGTCGGTTCAACCATCAGCGTACCCGCCACCGGGAACGACAtcgt       >  1:580246/1‑140 (MQ=255)                                                                                                              ccAGACACCGGCTTTCACCGGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTCAGTCGGTTCAACCATCAGCGTACCCGCCACCGGGAACGACAtcgtc      >  2:699790/1‑136 (MQ=255)                                                                                                              ccAGACACCGGCTTTCACCGGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTCAGTCGGTTCAACCATCAGCGTACCCGCCACCGGGAACGACAtcgtc      <  1:699790/136‑1 (MQ=255)                                                                                                              ccAGACACCGGCTTTCACCGGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTCAGTCGGTTCAACCATCAGCGTACCCGCCACCGGGAACGACATCGTCGgcg  >  2:693380/1‑140 (MQ=255)                                                                                                                                 |                                                                                                                         CCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGCGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTCAGTCGGTTCAACCATCAGCGTACCCGCCACCGGGAACGACATCGTCGGCG  >  NC_000913/3046263‑3046511

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