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	94.4%	87.4 / ‑3.3	36	T30T (ACG→ACA)	SAUSA300_RS13770	peptidase M23
Reads supporting (aligned to +/- strand):  ref base G (2/0);  new base A (29/5);  total (31/5)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.11e-01

GGAATATAAAAAGATACTAATTCGTTTATTAATTGCTTTTGCAGTACTTTTCTCAGCAGATTTCACTTATCAATCCGTCGAACAAACGCATCAATCGCATGCCGCAGTTAATTATTATAGTAAAAACCAATGTACATGGTGGGCATTTAAACGTCGCGCACAAGTCGGTAAACCTGTTTCTAAT  >  NC_007793/2679373‑2679556                                                                                        |                                                                                                 ggAATATAAAAAGATACTAATTCGTTTATTAATTGCTTTTGCAGTACTTTTCTCAGCAGATTTCACTTATCAATCCGTCGAACAAACACATCAATCGCAt                                                                                      >  1:2461127/1‑100 (MQ=255)           aaGATACTAATTCGTTTATTAATTGCTTTTGCAGTACTTTTCTCAGCAGATTTCACTTATCAATCCGACGAACAAACACATCAATCGCATGCCGCAGTTa                                                                            >  1:3821094/1‑100 (MQ=255)                aCTAATTCGTTTATTAATTGCTTTTGCAGTACTTTTCTCAGCAGATTTCACTTATCAATCCGTCGAACAAACACATCAATCGCATGCCGCAGTTAattat                                                                       >  1:1724222/1‑100 (MQ=255)                       agTTTATTAATTGCTTTTGCAGTACTTTTCTCAGCAGATTTCACTTATCAATCCGTCGAACAAACACATCAATCGCATGCCGCAGTTAATTATTATAGTa                                                                >  1:1303664/2‑100 (MQ=255)                                   gCTTTTGCAGTACTTTTCTCAGCAGATTTCACTTATCAATCCGTCGAACAAACACATCAATCGCATGCCGCAGTTAATTATTATAGTAAAAACCAATGTa                                                    <  1:3106417/100‑1 (MQ=255)                                      tttGCAGTACTTTTCTCAGCAGATTTCACTTATCAATCCGTCGACAAAACACATCAATCGCATGCCGCAGTTAATTATTATAGTAAAAACCAATGTACAt                                                 >  1:3926681/1‑100 (MQ=255)                                      cttGCAGTACTTTTCTCAGCAGATTTCACTTATCAATCCGTCGAACAAACACATCAATCGCATGCCGCAGTTAATTATTATAGTAAAAACCAATGTACAt                                                 >  1:2972978/2‑100 (MQ=255)                                       ttGCAGTACTTTTCTCAGCAGATTTCACTTATCAATCCGTCGAACAAACGCATCAATCGCATGCCGCAGTTAATTATTATAGTAAAAACCAATGTACAtg                                                >  1:2886841/1‑100 (MQ=255)                                       ttGCAGTACTTTTCTCAGCAGATTTCACTTATCAATCCGTCGAACAAACACATCAATCGCATGCCGCAGTTAATTATTATAGTAAAAACCAATGTACAtg                                                >  1:2149442/1‑100 (MQ=255)                                       ttGCAGTACTTTTCTCAGCAGATTTCACTTATCAATCCGTCGAACAAACACATCAATCGCATGCCGCAGTTAATTATTATAGTAAAAACCAATGTACAtg                                                >  1:2214652/1‑100 (MQ=255)                                       ttGCAGTACTTTTCTCAGCAGATTTCACTTATCAATCCGTCGAACAAACACATCAATCGCATGCCGCAGTTAATTATTATAGTAAAAACCAATGTACAtg                                                >  1:806509/1‑100 (MQ=255)                                        tGCAGCACTTTTCTCAGCAGATTTCACTTATCAATCCGTCGAACAAACACATCAATCGCATGCCGCAGTTAATTATTATAGTAAAAACCAATGTACAtgg                                               >  1:2869309/1‑100 (MQ=255)                                           aGTACTTTTCTCAGCAGATTTCACTTATCAATCCGTCGAACAAACACATCAATCGCATGCCGCAGTTAATTATTATAGTAAAAACCAATGTACAtggtgg                                            <  1:3098965/100‑1 (MQ=255)                                              aCTTTTCTCAGCAGATTTCACTTATCAATCCGTCGAACAAACACATCAATCGCATGCCGCAGTTAATTATTATAGTAAAAACCAATGTACATGGTGGGCa                                         >  1:4048265/1‑100 (MQ=255)                                                   tctcAGCAGATTTCACTTATCAATCCGTCGAACAAACACATCAATCGCATGCCGCAGTTAATTATTATAGTAAAAACCAATGTACATGGTGGGCATTTaa                                    >  1:147418/1‑100 (MQ=255)                                                    ctcAGCAGATTTCACTTATCAATCCGTCGAACAAACACATCAATCGCATGCCGCAGTTAATTATTATAGTAAAAACCAATGTACATGGTGGGCATTTaaa                                   >  1:1821950/1‑100 (MQ=255)                                                     tcAGCAGATTTCACTTATCAATCCGTCGAACAAACACATCAATCGCATGCCGCAGTTAATTATTATAGTAAAAACCAATGTACATGGTGGGCATTTAAAc                                  >  1:1250419/1‑100 (MQ=255)                                                      cagcagATTTCACTTATCAATCCGTCGAACAAACACATCAATCGCATGCCGCAGTTAATTATTATAGTAAAAACCAATGTACATGGTGGGCATTTAAACg                                 >  1:3046069/1‑100 (MQ=255)                                                      cagcagATTTCACTTATCAATCCGTCGAACAAACACATCAATCGCATGCCGCAGTTAATTATTATAGTAAAAACCAATGTACATGGTGGGCATTTAAACg                                 >  1:2498636/1‑100 (MQ=255)                                                      cagcagATTTCACTTATCAATCCGTCGAACAAACACATCAATCGCATGCCGCAGTTAATTATTATAGTAAAAACCAATGTACATGGTGGGCATTTAAACg                                 >  1:4548364/1‑100 (MQ=255)                                                      cagcagATTTCACTTATCAATCCGTCGAACAAACACATCAATCGCATGCCGCAGTTAATTATTATAGTAAAAACCAATGTACATGGTGGGCATTTAAACg                                 >  1:3602907/1‑100 (MQ=255)                                                      cagcagATTTCACTTATCAATCCGTCGAACAAACACATCAATCGCATGCCGCAGTTAATTATTATAGTAAAAACCAATGTACATGGTGGGCATTTAAACg                                 >  1:2052984/1‑100 (MQ=255)                                                       cgcagATTTCACTTATCAATCCGTCGAACAAACACATCAATCGCATGCCGCAGTTAATTATTATAGTAAAAACCAATGTACATGGTGGGCATTTAAACGt                                >  1:2345924/2‑100 (MQ=255)                                                               tCACTTATCAATCCGTCGAACAAACACGTCAATCGCATGCCGCAGTTAATTATTATAGTAAAAACCAATGTACATGGTGGGCATTTAAACGTCGTGcaca                        >  1:3611254/1‑100 (MQ=255)                                                                cACTTATCAATCCGTCGAACAAACACATCAATCGCATGCCGCAGTTAATTATTATAGTAAAAACCAATGTACATGGTGGGCATTTAAACGTCGCGCACaa                       >  1:4510844/1‑100 (MQ=255)                                                                 aCTTATCAATCCGTCGAACAAACACATCAATCGCATGCCGCAGTTAATTATTATAGTAAAAACCAATGTACATGGTGGGCATTTAAACGTCGCGCACAAg                      <  1:3868892/100‑1 (MQ=255)                                                                    tATCAATCCGTCGAACAAACACATCAATCGCATGCCGCAGTTAATTATTATAGTAAAAACCAATGTACATGGTGGGCATTTAAACGTCGCGcacagatcg                   >  1:4406654/1‑95 (MQ=255)                                                                    tATCAATCCGTCGAACAAACACATCAATCGCATGCCGCAGTTAATTATTATAGTAAAAACCAATGTACATGGTGGGCATTTAAACGTCGCGCACAAgatc                   >  1:1885215/1‑97 (MQ=255)                                                                    tATCAATCCGTCGAACAAACACATCAATCGCATGCCGCAGTTAATTATTATAGTAAAAACCAATGTACATGGTGGGCATTTAAACGTCGCGCACAAgatc                   >  1:4007266/1‑97 (MQ=255)                                                                    tATCAATCCGTCGAACAAACACATCAATCGCATGCCGCAGTTAATTATTATAGTAAAAACCAATGTACATGGTGGGCATTTAAACGTCGCGCACAAGTCg                   >  1:1634856/1‑100 (MQ=255)                                                                     aTCAATCCGTCGAACAAACGCATCAATCGCATGCCGCAGTTAATTATTATAGTAAAAACCAATGTACATGGTGGGCATTTAAACGTCGCGCACAAgatcg                  >  1:4228755/1‑96 (MQ=255)                                                                       ccaTCCGTCGAACAAACACATCAATCGCATGCCGCAGTTAATTATTATAGTAAAAACCAATGTACATGGTGGGCATTTAAACGTCGCGCACAAGTCGGTa                >  1:4829911/3‑100 (MQ=255)                                                                       cAATCCGTCGAACAAACACATCAATCGCATGCCGCAGTTAATTATTATAGTAAAAACCAATGTACATGGTGGGCATTTAAACGCCGCGCACAGATCggaa                >  1:3266114/1‑98 (MQ=255)                                                                                  aCAAACACATCAATCGCATGCCGCAGTTAATTATTATAGTAAAAACCAATGTACATGGTGGGCATTTAAACGTCGCGCACAAGTCGGTAAACCTGTTTCt     >  1:3585330/1‑100 (MQ=255)                                                                                     aaCACATCAATCGCATGCCGCAGTTAATTATTATAGTAAAAACCAATGTACATGGTGGGCATTTAAACGTCGCGCACAAGTCGGTAAACCTGTTTCTAAt  <  1:470866/100‑1 (MQ=255)                                                                                     aaCACATCAATCGCATGCCGCAGTTAATTATTATAGTAAAAACCAATGTACATGGTGGGCATTTAAACGTCGCGCACAAGTCGGTAAACCTGTTTCTAAt  <  1:5289261/100‑1 (MQ=255)                                                                                        |                                                                                                 GGAATATAAAAAGATACTAATTCGTTTATTAATTGCTTTTGCAGTACTTTTCTCAGCAGATTTCACTTATCAATCCGTCGAACAAACGCATCAATCGCATGCCGCAGTTAATTATTATAGTAAAAACCAATGTACATGGTGGGCATTTAAACGTCGCGCACAAGTCGGTAAACCTGTTTCTAAT  >  NC_007793/2679373‑2679556

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