Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,309,390	A→G	G593G (GGT→GGC)	pnp ←	polynucleotide phosphorylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,309,390	0	A	G	100.0%	84.7 / NA	26	G593G (GGT→GGC)	pnp	polynucleotide phosphorylase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (16/10);  total (16/10)

ACGGGTCACTTTACCAGTGTAGACGCGGCCCACTTCGATTTCTGCAGTGATCTCTTCGATACGACGAATAGCATGTTTCGCTTTCTCGCCGTCGGTCGCTGCAATCTTCACAGTACCGTCATCTTCGATTTCGATGGTGGTGCCAGTTTCTTCGGTCAGAGCACGGATTACAGAACCGCCTTTACCGATAACATCTTTGATCTTGTCCGGGTTGATCTTGATGGTATGGATACGCGGTGCGA  >  NC_000913/3309276‑3309517                                                                                                                   |                                                                                                                                aCGGGTCACTTTACCAGTGTAGACGCTGCCCACTTCGATTTCTGCAGTGATCTCTTCGATACGACGAATAGCATGTTTCGCTTTCTCGCCGTCGGTCGCTGCAATCTTCACAGTGCCGTCATCttcgatttcg                                                                                                               >  1:234180/1‑133 (MQ=255) aCGGGTCACTTTACCAGTGTAGACGCGGCCCACTTCGATTTCTGCAGTGATCTCTTCGATACGACGAATAGCATGTTTCGCTTTCTCGCCGTCGGTCGCTGCAATGTTCACAGTGCCGTCATCttcgattttg                                                                                                               >  1:619456/1‑131 (MQ=255) aCGGGTCACTTTACCAGTGTAGACGCGGCCCACTTCGATTTCTGCAGTGATCTCTTCGATACGACGAATAGCATGTTTCGCTTTCTCGCCGTCGGTCGCTGCAATCTTCACAGTGCCGTCATCttcaatttcg                                                                                                               <  2:619456/133‑1 (MQ=255)   gggTCACTTTACCAGTGTAGACGCGGCCCACTTCGATTTCTGCAGTGATCTCTTCGATACGACGAATAGCATGTTTCGCTTTCTCGCCGTCGGTCGCTGCAATCTTCACAGTGCCGTCATCttcgatttcgat                                                                                                             >  2:570020/1‑133 (MQ=255)     cacggattACCAGTGTAGACGCGGCCCACTTCGATTTCTGCAGTGATCTCTTCGATACGACGAATAGCATGTTTCGCTTTCTCGCCGTCGGTCGATGCAATCTTCACAGTGCCGTCATCTTCGATTTCGAtgg                                                                                                           <  2:15543/127‑1 (MQ=255)                    tAGACGCGGCCCACTTCGATTTCTGCAGTGATCTCTTCGATACGACGAATAGCATGTTTCGCTTTCTCGCCGTCGGTCGCTGCAATCTTCACAGTGCCGTCATCTTCGATTTCGATGGTGGTGCCAGTttct                                                                                             >  2:491626/1‑132 (MQ=255)                    tAGACGCGGCCCACTTCGATTTCTGCAGTGATCTCTTCGATACGACGAATAGCATGTTTCGCTTTCTCGCCGTCGGTCGCTGCAATCTTCACAGTGCCGTCATCTTCGATTTCGATGGTGGTGCCAGTttct                                                                                             <  1:491626/132‑1 (MQ=255)                              ccACTTCGGTTTCTGTAGTGATCTCTTCGATACGACGAATAGCATGTTTCGCTTTCTCGCCGTCGGTCGCTGCAATCTTCACAGTGCCGTCATCTTCGATTTCGATGGTGGTGCCAGTTTCTTCGGTCAGAGc                                                                                  <  1:231591/133‑1 (MQ=255)                              ccACTTCGCTTTCTGAAGGGATCTCTTCGACACGACGAATAGCATGTTTCGCTTTCTCGCCGTCGGTCGCTGCAATCTTCACAGTGCCGTCATCTTCGATTTCGATGGTGGTGCCAGTTTCTTCGGTCAGAGc                                                                                  <  1:169986/133‑1 (MQ=255)                                       tttCTGCAGTGATCTCTTCGATACGACGAATAGCATGTTTCGCTTTCTCGCCGTCGGTCGCTGCAATCTTCACAGTGCCGTCATCTTCGATTTCGATGGTGGAGCCAGTTGCTTAGGGCAGAGCGCGGATTAc                                                                         >  1:164405/1‑133 (MQ=255)                                                  aTCTCTTCGATACGACGAATAGCATGTTTCGCTTTCTCGCCGTCGGTCGCTGCAATCTTCACAGTGCCGTCATCTTCGATTTCGATGGTGGTGCCAGTTTCTTCGGTCAGAGCACGGTTTACAGAACCGCCtt                                                              >  1:247516/1‑133 (MQ=255)                                                     tctTCGATACGACGAATAGCATGTTTCGCTTTCTCGCCGTCGGTCGCTGCAATCTTCACAGTGCCGTCATCTTCGATTTCGATGGTGGTGCCAGTTTCTTCGGTCAGAGCACGGATTACAGaa                                                                     >  2:417628/1‑123 (MQ=255)                                                     tctTCGATACGACGAATAGCATGTTTCGCTTTCTCGCCGTCGGTCGCTGCAATCTTCACAGTGCCGTCATCTTCGATTTCGATGGTGGTGCCAGTTTCTTCGGTCAGAGCACGGATTACAGaa                                                                     <  1:417628/123‑1 (MQ=255)                                                         cGAAACGACGTATAGTATGTTTCGCTTTTTCGTCGTCGGTCACTGCAATCTTCACAGTGCCGTCATCTTCGATTTCGATGGTGGTGCCAGTTTCTTCGGTCAGAGCACGGATTACAGAACCGCCTTTACCGAt                                                       <  2:247516/133‑1 (MQ=255)                                                          gATACGACGAATAGCATGTTTCGCTTTCTCGCCGTCGGTCGCTGCAATCTTCACAGTGCCGTCATCTTCGATTTCGATGGTGGTGCCAGTTTCTTCGGTCAGAGCACGGAATACAGAACCGCCTCTACCGATa                                                      >  1:266997/1‑133 (MQ=255)                                                            tACGACGAATAGCATGTTTCGCTTTCTCGCCGTCGGTCGCTGCAATCTTCACAGTGCCGTCATCTTCGATTTCGATGGTGGTGCCAGTTTCTTCGGTCAGAGCACGGATTACAGAACCGCCTTTACCGATAAc                                                    >  2:41836/1‑133 (MQ=255)                                                            tACGACGAATAGCATGTTTCGCTTTCTCGCCGTCGGTCGCTGCAATCTTCACAGTGCCGTCATCTTCGATTTCGATGGTGGTGCCAGTTTCTTCGGTCAGAGCACGGATTACAGAACCGCCTTTACCGATAAc                                                    <  1:41836/133‑1 (MQ=255)                                                               gacgaATAGCATGTTTCGCTTTCTCGCCGTCGGTCGCTGCAATCTTCACAGTGCCGTCATCTTCGATTTCGATGGTGGTGCCAGTTTCTTCGGTCAGAACACGGATTACAGAACCGCCTTTACCGATAACATc                                                 >  2:502242/1‑133 (MQ=255)                                                                 cgaATAGCATGTTTCGCTTTCTCGCCGTCGGTCGCTGCAATCTTCACAGTGCCGTCATCTTCGATTTCGATGGTGGTGCCAGTTTCTTCGGTCAGAGCACGGATTACAGAACCGCCTTTACCGATAACATCtt                                               >  1:475424/1‑133 (MQ=255)                                                                             ttCGCTTTCTCGCCGTCGGTCGCTGCAATCTTCACAGTGCCGTCATCTTCGATTTCGATGGTGGTGCCAGTTTCTTCGGTCAGAGCACGGATTACAGAACCGCCTTTACCGATAACATCTTTGATCTTGTCCg                                   >  1:481872/1‑133 (MQ=255)                                                                             ttCGCTTTCTCGCCGTCGGTCGCTGCAATCTTCACAGTGCCGTCATCTTCGATTTCGATGGTGGTGCCAGTTTCTTCGGTCAGAGCACGGATTACAGAACCGCCTTTACCCATATCATCTTTGAGCTTGTCCg                                   >  1:291991/1‑133 (MQ=255)                                                                                 cTTTCTCGCCGTCGGTCGCTGCAATCTTCACAGTGCCGTCATCTTCGATTTCGATGGTGGTGCCAGTTTCTTCGGTCAGAGCACGGATTACAGAACCGCCTTTACAGATAACATCTTTGATCTTGTCCGGGtt                               >  1:242995/1‑133 (MQ=255)                                                                                     ctcGCCGTCGGTCGCTGCAATCTTCACAGTGCCGTCATCTTCGATTTCGATGGTGGTGCCAGTTTCTTCGGTCAGAGCACGGATTACAGAACCGCCTTTACCGATAACATATTTGATCTTGTCCGGGTTGATc                           >  1:205279/1‑133 (MQ=255)                                                                                       cGCCGTCGGTCGCTGCAATCTTCACAGTGCCGTCATCTTCGATTTCGATGGTGGTGCCAGTTTCTTCGGTCAGAGCACGGATTACAGAACCGCCTTTACCGATAACATCTTTGATCTTGTCCGGGTTGATCtt                         <  2:475424/133‑1 (MQ=255)                                                                                         ccGTCGGTCGTTGCAATCTTCACAGTGCCGTCATCTTCGATTTCGATGGTGGTGCCAGTTTCTTCGGTCAGAGCACGGATTACAGAACCGCCTTTACCGATAACATCTTTGATCTTGTCCGGGTTGATCTTGa                       <  2:603156/133‑1 (MQ=255)                                                                                                tcgCTGCAATCTTCACAGTGCCGTCATCTTCGATTTCGATGGTGGTGCCAGTTTCTTCGGTCAGAGCACGGATTACAGAACCGCCTTTACCGATAACATCTTTGATCTTGTCCGGGTTGATCTTGATGGTATg                >  2:117326/1‑133 (MQ=255)                                                                                                       aaTCTTCACACTGCCATCATCTTCAACTCCGATGGTGGTGCAAGTTTCTTCGGTCAGAGCACGGATTACAGAACCGCCTTTACCGATAACATCTTTGATCTTGTCCGGGTTGATCTTGATGGTATGGATAcgc         <  2:266997/133‑1 (MQ=255)                                                                                                              acaGTGCCGTCATCTTCGATTTCGATGGTGGTGCCAGTTTCTTCGGTCAGAGCACGGATTACAGAACCGCCTTTACCGATAACATCTTTGATCTTGTCCGGGTTGATCTTGATGGTATGGATACGCGGTGCGa  >  1:56728/1‑133 (MQ=255)                                                                                                                   |                                                                                                                                ACGGGTCACTTTACCAGTGTAGACGCGGCCCACTTCGATTTCTGCAGTGATCTCTTCGATACGACGAATAGCATGTTTCGCTTTCTCGCCGTCGGTCGCTGCAATCTTCACAGTACCGTCATCTTCGATTTCGATGGTGGTGCCAGTTTCTTCGGTCAGAGCACGGATTACAGAACCGCCTTTACCGATAACATCTTTGATCTTGTCCGGGTTGATCTTGATGGTATGGATACGCGGTGCGA  >  NC_000913/3309276‑3309517

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