Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,258,693	(C)5→6	coding (957/1365 nt)	tdcG ←	L‑serine deaminase III

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,258,688	1	.	C	100.0%	122.9 / NA	34	A321G (GCT→GGT)	tdcG	L‑serine deaminase III
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base C (17/17);  total (17/17)

TACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCAATCTCCCCCTGACAGCCGACTTCCGCGCCAGAGATGGAGGCGTTCATTTTATACAGCGCGCCAATAG‑CCCCCGCGGCCAGAAAATAGCGGGCAATTGACCGCTCGTTTACCGGACGACGGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCA  >  NC_000913/3258555‑3258819                                                                                                                                       |                                                                                                                                    tACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTACCCCTGCCGCCGCCATTGAACAGGCCACGCCAATCTCCCCCTGACAGCCGACTTCCGCGCCAGAGATGGAGGCGTTCATTTTATACAGCGCGCCAATAG‑cccccc                                                                                                                               <  1:453337/140‑2 (MQ=255)     tGCGCCGGACTGCCGCCCAGTAGTTCAGTTACCCCTGCCGCCGCCATTGAACAGGCCACGCCAATCTCCCCCTGACAGCCGACTTCCGCGCCAGAGATGGAGGCGTTCATTTTATACAGCGCGCCAATAG‑ccccccgcgg                                                                                                                           <  2:417157/140‑6 (MQ=255)     tGCGCCGGACTGCCGCCCAGGAGTTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCAATCTCCCCCTGACAGCCGACTTCCGCGCCAGAGATGGAGGCGTTCATTTTATACAGCGCGCCAATAG‑ccccccgcgg                                                                                                                           <  2:402082/140‑6 (MQ=255)       cgcCGGACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCAATCTCCCCCTGACAGCCGACTTCCGCGCCAGAGATGGAGGCGTTCATTTTATACAGCGCGCCAATAGCCCCCCGCGGcc                                                                                                                          >  2:315422/1‑140 (MQ=255)         ccGGACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCAATCTCCCCCTGACAGCCGACTTCCGCGCCAGAGATGGAGGCGTTCATTTTATACAGCGCGCCAATAGCCCCCCGCGGCCAg                                                                                                                        >  1:420511/1‑140 (MQ=255)                          agtTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCAATCTCCCCCTGACAGCCGACTTCCGCGCCAGAGATGGAGGCGTTCATTTTATACAGCGCGCCAATAGCCCCCCGCGGCCAGAAAATAGCGGGCAATTg                                                                                                       <  2:81494/140‑1 (MQ=255)                            ttCAGTTAGCCCTGCCGCCGCCATTGAACAGGCCACGCCAATCTCCCCCTGACAGCCGACTTCCGCGCCAGAGATGGAGGCGTTCATTTTATACAGCGCGCCAATAGCCCCCCGCGGCCAGAAAATAGCGGGCAATTGAc                                                                                                     <  1:365450/140‑1 (MQ=255)                                         gccgccgccATTGAACAGGCCACGCCAATCTCCCCCTGACAGCCGACTTCCGCGCCAGAGATGGAGGCGTTCATTTTATACAGCGCGCCAATAGCCCCCCGCGGCCAGAAAATAGCGGGCAATTGACCGCTCGTTTACCg                                                                                        >  2:301111/1‑140 (MQ=255)                                           cgccgccATTGAACAGGCCACGCCAATCTCCCCCTGACAGCCGACTTCCGCGCCAGAGATGGAGGCGTTCATTTTATACAGCGCGCCAATAGCCCCCCGCGGCCAGAAAATAGCGGGCAATTGACCGCTCGTTTACCGga                                                                                      <  2:147482/140‑1 (MQ=255)                                                       aCAGGCCACGCCAATCTCCCCCTGACAGCCGACTTCCGCGCCAGAGATGGAGGCGTTCATTTTATACAGCGCGCCAATAGCCCCCCGCGGCCAGAAAATAGCGGGCAATTGACCGCTCGTTTACCGGACGACGGAACTTa                                                                          <  1:56403/140‑1 (MQ=255)                                                        cAGGCCACGCCAATCTCCCCCTGACAGCCGACTTCCGCGCCAGAGATGGAGGCGTTCATTTTATACAGCGCGCCAATAGCCCCCCGCGGCCAGAAAATAGCGGGCAATTGACCGCTCGTTTACCGGACGACGGAACTTAt                                                                         <  2:251367/140‑1 (MQ=255)                                                           gCCACGCCAATCTCCCCCTGACAGCCGACTTCCGCGCCAGAGATGGAGGCGTTCATTTTATACAGCGCGCCAATAGCCCCCCGCGGCCAGAAAATAGCGGGCAATTGACCGCTCGTTTACCGGACGACGGAACTTATCGt                                                                      >  2:226168/1‑140 (MQ=255)                                                             cACGCCAATCTCCCCCTGACAGCCGACTTCCGCGCCAGAGATGGAGGCGTTCATTTTATACAGCGCGCCAATAGCCCCCCGCGGCCAGAAAATAGCGGGCAATTGACCGCTCGTTTACCGGACGACGGAACTTATCGtaa                                                                    <  1:226168/140‑1 (MQ=255)                                                             cACGCCAATCTCCCCCTGACAGCCGACTTCCGCGCCAGAGATGGAGGCGTTCATTTTATACAGCGCGCCAATAGCCCCCCGCGGCCAGAAAATAGCGGGCAATTGACCGCTCGTTTACCGGACGACGGAACTTATCGtaa                                                                    <  2:457458/140‑1 (MQ=255)                                                                gCCAATCTCCCCCTGACAGCCGACTTCCGCGCCAGAGATGGAGGCGTTCATTTTATACAGCGCGCCAATAGCCCCCCGCGGCCAGAAAATAGCGGGCAATTGACCGCTCGTTTACCGGACGACGGAACTTATCGtaataa                                                                 >  1:220220/1‑140 (MQ=255)                                                                     tctcCCCCTGACAGCCGACTTCCGCGCCAGAGATGGAGGCGTTCATTTTATACAGCGCGCCAATAGCCCCCCGCGGCCAGAAAATAGCGGGCAATTGACCGCTCGTTTACCGGACGACGGAACTTATCGTAATAAGCCAg                                                            >  1:159571/1‑140 (MQ=255)                                                                          cccTGACAGCCGACTTCCGCGCCAGAGATGGAGGCGTTCATTTTATACAGCGCGCCAATAGCCCCCCGCGGCCAGAAAATAGCGGGCAATTGACCGCTCGTTTACCGGACGACGGAACTTATCGTAATAAGCCAGTACTg                                                       >  1:35817/1‑140 (MQ=255)                                                                             tGACAGCCGACTTCCGCGCCAGAGATGGAGGCGTTCATTTTATACAGCGCGCCAATAGCCCCCCGCGGCCAGAAAATAGCGGGCAATTGACCGCTCGTTTACCGGACGACGGAACTTATCGTAATAAGCCAGTACTGCCg                                                    >  2:142564/1‑140 (MQ=255)                                                                                 aGCCGACTTCCGCGCCAGAGATGGAGGCGTTCATTTTATACAGCGCGCCAATAGCCCCCCGCGGCCAGAAAATAGCGGGCAATTGACCGCTCGTTTACCGGACGACGGAACTTATCGTAATAAGCCAGt                                                           <  2:228461/129‑1 (MQ=255)                                                                                 aGCCGACTTCCGCGCCAGAGATGGAGGCGTTCATTTTATACAGCGCGCCAATAGCCCCCCGCGGCCAGAAAATAGCGGGCAATTGACCGCTCGTTTACCGGACGACGGAACTTATCGTAATAAGCCAGt                                                           >  1:228461/1‑129 (MQ=255)                                                                                   ccGACTTCCGCGCCAGAGATGGAGGCGTTCATTTTATACAGCGCGCCAATAGCCCCCCGCGGCCAGAAAATAGCGGGCAATTGACCGCTCGTTTACCGGACGACGGAACTTATCGTAATAAGCCAGTACTGCCGGaataa                                              >  2:81910/1‑140 (MQ=255)                                                                                         tCCGCGCCAGAGATGGAGGCGTTCATTTTATACAGCGCGCCAATAGCCCCCCGCGGCCAGAAAATAGCGGGCAATTGACCGCTCGTTTACCGGACGACGGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCcgc                                        >  2:462659/1‑140 (MQ=255)                                                                                          ccGCGCCAGAGATGGAGGCGTTCATTTTATACAGCGCGCCAATAGCCCCCCGCGGCCAGAAAATAGCGGGCAATTGACCGCTCGTTTACCGGACGACGGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCcgca                                       <  1:142564/140‑1 (MQ=255)                                                                                          ccGCGCCAGAGATGGAGGCGTTCATTTTATACAGCGCGCCAATAGCCCCCCGCGGCCAGAAAATAGCGGGCAATTGACCGCTCGTTTACCGGACGACGGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCcgca                                       <  2:255226/140‑1 (MQ=255)                                                                                            gcgcCAGAGATGGAGGCGTTCATTTTATACAGCGCGCCAATAGCCCCCCGCGGCCAGAAAATAGCGGGCAATTGACCGCTCGTTTACCGGACGACGGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCcgcacg                                     <  1:431706/140‑1 (MQ=255)                                                                                                 agagATGGAGGCGTTCATTTTATACAGCGCGCCAATAGCCCCCCGCGGCCAGAAAATAGCGGGCAATTGACCGCTCGTTTACCGGACGACGGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCg                                <  2:91056/140‑1 (MQ=255)                                                                                                       ggAGGCGTTCATTTTATACAGCGCGCCAATAGCCCCCCGCGGCCAGAAAATAGCGGGCAATTGACCGCTCGTTTACCGGACGACGGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTc                          <  1:315422/140‑1 (MQ=255)                                                                                                         aGGCGTTCATTTTATACAGCGCGCCAATAGCCCCCCGCGGCCAGAAAATAGCGGGCAATTGACCGCTCGTTTACCGGACGACGGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCgg                        <  1:335402/140‑1 (MQ=255)                                                                                                          ggCGTTCATTTTATACAGCGCGCCAATAGCCCCCCGCGGCCAGAAAATAGCGGGCAATTGACCGCTCGTTTACCGGACGACGGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGt                       >  2:407003/1‑140 (MQ=255)                                                                                                              ttCATTTTATACAGCGCGCCAATAGCCCCCCGCGGCCAGAAAATAGCGGGCAATTGACCGCTCGTTTACCGGACGACGGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCg                   >  1:485955/1‑140 (MQ=255)                                                                                                                      ataCAGCGCGCCAATAGCCCCCCGCGGCCAGAAAATAGCGGGCAATTGACCGCTCGTTTACCGGACGACGGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCAcg           <  1:301111/140‑1 (MQ=255)                                                                                                                      ataCAGCGCGCCAATAGCCCCCCGCGGCCAGAAAATAGCGGGCAATTGACCGCTCGTTTACCGGACGACGGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCAcg           <  2:220220/140‑1 (MQ=255)                                                                                                                            cgcgcCAATAGCCCCCCGCGGCCAGAAAATAGCGGGCAATTGACCGCTCGTTTACCGGACGACGGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGcgcccg     >  2:332964/1‑140 (MQ=255)                                                                                                                            cgcgcCAATAGCCCCCCGCGGCCAGAAAATAGCGGGCAATTGACCGCTCGTTTACCGGACGACGGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGcgcccg     >  2:252591/1‑140 (MQ=255)                                                                                                                            cgcgcCAATAGCCCCCCGCGGCCAGAAAATAGCGGGCAATTGACCGCTCGTTTACCGGACGACGGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGcgcccg     >  1:340499/1‑140 (MQ=255)                                                                                                                            cgcgcCAATAGCCCCCCGCGGCCAGAAAATAGCGGGCAATTGACCGCTCGTTTACCGGACGACGGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGcgcccg     >  2:231589/1‑140 (MQ=255)                                                                                                                               gcCAATAGCCCCCCGCGGCCAGAAAATAGCGGGCAATTGACCGCTCGTTTACCGGACGACGGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCa  <  1:332964/140‑1 (MQ=255)                                                                                                                                       |                                                                                                                                    TACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCAATCTCCCCCTGACAGCCGACTTCCGCGCCAGAGATGGAGGCGTTCATTTTATACAGCGCGCCAATAG‑CCCCCGCGGCCAGAAAATAGCGGGCAATTGACCGCTCGTTTACCGGACGACGGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCA  >  NC_000913/3258555‑3258819

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: ‑