Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	1,331,771	A→G	T45A (ACA→GCA)	SAUSA300_RS06580 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	1,331,771	0	A	G	91.8%	113.4 / 1.1	50	T45A (ACA→GCA)	SAUSA300_RS06580	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (3/1);  new base G (42/3);  total (46/4)
Fisher's exact test for biased strand distribution p-value = 2.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.07e-01

TTAACGCTTACTAATAATAAAATATTAATTGGAAAAGTGATAGACTTTGATGATAAAGTTGATAATTTTGATGGTTATAATTCGATAGAAATTGATACAGGTAGAATTACATATGATATATCAGAAAATAAAATTAAAACTATTGTTTTACAATAATAGCATCCTTTCTACACAGATAAAGAGAAAGCGGTGC  >  NC_007793/1331675‑1331867                                                                                                 |                                                                                                 ttAACGCTTACTAATAATAAAATATTAATTGGAAAAGTGATAGACTTTGATGATAAAGTTGATAATTTTGATGGTTATAATTCGATAGAAATTGATGCAg                                                                                               >  1:1276665/1‑100 (MQ=255) atAACGCTTACTAATAATAAAATATTAATTGGAAAAGTGATAGACTTTGATGATAAAGTTGATAATTTTGATGGTTATAATTCGATAGAAATTGATGCAg                                                                                               >  1:1973019/2‑100 (MQ=255)  tAACGCTTACTAATAATAAAATATTAATTGGAAAAGTGATAGACTTTGATGATAAAGTTGATAATTTTGATGGTTATAATTCGATAGAAATTGATGCAgg                                                                                              >  1:4382370/1‑100 (MQ=255)  tAACGCTTACTAATAATAAAATATTAATTGGAAAAGTGATAGACTTTGATGATAAAGTTGATAATTTTGATGGTTATAATTCGATAGAAATTGATGCAgg                                                                                              >  1:1950671/1‑100 (MQ=255)  tAACGCTTACTAATAATAAAATATTAATTGGAAAAGCGATAGACTTTGATGATAAAGTTGATAATTTTGATGGTTATAATTCGATAGAAATCGATGCAga                                                                                              >  1:4646912/1‑99 (MQ=255)  gaaCGCTTACTAATAATAAAATATTAATTGGAAAAGTGATAGACTTTGATGATAAAGTTGATAATTTTGATGGTTATAATTCGATAGAAATTGATGCAgg                                                                                              >  1:3928372/2‑100 (MQ=255)  caaCGCTTACTAATAATAAAATATTAATTGGAAAAGTGATAGACTTTGATGATAAAGTTGATAATTTTGATGGTTATAATTCGATAGAAATTGATGCaag                                                                                              >  1:5508936/2‑98 (MQ=255)  caaCGCTTACTAATAATAAAATATTAATTGGAAAAGTGATAGACTTTGATGATAAAGTTGATAATTTTGATGGTTATAATTCGATAGAAATTGATGCAga                                                                                              >  1:2547659/2‑99 (MQ=255)  aaaCGCTTACTAATAATAAAATATTAATTGGAAAAGTGATAGACTTTGATGATAAAGTTGATAATTTTGATGGTTATAATTCGATAGAAATTGATGCAgg                                                                                              >  1:4702963/2‑100 (MQ=255)   cgatcTTACTAATAATAAAATATTAATTGGAAAAGTGATAGACTTTGATGATAAAGTTGATAATTTTGATGGTTATAATTCGATAGAAATTGATGCAGGt                                                                                             <  1:1941696/96‑1 (MQ=255)   caCGCTTACTAATAATAAAATATTAATTGGAAAAGTGATAGACTTTGATGATAAAGTTGATAATTTTGATGGTTATAATTCGATAGAAATTGATACAGGt                                                                                             <  1:2147025/99‑1 (MQ=255)   accGCTTACTAATAATAAAATATTAATTGGAAAAGTGATAGACTTTGATGATAAAGTTGATAATTTTGATGGTTATAATTCGATAGAAATTGATGCAGGt                                                                                             >  1:4326892/3‑100 (MQ=255)   accGCTTACTAATAATAAAATATTAATTGGAAAAGTGATAGACTTTGATGATAAAGTTGATAATTTTGATGGTTATAATTCGATAGAAATTGATGCAGGt                                                                                             >  1:3124930/3‑100 (MQ=255)   aaCGCTTACTAATAATAAAATGTTAATTGGAAAAGTGATAGACTTTGATGATAAAGTTGATAATTTTGATGGTTATAATTCGATAGAGATTGATGCAGGt                                                                                             >  1:4685268/1‑100 (MQ=255)   aaCGCTTACTAATAATAAAATATTAATTGGAAAAGTGATAGACTTTGATGATAAAGTTGATAATTTTGATGGTTATAATTCGATAGAAATTGATGCAggc                                                                                             >  1:3855914/1‑99 (MQ=255)   aaCGCTTACTAATAATAAAATATTAATTGGAAAAGTGATAGACTTTGATGATAAAGTTGATAATTTTGATGGTTATAATTCGATAGAAATTGATGCAgat                                                                                             >  1:698404/1‑98 (MQ=255)   aaCGCTTACTAATAATAAAATATTAATTGGAAAAGTGATAGACTTTGATGATAAAGTTGATAATTTTGATGGTTATAATTCGATAGAAATTGATGCAGGt                                                                                             >  1:4414819/1‑100 (MQ=255)   aaCGCTTACTAATAATAAAATATTAATTGGAAAAGTGATAGACTTTGATGATAAAGTTGATAATTTTGATGGTTATAATTCGATAGAAATTGATGCAGGt                                                                                             >  1:3862625/1‑100 (MQ=255)   aaCGCTTACTAATAATAAAATATTAATTGGAAAAGTGATAGACTTTGATGATAAAGTTGATAATTTTGATGGTTATAATTCGATAGAAATTGATGCAGGt                                                                                             >  1:1433946/1‑100 (MQ=255)   aaCGCTTACTAATAATAAAATATTAATTGGAAAAGTGATAGACTTTGATGATAAAGTTGATAATTTTGATGGTTATAATTCGATAGAAATTGATGCAGGt                                                                                             <  1:1597494/100‑1 (MQ=255)   aaCGCTTACTAATAATAAAATATTAATTGGAAAAGTGATAGACTTTGATGATAAAGTTGATAATTTTGATGGTTATAATTCGATAGAAATTGATGCAGGt                                                                                             >  1:4822773/1‑100 (MQ=255)   aaCGCTTACTAATAATAAAATATTAATTGGAAAAGTGATAGACTTTGATGATAAAGTTGATAATTTTGATGGTTATAATTCGATAGAAATTGATACAGGt                                                                                             >  1:3692225/1‑100 (MQ=255)   aaCGCTTACTAATAATAAAATATTAATTGGAAAAGTGATAGACTTTGATGATAAAGTTGATAATTTTGATGGTTATAATTCGATAGAAATTGATACAGGt                                                                                             >  1:383276/1‑100 (MQ=255)   aaCGCTTACTAATAATAAAATATTAATTGGAAAAGTGATAGACTTTGATGATAAAGTTGATAATTTTGATGGTTATAATTCGATAGAAATCGATGCAgat                                                                                             >  1:3182567/1‑98 (MQ=255)   aaCGCTTACTAATAATAAAATATTAATTGGAAAAGTGATAGACTTTGATGATAAAGTTGATAATTTTGATGGTTATAATTCGATAGAAATAGATGCAgat                                                                                             >  1:792179/1‑98 (MQ=255)    aCGCTTACTAATAATAAAATATTAATTGGAAAAGTGATAGACTTTGATGATAAAGTTGATAATTTTGATGGTTATAATTCGATAGAAATTGATGCAGGTa                                                                                            >  1:2979039/1‑100 (MQ=255)     cGCTTACTAATAATAAAATATTAATTGGAAAAGTGATAGACTTTGATGATAAAGTTGATAATTTTGATGGTTATAATTCGATAGAAATTGATGCAGGTAg                                                                                           >  1:300449/1‑100 (MQ=255)     accaggcTAATAATAAAATATTAATTGGAAAAGTGATAGACTTTGATGATAAAGTTGATAATTTTGATGGTTATAATTCGATAGAAATTGATACAGGTAg                                                                                           >  1:416974/7‑100 (MQ=255)      gCTTACTACTAATAAAATATTAATTGGAAAAGTGATAGACTTTGATGATAAAGTTGATAATTTTGATGGTTATAATTCGATAGAAATTGATGCAGGTAGa                                                                                          >  1:2552275/1‑100 (MQ=255)      gCTTACTAATAATAAAATATTAATTGGAAAAGTGATAGACTTTGATGATAAAGTTGATAATTTTGATGGTTATAATTCGATAGAAATTGATGCAgatcgg                                                                                          >  1:789795/1‑95 (MQ=255)         tACTAATAATAACATATTAATTGGAAAAGTGATAGACTTTGATGATAAAGTTGATACTTTTGATGGTTATAATTCGATAGCAATTGATGCAGGTAGatcg                                                                                       >  1:1032261/1‑97 (MQ=255)         tACTAATAATAAAATATTAATTGGAAAAGTGATAGACTTTGATGATAAAGTTGATAATTTTGATGGTTATAATTCGATAGAAATTGATGCAGGTAGatcg                                                                                       >  1:3865654/1‑97 (MQ=255)         tACTAATAATAAAATATTAATTGGAAAAGTGATAGACTTTGATGATAAAGTTGATAATTTTGATGGTTATAATTCGATAGAAATTGATGCAGGTAGAAtt                                                                                       >  1:1199092/1‑100 (MQ=255)         tACTAATAATAAAATATTAATTGGAAAAGTGATAGACTTTGATGATAAAGTTGATAATTTTGATGGTTATAATTCGATAGAAATTGATCCAGGCAGatcg                                                                                       >  1:4843207/1‑97 (MQ=255)         tACTAATAATAAAATATTAATTGGAAAAGGGATAGACTTTGATGATAAAGTTGATAATTTTGATGGTTATAATTCGATAGAAATTTATGCAGGTAGatcg                                                                                       >  1:1790623/1‑97 (MQ=255)          ccTAATAATAAAATATTAATTGGAAAAGTGATAGACTTTGATGATAAAGTTGATAATTTTGATGGTTATAATTCGATAGAAATTGATGCAGGTAGAATTa                                                                                      >  1:657825/2‑100 (MQ=255)          aCTAATAATAAAATATTAATTGGAAAAGTGATAGACTTTGATGGTAAAGTTGATAATTTTGATGGTTATAATTCGATAGAAATAGATGCAgatcggaaga                                                                                      >  1:3218120/1‑91 (MQ=255)                    aaaTATTAATTGGAAAAGTGATAGACTTTGATGATAAAGTTGATAATTTTGATGGTTATAATTCGATAGAAATTGATGCAGGTAGAATTACATATGatat                                                                            >  1:1471264/1‑100 (MQ=255)                                     gTGATAGACTTTGATGATAAAGTTGATAATTTTGATGGTTATAATTCGATAGAAATTGATGCAGGTAGAATTACATATGATATATCAGAAAATAAAATTa                                                           >  1:2646701/1‑100 (MQ=255)                                                      tAAAGTTGATAATTTTGATGGTTATAATTCGATAGAAATTGATGCAGGTAGAATTACATATGATATATCAGAAAATAAAATTAAAACTATTGTTTTACaa                                          >  1:898314/1‑100 (MQ=255)                                                             gATAATTTTGATGGTTATAATTCGATAGAAATTGATGCAGGAAGAATTACATATGATATATCAGAAAATAAAATTAAAACTATTGTTTTACAATAATAGc                                   >  1:2647471/1‑100 (MQ=255)                                                                      gctGGTTATAATTCGATAGAAATTGATGCAGGTAGAATTACATATGATATATCAGAAAATAAAATTAAAACTATTGTTTTACAATAATAGCATCCTTTCt                          >  1:2275363/3‑100 (MQ=255)                                                                      gATGGTTATAATTCGATAGAAATTGATGCAGGTAGAATTACATATGATATATCAGAAAATAAAATTAAAACTATTGTTTTACAATAATAGCATCCTTTCt                          >  1:3077972/1‑100 (MQ=255)                                                                      cctGGTTATAATTCGATAGAAATTGATGCAGGTAGAATTACATATGATATATCAGAAAATAAAATTAAAACTATTGTTTTACAATAATAGCATCCTTTCt                          >  1:4802813/3‑100 (MQ=255)                                                                         ggTTATAATTCGATAGAAATTGATGCAGGTAGAATTACATATGATATATCAGAAAATAAAATTAAAACTATTGTTTTACAATAATAGCATCCTTTCTaca                       >  1:490382/1‑100 (MQ=255)                                                                             ataATTCGATAGAAATTGATGCAGGTAGAATTACATATGATATATCAGAAAATAAAATTAAAACTATTGTTTTACAATAATAGCATCCTTTCTagggcta                   <  1:2639848/100‑7 (MQ=255)                                                                                             tGATGCAGGTAGAATTACATATGATATATCAGAAAATAAAATTAAAACTATTGTTTTACAATAATAGCATCCTTTCTACACAGATAAAGAGAAAGCggcg   >  1:3339907/1‑98 (MQ=255)                                                                                             tGATGCAGGTAGAATTACATATGATATATCAGAAAATAAAATTAAAACTATTGTTTTACAATAATAGCATCCTTTCTACACAGATAAAGAGAAAGCGGTg   >  1:5011254/1‑100 (MQ=255)                                                                                             tGATGCAGGTAGAATTACATATGATATATCAGAAAATAAAATTAAAACTATTGTTTTACAATAATAGCATCCTTTCTACACAGATAAAGAGAAAGCGGTg   >  1:2793363/1‑100 (MQ=255)                                                                                              gATGCAGGTAGAATTACATATGATATATCAGAAAATAAAATTAAAACTATTGTTTTACAATAATAGCATCCTTTCTACACAGATAAAGAGAAAGCGGTGc  >  1:4855968/1‑100 (MQ=255)                                                                                                 |                                                                                                 TTAACGCTTACTAATAATAAAATATTAATTGGAAAAGTGATAGACTTTGATGATAAAGTTGATAATTTTGATGGTTATAATTCGATAGAAATTGATACAGGTAGAATTACATATGATATATCAGAAAATAAAATTAAAACTATTGTTTTACAATAATAGCATCCTTTCTACACAGATAAAGAGAAAGCGGTGC  >  NC_007793/1331675‑1331867

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