Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	648,732	A→G	T41A (ACA→GCA)	SAUSA300_RS03070 →	phosphomevalonate kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	648,732	0	A	G	100.0%	133.1 / NA	50	T41A (ACA→GCA)	SAUSA300_RS03070	phosphomevalonate kinase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (49/1);  total (49/1)

ATATTGCTGGAGAATATGCTGTAACAGAACCAGGATATAAATCTGTACTTATTGCGTTAGATCGTTTTGTAACTGCTACTATTGAAGAAACAGACCAATATAAAGGTACCATTCATTCAAAAGCATTACATCATAACCCAGTTACATTTAGTAGAGATGAAGATAGTATTGTCATTTCAGATCCA  >  NC_007793/648643‑648827                                                                                          |                                                                                                ggcttGCTGGAGAATATGCTGTAACAGAACCAGGATATAAATCTGTACTTATTGCGTTAGATCGTTTTGTAACTGCTACTATTGAAGAAGCAGACCAata                                                                                       >  1:1544592/4‑100 (MQ=255)  gcctGCTGGAGAATATGCTGTAACAGAACCAGGATATAAATCTGTACTTATTGCGTTAGATCGTTTTGTAACTGCTACTATTGAAGAAGCAGACCAatat                                                                                      >  1:3626930/4‑100 (MQ=255)  catTGCTGGAGAATATGCTGTAACAGAACCAGGATATAAATCTGTACTTATTGCGTTAGATCGTTTTGTAACTGCTACTATTGAAGAAGCAGACCAatat                                                                                      >  1:3397787/2‑100 (MQ=255)   atTGCTGGAGAATATGCTGTAACAGAACCAGGATATAAATCTGTACTTATTGCGTTAGATCGTTTTGTAACTGCTACTATTGAAGAAGCAGACCAAtata                                                                                     >  1:6308580/1‑100 (MQ=255)     tGCTGGAGAATATGCTGTAACAGAACCAGGATATACATCTGTACTTATTGCGTTAGGTCGTTTTGTAACTGCTACTATTGAGGAAGCAGatcggaagagc                                                                                   >  1:6271957/1‑90 (MQ=255)     tGCTGGAGAATATGCTGTAACAGAACCAGGATATAAATCTGTACTTATTGCGTTAGATCGTTTTGTAACTGCTACTATTGAAGAAGCAGACCAATATaaa                                                                                   >  1:1823546/1‑100 (MQ=255)      gCTGGAGAATATGCTGTAACAGAACCAGGATATAAATCTGTACTTATTGCGTTAGATCGTTTTGTAACTGCTACTATTGAAGAAGCAGACCAATATAAAg                                                                                  >  1:6089863/1‑100 (MQ=255)      gCTGGAGAATATGCTGTAACAGAACCAGGATATAAATCTGTACTTATTGCGTTAGATCGTTTTGTAACTGCTACTATTGAAGAAGCAGACCAATATAAAg                                                                                  >  1:4037607/1‑100 (MQ=255)           agaATATGCTGTAACAGAACCAGGATATAAATCTGTACTTATTGCGTTAGATCGTTTTGTAACTGCTACTATTGAAGAAGCAGACCAATATAAAGGTAcc                                                                             >  1:2093464/1‑100 (MQ=255)            gaATATGCTGTAACAGAACCAGGATATAAATCTGTACTTATTGCGTTAGATCGTTTTGTAACTGCTACTATTGAAGAAGCAGACCAATATAAAGGCACca                                                                            >  1:4286467/1‑100 (MQ=255)               tatGCTGTAACAGAACCAGGATATAAATCTGTACTTATTGCGTTAGATCGTTTTGTAACTGCTACTATTGAAGAAGCAGACCAATATAAAGGTACcattc                                                                         >  1:305557/1‑100 (MQ=255)                atGCTGTAACAGAACCAGGATATAAATCTGTACTTATTGCGTTAGATCGTTTTGTAACTGCTACTATTGAAGAAGCAGACCAATATAAAGGTACcattca                                                                        >  1:1429427/1‑100 (MQ=255)                acgCTGTAACAGAACCAGGATATAAATCTGTACTTATTGCGTTAGATCGTTTTGTAACTGCTACTATTGAAGAAGCAGACCAATATAAAGGTACcattca                                                                        >  1:3676800/3‑100 (MQ=255)                 tGCTGTAACAGAACCAGGATATAAATCTGTACTTATTGCGTTAGATCGTTTTGTAACTGCTACTATTGAAGAAGCAGACCAATATAAAGGTACcattcat                                                                       >  1:653333/1‑100 (MQ=255)                     gTAACAGAACCAGGATATAAATCTGTACTTATTGCGTTAGATCGTTTTGTAACTGCTACTATTGAAGAAGCAGACCAATATAAAGGTACCATTCATTCaa                                                                   >  1:5523968/1‑100 (MQ=255)                      tAACAGAACCAGGATATAAATCTGTACTTATTGCGTTAGATCGTTTTGTAACTGCTACTATTGAAGAAGCAGACCAATATAAAGGTACCATTCATTCaaa                                                                  >  1:6258935/1‑100 (MQ=255)                       aaCAGAACCAGGATATAAATCTGTACTTATTGCGTTAGATCGTTTTGTAACTGCTACTATTGAAGAAGCAGACCAATATAAAGGTACCATTCATTCaaaa                                                                 >  1:4746486/1‑100 (MQ=255)                        aCAGAACCAGGATATAAATCTGTACTTATTGCGTTAGATCGTTTTGTAACTGCTACTATTGAAGAAGCAGACCAATATAAAGGTACCATTCATTCAAAAg                                                                >  1:6533289/1‑100 (MQ=255)                         cAGAACCAGGATATAAATCTGTACTTATTGCGTTAGATCGTTTTGTAACTGCTACTATTGAAGAAGCAGACCAATATAAAGGTACCATTCATTCAAAAGc                                                               >  1:2814555/1‑100 (MQ=255)                          aGAACCAGGATATAAATCTGTACTTATTGCGTTAGATCGTTTTGTAACTGCTACTATTGAAGAAGCAGACCAATATAAAGGTACCATTCATTCAAAAGCa                                                              >  1:3128682/1‑100 (MQ=255)                          aGAACCAGAATATAAATCTGTACTTATTGCGTTAGATCGTTTTGTAACTGCTACTATTGAAGAAGCAGACCAATATAAAGGCACCATTCATTCAAAAGCa                                                              >  1:1960663/1‑100 (MQ=255)                           gAACCAGGATATAAATCTGTACTTATTGCGTTAGATCGTTTTGTAACTGCTACTATTGAAGAAGCAGACCAATATAAAGGTACCATTCATTCAAAAGCAt                                                             >  1:3083622/1‑100 (MQ=255)                             cccAGGATATAAATCTGTACTTATTGCGTTAGATCGTTTTGTAACTGCTACTATTGAAGAAGCAGACCAATATAAAGGTACCATTCATTCAAAAGCATTa                                                           >  1:5260946/2‑100 (MQ=255)                              ccAGGATATAAATCTGTACTTATTGCGTTAGATCGTTTTGTAACTGCTACTATTGAAGAAGCAGACCAATATAAAGGTACCATTCATTCAAAAGCATTAc                                                          >  1:4239427/1‑100 (MQ=255)                              ccAGGATATAAATCTGTACTTATTGCGTTAGATCGTTTTGTAACTGCTACTATTGAAGAAGCAGACCAATATAAAGGTACCATTCATTCAAAAGCATTAc                                                          >  1:2226280/1‑100 (MQ=255)                                     ataAATCTGTACTTATTGCGTTAGATCGTTTTGTAACTGCTACTATTGAAGAAGCAGACCAATATAAAGGTACCATTCATTCAAAAGCATTACATCATaa                                                   >  1:3393934/1‑100 (MQ=255)                                       aaaTCTGTACTTATTGCGTTAGATCGTTTTGTAACTGCTGCTATTGAAGAAGCAGCCCAATATAAAGGTACCATTCATTCCAAAGCATTACATCATAAcc                                                 >  1:2195801/1‑100 (MQ=255)                                           cTGTACTTATTGCGTTAGATCGTTTTGTAACTGCTACTATTGAAGAAGCAGACCAATATAAAGGTACCATTCATTCAAAAGCATTACATCATAACCCAGt                                             >  1:6480046/1‑100 (MQ=255)                                           cTGTACTTATTGCGTTAGATCGTTTTGTAACTGCTACTATTGAAGAAGCAGACCAATATAAAGGTACCATTCATTCAAAAGCATTACATCATAACCCAGt                                             >  1:6048723/1‑100 (MQ=255)                                           cTGTACTTATTGCGTTAGATCGTTTTGTAACTGCTACTATTGAAGAAGCAGACCAATATAAAGGTACCATTCATTCAAAAGCATTACATCACAACCCAga                                             >  1:2441761/1‑99 (MQ=255)                                             gTACTTATTGCGTTAGATCGTTTTGTAACTGCTACTATTGAAGAAGCAGACCAATATAAAGGTACCATTCCTTCAAAAGCATTACATCATAACCCAGTTa                                           >  1:6322383/1‑100 (MQ=255)                                             gTACTTATTGCGTTAGATCGTTTTGTAACTGCTACTATTGAAGAAGCAGACCAATATAAAGGTACCATTCATTCAAAAGCATTACGTCATAACCCAGTTa                                           >  1:4938651/1‑100 (MQ=255)                                             gTACTTATTGCGTTAGATCGTTTTGTAACTGCTACTATTGAAGAAGCAGACCAATATAAAGGTACCATTCATTCAAAAGCATTACATCATAACCCAGTTa                                           >  1:893796/1‑100 (MQ=255)                                             gTACTTATTGCGTTAGATCGTTTTGTAACTGCTACTATTGAAGAAGCAGACCAATATAAAGGTACCATTCATTCAAAAGCATTACATCATAACCCAGTTa                                           >  1:5999610/1‑100 (MQ=255)                                             gTACTTATTGCGTTAGATCGTTTTGTAACTGCTACTATTGAAGAAGCAGACCAATATAAAGGTACCATTCATTCAAAAGCATTACATCATAACCCAGTTa                                           >  1:5047153/1‑100 (MQ=255)                                             gTACTTATTGCGTTAGATCGTTTTGTAACTGCTACTATTGAAGAAGCAGACCAATATAAAGGTACCATTCATTCAAAAGCATTACATCATAACCCAGTTa                                           >  1:4083840/1‑100 (MQ=255)                                             gTACTTATTGCGTTAGATCGTTTTGTAACTGCTACTATTGAAGAAGCAGACCAATATAAAGGTACCATTCATTCAAAAGCATTACATCATAACCCAGTTa                                           >  1:3906210/1‑100 (MQ=255)                                             gTACTTATTGCGTTAGATCGTTTTGTAACTGCTACTATTGAAGAAGCAGACCAATATAAAGGTACCATTCATTCAAAAGCATTACATCATAACCCAGTTa                                           >  1:3072924/1‑100 (MQ=255)                                             gTACTTATTGCGTTAGATCGTTTTGTAACTGCTACTATTGAAGAAGCAGACCAATATAAAGGTACCATTCATTCAAAAGCATTACATCATAACCCAGTTa                                           >  1:3064786/1‑100 (MQ=255)                                             gTACTTATTGCGTTAGATCGTTTTGTAACTGCTACTATTGAAGAAGCAGACCAATATAAAGGTACCATTCATTCAAAAGCATTACATCATAACCCAGTTa                                           >  1:2529325/1‑100 (MQ=255)                                               aCTTATTGCGTTAGATCGTTTTGTAACTGCTACTATTGAAGAAGCAGACCAATATAAAGGTACCATTCATTCAAAAGCATTACATCATAACCCAGTTACa                                         <  1:1868161/100‑1 (MQ=255)                                                  tATTGCGTTAGATCGTTTTGTAACTGCTACTATTGAAGAAGCAGACCAATATAAAGGTACCATTCATTCAAAAGCATTACATCATAACCCAGTTACAttt                                      >  1:4543092/1‑100 (MQ=255)                                                  tATTGCGTTAGATCGTTTTGTAACTGCTACTATTGAAGAAGCAGACCAATATAAAGGTACCATTCATTCAAAAGCATTACATCATAACCCAGTTACAttt                                      >  1:2158702/1‑100 (MQ=255)                                                   aTTGCGTTAGATCGTTTTGTAACTGCTACTATTGAAGAAGCAGACCAATATAAAGGTACCATTCATTCAAAAGCATTACATCATAACCCAGTTACAttga                                     >  1:655922/1‑98 (MQ=255)                                                   aTTGCGTTAGATCGTTTTGTAACTGCTACTATTGAAGAAGCAGACCAATATAAAGGTACCATTCATTCAAAAGCATTACATCATAACCCAGTTACATTta                                     >  1:100231/1‑100 (MQ=255)                                                         ctAGATCGTTTTGTAACTGCTACTATTGAAGAAGCAGACCAATATAAAGGTACCATTCATTCAAAAGCATTACATCATAACCCAGTTACATTTAGTagag                               >  1:2380991/2‑100 (MQ=255)                                                           aGATCGTTTTGTAACTGCTACTATTGAAGAAGCAGACCAATATAAAGGTACCATTCATTCAAAAGCATTACATCATAACCCAGTTGCATTTAGTAGAGAt                             >  1:4790347/1‑100 (MQ=255)                                                            gATCGTTTTGTAACTGCTACTATTGGAGAAGCAGACCAATATAAAGGTACCATTCATTCAAAAGCATTACATCATAACCCAGTTACATTTAGTAGAGATg                            >  1:5898798/1‑100 (MQ=255)                                                                           cctactaTTGAAGAAGCAGACCAATATAAAGGTACCATTCATTCAAAAGCATTACATCATAACCCAGTTACATTTAGTAGAGATGAAGATAGTATTGTCa             >  1:4322410/2‑100 (MQ=255)                                                                                      agaaGCAGACCAATATAAAGGTACCATTCATTCAAAAGCATTACATCATAACCCAGTTACATTTAGTAGAGATGAAGATAGTATTGTCATTTCAGATcgg  >  1:745923/1‑98 (MQ=255)                                                                                          |                                                                                                ATATTGCTGGAGAATATGCTGTAACAGAACCAGGATATAAATCTGTACTTATTGCGTTAGATCGTTTTGTAACTGCTACTATTGAAGAAACAGACCAATATAAAGGTACCATTCATTCAAAAGCATTACATCATAACCCAGTTACATTTAGTAGAGATGAAGATAGTATTGTCATTTCAGATCCA  >  NC_007793/648643‑648827

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