Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	1,923,107	0	A	T	78.6%	29.6 / 2.4	30	I33F (ATT→TTT)	SAUSA300_RS09510	calcium‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (5/1);  new base T (22/0);  total (29/1)
Fisher's exact test for biased strand distribution p-value = 2.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.51e-01

TTGAATACAGAGAACAACAAGAATCAAAACCAATCTGTTAAAAATTCTGAAAGACGCGGCATGTTAAAAGGATGCGGCGGTTGCCTTATTTCTTTTATTTTATTAATAATCTTATTATCAGCCTGTTCAATGATGTTTAGTAATAATGACAATTCCACTAATAATCAATCATCAAA  >  NC_007793/1923011‑1923186                                                                                                 |                                                                                ttGAATACAGAGAACAACAAGAATCAAAACCAATCTGTTAAAAATTCTGAAAGACGCGGCATGTTAAAAGGATGCGGCGGGTTACCTCTTTCtttttttt                                                                              >  1:3639963/1‑100 (MQ=255)   gAATACAGAGAACAACAAGAATCAAAACCAATCTGTTAAAAATTCTGAAAGACGCGGCATGTTAAAAGGATGCGGCGGTTTCCCTATttttttttttttt                                                                            >  1:3046520/1‑99 (MQ=255)     aTACAGAGAACAACAAGAATCAAAACCAATCTGTTAAAAATTCTGAAAGACGCGGCATGTTAAAAGGATGCGGCGGGTGGCTTAtttttttttttttttt                                                                          >  1:5388386/1‑100 (MQ=255)               gtgcaaGAATCAAAACCAATCTGTTAAAAATTCTGAAAGACGCGGCATGTTAAAAGGATGCGGCGGTTGCCTTTTTTTTTTTTTTTTTTTTATAATCttt                                                                >  1:549578/4‑99 (MQ=255)               caacaaGAATCAAAACCAATCTGTTAAAAATTCTGGAAGACGCGGCATGTTAAAAGGATGCGGCGGTTATCATATTTTTTTTCTTTTTTTTAAAATTtta                                                                >  1:5387561/1‑100 (MQ=255)                           aaaCCAATCTGTTAAAAATTCTGAAAGACGCGGCATGTTAAAAGGATGCGGCGGTTGCCTTATTTCTTTTATTTTATTAATAATCTTATTATCAGCCTGt                                                    <  1:4193130/100‑1 (MQ=255)                           aaaCCAATCTGTTAAAAATTCTGAAAGACGCGGCATGTTAAAAGGATGCGGCGGATTTCTTCTTTTTTTTCTTTTTTTCATAATTCTATTATCCACCTGt                                                    >  1:4234267/1‑100 (MQ=255)                                     gTTAAAAATTCTGAAAGACGCGGCATGTTAAAAGGATGCGGCGGATTCCTTATTTTTTTTTTTTTTTTAATAAACTTTTTTTCAACCTGTTCTATGAAGt                                          >  1:4149704/1‑100 (MQ=255)                                         aaaaTTCTGAAAGACGCGGCATGTTAAAAGGATGCGGCGGTTTCCTTTTTTCTTTTTATTTATTAATAATCTTATTATTAGCCTGTTTAATGTTGTTTAg                                      >  1:5444543/1‑100 (MQ=255)                                           aaTTCTGAAAGACGCGGCATGTTAAAAGGATGCGGCGGTTGCCTTATTTCTTTTTTTTTTTTAATAAACCTTTTTTTAACCTGCTCCATGATGTTTTGta                                    >  1:5660833/1‑100 (MQ=255)                                             cccTGAAAGACGCGGCATGTTAAAAGGATGCGGCGGTTGCCTTTTTTTTTTTTTTTTTTTAATAATCTTTTTTTCAGCCCGTTCAATTATGTTTTGtaat                                  >  1:4444223/3‑100 (MQ=255)                                               cTGAAAGACGCGGCATGTTAAAAGGATGCGGCGGTTGCCTTTTTTCTTTTTTTTTTTTTAAAATCTTTTTTTCAGGCCGGTTTATGATGGTTAGtaataa                                >  1:3760350/1‑100 (MQ=255)                                                 gAAAGACGCGGCATGTTAAAAGGATGCGGCGGTTTCCTTTTTTTTTTTATTTTTTTTATAATCTTTTTTTTCGCCTGTTCAATTATTTTTTTTtataatt                              >  1:3770324/1‑99 (MQ=255)                                                 gAAAGACGCGGCATGTTAAAAGGATGCGGCGGTTGCCTTATTTTTTTTTTTTTTTTAATAATCTTTTTATCAGCCCGTTCAATGATGTTTTGTAATAATg                              >  1:4630972/1‑100 (MQ=255)                                                 gAAAGACGCGGCATGTTAAAAGGATGCGGCGGTTGCCTTATTTTTTTTTTTTTATTAATAATCTTATTTTCAGCCTGTTCAATGATGTTTAGTaataatt                              >  1:3640134/1‑99 (MQ=255)                                                 gAAAGACGCGGCATGTTAAAAGGATGCGGCGGTTGCCCTATTTCTTTTTTTTTATTTAAAATCTTATTATTATCCTGGTCAATGTTGGTTTGAAATAATg                              >  1:5612177/1‑100 (MQ=255)                                                 gAAAGACGCGGCATGTTAAAAGGATGCGGCGGGTGCCTCTTTTTTTTTTTTTTATTAATAATATTATTAACAGCCCGGTCAAAGATGTTTTTTaataata                              >  1:5262961/1‑99 (MQ=255)                                                 gAAAGACGCGGCATGTTAAAAGGATGCGGCGGCTGCCTTATTTTTTTTTTTTTTTTAATAATCTTATTTTCAGCCCGTTCTATGATGTTTTGTaataatt                              >  1:585226/1‑99 (MQ=255)                                                  aaaGACGCGGCATGTTAAAAGGATGCGGCGGTTTCCTTTTTTCTTTTATTTTTTTAATAATCTTTTTTTTAACCTGTTTAATGATGTTTAGTAATAATGa                             >  1:4704156/1‑100 (MQ=255)                                                       cacgGCATGTTAAAAGGATGCGGCGGTTGTCTTATTTTTTTTTTTTTTTTTATAAACTTATTATCAGCCTCTTCAATTATGTTTTGTAATAATGACAAtt                        >  1:5327098/3‑100 (MQ=255)                                                        ccgGCATGCTAAAAGGATGCGGCGGATACCTTTTTTTTTTTTTTTTTTTAATAATCTTATTTTCAACCTCGCTAATCATTTTTTGAAAAAATTACAAttt                       >  1:5634402/2‑99 (MQ=255)                                                        acgGCATGTTAAAAGGATGCGGCGGTTGCCCTTTTTTTTTTTTTTTTTTTATAATCTTATTTTCAGCCCGTTCACTGATGTTTTGTAATAATGACCAttt                       >  1:634114/2‑99 (MQ=255)                                                         cgGCATGTTAAAAGGATGCGGCGGTTTCCTTTTTTTTTTTTTTTTTTTTATAATCTTTTTTTCAACCTGTTCCATTATGTTTTGTAATAATGACAATTcc                      >  1:4019894/1‑100 (MQ=255)                                                         cgGCATGTTAAAAGGATGCGGCGGTTTCCTTATTTCTTTTTTTTTTTTAATAATCTTATTATCAGCCTGTTCAATGATGTTTAGTAATAATGACAATTcc                      >  1:4714326/1‑100 (MQ=255)                                                         cgGCATGTTAAAAGGATGCGGCGGTTGTCCTATTTCTTTTTTTTTTTTAATAATCTTATTTTTAACCTGTTTAATGATGTTTAGTAAAAATAACAAtttc                      >  1:3067240/1‑98 (MQ=255)                                                         cgGCATGTTAAAAGGATGCGGCGGTTGGCTTTTTTTTTTTTTTTTATTTATAATCTTTTTATTAGCCTGGTGAAAGATGTTTAGTAATAATGACAAttac                      >  1:5376431/1‑98 (MQ=255)                                                          ggCATGTTAAAAGGATGCGGCGGTTTCCTTTTTTCTTTTTTTTTTTTTATAATCTTTTTTTCAGCCTGTTCCATGATGTTTTGTAATAATTACAATTccc                     >  1:4806507/1‑99 (MQ=255)                                                          ggCATGTTAAAAGGATGCGGCGGTTTCCTTATTTCTTTTATTTTATTAATAATCTTATTATTAGCCTCTTCCATGATGTTTATTAATAATTACAATTccc                     >  1:5665323/1‑99 (MQ=255)                                                          ggCATGTTAAAAGGATGCGGCGGTTTCCCTATTTATTTTTTTTTTTTACTAATATTATTATCAGCCTGTTTAATTTTGTTTTGTAATAATGACAATTTCa                     >  1:5478271/1‑100 (MQ=255)                                                          ggCATGTTAAAAGGATGCGGCGGTTGCCTTATTTTTTTTTTTTTTTTTATAATATTTTTTTCAGCCCGTTCAATGATGTTTAGTAATAATGACAAATCCa                     >  1:1569605/1‑100 (MQ=255)                                                           gCATGTTAAAGGGATGCGGCGGTTTGCCTCTTTTTTTTATTTTATTAAAAATCTTATTTTTAGCCTGTTCAATTATGTTTTGTAATAATGACAATTTCAc                    >  1:2729767/1‑100 (MQ=255)                                                                  aaaaGGATGCGGCGGTTTCCTTTTTTCTTTTTTTTTTTTAATAATCTTTTTATCAGCCTGTTCAATGATGTTTTGTAATAATTACAATTTCCCTaataat             >  1:3793352/1‑100 (MQ=255)                                                                             tcggTTGCCTTATTTCTTTTATTTTATTAATAATCTTATTATCAGCCTGTTCAATGATGTTTAGTAATAATGACAATTCCACTAATAATCAATCATCaaa  >  1:3291367/2‑100 (MQ=255)                                                                                                 |                                                                                TTGAATACAGAGAACAACAAGAATCAAAACCAATCTGTTAAAAATTCTGAAAGACGCGGCATGTTAAAAGGATGCGGCGGTTGCCTTATTTCTTTTATTTTATTAATAATCTTATTATCAGCCTGTTCAATGATGTTTAGTAATAATGACAATTCCACTAATAATCAATCATCAAA  >  NC_007793/1923011‑1923186

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