Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	61,019	T→C	I64V (ATA→GTA)	SAUSA300_RS00265 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	61,019	0	T	C	82.6%	37.3 / 2.7	24	I64V (ATA→GTA)	SAUSA300_RS00265	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (0/4);  new base C (4/15);  total (4/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

AATGCTAAATATTAATAATTGATGATCCAAAGGTTTATCAAAGACTTCTTTAGAAAAAATAAATGCCAATCTACCTATAAAAAGACTAACCCAAGTCATAATATGATTAATGAATGCTTTTTTGCCATGTTCTGATGTTGGTTCATCAGCAAGTATCCATACAAATAAAGGTAAAA  >  NC_007793/60942‑61117                                                                              |                                                                                                   aataacaaaTATTAATAATTGATGATCAAACGGTTTATCAAAGACTCCTTTTGAAAAAATAAAAATTCCTCTGCCTACACAAAAACTAATCCATGTCata                                                                              <  1:1696896/94‑1 (MQ=255)                                                           tgctcttccgaTCTGCCTACACAAAAACTAATCCATGTCATAATATGATTCCCTAATGCTTTTTTACCATGTCCTGATGTTGGCTCATCAGCAAGTATcc                    <  1:1342785/90‑1 (MQ=255)                                                           tgctcttccgaTCTGCCTACACAAAAACTAATCCATGTCATAATATGATTCCCTAATGCTTTCTTACCATGTCCTGATGTTGGCTCATCAGCAAGTATcc                    <  1:1447598/90‑1 (MQ=255)                                                            tcttccgatctCTGCCTACACAAAAACTAATCCATGTCATAATATGATTCCCTAATGCTTTTTTACCATGTCCTGATGTTGGCTCATCAGCAAGTATCCa                   <  1:2455282/90‑1 (MQ=255)                                                            gctcatccgaTCTGCCTATACAAAAACTAATCCATGTCATAATATGATTCCCTAATGCTTTTTTACCATGTCCTGATGTTGGCTCATCAGCAAGTATCCa                   <  1:220221/91‑1 (MQ=255)                                                             ttccgatcttCTGCCTACACAAAAACTAATCCATGTCATAATATGATTCCCTAATGCTTTTTTACCATGTCCTGATGTTGGCTCATCAGCAAGTATCCAt                  <  1:786688/91‑1 (MQ=255)                                                              tcttccgaTCTGCCTATACAAAAACTAATCCATGTCATAATATGATTCCCTAATGCTTTTTTACCATGTCCTGATGTTGGCTCATCAGCAAGTATCCATa                 <  1:1142155/93‑1 (MQ=255)                                                              tccgatcttCTGCCTACACAAAAACTAATCCATGTCATAATATGATTCCCTAATGCTTTTTTACCATGTCCTGATGTTGGCTCATCAGCAAGTATCCAtg                 <  1:4467611/92‑2 (MQ=255)                                                               cttccgaTCTGCCTACACAAAAACTAATCCATGTCATAATATGATTCCCTAATGCTTTTTTACCATGTCCTGATGTTGGCTCATCAGCAAGTATCCATAc                <  1:3360591/94‑1 (MQ=255)                                                                ttccggtCTGCCTATACAAAAACTAATCCATGTCAGAATATGATTTCCTAAAGCTTTTTTACCATGTCCTGATGTTGGCTCATCAGCAAGTATCCATACa               <  1:4409056/94‑1 (MQ=255)                                                                ttccgaTCTGCCTAGACAAAAACTAATCCATGTCATAATATGATTCCCTAATGCTTTTTTACCATGTCCTGATGTTGGCTCATCAGCAAGTATCCATACa               <  1:4724217/95‑1 (MQ=255)                                                                ttccgaTCTGCCTACACAAAAACTAATCCATGTCATAATATGATTCCCTAATGCTTTTTTACCATGTCCTGATGTTGGCTCATCAGCAAGTATCCATACa               <  1:2485526/95‑1 (MQ=255)                                                                ttccgaTCTGCCTACACAAAAACTAATCCATGTCATAATATGATTCCCTAATGCTTTTTTACCATGTCCTGATGTTGGCTCATCAGCAAGTATCCATACa               <  1:5716041/95‑1 (MQ=255)                                                                ttccgaTCTGCCTACACAAAAACTAATCCATGTCATAATATGATTCCCTAATGCTTTTTTACCATGTCCTGATGTTGGCTCATCAGCAAGTATCCATACa               <  1:5224148/95‑1 (MQ=255)                                                                ccgatctCTGCCTACACAAAAACTAATCCATGTCATAATATGATTCCCTAATGCTTTTTTACCATGTCCTGATGTTGGCTCATCAGCAAGTATCCATACa               <  1:2334305/94‑1 (MQ=255)                                                                  gatctCTGCCTACACAAAAACTAATCCATGTCATAATATGATTCCCTAATGCTTTTTTACCATGTCCTGATGTTGGCTCATCAGCAAGTATCCATACaaa             <  1:2776804/96‑1 (MQ=255)                                                                     aTCTGCCTACACAAAAACTAATCCATGTCATAATATGATTCCCTAATGCTTTTTTACCATGTCCTGATGTTGGCTCATCAGCAAGTATCCATACAAATaa          <  1:3021366/100‑1 (MQ=255)                                                                         gccTATACAAAAACTAATCCATGTCATAATATGATTCCCTAATGCTTTTTTACCATGTCCTGATGTTGGCTCATCAGCAAGTATCCATACAAATAAAGGt      <  1:49580/99‑1 (MQ=255)                                                                         gccTACACAAAAACTAATCCATGTCATAATATGATTCCCTAATGCTTTTTTACCATGTCCTGATGTTGGCTCATCAGCAAGTATCCATACAAATAAAGGt      >  1:4033749/2‑100 (MQ=255)                                                                         gccTACACAAAAACTAATCCATGTCATAATATGATTCCCTAATGCTTTTTTACCATGTCCTGATGTTGGCTCATCAGCAAGTATCCATACAAATAAAGGt      >  1:3785813/2‑100 (MQ=255)                                                                         gccTACACAAAAACTAATCCATGTCATAATATGATTCCCTAATGCTTTTTTACCATGTCCTGATGTTGGCTCATCAGCAAGTATCCATACAAATAAAGGt      >  1:3628319/2‑100 (MQ=255)                                                                         gccTACACAAAAACTAATCCATGTCATAATATGATTCCCTAATGCTTTTTTACCATGTCCTGATGTTGGCTCATCAGCAAGTATCCATACAAATAAAGGt      <  1:4712518/99‑1 (MQ=255)                                                                         gccTACACAAAAACTAATCCATGTCATAATATGATTCCCTAATGCTTTTTTACCATGTCCTGATGTTGGCTCATCAGCAAGTATCCATACAAATAAAGGt      >  1:5478884/2‑100 (MQ=255)                                                                         gccTACACAAAAACTAATCCATGTCATAATATGAGTCCCTAATGCTTTTTTACCATGTCCTGATGTTGGGTCATCAGCAAGTATCCATACAAATAAAGGt      <  1:3935916/99‑1 (MQ=255)                                                                             ataCAAAAACTAATCCATGTCATAATATGATTCCCTAAGGCGTTTTTACCATGTCCTGATGTTGGCTCATCAGCAAGTATCCATACAAATAAAGGTagca  <  1:2640796/100‑4 (MQ=255)                                                                              |                                                                                                   AATGCTAAATATTAATAATTGATGATCCAAAGGTTTATCAAAGACTTCTTTAGAAAAAATAAATGCCAATCTACCTATAAAAAGACTAACCCAAGTCATAATATGATTAATGAATGCTTTTTTGCCATGTTCTGATGTTGGTTCATCAGCAAGTATCCATACAAATAAAGGTAAAA  >  NC_007793/60942‑61117

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