Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	2,614,129	G→A	P131L (CCA→CTA)	SAUSA300_RS13460 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	2,614,129	0	G	A	96.7%	80.1 / ‑4.3	30	P131L (CCA→CTA)	SAUSA300_RS13460	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (0/29);  minor base T (0/1);  total (0/30)
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.37e-01

AGTTTTCGATTTCTTTTCTTAATTTATCATCAGCTATCGGTTTGGTTGGGATAATTTTATTATGTTCCATTTTTACTGGATATTCTTTTTCTTTACTACGTGAGTATCCTTTTTTATCTTCCCATAATTCCCTCACAATAAATTCACCTTTTGTCGTTCTGGTATTGCGGTCT  >  NC_007793/2614052‑2614224                                                                              |                                                                                                ctcttccGATCTCTTTTCTTAATTTATCATCAGCTATCGGTTTGGTTGGGATAATTTTATTATGTTCCATTTTTACTAGATATTCTTTTTCTTTACTACg                                                                           <  1:3742760/94‑1 (MQ=37)  gctcttccgatctcttcttAATTTATCATCAGCTATCGGTTTGGTTGGGATAATTTTATTATGTTCCATTTTTACTAGATATTCTTTTTCTTTACTACGt                                                                          <  1:34733/90‑1 (MQ=37)   ttttCGATTTCTTTTCTTAATTTATCATCAGCTATCGGTTTGGTTGGGATAATTTTATTATGTTCCATTTTTACTAGATATTCTTTTTCTTTACTACGTg                                                                         <  1:1093500/100‑1 (MQ=255)   ttttCGATTTCTTTTCTTAATTTATCATCAGCTATCGGTGTGGTTGGGATAATTTTATTATGTTCCATTTTTACTAGATATTCTTTTTCTTTACTACGTg                                                                         <  1:5493670/100‑1 (MQ=255)   ttttCGATTTCCTTTCTTAATTTATCATCAGCTATCGGTTTGGTTGGGATAATTTTATTATGTTCCATTTTTACTAGATATTCTTTTTCTTTACTACGTg                                                                         <  1:6052599/100‑1 (MQ=255)   ttccgatcttcttttcttAATTTATCATCAGCTATCGGTTTGGTTGGGATAATTTTATTATGTTCCATTTTTACTAGATATTCTTTTTCTTTACTACggg                                                                         <  1:3721134/92‑3 (MQ=38)   ttccgatcttcttttcttAATTTATCATCAGCTATCGGTTTGGTTGGGATAATTTTATTATGTTCCATTTTTACTAGATATTCTTTTTCTTTACTACGTg                                                                         <  1:4166849/92‑1 (MQ=38)    tttCGATTTCTTTTCTTAATTTATCATCAGCTATCGGTTTGGTTGGGATAATTTTATTATGTTCCATTTTTACTAGATATTCTTTTTCTTTACTACGTgg                                                                        <  1:4474317/100‑2 (MQ=255)      tCGATTTCTTTTCTTAATTTATCATCAGCTATCGGTTTGGTTGGGATAATTTTATTATGTTCCATTTTTACTAGATATTCTTTTTCTTTACTACGTGAGt                                                                      <  1:2278394/100‑1 (MQ=255)      tCGATTTCTTTTCTCAATTTATCATCAGCTATCGGTTTGGTTGGGATAATTTTATTATGTTCCATTTTTACTAGATATTCTTTTTCTTTACTACGTGAGt                                                                      <  1:4995314/100‑1 (MQ=255)      tCGATTTCGTTTCTTAATTTATCATCAGCTATCGGTTTGGTTGGGATAATTTTATTATGTTCCATATTTACTAGATATTCTTTTTCTTTACTACGTGAGt                                                                      <  1:6503691/100‑1 (MQ=255)       cGATCTCCTTTCTTAATTTATCATCAGCTATCGGTTTGGTTGGGATAATTTTATTATGTTCCATTTTTACTAGATATTCTTTTTCTTTACTACGTGAGtc                                                                     <  1:2339627/100‑2 (MQ=255)                       tttATCATCAGCTATCGGTTTGGTTGGGATAATTTTATTATGTTCCATTTTTACTAGATATTCTTTTTCTTTACTACGTGAGTATCCTTTTTTaccccat                                                     <  1:2312977/100‑7 (MQ=255)                           tcatcaGCTATCGGTTTGGTTGGGATAATTTTATTATGTTCCATTTTTACTAGATATTCTTTTTCTTTACTACGTGAGTATCCTTTTTTATCTTCCCATa                                                 <  1:5162251/100‑1 (MQ=255)                                  cTATCGGTTTGGTTGGGATAATTTTATTATGTTCCGTTTTTACTAGATATTCTTTTTCTTTACTACGTGAGTATCCTTTTTTATCTTCCCATAATTCCCt                                          <  1:354322/100‑1 (MQ=255)                                   tATCGGTTTGGTTGGGATAATTTTATTATGTTCCATTTTTACTAGATATTCTTTTTCTTTACTACGTGAGTATCCTTTTTTATCTTCCCATAATTCCCTc                                         <  1:60696/100‑1 (MQ=255)                                      cGGTTTGGTTGGGATAATTTTATTATGTTCCATTTTTACTAGATATTCTTTTTCTTTACTACGTGAGTATCCTTTTTTATCTTCCCATAATTCCCTcaca                                      <  1:21120/100‑1 (MQ=255)                                      cGGTTTGGTTGGGACAATTTTATTATGTTCCATTTTTACTAGATATTCTTTTTCTTTACTACGTGAGTATCCTTTTTTATCTTCCCATAATTCCCTcaca                                      <  1:4842374/100‑1 (MQ=255)                                       ggTTTGGTTGGGATAATTTTATTATGTTCCATTTTTACTAGATATTCTTTTTCTTTACTACGTGAGTATCCTTTTTTATCTTCCCATAATTCCCTCACaa                                     <  1:1316317/100‑1 (MQ=255)                                        gTTTGGTTGGGATAATTTTATTATGTTCCATTTTTACTAGATATTCTTTTTCTTTACTACGTGAGTATCCTTTTTTATCTTCCCATAATTCCCTCACAAt                                    <  1:4869231/100‑1 (MQ=255)                                              ttggGATAATTTTATTATGTTCCATTTTTACTAGATATTCTTTTTCTTTACTACGTGAGTATCCTTTTTTATCTTCCCATAATTCCCTCACAATAAATTc                              <  1:5586902/100‑1 (MQ=255)                                               tggGATAATTTTATTATGTTCCATTTTTACTAGATATTCTTTTTCTTTACTACGTGAGTATCCTTTTTTATCTTCCCATAATTCCCTCACAATAAATTCa                             <  1:4517102/100‑1 (MQ=255)                                                      aTTTTATTATGTTCCATTTTTACTAGATATTCTTTTTCTTTACTACGTGAGTATCCTTTTTTATCTTCCCATAATTCCCTCACAATAAATTCACCTTTTg                      <  1:2653735/100‑1 (MQ=255)                                                              atGTTCCATTTTTACTAGATATTCTTTTTCTTTACTACGTGAGTATCCTTTTTTATCTTCCCATAATTCCCTCACAATAAATTCACCTTTTGTCGTTCTg              <  1:4093115/100‑1 (MQ=255)                                                              atGTTCCATTTTTACTAGATATTCTTTTTCTTTACTACGTGAGTATCCTTTTTTATCTTCCCATAATTCCCTCACAATAAATTCACCTTTTGTCGTTCTg              <  1:4012911/100‑1 (MQ=255)                                                              atGTTCCATTTTTACTAGATATTCTTTTTCTTTACTACGTGAGTATCCTTTTTTATCTTCCCATAATTCCCTCACAATAAATTCACCTTTTGTCGTTCTg              <  1:3900851/100‑1 (MQ=255)                                                                gTTCCATTTTTACTAGATATTCTTTTTCTTTACTACGTGAGTATCCTTTTTTATCGTCCCATAATTCCCTCACAATAAATTCACCTTTTGTCGTTCTGGt            <  1:1508472/100‑1 (MQ=255)                                                                       ttttACTAGATATTCTTTTTCTTTACTACGTGAGTATCCTTTTTTATCTTCCCATAATTCCCTCACAATAAATTCACCTTTTGTCGTTCTgggcgtgtta     <  1:1751277/100‑9 (MQ=38)                                                                          tACTTGATATTCTTTTTCTTTACTACGTGAGTATCCTTTTTTATCTTCCCATAATTCCCTCACAATAAATTCACCTTTTGTCGTTCTGGTATTGCGGTCt  <  1:3185084/100‑1 (MQ=255)                                                                          tACTAGATATTCTTCTTCTTTACTACGTGAGTATCCTTTTTTATCTTCCCATAATTCCCTCACAATAAATTCACCTTTTGTCGTTCTGGTATTGCggcgt  <  1:4801373/100‑4 (MQ=255)                                                                              |                                                                                                AGTTTTCGATTTCTTTTCTTAATTTATCATCAGCTATCGGTTTGGTTGGGATAATTTTATTATGTTCCATTTTTACTGGATATTCTTTTTCTTTACTACGTGAGTATCCTTTTTTATCTTCCCATAATTCCCTCACAATAAATTCACCTTTTGTCGTTCTGGTATTGCGGTCT  >  NC_007793/2614052‑2614224

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: ‑