Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	2,113,184	0	A	G	70.4%	20.6 / 14.5	28	F61L (TTT→CTT)	SAUSA300_RS10710	acetyltransferase
Reads supporting (aligned to +/- strand):  ref base A (1/7);  new base G (9/10);  total (10/18)
Fisher's exact test for biased strand distribution p-value = 1.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.59e-01

TCTTCTATATATAGCTACTTCCTTGTTCGCTCTGCTCTCAATCTGTAGTTCGCTAAGTCTAAAACGGTAAAATTCGTATCTTCCAAACAATTCATTTTTAAGCGTGCGCCACATGTTCTCCAGCTCTTCGTTACGCTCTCTTAACTTAGCTATATCCCCAATAAGCTCGTCACGTTGCTTCTT  >  NC_007793/2113098‑2113280                                                                                       |                                                                                                 tcttctGTATATAGCTACTTCCTTGTTCGCTCTGCTCTCAATCTGTAGTTCGCTAAGTCTAAAACGGTAAAATTCGTATCTTCCAAACAATTCATTTTTa                                                                                     <  1:2643049/100‑1 (MQ=255)       atatataGCTACTTCCTTGTTCGCTCTGCTCTCAATCTGTAGTTCGCTAAGTCTAAAACGGTAAAATTCGTATCTTCCAAACAATTCATTTTTAAGCGTg                                                                               <  1:3906685/100‑1 (MQ=255)             aGCTACTTCCTTGTTCGCTCTGCTCTCAATCTGTAGTTCGCTAAGTCTAAAACGGTAAAATTCGTATCTTCCAAACAATTCATTTTTAAGCGTGCGCcac                                                                         <  1:5339004/100‑1 (MQ=255)                          ttCGCTTTGCTCTCAGGGTGTAGTTCTCTAAATTTAAAACAGTAATGTTCGTATCTTCCAAGCAATTCATTTTTGACTGTGCGCCACATGTTCTCCAGcg                                                            <  1:1039127/100‑2 (MQ=255)                              ctctgctctCAATCTGTAGTTCGCTAAGTCTAAAACGGTAAAATTCGTATCTTCCAAACAATTCATTTTTAAGCGTGCGCCACATGTTCTCCAGccatcg                                                        <  1:362607/100‑6 (MQ=255)                                   ctctCAATCTGTAGTTCGCTAAGTCTAAAACGGTAAAATTCGTATCTTCCAAACAATTCATTTTTAAGCGTGCGCCACATGTTCTCCAGCTCTTCGTTag                                                   <  1:5514309/100‑2 (MQ=255)                                     ctcagggtGTAGTTCTCTAAATTTAAAACAGTAATGTTCGTATCTTCCAAGCAATTCATTTTTGACTGTGCGCCACATGTTCTCCAGCTCTTCGTTACGc                                                 >  1:5764017/8‑100 (MQ=255)                                         tggtGTAGTTCTCTAAATTTAAAACAGTAATGTTCGTATCTTCCAAGCAATTCATTTTTGACTGTGCGCCACATGTTCTCCAGCTCTTCGTTACGctatc                                             >  1:4858793/4‑100 (MQ=255)                                          ggtGTAGTTCTCTAAATTTAAAACAGTAATGTTCGTATCTTCCAAGCAATTCATTTTTGACTGTGCGCCGCATGTTCTCCAGCTCTTCGTTACGctatga                                            <  1:786811/98‑5 (MQ=255)                                          ggtGTAGTTCTCTAAATTTAAAACAGTAATGTTCGTATCTTCCAAGCAATTCATTTTTGACTGTGCGCCACATGTTCTCCAGCTCTTCGTTACGctatga                                            <  1:5662437/98‑5 (MQ=255)                                           gtGTAGTTCTCTAAATTTAAAACAGTAATGTTCGTATCTTCCAAGCAATTCATTTTTGACTGTGCGCCACATGTTCTCCAGCTCTTCGTTACGCTATCtt                                           <  1:2964632/99‑1 (MQ=255)                                               aGTTCTCTAAATTTAAAACAGTAATGTTCGTATCTTCCAAGCAATTCATTTTTGACTGTGCGCCACATGTTCTCCAGCTCTTCGTTACGCTATCTTAACt                                       >  1:3038638/1‑100 (MQ=255)                                                 ttCGCTAAGTCTAAAACGGTAAAATTCGTATCTTCCAAACAATTCACTTTTAAGCGTGCGCCACATGTTCTCCAGCTCTTCGTTACGCTCTCTTAACTTa                                     <  1:5226265/100‑1 (MQ=255)                                                     cTAAATTTAAAACAGTAATGTTCGTATCTTCCAAGCAATTCATTTTTGACTGTGCGCCACATGTTCTCCAGCTCTTCGTTACGCTATCTTAACTTAGCta                                 >  1:3455713/1‑100 (MQ=255)                                                     cTAAATTTAAAACAGTAATGTTCGTATCTTCCAAGCAATTCATTTTTGACTGTGCGCCACATGTTCTCCAGCTCTTCGTTACGCTATCTTAACTTAGCta                                 >  1:460325/1‑100 (MQ=255)                                                      taaatttAAAACAGTAATGTTCGTATCTTCCAAGCAATTCATTTTTGACTGTGCGCCACATGTTCTCCAGCTCTTCGTTACGCTATCTTAACTTAGCtat                                <  1:2256481/94‑1 (MQ=255)                                                       aaatttAAAACAGTAATGTTCGTATCTTCCAAGCAATTCATTTTTGACTGTGCGCCACATGTTCTCCAGCTCTTCGTTACGCTATCTTAACTTAGCtatg                               <  1:205003/95‑2 (MQ=255)                                                            tAAAACAGTAATGTTCGTATCTTCCAAGCAATTCATTTTTGACTGTGCGCCACATGTTCTCCAGCTCTTCGTTACGCTATCTTAACTTAGCTATATccct                          <  1:4565642/100‑2 (MQ=255)                                                                aCGGTAAAATTCGTATCTTCCAAACAATTCATTTTTAAGCGTGCGCCACATGTTCTCCAGCTCTTCGTTACGCTCTCTTAACTTAGCTATATCCCCAATa                      <  1:318589/100‑1 (MQ=255)                                                                 cagtaatgttCGTATCTTCCAAGCAATTCATTTTTGACTGTGCGCCACATGTTCTCCAGCTCTTCGTTACGCTATCTTAACTTAGCTATATCCCCAATaa                     >  1:4341592/9‑100 (MQ=255)                                                                    taatgttCGTATCTTCCAAGCAATTCATTTTTGACTGTGCGCCACATGTTCTCCAGCTCTTCGTTACGCTATCTTAACTTAGCTATATCCCCAATAGGCt                  <  1:2077472/95‑1 (MQ=255)                                                                        cttccgaTCTTCCAAGCAATTCATTTTTGACTGTGCGCCACATGTTCTCCAGCTCTTCGTTACGCTATCTTAACTTAGCTATATCCCCAATAAGCTCGTc              <  1:4321560/94‑1 (MQ=255)                                                                               tCTTCCGATCTATTCATTTTTAAGCGTGCGCCACATGTTCTCCAGCTCTTCGTTACGCTCTCTTAACTTAGCTATATCCCCAATAAGCTCGTCACGTTGc       <  1:3943430/100‑1 (MQ=255)                                                                                 ttCCAAGCAATTCATTTTTGACTGTGCGCCACATGTTCTCCAGCTCTTCGTTACGCTATCTTAACTTAGCTATATCCCCAATAAGCTCGTCACGTTGctt     <  1:880033/100‑1 (MQ=255)                                                                                  tCCAAGCAATTCATTTTTGACTGTGCGCCACATGTTCTCCAGCTCTTCGTTACGCTATCTTAACTTAGCTATATCCCCAATAAGCTCGTCACGTTGcttc    >  1:2891234/1‑100 (MQ=255)                                                                                  cccAAACAATTCATTTTTAAGCGTGCGCCACATGTTCTCCAGCTCTTCGTTACGCTCTCTTAACTTAGCTATATCCCCAATAAGCTCGTCACGTTGcttc    >  1:983401/2‑100 (MQ=255)                                                                                   ccAAGCAATTCATTTTTGACTGTGCGCCACATGTTCTCCAGCTCTTCGTTACGCTATCTTAACTTAGCTATATCCCCAATAAGCTCGTCACGTTGcttct   >  1:598432/1‑100 (MQ=255)                                                                                    cAAGCAATTCATTTTTGACTGTGCGCCACATGTTCTCCAGCTCTTCGTTACGCTATCTTAACTTAGCTATATCCCCAATAAGCTCGTCACGTTGcttctt  >  1:745713/1‑100 (MQ=255)                                                                                       |                                                                                                 TCTTCTATATATAGCTACTTCCTTGTTCGCTCTGCTCTCAATCTGTAGTTCGCTAAGTCTAAAACGGTAAAATTCGTATCTTCCAAACAATTCATTTTTAAGCGTGCGCCACATGTTCTCCAGCTCTTCGTTACGCTCTCTTAACTTAGCTATATCCCCAATAAGCTCGTCACGTTGCTTCTT  >  NC_007793/2113098‑2113280

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑