Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	75,750	G→A	E62E (GAG→GAA)	SAUSA300_RS00340 →	universal stress protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	75,750	0	G	A	92.9%	71.8 / ‑0.6	28	coding (186/432 nt)	SAUSA300_RS00340	universal stress protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (19/7);  minor base . (1/1);  total (20/8)
Fisher's exact test for biased strand distribution p-value = 4.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.04e-01

TGTAATTTATCACGTTGTAACAGAAAATGATCTTGCAATTAGTGTCAAGTATTATAATAAGCATAAAGAAGATATTATTAGAGAAAAAGAGAAAAAACTCACTCCATTTTTACGTGAATTAGAAAAAAGAGATATTCAATATAAAATAGATGTAGATTTTGGGCATATTAAAGATAC  >  NC_007793/75660‑75836                                                                                           |                                                                                       tgtAATTTATCACGTTGTAACAGAAAATGATCTTGCAATTAGTGTCAAATATTATAATAAGCATAAAGAAGATATTATTAGAGAAAAAGAAAAAAAactc                                                                               >  1:1102487/1‑100 (MQ=255)                 gTAACAGAAAATGATCTTGCAATTAGTGTCAAATATTATAATAAGCATAAAGAAGATATTATTAGAGAAAAAGAAAAAAAACTCACTCCATTTTTACGTg                                                               >  1:4203641/1‑100 (MQ=255)                   aaCAGAAAATGATCTTGCAATTAGTGTCAAATATTATAATAAGCATAAAGAAGATATTATTAGAGAAAAAAAAAAAAAAATCAATCCCTTTTTaacggaa                                                             >  1:4676626/1‑94 (MQ=255)                              aTCTTGCAATTAGTGTCAAATATTATAATAAGCATAAAGAAGATATTATTAGAGAAAAAGAAAAAAAACTCACTCCATTTTTACGTGAATTAGAAAAAag                                                  <  1:2754868/100‑1 (MQ=255)                                  tGCAATTAGTGTCAAATATTATAATAAGCATAAAGAAGATATTATTAGAGAAAAAGAAAAAAAACTCACTCCATTTTTACGTGAATTATAAAAAAGAGat                                              >  1:3021950/1‑100 (MQ=255)                                  tGCAATTAGTGTCAAATATTATAATAAGCATAAAGAAGATATTATTAGAGAAAAAGAAAAAAAACTCACTCCATTTTTACGTGAATTAGAAAAAAGAGat                                              >  1:4707936/1‑100 (MQ=255)                                  tGCAATTAGTGTCAAATATTATAATAAGCATAAAGAAGATATTATTAGAGAAAAAGAAAAAAAACTCACTCCATTTTTACGTGAATTAGAAAAAAGAGat                                              >  1:4467156/1‑100 (MQ=255)                                  tGCAATTAGTGTCAAATATTATAATAAGCATAAAGAAGATATTATTAGAGAAAAAGAAAAAAAACTCACTCCATTTTTACGTGAATTAGAAAAAAGAGat                                              >  1:433408/1‑100 (MQ=255)                                  tGCAATTAGTGTCAAATATTATAATAAGCATAAAGAAGATATTATTAGAGAAAAAGAAAAAAAACTCACTCCATTTTTACGTGAATTAGAAAAAAGAGat                                              >  1:4244415/1‑100 (MQ=255)                                  tGCAATTAGTGTCAAATATTATAATAAGCATAAAGAAGATATTATTAGAGAAAAAGAAAAAAAACTCACTCCATTTTTACGTGAATTAGAAAAAAGAGat                                              >  1:4780665/1‑100 (MQ=255)                                  tGCAATTAGTGTCAAATATTATAATAAGCATAAAGAAGATATTATTAGAGAAAAAGAAAAAAAACTCACTCCATTTTTACGTGAATTAGAAAAAAGAGat                                              >  1:3761745/1‑100 (MQ=255)                                  tGCAATTAGTGTCAAATATTATAATAAGCATAAAGAAGATATTATTAGAGAAAAAGAAAAAAAACTCACTCCATTTTTACGTGAATTAGAAAAAAGAGat                                              >  1:3673030/1‑100 (MQ=255)                                  tGCAATTAGTGTCAAATATTATAATAAGCATAAAGAAGATATTATTAGAGAAAAAGAAAAAAAACTCACTCCATTTTTACGTGAATTAGAAAAAAGAGat                                              >  1:3251474/1‑100 (MQ=255)                                  tGCAATTAGTGTCAAATATTATAATAAGCATAAAGAAGATATTATTAGAGAAAAAGAAAAAAAACTCACTCCATTTTTACGTGAATTAGAAAAAAGAGat                                              >  1:596209/1‑100 (MQ=255)                                  tGCAATTAGTGTCAAATATTATAATAAGCATAAAGAAGATATTATTAGAGAAAAAGAAAAAAAACTCACTCCATTTTTACGTGAATTAGAAAAAAGAGat                                              >  1:2681230/1‑100 (MQ=255)                                  tGCAATTAGTGTCAAATATTATAATAAGCATAAAGAAGATATTATTAGAGAAAAAGA‑AAAAAACTCACTCCATTTTTACGTGAATTAGAGAAAAGAGatc                                             >  1:346909/1‑99 (MQ=255)                                   gCAATTAGTGTCAAATATTATAATAAGCATAAAGAAGATATTATTAGAGAAAAAGAAAAAAAACTCACTCCATTTTTACGTGAATTAGAAAAAAGAGata                                             >  1:3142098/1‑100 (MQ=255)                                                      ataataAGCATAAAGAAGATATTATTAGAGAAAAAGAAAAAAAACTCACTCCATTTTTACGTGAATTAGAAAAAAGAGATATTCAATATAAAATAGACGt                          <  1:2712756/100‑1 (MQ=255)                                                       taataaGCATAAAGAAGATATTATTAGAGAAAAAGAAAAAAAACTCACTCCATTTTTACGTGAATTAGAAAAAGGAGATATTCAATATAAAATAGATGTa                         <  1:5649972/100‑1 (MQ=255)                                                        aataaGCATAAAGAAGATATTATTAGAGAAAAAGAAAAAAAACTCACTCCATTTTTACGTGAATTAGAAAAAAGAGATATTCAATATAAAATGGATGTAg                        >  1:172848/1‑100 (MQ=255)                                                        aataaGCATAAAGAAGATATTATTAGAGAAAAAGAAAAAAAACTCACTCCATTTTTACGTGAATTAGAAAAAAGAGATATTCAATATAAAATAGacgaga                        >  1:4354795/1‑95 (MQ=255)                                                        aataaGCATAAAGAAGATATTATTAGAGAAAAAGAAAAAAAACTCACTCCATTTTTACGTGAATTAGAAAAAAGAGATATTCAATATAAAATAGATGAAg                        >  1:3825382/1‑100 (MQ=255)                                                           aaGCATAAAGAAGATATTATTAGAGAAAAAGAAAAAAAACTCACTCCATTTTTACGTGAATTAGAAAAAAGAGATATTCAATATTAAATAGAtcggaaga                     >  1:1554664/1‑93 (MQ=255)                                                                     aagaTATTATTAGAGAAAAAGAAAAAAAACTCACTCCATTTTTACGTGAATTAGAAAAAAGAGATATTCAATATAAAATAGATGTAGATTTTGGGCataa           <  1:5473829/100‑2 (MQ=255)                                                                      agaTATTATTAGAGAAAAAGAAAAAAAACTCACTCCATTTTTACGTGAATTAGAAAAAAGAGATATTCAATATAAAATAGATGTAGATTTTGGGCataag          <  1:2790010/100‑3 (MQ=255)                                                                       gaTATAATAAGAGAAAAAGA‑AAAAAACTCACTCCATTTTTACGTGAATTAGAAAAAAGAGATATTCAATATAAAATAGATGTAGATTTTGGGCATATTaa        <  1:692108/100‑1 (MQ=255)                                                                        atatTATTAGAGAAAAAGAAAAAAAACTCACTCCATTTTTACGTGAATTAGAAAAAAGAGATATTCAATATAAAATAGATGTAGATTTTGGGCATATTaa        <  1:4206479/100‑1 (MQ=255)                                                                              ttaGAGAAAAAGAAAAAAAACTCACTCCATTTTTACGTGAATTAGAAAAAAGAGATATTCAATATAAAATAGATGTAGATTTTGGGCATATTAAAGATAc  <  1:3772011/100‑1 (MQ=255)                                                                                           |                                                                                       TGTAATTTATCACGTTGTAACAGAAAATGATCTTGCAATTAGTGTCAAGTATTATAATAAGCATAAAGAAGATATTATTAGAGAAAAAGAGAAAAAACTCACTCCATTTTTACGTGAATTAGAAAAAAGAGATATTCAATATAAAATAGATGTAGATTTTGGGCATATTAAAGATAC  >  NC_007793/75660‑75836

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