Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	521,750	A→G	R141R (AGA→AGG)	SAUSA300_RS02475 →	DNA polymerase III subunit delta'

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	521,750	0	A	G	100.0%	109.5 / NA	40	R141R (AGA→AGG)	SAUSA300_RS02475	DNA polymerase III subunit delta'
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (39/1);  total (39/1)

TGAAATTTCTTGAAGAACCACCGGACAATACGATTGCTATTTTATTGTCTACAAAACCTGAGCAAATTTTAGACACAATCCATTCAAGATGTCAGCATGTATATTTCAAGCCTATTGATAAAGAAAAGTTTATAAATAGATTAGTTGAACAAAACATGTCTAAGCCAGTAGCTGAAATGATTAG  >  NC_007793/521662‑521845                                                                                         |                                                                                                tGAAATTTCTTGAAGAACCACCGGACAATACGATTGCTATTTTATTGTCTACAAAACCTGAGCAAATTTTAGACACAATCCATTCAAGGTGTCAGCATGt                                                                                      >  1:690818/1‑100 (MQ=255) tGAAATTTCTTGAAGAACCACCGGACAATACGATTGCTATTTTATTGTCTACAAAACCTGAGCAAATTTTAGACACAATCCATTCAAGGTGTCAGCATGt                                                                                      >  1:269998/1‑100 (MQ=255)  gAAATTTCTTGAAGAACCACCGGACAATACGATTGCTATTTTATTGTCTACAAAACCTGAGCAAATTTTAGACACAATCCATTCAAGGTGTCAGCATGta                                                                                     >  1:713798/1‑100 (MQ=255)        gcTTGAAGAACCACCGGACAATACGATTGCTATTTTATTGTCTACAAAACCTGAGCAAATTTTAGACACAATCCATTCAAGGTGTCAGCATGTATATTTc                                                                               >  1:2154198/2‑100 (MQ=255)           tGAAGAACCACCGGACAATACGATTGCTATTTTATTGTCTACAAAACCTGAGCAAATTTTAGACACAATCCATTCAAGGTGTCAGCATGTATATTTCAAg                                                                            >  1:4083178/1‑100 (MQ=255)           tGAAGAACCACCGGACAATACGATTGCTATTTTATTGTCTACAAAACCTGAGCAAATTTTAGACACAATCCATTCAAGGTGTCAGCATGTATATTTCAAg                                                                            >  1:1183252/1‑100 (MQ=255)           tGAAGAACCACCGGACAATACGATTGCTATTTTATTGTCTACAAAACCTGAGCAAATTTTAGACACAATCCATTCAAGGTGTCAGCATGTATATTTCAAg                                                                            >  1:1715869/1‑100 (MQ=255)              agaaCCACCGGACAATACGATTGCTATTTTATTGTCTACAAAACCTGAGCAAATTTTAGACACAATCCATTCAAGGTGTCAGCATGTATATTTCAAGCCt                                                                         >  1:4754455/1‑100 (MQ=255)               gaaCCACCGGACAATACGATTGCTATTTTATTGTCTACAAAACCTGAGCAAATTTTAGACACAATCCATTCAAGGTGTCAGCATGTATATTTCAAGCCTa                                                                        >  1:1438439/1‑100 (MQ=255)                gaccaccGGACAATACGATTGCTATTTTATTGTCTACAAAACCTGAGCAAATTTTAGACACAATCCATTCAAGGTGTCAGCATGTATATTTCAAGCCTAt                                                                       >  1:1113532/2‑100 (MQ=255)                 accaccGGACAATACGATTGCTATTTTATTGTCTACAAAACCTGAGCAAATTTTAGACACAATCCATTCAAGGTGTCAGCATGTATATTTCAAGCCTAtt                                                                      >  1:3495367/1‑100 (MQ=255)                 accaccGGACAATACGATTGCTATTTTATTGTCTACAAAACCTGAGCAAATTTTAGACACAATCCATTCAAGGTGTCAGCATGTATATTTCAAGCCTAtt                                                                      >  1:4748293/1‑100 (MQ=255)                  ccaccGGACAATACGATTGCTATTTTATTGTCTACAAAACCTGAGCAAATTTTAGACACAATCCATTCAAGGTGTCAGCATGTATATTTCAAGCCTAtta                                                                     <  1:4925655/100‑2 (MQ=255)                                          ttATTGTCTACAAAACCTGAGCAAATTTTAGACACAATCCATTCAAGGTGTCAGCATGTATATTTCAAGCCTATTGATAAAGAAAAGTTTATAAATAGAt                                             >  1:4575653/1‑100 (MQ=255)                                                   aCAAAACCTGAGCAAATTTTAGACACAATCCATTCAAGGTGTCAGCATGTATATTTCAAGCCTATTGATAAAGAAAAGTTTATAAATAGagatcggaaga                                    >  1:2445022/1‑90 (MQ=255)                                                   aCAAAACCTGAGCAAATTTTAGACACAATCCATTCAAGGTGTCAGCATGTATATTTCAAGCCTATTGATAAAGAAAAGTTTATAAATAGAtcggaagagc                                    >  1:622627/1‑91 (MQ=255)                                                   aCAAAACCTGAGCAAATTTTAGACACAATCCATTCAAGGTGTCAGCATGTATATTTCAAGCCTATTGATAAAGAAAAGTTTATAAATAGATTAGTTGAAc                                    >  1:4791908/1‑100 (MQ=255)                                                   aCAAAACCTGAGCAAATTTTAGACACAATCCATTCAAGGTGTCAGCATGTATATTTCAAGCCTATTGATAAAGAAAAGTTTATAAATAGATTAGTTGAAc                                    >  1:3193360/1‑100 (MQ=255)                                                   aCAAAACCTGAGCAAATTTTAGACACAATCCATTCAAGGTGTCAGCATGTATATTTCAAGCCTATTGATAAAGAAAAGTTTATAAATAGATTAGTTGAAc                                    >  1:5470871/1‑100 (MQ=255)                                                   aCAAAACCTGAGCAAATTTTAGACACAATCCATTCAAGGTGTCAGCATGTATATTTCAAGCCTATTGATAAAGAAAAGTTTATAAATAGATTAGTTGAAc                                    >  1:4586666/1‑100 (MQ=255)                                                   aCAAAACCTGAGCAAATTTTAGACACAATCCATTCAAGGTGTCAGCATGTATATTTCAAGCCTATTGATAAAGAAAAGTTTATAAATAGATTAGTTGAAc                                    >  1:5695985/1‑100 (MQ=255)                                                   aCAAAACCTGAGCAAATTTTAGACACAATCCATTCAAGGTGTCAGCATGTATATTTCAAGCCTATTGATAAAGAAAAGTTTATAAATAGATTAGTTGAAc                                    >  1:4131739/1‑100 (MQ=255)                                                   aCAAAACCTGAGCAAATTTTAGACACAATCCATTCAAGGTGTCAGCATGTATATTTCAAGCCTATTGATAAAGAAAAGTTTATAAATAGATTAGTTGAAc                                    >  1:724555/1‑100 (MQ=255)                                                   aCAAAACCTGAGCAAATTTTAGACACAATCCATTCAAGGTGTCAGCATGTATATTTCAAGCCTATTGATAAAGAAAAGTTTATAAATAGATTAGTTGAAc                                    >  1:3751381/1‑100 (MQ=255)                                                   aCAAAACCTGAGCAAATTTTAGACACAATCCATTCAAGGTGTCAGCATGTATATTTCAAGCCTATTGATAAAGAAAAGTTTATAAATAGATTAGTTGAAc                                    >  1:1223027/1‑100 (MQ=255)                                                   aCAAAACCTGAGCAAATTTTAGACACAATCCATTCAAGGTGTCAGCATGTATATTTCAAGCCTATTGATAAAGAAAAGTTTATAAATAGATTAGTTGAAc                                    >  1:2843991/1‑100 (MQ=255)                                                   aCAAAACCTGAGCAAATTTTAGACACAATCCATTCAAGGTGTCAGCATGTATATTTCAAGCCTATTGATAAAGAAAAGTTTATAAATAGATTAGTTGAAc                                    >  1:2154789/1‑100 (MQ=255)                                                   aCAAAACCTGAGCAAATTTTAGACACAATCCATTCAAGGTGTCAGCATGTATATTTCAAGCCTATTGATAAAGAAAAGTTTATAAATAGATTAGTTGAAc                                    >  1:201156/1‑100 (MQ=255)                                                   aCAAAACCTGAGCAAATTTTAGACACAATCCATTCAAGGTGTCAGCATGTATATTTCAAGCCTATTGATAAAGAAAAGTTTATAAATAGATTAGTTGAAc                                    >  1:1861352/1‑100 (MQ=255)                                                   aCAAAACCTGAGCAAATTTTAGACACAATCCATTCAAGGTGTCAGCATGTATATTTCAAGCCTATTGATAAAGAAAAGTTTATAAATAGATTAGTTGAAc                                    >  1:1691066/1‑100 (MQ=255)                                                   aCAAAACCTGAGCAAATTTTAGACACAATCCATTCAAGGTGTCAGCATGTATATTTCAAGCCTATTGATAAAGAAAAGTTTATAAATAGATTAGTTGAAc                                    >  1:1328990/1‑100 (MQ=255)                                                   aCAAAACCTGAGCAAATTTTAGACACAATCCATTCAAGGTGTCAGCATATATATTTCAAGCCTATTGATAAAGAAAAGTTTATAAATAGATTAGTTGAAc                                    >  1:5645551/1‑100 (MQ=255)                                                       aaCCTGAGCAAATTTTAGACACAATCCATTCAAGGTGTCAGCATGTATATTTCAAGCCTATTGATAAAGAAAAGTTTATAAATAGATTAGTTGAACaaaa                                >  1:4689853/1‑100 (MQ=255)                                                          cTGAGCAAATTTTAGACACAATCCATTCAAGGTGTCAGTATGTATATTTCAAGCCTATCGATAAAGAAAAGTTTATAAATAGATTAGTTGAACAAAACAt                             >  1:2596617/1‑100 (MQ=255)                                                          cTGAGCAAATTTTAGACACAATCCATTCAAGGTGTCAGCATGTATATTTCAAGCCTATTGATAAAGAAAAGTTTATAAATAGATTAGTTGAACAAAACAt                             >  1:1936713/1‑100 (MQ=255)                                                                      aaGACACAATCCATTCAAGGTGTCAGCATGTATATTTCAAGCCTATTGATAAAGAAAAGTTTATAAATAGATTAGTTGAACAAAACATGTCTAAGCCAGt                 >  1:2378511/2‑100 (MQ=255)                                                                               tCCATTCAAGGTGTCAGCATGTATATTTCAAGCCTATTGATAAAGAAAAGTTTATAAATAGATTAGTTGAACAAAACATGTCTAAGCCAGTAGCTGAAAt        >  1:1113448/1‑100 (MQ=255)                                                                                ccATTCAAGGTGTCAGCATGTATATTTCAAGCCTATTGATAAAGAAAAGTTTATAAATAGATTAGTTGAACAAAACATGTCTAAGCCAGTAGCTGAAATg       >  1:4213776/1‑100 (MQ=255)                                                                                    tCAAGGTGTCAGCATGTATATTTCAAGCCTATTGATAAAGAAAAGTTTATAAATAGATTAGTTGAACAAAACATGTCTAAGCCAGTAGCTGAAATGATTa   >  1:1381160/1‑100 (MQ=255)                                                                                     cAAGGTGTCAGCATGTATATTTCAAGCCTATTGATAAAGAAAAGTTTATAAATAGATTAGTTGAACAAAACATGTCTAAGCCAGTAGCTGAAATGATTAg  >  1:1017349/1‑100 (MQ=255)                                                                                         |                                                                                                TGAAATTTCTTGAAGAACCACCGGACAATACGATTGCTATTTTATTGTCTACAAAACCTGAGCAAATTTTAGACACAATCCATTCAAGATGTCAGCATGTATATTTCAAGCCTATTGATAAAGAAAAGTTTATAAATAGATTAGTTGAACAAAACATGTCTAAGCCAGTAGCTGAAATGATTAG  >  NC_007793/521662‑521845

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