Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	1,977,282	G→A	A140V (GCA→GTA)	SAUSA300_RS09815 ←	DNA repair exonuclease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	1,977,282	0	G	A	97.2%	95.8 / ‑5.0	36	A140V (GCA→GTA)	SAUSA300_RS09815	DNA repair exonuclease
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (8/27);  total (8/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.32e-01

AACTGAAGATTTACTATACGTACCATGCAAGACACCAATATGTATGCCTTTTTGGCCTTGACTTGATGGATATTCATCAATCTTGTTCTCATAACTTGCTCTATTTTCATAACTAAATCCGTGAATATAAATTGTTTCACCAGATTTAGTAATTGCTTCATACGTCTCAACTTTATTTGAAAATACAGAA  >  NC_007793/1977185‑1977374                                                                                                  |                                                                                             aaCTGAAGATTTACTATACGTACCATGCAAGACACCAATATGTATGCCTTTTTGGCCTTGACTAGATGGATATTCATCAATCTTGTTCTCATAACTTAct                                                                                            >  1:1898806/1‑100 (MQ=255)       aGATTTACTATACGTACCATGCAAGACACCAATATGTATGCCTTTTTGGCCTTGACTTGATGGATATTCATCAATCTTGTTCTCATAACTTACTCTAttg                                                                                      <  1:992616/100‑2 (MQ=255)            tACTATACGTACCATGCAAGACACCAATATGTATGCCTTTTTGGCCTTGACTTGATGGATATTCATCAATCTTGTTCTCATAACTTACTCTATTTTCATa                                                                                 >  1:2820954/1‑100 (MQ=255)            tACTATACGTACCATGCAAGACACCAATATGTATGCCTTTTTGGCCTTGACTTGATGGATATTCATCAATCTTGTTCTCATAACTTACTCTATTTTCATa                                                                                 <  1:4946143/100‑1 (MQ=255)                        cATGCAAGACACCAATATGTATGCCTTTTTGGCCTTGACTTGATGGATATTCATCAATCTTGTTCTCATAACTTACTCTATTTTCATAACTAAATCCGtt                                                                     <  1:5034305/100‑2 (MQ=255)                           gCAAGACACCAATATGTATGCCTTTTTGGCCTTGACTTGATGGATATTCATCAATCTTGTTCTCATAACTTACTCTATTTTCATAACTAAATCCGTGAat                                                                  >  1:996259/1‑100 (MQ=255)                                   gcAATATGTATGCCTTTTTGGCCTTGACTTGATGGATATTCATCAATCTTGTTCTCATAACTTACTCTATTTTCATAACTAAATCCGTGAATATAAATTg                                                          >  1:742488/2‑100 (MQ=255)                                        atgtatgCCTTTTTGGCCTTGACTTGATGGATATTCATCAATCTTGTTCTCATAACTTACTCTATTTTCATAACTAAATCCGTGAATATAAATTGTTTCa                                                     <  1:1241203/100‑1 (MQ=255)                                          gtatgCCTTTTTGGCCTTGACTTGATGGATATTCATCAATCTTGTTCTCATAACTTACTCTATTTTCATAACTAAATCCGTGAATATAAATTGTTTCAcc                                                   <  1:1786993/100‑1 (MQ=255)                                          gtatgCCTTTTTGGCCTTGACTTGATGGATATTCATCAATCTTGTTCTCATAACTTACTCTATTTTCATAACTAAATCCGTGAATATAAATTGTTTCAcc                                                   <  1:3334236/100‑1 (MQ=255)                                          gtatgCCTTTTTGGCCTTGACTTGATGGATATTCATCAATCTTGTTCTCATAACTTACTCTATTTTCATAACTAAATCCGTGAATATAAATTGTTTCAcc                                                   <  1:5632754/100‑1 (MQ=255)                                          gtatgCCTTTTTGGCCTTGACTTGATGGATATTCATCAATCTTGTTCTCATAACTTACTCTATTTTCATAACTAAATCCGTGAATATAAATTGTTTCAcc                                                   <  1:3576309/100‑1 (MQ=255)                                          gtatgCCTTTTTGGCCTTGACTTGATGGATATTCATCAATCTTGTTCTCATAACTTACTCTATTTTCATAACTAAATCCGTGAATATAAATTGTTTCAcc                                                   <  1:3672009/100‑1 (MQ=255)                                           tatgCCTTTTTGGCCTTGACTTGATGGATATTCATCAATCTTGTTCTCATAACTTACTCTATTTTCATAACTAAATCCGTGAATATAAATTGTTTCACCa                                                  <  1:655621/100‑1 (MQ=255)                                            ttccgatcttGGCCTTGACTTGATGGATATTCATCAATCTTGTTCTCATAACTTACTCTATTTTCATAACTAAATCCGTGAATATAAATTGTTTCACCAg                                                 <  1:130195/92‑1 (MQ=255)                                            atgCCTTTTTGGCCTTGACTTGATGGATATTCATCAATCTTGTTCTCATAACTTGCTCTATTTTCATAACTAAATCCGTGAATATAAATTGTTTCACCAg                                                 <  1:13918/100‑1 (MQ=255)                                            atgCCTTTTTGGCCTTGACTTGATGGATATTCATCAATCTTGTTCTCATAACTTACTCTATTTTCATAACTAAATCCGTGAATATAAATTGTTTCACCat                                                 <  1:5742698/100‑2 (MQ=255)                                                            gacttgaTGGATATTCATCAATCTTGTTCTCATAACTTACTCTATTTTCATAACTAAATCCGTGAATATAAATTGTTTCACCAGATTTAGTAATTGCtgc                                 <  1:4637693/100‑3 (MQ=255)                                                            gacttgaTGGATATTCATCAATCTTGTTCTCATAACTTACTCTATTTTCATAACTAAATCCGTGAATATAAATTGTTTCACCAGATTTAGTAATTGCTTc                                 >  1:3641597/1‑100 (MQ=255)                                                             acttgaTGGATATTCATCAATCTTGTTCTCATGACTTACTCTATTTTCATAACTAAATCCGTGAATATAAATTGTTTCACCAGATTTAGTAATTGCtaag                                <  1:4998792/100‑4 (MQ=255)                                                              cttgaTGGATATTCATCAATCTTGTTCTCATAACTTACTCTATTTTCATAACTAAATCCGTGAATATAAATTGTTTCACCAGATTTAGTAATTGCTTCAt                               <  1:605793/100‑1 (MQ=255)                                                              cctgaTGGATATTCATCAATCTTGTTCTCATAACTTACTCTATTTTCATAACTAAATCCGTGAATATAAATTGTTTCACCAGATTTAGTAATTGCTTCAt                               <  1:1803752/98‑1 (MQ=255)                                                                  aTGGATATTCATCAATCTTGTTCTCATAACTTACTCTATTTTCATAACTAATTCCGTGAATATAAATTGTTTCACCAGATTTAGTAATTGCTTCATACGt                           <  1:3840159/100‑1 (MQ=255)                                                                   tGGATATTCATCAATCTTGTTCTCATAACTTACTCTATTTTCATAACTAAATCCGTGAATATAAATTGTTTCACCAGATTTAGTAATTGCTTCATCCggg                          <  1:613656/100‑3 (MQ=255)                                                                     gATATTCATCAATCTTGTTCTCATAACTTACTCTATTTTCATAACTAAATCCGTGAGTATAAATTGTTTCACCAGATTTAGTAATTGCTTCATACGtctc                        <  1:3340731/100‑1 (MQ=255)                                                                     gATATTCATCAATCTTGTTCTCATAACTTACTCTATTTTCATAACTAAATCCGTGAATATAAATTGTTTCACCAGATTTAGTAATTGCTTCATACGtctc                        <  1:2263651/100‑1 (MQ=255)                                                                      atatTCATCAATCTTGTTCTCATAACTTACTCTATTTTCATAACTAAATCCGTGAATATAAATTGTTTCACCAGATTTAGTAATTGCTTCATACGTCTCa                       <  1:5033288/100‑1 (MQ=255)                                                                      atatTCATCAATCTTGTTCTCATAACTTACTCTATTTTCATAACTAAATCCGTGAATATAAATTGTTTCACCAGATTTAGTAATTGCTTCATACGTCTCa                       <  1:4811248/100‑1 (MQ=255)                                                                      atatTCATCAATCTTGTTCTCATAACTTACTCTATTTTCATAACTAAATCCGTGAATATAAATTGTTTCACCAGATTTAGTAATTGCTTCATACGTCTCa                       <  1:4600476/100‑1 (MQ=255)                                                                       gctcttccgaTCTTGTTCTCATAACTTACTCTATTTTCATAACTAAATCCGTGAATATAAATTGTTTCACCAGATTTAGTAATTGCTTCATACGTCTCaa                      <  1:909376/91‑1 (MQ=255)                                                                          tcatcaATCTTGTTCTCATAACTTACTCTATTTTCATAACTAAATCCGTGAATATAAATTGTTTCACCAGATTTAGTAATTGCTTCATACGTCTCAAcat                   <  1:4948317/100‑3 (MQ=255)                                                                           ttccgaTCTTGTTCTCATAACTTACTCTATTTTCATAACTAAATCCGTGAATATAAATTGTTTCACCAGATTTAGTAATTGCTTCATACGTCTCAACttt                  <  1:4229164/95‑1 (MQ=255)                                                                                    tGTTCTCATAACTTACTCTATTTTCATAACTAAATCCGTGAATATAAATTGTTTCACCAGATTTAGTAATTGCTTCATACGTCTCAACTTTATTTGaaaa         >  1:3425007/1‑100 (MQ=255)                                                                                          cATAACTTACTCTATTTTCATAACTAAATCCGTGAATATAAATTGTTTCACCAGATTTAGTAATTGCTTCATACGTCTCAACTTTATTTGAAAATACAGa   >  1:5682462/1‑100 (MQ=255)                                                                                          cATAACTTACTCTATTTTCATAACTAAATCCGTGAATATAAATTGTTTCACCAGATTTAGTAATTGCTTCATACGTCTCAACTTTATTTGAAAATACAGa   >  1:4195771/1‑100 (MQ=255)                                                                                           ttAACTTACTCTATTTTCATAACTAAATCCGTGAATATAAATTGTTTCACCAGATTTAGTAATTGCTTCATACGTCTCAACTTTATTTGAAAATACAGaa  <  1:2592201/99‑1 (MQ=255)                                                                                                  |                                                                                             AACTGAAGATTTACTATACGTACCATGCAAGACACCAATATGTATGCCTTTTTGGCCTTGACTTGATGGATATTCATCAATCTTGTTCTCATAACTTGCTCTATTTTCATAACTAAATCCGTGAATATAAATTGTTTCACCAGATTTAGTAATTGCTTCATACGTCTCAACTTTATTTGAAAATACAGAA  >  NC_007793/1977185‑1977374

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