Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	425,833	A→T	E106D (GAA→GAT)	SAUSA300_RS02000 →	phosphoglycerate mutase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	425,833	0	A	T	100.0%	106.2 / NA	36	E106D (GAA→GAT)	SAUSA300_RS02000	phosphoglycerate mutase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (25/11);  total (25/11)

GCATCAACGACTACCCCTCTATTAAGAGAACGTTCACTTGGCGTATTTGAGGGTGAATATAAAGATGAAATCAGTGCGAATCCGAAATATGAAAAATATTTCAATGATCCAAACTTTAAAGACTTTCGTCATAGTTTTTCAC‑AAAAAGCGCCTGAAGGAGAAAGTTATGAAGATGTATATCAACGCGT  >  NC_007793/425741‑425928                                                                                             |                                                                                                 ccATCAACGACTACCCCTCTATTAAGAGAACGTTCACTTGGGGTATTTGAGGGTGAATATAAAGATGAAATCAGTGCGAATCCGAAATATGATAAATAtt                                                                                           >  1:5966118/2‑100 (MQ=255)   aTCAACGACTACCCCTCTATTAAGAGAACGTTCACTTGGCGTATTTGAGGGTGAATATAAAGATGAAATCAGTGCGAATCCGAAATATGATAAATATTTc                                                                                         >  1:797794/1‑100 (MQ=255)    tgcaggACTACCCCTCTATTAAGAGAACGTTCACTTTGCGTATTTGAGGGTGAATATAAAGATGAAATCAGTGCGAATCCGAAATATGATAAATATTTCa                                                                                        >  1:1312525/6‑100 (MQ=255)     taaCGACTACCCCTCTATTAAGAGAACGTTCACTTGGCGTATTTGAGGGTGAATATAAAGATGAAATCAGTGCGAATCCGAAATATGATAAATATTTCaa                                                                                       >  1:2909826/2‑100 (MQ=255)      aaCGACTACCCCTCTATTAAGAGAACGTTCACTTGGCGTATTTGAGGGTGAATATAAAGATGAAATCAGTGCGAATCCGAAATATGATAAATATTTCAAt                                                                                      >  1:2694099/1‑100 (MQ=255)          ccTACCCCTCTATTAAGAGAACGTTCACTTGGCGTATTTGAGGGTGAATATAAAGATGAAATCAGTGCGAATCCGAAATATGATAAATATTTCAATGATc                                                                                  >  1:983258/2‑100 (MQ=255)             cccccTCTATTAAGAGAACGTTCACTTGGCGTATTTGAGGGTGAATATAAAGATGAAATCAGTGCGAATCCGAAATATGATAAATATTTCAATGATCCaa                                                                               >  1:1401076/2‑100 (MQ=255)             aCCCCTCTATTAAGAGAACGTTCACTTGGCGTATTTGAGGGTGAATATAAAGATGAAATCAGTGCGAATCCGAAATATGATAAATATTTCAATGATCCaa                                                                               >  1:3859261/1‑100 (MQ=255)             aCCCCTCTATTAAGAGAACGTTCACTTGGCGTATTTGAGGGTGAATATAAAGATGAAATCAGTGCGAATCCGAAATATGATAAATATTTCAATGATCCaa                                                                               >  1:1584884/1‑100 (MQ=255)             aCCCCTCTATTAAGAGAACGTTCACTTGGCGTATTTGAGGGTGAATATAAAGATGAAATCAGTGCGAATCCGAAATATGATAAATATTTCAATGATCCaa                                                                               >  1:5700507/1‑100 (MQ=255)              ccccTCTATTAAGAGAACGTTCACTTGGCGTATTTGAGGGTGAATATAAAGACGAAATCAGTGCGAATCCGAAATATGATAAATATTTCAATGATCCaaa                                                                              >  1:4598018/1‑100 (MQ=255)                     aTTAAGAGAACGTTCACTTGGCGTATTTGAGGGTGAATATAAAGATGAAATCAGTGCGAATCCGAAATATGATAAATATTTCAATGATCCAAACTTTaaa                                                                       >  1:1574245/1‑100 (MQ=255)                     aTTAAGAGAACGTTCACTTGGCGTATTTGAGGGTGAATATAAAGATGAAATCAGTGCGAATCCGAAATATGATAAATATTTCAATGATCCAAACTTTaaa                                                                       >  1:2992392/1‑100 (MQ=255)                          gagaACGTTCACTTGGCGTATTTGAGGGTGAATATAAAGATGAAATCAGTGCGAATCCGAAATATGATAAATATTTCAATGATCCAAACTTTAAAGACtt                                                                  <  1:2445628/100‑1 (MQ=255)                                       tGGCGTATTTGAGGGTGAATATAAAGATGAAATCAGTGCGAATCCGAAATATGATAAATATTTCAATGATCCAAACTTTAAAGACTTTCGTCATAGtttt                                                     >  1:5481271/1‑100 (MQ=255)                                       tGGCGTATTTGAGGGTGAATATAAAGATGAAATCAGTGCGAATCCGAAATATGATAAATATTTCAATGATCCAAACTTTAAAGACTTTCGTCATAGtttt                                                     >  1:5398361/1‑100 (MQ=255)                                       tGGCGTATTTGAGGGTGAATATAAAGATGAAATCAGTGCGAATCCGAAATATGATAAATATTTCAATGATCCAAACTTTAAAGACTTTCGTCATAGtttt                                                     >  1:5937638/1‑100 (MQ=255)                                       tGGCGTATTTGAGGGTGAATATAAAGATGAAATCAGTGCGAATCCGAAATATGATAAATATTTCAATGATCCAAACTTTAAAGACTTTCGTCATAGtttt                                                     >  1:330575/1‑100 (MQ=255)                                       tGGCGTATTTGAGGGTGAATATAAAGATGAAATCAGTGCGAATCCGAAATATGATAAATATTTCAATGATCCAAACTTTAAAGACTTTCGTCATAGtttt                                                     >  1:363390/1‑100 (MQ=255)                                       tGGCGTATTTGAGGGTGAATATAAAGATGAAATCAGTGCGAATCCGAAATATGATAAATATTTCAATGATCCAAACTTTAAAGACTTTCGTCATAGtttt                                                     >  1:4498407/1‑100 (MQ=255)                                             aTTTGAGGGTGAATATAAAGATGAAATCAGTGCGAATCCGAAATATGATAAATATTTCAATGATCCAAACTTTAAAGACTTTCGTCATAGTTTTTCAC‑aa                                              >  1:1317142/1‑100 (MQ=255)                                               ttGAGGGTGAATATAAAGATGAAATCAGTGCGAATCCGAAATATGATAAATATTTCAATGATCCAAACTTTAAAGACTTTCGTCATAGTTTTTCAC‑aaaa                                            >  1:3485526/1‑100 (MQ=255)                                                  ggggTGAATATAAAGATGAAATCAGTGCGAATCCGAAATATGATAAATATTTCAATGATCCAAACTTTAAAGACTTTCGTCATAGTTTTTCAC‑AAAAAgc                                         >  1:2183178/2‑100 (MQ=255)                                                    ggTGAATATAAAGATGAAATCAGTGCGAATCCGAAATATGATAAATATTTCAATGATCCAAACTTTAAAGACTTTCGTCATAGTTTTTCAC‑AAAAAgcac                                       >  1:2030216/1‑98 (MQ=255)                                                       gAATATAAAGATGAAATCAGTGCGAATCCGAAATATGATAAATATTTCAATGATCCAAACTTTAAAGACTTTCGTCATAGTTTTTCAC‑AAAAAGCGCCTg                                    >  1:1144725/1‑100 (MQ=255)                                                         atatAAAGATGAAATCAGTGCGAATCCGAAATATGATAAATATTTCAATGATCCAAACTTTAAAGACTTTCGTCATAGTTTTTCAC‑AAAAAGCGCCTGaa                                  >  1:1817723/1‑100 (MQ=255)                                                           ataAAGATGAAATCAGTGCGAATCCGAAATATGATAAATATTTCAATGATCCAAACTTTAAAGACTTTCGTCATAGTTTTTCAC‑AAAAAGCGCCTGAAgg                                <  1:6058610/100‑1 (MQ=255)                                                            taAAGATGAAATCAGTGCGAATCCGAAATATGATAAATATTTCAATGATCCAAACTTTAAAGACTTTCGTCATAGTTTTTCACAAAAAAGCGCCTGAAgg                                <  1:6074161/100‑1 (MQ=255)                                                             aaaGATGGAATCAGTGCGAATCCGAAATATGATAAATATTTCAATGATCCAAACTTTAAAGACTTTCGTCATAGTTTTTCAC‑AAAAAGCGCCTGAAGgag                              <  1:3840313/100‑1 (MQ=255)                                                                       tCAGTGCGAATCCGAAATATGATAAATATTTCAATGATCCAAACTTTAAAGACTTTCGTCATAGTTTTTCAC‑AAAAAGCGCCTGAAGGAGAAAGTTATGa                    <  1:2941419/100‑1 (MQ=255)                                                                       tCAGTGCGAATCCGAAATATGATAAATATTTCAATGATCCAAACTTTAAAGACTTTCGTCATAGTTTTTCAC‑AAAAAGCGCCTGAAGGAGAAAGTTATGa                    <  1:2882733/100‑1 (MQ=255)                                                                       tCAGTGCGAATCCGAAATATGATAAATATTTCAATGATCCAAACTTTAAAGACTTTCGTCATAGTTTTTCAC‑AAAAAGCGCCTGAAGGAGAAAGTTATGa                    <  1:2582419/100‑1 (MQ=255)                                                                       tCAGTACGAATCCGAAATATGATAAATATTTCAATGATCCAAACTTTAAAGACTTTCGTCATAGTTTTTCAC‑AAAAAGCGCCTGAAGGAGAAAGTTATGa                    <  1:2268605/100‑1 (MQ=255)                                                                                         atGATAAATATTTCAATGATCCAAACTTTAAAGACTTTCGTCATAGTTTTTCAC‑AAAAAGCGCCTGAAGGAGAAAGTTATGAAGATGTATATCAACGCGt  <  1:3702061/100‑1 (MQ=255)                                                                                         atGATAAATATTTCAATGATCCAAACTTTAAAGACTTTCGTCATAGTTTTTCAC‑AAAAAGCGCCTGAAGGAGAAAGTTATGAAGATGTATATCAACGCGt  <  1:3467284/100‑1 (MQ=255)                                                                                         atGATAAATATTTCAATGATCCAAACTTTAAAGACTTTCGTCATAGTTTTTCAC‑AAAAAGCGCCTGAAGGAGAAAGTTATGAAGATGTATATCAACGCGt  <  1:23918/100‑1 (MQ=255)                                                                                             |                                                                                                 GCATCAACGACTACCCCTCTATTAAGAGAACGTTCACTTGGCGTATTTGAGGGTGAATATAAAGATGAAATCAGTGCGAATCCGAAATATGAAAAATATTTCAATGATCCAAACTTTAAAGACTTTCGTCATAGTTTTTCAC‑AAAAAGCGCCTGAAGGAGAAAGTTATGAAGATGTATATCAACGCGT  >  NC_007793/425741‑425928

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