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%	147.4 / NA	52	T41A (ACA→GCA)	SAUSA300_RS03070	phosphomevalonate kinase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (52/0);  total (52/0)

ATATTGCTGGAGAATATGCTGTAACAGAACCAGGATATAAATCTGTACTTATTGCGTTAGATCGTTTTGTAACTGCTACTATTGAAGAAACAGACCAATATAAAGGTACCATTCATTCAAAAGCATTACATCATAACCCAGTTACATTTAGTAGAGATGAAGATAGTATTGTCATTTCAGATCCACATG  >  NC_007793/648643‑648831                                                                                          |                                                                                                    ctatTGCTGGAGAATATGCTGTAACAGAACCAGGATATAAATCTGTACTTATTGCGTTAGATCGTTTTGTAACTGCTACTATTGAAGAAGCAGACCAata                                                                                           >  1:476313/2‑100 (MQ=255)  tatTGCTGGAGAATATGCTGTAACAGAACCAGGATATAAATCTGTACTTATTGCGTTAGATCGTTTTGTAACTGCTACTATTGAAGAAGCAGACCAatat                                                                                          >  1:4406179/1‑100 (MQ=255)   atTGCTGGAGAATATGCTGTAACAGAACCAGGATATAAATCTGTACTTATTGCGTTAGATCGTTTTGTAACTGCTACTATTGAAGAAGCAGACCAAtata                                                                                         >  1:4309090/1‑100 (MQ=255)    gtGCTGGAGAATATGCTGTAACAGAACCAGGATATAAATCTGTACTTATTGCGTTAGATCGTTTTGTAACTGCTACTATTGAAGAAGCAGACCAATATaa                                                                                        >  1:2472364/2‑100 (MQ=255)    ctGCTGGAGAATATGCTGTAACAGAACCAGGATATAAATCTGTACTTATTGCGTTAGATCGTTTTGTAACTGCTACTATTGAAGAAGCAGACCAATATaa                                                                                        >  1:3403996/2‑100 (MQ=255)     ggCTGGAGAATATGCTGTAACAGAACCAGGATATAAATCTGTACTTATTGCGTTAGATCGTTTTGTAACTGCTACTATTGAAGAAGCAGACCAATATaaa                                                                                       >  1:885390/2‑100 (MQ=255)      gCTGGAGAATATGCTGTAACAGAACCAGGATATAAATCTGTACTTATTGCGTTAGATCGTTTTGTAACTGCTACTATTGAAGAAGCAGACCAATATAAAg                                                                                      >  1:3959578/1‑100 (MQ=255)      gCTGGAGAATATGCTGTAACAGAACCAGGATATAAATCTGTACTTATTGCGTTAGATCGTTTTGTAACTGCTACTATTGAAGAAGCAGACCAATATAAAg                                                                                      >  1:3263859/1‑100 (MQ=255)      gCTGGAGAATATGCTGTAACAGAACCAGGATATAAATCTGTACTTATTGCGTTAGATCGTTTTGTAACTGCTACTATTGAAGAAGCAGACCAATATAAAg                                                                                      >  1:2537097/1‑100 (MQ=255)              ctatGCTGTAACAGAACCAGGATATAAATCTGTACTTATTGCGTTAGATCGTTTTGTAACTGCTACTATTGAAGAAGCAGACCAATATAAAGGTACcatt                                                                              >  1:674539/2‑100 (MQ=255)                atGCTGTAACAGAACCAGGATATGAATCTGTACTTATTGCGTTAGATCGTTTTGTAACTGCTACTATTGAAGAAGCAGACCAATATAAAGGTACcattca                                                                            >  1:1745118/1‑100 (MQ=255)                atGCTGTAACAGAACCAGGATATAAATCTGTACTTATTGCGTTAGATCGTTTTGTAACTGCTACTATTGAAGAAGCAGACCAATATAAAGGTACcattca                                                                            >  1:3061416/1‑100 (MQ=255)                 tGCTGTAACAGAACCAGGATATAAATCTGTACTTATTGCGTTAGATCGTTTTGTAACTGCTACTATTGAAGAAGCAGACCAATATAAAGGTACcattcat                                                                           >  1:2200234/1‑100 (MQ=255)                     gTAACAGAACCAGGATATAAATCTGTACTTATTGCGTTAGATCGTTTTGTAACTGCTACTATTGAAGAAGCAGACCAATATAAAGGTACCATTCATTCaa                                                                       >  1:2632748/1‑100 (MQ=255)                     gTAACAGAACCAGGATATAAATCTGTACTTATTGCGTTAGATCGTTTTGTAACTGCTACTATTGAAGAAGCAGACCAATATAAAGGTACCATTCATTCaa                                                                       >  1:4502555/1‑100 (MQ=255)                     gTAACAGAACCAGGATATAAATCTGTACTTATTGCGTTAGATCGTTTTGTAACTGCTACTATTGAAGAAGCAGACCAATATAAAGGTACCATTCATTCaa                                                                       >  1:3269453/1‑100 (MQ=255)                     gTAACAGAACCAGGATATAAATCTGTACTTATTGCGTTAGATCGTTTTGTAACTGCTACTATTGAAGAAGCAGACCAATATAAAGGTACCATTCATTCaa                                                                       >  1:1505458/1‑100 (MQ=255)                     gTAACAGAACCAGGATATAAATCTGTACTTATTGCGTTAGATCGTTTTGTAACTGCTACTATTGAAGAAGCAGACCAATATAAAGGTACCATTCATTCaa                                                                       >  1:1572560/1‑100 (MQ=255)                     gTAACAGAACCAGGATATAAATCTGTACTTATTGCGTTAGATCGCTTTGTAACTGCTACTATTGAAGAAGCAGACCAATATAAAGGTACCATTCATTCaa                                                                       >  1:2276529/1‑100 (MQ=255)                      aaaCAGAACCAGGATATAAATCTGTACTTATTGCGTTAGATCGTTTTGTAACTGCTACTATTGAAGAAGCAGACCAATATAAAGGTACCATTCATTCaaa                                                                      >  1:4962946/2‑100 (MQ=255)                       aaCAGAACCAGGATATAAATCTGTACTTATTGCGTTAGATCGTTTTGTAACTGCTACTATTGAAGAAGCAGACCAATATAAAGGTACCATTCATTCaaaa                                                                     >  1:422187/1‑100 (MQ=255)                          aGAACCAGGATATAAATCTGTACTTATTGCGTTAGATCGTTTTGTAACTGCTACTATTGAAGAAGCAGACCAATATAAAGGTACCATTCATTCAAAAGCa                                                                  >  1:5163192/1‑100 (MQ=255)                          aGAACCAGGATATAAATCTGTACTTATTGCGTTAGATCGTATTGTAACTGCTACTATTGAAGAAGCAGACCAATATAAAGGTACCATTCATTCAACAGCa                                                                  >  1:5429549/1‑100 (MQ=255)                           taaCCAGGATATAAATCTGTACTTATTGCGTTAGATCGTTTTGTAACTGCTACTATTGAAGAAGCAGACCAATATAAAGGTACCATTCATTCAAAAGCAt                                                                 >  1:1340515/2‑100 (MQ=255)                           gAACCAGGATATAAATCTGTACTTATTGCGTTAGATCGTTTTGTAACTGCTACTATTGAAGAAGCAGACCAATATAAAGGTACCATTCATTCAAAAGCAt                                                                 >  1:2736318/1‑100 (MQ=255)                           gAACCAGGATATAAATCTGTACTTATTGCGTTAGATCGTTTTGTAACTGCTACTATTGAAGAAGCAGACCAATATAAAGGTACCATTCATTCAAAAGCAt                                                                 >  1:2519075/1‑100 (MQ=255)                            aaCCAGGATATAAATCTGTACTTATTGCGTTAGATCGTTTTGTAACTGCTACTATTGAAGAAGCAGACCAATATAAAGGTACCATTCATTCAAAAGCAtt                                                                >  1:530357/1‑100 (MQ=255)                            aaCCAGGATATAAATCTGTACTTATTGCGTTAGATCGTTTTGTAACTGCTACTATTGAAGAAGCAGACCAATATAAAGGTACCATTCATTCAAAAGCAtt                                                                >  1:1057204/1‑100 (MQ=255)                            aaCCAGGATATAAATCTGTACTTATTGCGTTAGATCGTTTTGTAACTGCTACTATTGAAGAAGCAGACCAATATAAAGGTACCATTCATTCAAAAGCAtt                                                                >  1:1781742/1‑100 (MQ=255)                             aCCAGGATATAAATCTGTACTTATTGCGTTAGATCGTTTTGTAACTGCTACTATTGAAGAAGCAGACCAATATAAAGGTACCATTCACTCAAAAGCATTa                                                               >  1:351555/1‑100 (MQ=255)                             aCCAGGATATAAATCTGTACTTATTGCGTTAGATCGTTTTGTAACTGCTACTATCGAAGAAGCAGACCAATATAAAGGTACCATTCATTCAAAAGCATTa                                                               >  1:5416075/1‑100 (MQ=255)                              gcAGGATATAAATCTGTACTTATTGCGTTAGATCGTTTTGTAACTGCTACTATTGAAGAAGCAGACCAATATAAAGGTACCATTCATTCAAAAGCATTAc                                                              >  1:310566/2‑100 (MQ=255)                              ccAGGATATAAATCTGTACTTATTGCGTTAGATCGTTTTGTAACTGCTACTATTGAAGAAGCAGACCAATATAAAGGTACCATTCATTCAAAAGCATTAc                                                              >  1:2442933/1‑100 (MQ=255)                               cAGGATATAAATCTGTACTTATTGCGTTAGATCGTTTTGTAACTGCTACTATTGAAGAAGCAGACCAATATAAAGGTACCATTCATTCAAAAGCATTAca                                                             >  1:3360036/1‑100 (MQ=255)                                aGGATATAAATCTGTACTTATTGCGTTAGATCGTTTTGTAACTGCTACTATTGAAGAAGCAGACCAATATAAAGGTACCATTCATTCAAAAGCATTAcat                                                            >  1:4170023/1‑100 (MQ=255)                                aGGATATAAATCTGTACTTATTGCGTTAGATCGTTTTGTAACTGCTACTATTGAAGAAGCAGACCAATATAAAGGTACCATTCATTCAAAAGCATTAcat                                                            >  1:2689270/1‑100 (MQ=255)                                          tCTGTACTTATTGCGTTAGATCGTTTTGTAACTGCTACTATTGAAGAAGCAGACCAATATAAAGGTACCATTCATTCAAAAGCATTACATCATAACCCAg                                                  >  1:4180069/1‑100 (MQ=255)                                           cTGTACTTATTGCGTTAGATCGTTTTGTAACTGCTACTATTGAAGCAGCAGACCAATATAAAGGTACCATTCATTCAAAAGCATTACATCATAACCCAGt                                                 >  1:4926967/1‑100 (MQ=255)                                           cTGTACTTATTGCGTTAGATCGTTTTGTAACTGCTACTATTGAAGAAGCAGACCAATATAAAGGTACCATTCATTCAAAAGCATTACATCATAACCCAGt                                                 >  1:1111773/1‑100 (MQ=255)                                           cTGTACTTATTGCGTTAGATCGTTTTGTAACTGCTACTATTGAAGAAGCAGACCAATATAAAGGTACCATTCATTCAAAAGCATTACATCATAACCCAGt                                                 >  1:5320222/1‑100 (MQ=255)                                           cTGTACTTATTGCGTTAGATCGTTTTGTAACTGCTACTATTGAAGAAGCAGACCAATATAAAGGTACCATTCATTCAAAAGCATTACATCATAACCCAGt                                                 >  1:4676745/1‑100 (MQ=255)                                           cTGTACTTATTGCGTTAGATCGTTTTGTAACTGCTACTATTGAAGAAGCAGACCAATATAAAGGTACCATTCATTCAAAAGCATCACATCATAACCCAGt                                                 >  1:3996979/1‑100 (MQ=255)                                             gTACTTATTGCGTTAGATCGTTTTGTAACTGCTACTATTGAAGAAGCAGACCAATATAAAGGTACCATTCATTCAAAAGCATTACATCATAACCCAGTTa                                               >  1:3331305/1‑100 (MQ=255)                                             gTACTTATTGCGTTAGATCGTTTTGTAACTGCTACTATTGAAGAAGCAGACCAATATAAAGGTACCATTCATTCAAAAGCATTACATCATAACCCAGTTa                                               >  1:4649924/1‑100 (MQ=255)                                             gTACTTATTGCGTTAGATCGTTTTGTAACTGCTACTATTGAAGAAGCAGACCAATATAAAGGTACCATTCATTCAAAAGCATTACATCATAACCCAGTTa                                               >  1:2994074/1‑100 (MQ=255)                                             gTACTTATTGCGTTAGATCGTTTTGTAACTGCTACTATTGAAGAAGCAGACCAATATAAAGGTACCATTCATTCAAAAGCATTACATCATAACCCAGTTa                                               >  1:3552166/1‑100 (MQ=255)                                                cTTATTGCGTTAGATCGTTTTGTAACTGCTACTATTGAAGAAGCAGACCAATATAAAGGTACCATTCATTCAAAAGCATTACATCATAACCCAGTTACAt                                            >  1:5529565/1‑100 (MQ=255)                                                  tATTGCGTTAGATCGTTTTGTAACTGCTACTATTGAAGAAGCAGACCAATATAAAGGTACCATTCATTCAAAAGCATTACATCATAACCCAGTTACAttt                                          >  1:2301572/1‑100 (MQ=255)                                                  tATTGCGTTAGATCGTTTTGTAACTGCTACTATTGAAGAAGCAGACCAATATAAAGGTACCATTCATTCAAAAGCATTACATCATAACCCAGTTACAttt                                          >  1:5688148/1‑100 (MQ=255)                                                  tATTGCGTTAGATCGTTTTGTAACTGCTACTATTGAAGAAGCAGACCAATATAAAGGTACCATTCATTCAAAAGCATTACATCATAACCCAGTTACAttc                                          >  1:1485430/1‑99 (MQ=255)                                                                              actaTTGAAGAAGCAGACCAATATAAAGGTACCATTCATTCAAAAGCATTACATCATAACCCAGTTACATTTAGTAGAGATGAAGATAGTATTGTCAttt              >  1:207537/1‑100 (MQ=255)                                                                              actaTTGAAGAAGCAGACCAATATAAAGGTACCATTCATTCAAAAGCATTACATCATAACCCAGTTACATTTAGTAGAGATGAAGATAGTATTGTCAttt              >  1:4309313/1‑100 (MQ=255)                                                                                          aCAGACCAATATAAAGGTACCATTCATTCAAAAGCATTACATCATAACCCAGTTACATTTAGTAGAGATGAAGATAGTATTGTCATTTCAGATCCACATg  >  1:1530129/1‑100 (MQ=255)                                                                                          |                                                                                                    ATATTGCTGGAGAATATGCTGTAACAGAACCAGGATATAAATCTGTACTTATTGCGTTAGATCGTTTTGTAACTGCTACTATTGAAGAAACAGACCAATATAAAGGTACCATTCATTCAAAAGCATTACATCATAACCCAGTTACATTTAGTAGAGATGAAGATAGTATTGTCATTTCAGATCCACATG  >  NC_007793/648643‑648831

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