Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	726,738	A→G	R174R (CGT→CGC)	SAUSA300_RS03490 ←	peptidase M23

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	726,738	0	A	G	97.4%	94.5 / ‑4.6	39	R174R (CGT→CGC)	SAUSA300_RS03490	peptidase M23
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (10/28);  total (11/28)
Fisher's exact test for biased strand distribution p-value = 2.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.23e-01

ATGGTGTGTTGTTTACTTGGTAACCATCATTAGCTGCGTTACCAGCCCAATACTTAGCGTCTGACCAATATGTGCTAATTGGACTACCAGCTTGAGCACGACGGTCAAATACGTACCATGTACATTGACCAGCAGTGTATAAATTTTGGTGATTAAAAGATGATGCATTGCCATTGCTACCTGTTGTAGC  >  NC_007793/726641‑726830                                                                                                  |                                                                                             aTGGCGTGTTGTTTACTTGGTAACCATCATTAGCTGCGTTACCAGCCCAATACTTAGCGTCTGACCAATATGTGCTAATTGGACTACCAGCTTGAGcgcg                                                                                            <  1:1260829/100‑1 (MQ=255)       tgttgtTTACTTGGTAACCATCATTAGCTGCGTTACCAGCCCAATACTTAGCGTCTGACCAATATGTGCTAATTGGACTACCAGCTTGAGCGCGACGGtg                                                                                      <  1:4846837/100‑2 (MQ=255)        gttgttTACTTGGTAACCATCATTAGCTGCGTTACCAGCCCAATACTTAGCGTCTGACCAATATGTGCTAATTGGACTACCAGCTTGAGCGCGACGGTCa                                                                                     <  1:3662320/100‑1 (MQ=255)                       cccATCATTAGCTGCGTTACCAGCCCAATACTTAGCGTCTGACCAATATGTGCTAATTGGACTACCAGCTTGAGCGCGACGGTCAAATACGTACCATGTa                                                                      >  1:4776295/2‑100 (MQ=255)                               tAGCTGCGTTACCAGCCCAATACTTAGCGTCTGACCAATATGTGCTAATTGGACTACCAGCTTGAGCGCGACGGTCAAATACGTACCATGTACATTGAcc                                                              >  1:5342056/1‑100 (MQ=255)                                        tACCAGCCCAATACTTAGCGTCTGACCAATATGTGCTAATTGGACTACCAGCTTGAGCGCGACGGTCAAATACGTACCATGTACATTGACCAGCAGTGta                                                     <  1:3366183/100‑1 (MQ=255)                                         aCCAGCCCAATACTTAGCGTCTGACCAATATGTGCTAATTGGACTACCAGCTTGAGCGCGACGGTCAAATACGTACCATGTACATTGACCAGCAGTGtat                                                    >  1:948251/1‑100 (MQ=255)                                          ccAGCCCAATACTTAGCGTCTGACCAATATGTGCTAATTGGACTACCAGCTCGAGCGCGACGGTCAAATACGTACCATGTACATTGACCAGCAGTGtata                                                   <  1:3777547/100‑1 (MQ=255)                                                       tAGCGTCTGACCAATATGTGCTAATTGGACTACCAGCTTGAGCGCGACGGTCAAATACGTACCATGTACATTGACCAGCAGTGTATAAATTTTGGTggca                                      <  1:191325/100‑4 (MQ=255)                                                               gACCAATATGTGCTAATTGGACTACCAGCTTGAGCGCGACGGTCAAATACGTACCATGTACATTGACCAGCAGTGTATAAATTTTGGTGATTAAAAgatg                              >  1:1619708/1‑100 (MQ=255)                                                               gACCAATATGTGCTAATTGGACTACCAGCTTGAGCGCGACGGTCAAATACGTACCATGTACATTGACCAGCAGTGTATAAATTTTGGTGATTAAAAgatg                              <  1:3280634/100‑1 (MQ=255)                                                                             aaTTGGACTACCAGCTTGAGCGCGACGGTCAAATACGTACCATGTACATTGACCAGCAGTGTATAAATTTTGGTGATTAAAAGATGATGcattgccattg                >  1:394857/1‑100 (MQ=255)                                                                                 gatcTACCAGCTTGAGCGCGACGGTCAAATACGTACCATGTACATTGACCAGCAGTGTATAAATTTTGGTGATTAAAAGATGATGCATTGCCATTGCTAc            <  1:2823798/97‑1 (MQ=255)                                                                                  gACTACCAGCTTGAGCACGACGGTCAAATACGTACCCTGTACATTGACCAGCAGTGTATAAATTTTGGTGATTAAAAGATGATGCATTGCCATTGCTAcc           >  1:5542681/1‑100 (MQ=255)                                                                                  atcTACCAGCTTGAGCGCGACGGTCAAATACGTACCATGTACATTGACCAGCAGTGTATAAATTTTGGTGATTAAAAGATGATGCATTGCCATTGCTAcc           <  1:2702103/98‑1 (MQ=255)                                                                                     tACCAGTTTGAGCGCGACGGTCAAATACGTACCATGTACATTGACCAGCAGTGTATAAATTTTGGTGATTAAAAGATGATGCATTGCCATTGCTACCtga        <  1:446552/100‑2 (MQ=255)                                                                                      aCCTGCTTGAGCGCGACGGTCAAATACGTACCATGTACATTGACCAGCAGTGTATAAATTTTGGTGATTAAAAGATGATGCATTGCCATTGCTACCtgtt       <  1:1192520/100‑1 (MQ=255)                                                                                      aCCAGGTTGAGCGCGACGGTCAAATACGTACCATGTACATTGACCAGCAGTGTATAAATTTTGGTGATTAAAAGATGATGCATTGCCATTGCTACCtgtt       <  1:1710283/100‑1 (MQ=255)                                                                                      aCCAGCTTGAGCGCGACGGTCAAATACGTACCATGTACATTGACCAGCAGTGTATAAATTTTGGTGATTAAAAGATGATGCATTGCCATTGCTACCtgtt       >  1:4128499/1‑100 (MQ=255)                                                                                      aCCAGCTTGAGCGCGACGGTCAAATACGTACCATGTACATTGACCAGCAGTGTATAAATTTTGGTGATTAAAAGATGATGCATTGCCATTGCTACCtgtt       >  1:5244797/1‑100 (MQ=255)                                                                                      aCCAGCTTGAGCGCGACGGTCAAATACGTACCATGTACATTGACCAGCAGTGTATAAATTTTGGTGATTAAAAGATGATGCATTGCCATTGCTACCtgtt       <  1:469966/100‑1 (MQ=255)                                                                                      aCCAGCTTGAGCGCGACGGTCAAATACGTACCATGTACATTGACCAGCAGTGTATAAATTTTGGTGATTAAAAGATGATGCATTGCCATTGCTACCtgtt       >  1:4824030/1‑100 (MQ=255)                                                                                      aCCAGCTTGAGCGCGACGGTCAAATACGTACCATGTACATTGACCAGCAGTGTATAAATTTTGGTGATTAAAAGATGATGCATTGCCATTGCTACCtgtt       >  1:2016233/1‑100 (MQ=255)                                                                                      aCCAGCTTGAGCGCGACGGTCAAATACGTACCATGTACATTGACCAGCAGTGTATAAATTTTGGTGATTAAAAGATGATGCATTGCCATTGCTACCtgtc       >  1:2738674/1‑99 (MQ=255)                                                                                       tcAGCTTGAGCGCGACGGTCAAATACGTACCATGTACATTGACCAGCAGTGTATAAATTTTGGTGATTAAAAGATGATGCATTGCCATTGCTACCtgttg      <  1:4666469/99‑1 (MQ=255)                                                                                       ccAGCTTGAGCGCGACGGTCAAATACGTACCATGTACATTGACCAGCAGTGTATAAATTTTGGTGATTAAAAGATGATGCATTGCCATTGCTACCtgttg      <  1:4324789/100‑1 (MQ=255)                                                                                        taGCTTGAGCGCGACGGTCAAATACGTACCATGTACATTGACCAGCAGTGTATAAATTTTGGTGATTAAAAGATGATGCATTGCCATTGCTACCtgttgt     <  1:630548/99‑1 (MQ=255)                                                                                        cggtttGAGCGCGACGGTCAAATACGTACCATGTACATTGACCAGCAGTGTATAAATTTTGGTGATTAAAAGATGATGCATTGCCATTGCTACCtgttgt     <  1:375978/96‑1 (MQ=255)                                                                                        cAGTTTGAGCGCGACGGTCAAATACGTACCATGTACATTGACCAGCAGTGTATAAATTTTGGTGATTAAAAGATGATGCATTGCCATTGCTACCtgttgt     <  1:2504749/100‑1 (MQ=255)                                                                                        cAGGTTGAGCGCGACGGTCAAATACGTACCATGTACATTGACCAGCAGTGTATAAATTTTGGTGATTAAAAGATGATGCATTGCCATTGCTGCCtgttgt     <  1:3351338/100‑1 (MQ=255)                                                                                        cAGCTTGAGCGCGACGGTCAAATACGTACCATGTACATTGACCAGCAGTGTATAAATTTTGGTGATTAAAAGATGATGCATTGCCATTGCTACCtgttgt     <  1:1117158/100‑1 (MQ=255)                                                                                        cAGCTTGAGCGCGACGGTCAAATACGTACCATGTACATTGACCAGCAGTGTATAAATTTTGGTGATTAAAAGATGATGCATTGCCATTGCTACCtgttgt     <  1:3301966/100‑1 (MQ=255)                                                                                        cAGCTTGAGCGCGACGGTCAAATACGTACCATGTACATTGACCAGCAGTGTATAAATTTTGGTGATTAAAAGATGATGCATTGCCATTGCTACCtgttgt     <  1:2017642/100‑1 (MQ=255)                                                                                        cAGCTTGAGCGCGACGGTCAAATACGTACCATGTACATTGACCAGCAGTGTATAAATTTTGGTGATTAAAAGATGATGCATTGCCATTGCTACCtgttgt     <  1:88028/100‑1 (MQ=255)                                                                                        cAGCTTGAGCGCGACGGTCAAATACGTACCATGTACATTGACCAGCAGTGTATAAATTTTGGTGATTAAAAGATGATGCATTGCCATTGCTACCtgttgt     <  1:1483668/100‑1 (MQ=255)                                                                                        cAGCTTGAGCGCGACGGTCAAATACGTACCATGTACATTGACCAGCAGTGTATAAATTTTGGTGATTAAAAGATGATGCATTGCCATTGCTACCtgttgt     <  1:95650/100‑1 (MQ=255)                                                                                         ggCTTGAGCGCGACGGTCAAATACGTACCATGTACATTGACCAGCAGTGTATAAATTTTGGTGATTAAAAGATGATGCATTGCCATTGCTACCTGTTGTa    <  1:4095241/99‑1 (MQ=255)                                                                                          gCTTGAGCGCGACGGTCAAATACGTACCATGTACATTGACCAGCAGTGTATAAATTTTGGTGATTAAAAGATGATGCATTGCCATTGCTACCTGTTGTac   <  1:55393/100‑2 (MQ=255)                                                                                          gCTTGAGCGCGACGGTCAAATACGTACCATGTACATTGACCAGCAGTGTATAAATTTTGGTGATTAAAAGATGATGCATTGCCATTGCTACCTGTTGTAg   <  1:1564804/100‑1 (MQ=255)                                                                                           cTTGAGCGCGACGGTCAAATACGTACCATGTACATTGACCAGCAGTGTATAAATTTTGGTGATTAAAAGATGATGCATTGCCATTGCTACCTGTTGTAgg  <  1:3266076/100‑2 (MQ=255)                                                                                                  |                                                                                             ATGGTGTGTTGTTTACTTGGTAACCATCATTAGCTGCGTTACCAGCCCAATACTTAGCGTCTGACCAATATGTGCTAATTGGACTACCAGCTTGAGCACGACGGTCAAATACGTACCATGTACATTGACCAGCAGTGTATAAATTTTGGTGATTAAAAGATGATGCATTGCCATTGCTACCTGTTGTAGC  >  NC_007793/726641‑726830

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