Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	2,679,460	G→A	T30T (ACG→ACA)	SAUSA300_RS13770 →	peptidase M23

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	2,679,460	0	G	A	97.4%	97.5 / ‑5.6	38	T30T (ACG→ACA)	SAUSA300_RS13770	peptidase M23
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (31/6);  total (31/7)
Fisher's exact test for biased strand distribution p-value = 1.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

AGATAAAAATGGAATATAAAAAGATACTAATTCGTTTATTAATTGCTTTTGCAGTACTTTTCTCAGCAGATTTCACTTATCAATCCGTCGAACAAACGCATCAATCGCATGCCGCAGTTAATTATTATAGTAAAAACCAATGTACATGGTGGGCATTTAAACGTCGCGCACAAGTCGGTAAACCTGTTTCTAAT  >  NC_007793/2679363‑2679556                                                                                                  |                                                                                                 agaTAAAAATGGAATATAAAAAGATACTAATTCGTTTATTAATTGCTTTTGCAGTACTTTTCTCAGCAGATTTCACTTATCAATCCGTCGAACAAAcaca                                                                                                >  1:2268315/1‑100 (MQ=255)    tAAAAATGGAATATAAAAAGATACTAATGCGTTTATTAATTGCTTTTGCAGTACTTTTCTCAGCAGATTTCACTTATCAATCCGTCGAACAAACACATCa                                                                                             <  1:4460701/100‑1 (MQ=255)              atatAAAAAGATACTAATTCGTTTATTAATTGCTTTTGCAGTACTTTTCTCAGCAGATTTCACTTATCAATCCGTCGAACAAACACATCAATCGCATGcc                                                                                   >  1:3159550/1‑100 (MQ=255)                ataAAAAGATACTAATTCGTTTATTAATTGCTTTTGCAGTACTTTTCTCAGCAGATTTCACTTATCAATCCGTCGAACAAACACATCAATCGCATGCCGc                                                                                 <  1:2700985/100‑1 (MQ=255)                     aaGATACTAATTCGTTTATTAATTGCTTTTGCAGTACTTTTCTCAGCAGATTTCACTTATCAATCCGTCGAACAAACACATCAATCGCATGCCGCAgatc                                                                            >  1:6155492/1‑97 (MQ=255)                      aGATACTAATTCGTTTATTAATTGCTTTTGCAGTACTTTTCTCAGCAGATTTCACTTATCAATCCGTCGAACAAACACATCAATCGCACGcagcagatcg                                                                           >  1:1594044/1‑96 (MQ=255)                           cTAATTCGTTTATTAATTGCTTTTGCAGTACTTTTCTCAGCAGATTTCACTTATCAATCCGTCGAACAAACACATCAATCGCATGCCGCAGTTAattatt                                                                      >  1:6577715/1‑100 (MQ=255)                             aaTTCGTTTATTAATTGCTTTTGCAGTACTTTTCTCAGCAGATTTCACTTATCAATCCGTCGAACAAACACATCAATCGCATGCCGCAGTTAATtattat                                                                    >  1:108307/1‑100 (MQ=255)                             aaTTCGTTTATTAATTGCTTTTGCAGTACTTTTCTCAGCAGATTTCACTTATCAATCCGTCGAACAAACACATCAATCGCATGCCGCAGTTAATtattat                                                                    >  1:6368847/1‑100 (MQ=255)                             aaTTCGTTTATTAATTGCTTTTGCAGTACTTTTCTCAGCAGATTTCACTTATCAATCCGTCGAACAAACACATCAATCGCATGCCGCAGTTAATtattat                                                                    >  1:2500492/1‑100 (MQ=255)                                                tttGCAGTACTTTTCTCAGCAGATTTCACTTATCAATCCGTCGAACAAACACATCAATCGCATGCCGCAGTTAATTATTATAGTAAAAACCAATGTACAt                                                 >  1:1611340/1‑100 (MQ=255)                                                 ttGCAGTACTTTTCTCAGCAGATTTCACTTATCAATCCGTCGAACAAACACATCAATCGCATGCTGCAGTTAATTATTATAGTAAAAACCAATGTACAtg                                                >  1:542942/1‑100 (MQ=255)                                                 ttGCAGTACTTTTCTCAGCAGATTTCACTTATCAATCCGTCGAACAAACACATCAATCGCATGCCGCAGTTAATTATTATAGTAAAAACCAATGTACAtg                                                >  1:6318701/1‑100 (MQ=255)                                                 ctGCAGTACTTTTCTCAGCAGATTTCACTTATCAATCCGTCGAACAAACACATCAATCGCATGCCGCAGTTAATTATTATAGTAAAAACCAATGTACAtg                                                >  1:1200996/2‑100 (MQ=255)                                                 ctGCAGTACTTTTCTCAGCAGATTTCACTTATCAATCCGTCGAACAAACACATCAATCGCATGCCGCAGTTAATTATTATAGTAAAAACCAATGTACAtg                                                >  1:5167497/2‑100 (MQ=255)                                                  tGCAGTACTTTTCTCAGCAGATTTCACTTATCAATCCGTCTAACAAACACATCAATCGCATGCCGCAGTTAATTATTATAGTAAAAACCAATGTACAtgg                                               >  1:2992577/1‑100 (MQ=255)                                                  tGCAGTACTTTTCTCAGCAGATTTCACTTATCAATCCGTCGAACAAACACATCAATCGCATGCCGCAGTTAATTATTATAGTAAAAACCAATGTAGAtcg                                               >  1:1464258/1‑98 (MQ=255)                                                  tGCAGTACTTTTCTCAGCAGATTTCACTTATCAATCCGTCGAACAAACACATCAATCGCATGACGCAGTTAATTATTATAGTAAAAACCAATGTACAtgg                                               >  1:6256480/1‑100 (MQ=255)                                                   gCAGTACTTTTCTCAGCAGATTTCACTTATCAATCCGTCGAACAAACACATCAATCGCATGCCGCAGTTAATTATTATAGTAAAAACCAATGTACAtggt                                              >  1:3722424/1‑100 (MQ=255)                                                       tACTTTTCTCAGCAGATTTCACTTATCAATCCGTCGAACAAACACATCAATCGCATGCCGCAGTTAATTATTATAGTAAAAACCAATGTACATGGTgggt                                          <  1:4316614/100‑2 (MQ=255)                                                       cttccgatctagcagATTTCACTTATCAATCCGTCGAACAAACGCATCAATCGCATGCCGCAGTTAATTGTTATAGTAAAAACCAATGTACATGGTGGGc                                          <  1:6591231/90‑1 (MQ=255)                                                            ttCTCAGCAGATTTCACTTATCAATCCGTCGAACAAACACATCAATCGCATGCCGCAGTTAATTATTATAGTAAAAACCAATGTACATGGTGGGCATTTa                                     >  1:3115909/1‑100 (MQ=255)                                                                cagcagATTTCACTTATCAATCCGTCGAACAAACACATCAATCGCATGCCGCAGTTAATTATTATAGTAAAAACCAATGTACATGGTGGGCATTTAAACg                                 >  1:6356474/1‑100 (MQ=255)                                                                cagcagATTTCACTTATCAATCCGTCGAACAAACACATCAATCGCATGCCGCAGTTAATTATTATAGTAAAAACCAATGTACATGGTGGGCATTTAAACg                                 >  1:5650638/1‑100 (MQ=255)                                                                cagcagATTTCACTTATCAATCCGTCGAACAAACACATCAATCGCATGCCGCAGTTAATTATTATAGTAAAAACCAATGTACATGGTGGGCATTTAAACg                                 >  1:55340/1‑100 (MQ=255)                                                                cagcagATTTCACTTATCAATCCGTCGAACAAACACATCAATCGCATGCCGCAGTTAATTATTATAGTAAAAACCAATGTACATGGCGGGCATTTAAACg                                 >  1:954533/1‑100 (MQ=255)                                                                 ggcagATTTCACTTATCAATCCGTCGAACAAACACATCAATCGCATGCCGCAGTTAATTATTATAGTAAAAACCAATGTACATGGTGGGCATTTAAACGt                                >  1:5195233/2‑100 (MQ=255)                                                                        ttCACTTATCAATCCGTCGAACAAACACATCAATCGCATGCCGCAGTTAATTATTATAGTAAAAACCAATGTACATGGTGGGCATTTAAACGTCGCGcac                         <  1:4370144/100‑1 (MQ=255)                                                                        tccACTTATCAATCCGTCGAACAAACACATCAATCGCATGCCGCAGTTAATTATTATAGTAAAAACCAATGTACATGGTGGGCATTTAAACGTCGCGcac                         >  1:4132331/3‑100 (MQ=255)                                                                         tCACTTATCAATCCGTCGAACAAACACATCAATCGCATGCCGCAGTTAATTATTATAGTAAAAACCAATGTACATGGTGGGCATTTAAACGTCCCGcaca                        >  1:2398291/1‑100 (MQ=255)                                                                          cACTTATCAATCCGTCGAACAAACACATCAATCGCATGCCGCAGTTAATTATTATAGTAAAAACCAATGTACATGGTGGGCATTTAAACGTCGCGcacta                       >  1:196865/1‑98 (MQ=255)                                                                          cACTTATCAATCCGTCGAACAAACACATCAATCGCATGCCGCAGTTAATTATTATAGTAAAAACCAATGTACATGGTGGGCATTTAAACGTCGCGcacat                       >  1:4460712/1‑99 (MQ=255)                                                                            gttATCAATCCGTCGAACAAACACATCAATCGCATGCCGCAGTTAATTATTATAGTAAAAACCAATGTACATGGTGGGCATTTAAACGTCGCGCACAAGt                     >  1:3555456/2‑100 (MQ=255)                                                                              tATCAATCCGTCGAACAAACACATCAATCGCATGCCGCAGTTAATTATTATAGTAAAAACCAATGTACATGGTGGGCATTTAAACGTCGCGCACAAgatc                   >  1:5029075/1‑97 (MQ=255)                                                                              tATCAATCCGTCGAACAAACACATCAATCGCATGCCGCAATTAATTATTATAGTAAAAACCAATGTACATGGTGGGCATTTAAACGTCGCGCACaaagat                   >  1:4565729/1‑96 (MQ=255)                                                                                 cAATCCGTCGAACAAACACATCAATCGCATGCCGCAGTTAATTATTATAGTAAAAACCAATGTACATGGTGGGCATTTAAACGTCGCGCACAAGTCGGTa                <  1:6540737/100‑1 (MQ=255)                                                                                   aTCCGTCGAACAAACACATCAATCGCATGCCGCAGTTAATTATTATAGTAAAAACCAATGTACATGGTGGGCATTTAAACGTCGCGCACAGATCggaaga              >  1:3989690/1‑96 (MQ=255)                                                                                         agAACAAACACATCAATCGCATGCCGCAGTTAATTATTATAGTAAAAACCAATGTACATGGTGGGCATTTAAACGTCGCGCACAAGTCGGTAAACCTGtt        >  1:6009324/2‑100 (MQ=255)                                                                                               aaCACATCAATCGCATGCCGCAGTTAATTATTATAGTAAAAACCAATGTACATGGTGGGCATTTAAACGTCGCGCACAAGTCGGTAAGCCTGTTTCTAAt  <  1:4072581/100‑1 (MQ=255)                                                                                                  |                                                                                                 AGATAAAAATGGAATATAAAAAGATACTAATTCGTTTATTAATTGCTTTTGCAGTACTTTTCTCAGCAGATTTCACTTATCAATCCGTCGAACAAACGCATCAATCGCATGCCGCAGTTAATTATTATAGTAAAAACCAATGTACATGGTGGGCATTTAAACGTCGCGCACAAGTCGGTAAACCTGTTTCTAAT  >  NC_007793/2679363‑2679556

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