Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0077932,249,0090AT73.7% 18.0 / 8.6 19intergenic (‑396/+3)SAUSA300_RS11480/SAUSA300_RS11485hypothetical protein/ATP‑grasp domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base A (2/3);  new base T (0/14);  total (2/17)
Fisher's exact test for biased strand distribution p-value = 5.85e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.05e-02

TCGAAAAAAGCTTGTTACAAACGCATATTCAGTCAACTACTGGCAACACAACATTGTAGCGCCTTATGTTCCGAAC‑TC‑TAATAGTGTTCAAGTTCAGCTAAAAATGCGCGTTCATTCTCCAGTACTGTTTCTTTTGAATCACCATGCCAAATACAGCCAAACCTTGA  >  NC_007793/2248935‑2249101
                                                                          |                                                                                              
tCGAAAAAAGCTTGTTACAAACGCATATTCAGTCAACTACTGGCAACACAACATTGTAGCGCCTTATGTTCCGAAC‑TC‑TAATAGTGTTCAAGTTCAGCTa                                                                     <  1:93688/100‑1 (MQ=255)
                       cATATTCAGTCAACTACTGGCAACACAACATTGTAGCGCCTTATGTTCCGAAC‑TC‑TAATAGTGTTCAAGTTCAGCTAAAAATGCGCGTTCATTCTCCacg                                              <  1:3536060/100‑3 (MQ=255)
                                      actGGCAACACAACATTGTAGCGCCTTATGTTCCGAAC‑TC‑TAATAGTGTTCAAGTTCAGCTAAAAATGCGCGTTCATTCTCCAGTACTGTTTCTTTTGaa                               >  1:187323/1‑100 (MQ=255)
                                             acacaacaTTGTAGCGCCTTATGTTCCGAAC‑TC‑TAATAGTGTTCAAGTTCAGCTAAAAATGCGCGTTCATTCTCCAGTACTGTTTCTTTTGAATCACCAt                        >  1:1237608/1‑100 (MQ=255)
                                                          acgtgtgctcttCCGATC‑TC‑TTTTAGTGTTCAAGTTCAGCTAAAAATGCGCGTTCATTCTCCAGTACTGTTTCTTTTGAATCACCATGCCAAATACAGcc           <  1:335663/90‑1 (MQ=255)
                                                          acgtgtgctcttCCGATC‑TC‑ATTTAGTGTTCAAGTTCAGCTAAAAATGCGCGTTCATTCTCCAGTACTGTTTCTTTTGAATCACCATGCCAAATACAGcc           <  1:2793339/90‑1 (MQ=255)
                                                           cgtgtgctcttCCGATC‑TC‑TAATGGTGTTCAAGTTCAGCTAAAAATGCGCGTTCATTCTCCAGTACTGTTTCTTTTGAATCACCATGCCAAATACAGCCa          <  1:706430/91‑1 (MQ=255)
                                                           cgtgtgctcttCCGATC‑TC‑TAATAGTGTTCAAGTTCAGCTAAAAATGCGCGTTCATTCTCCAGTACTGTTTCTTTTGAATCACCATGCCAAATACAGCCa          <  1:2979958/91‑1 (MQ=255)
                                                           cgtgtgctcttCCGATC‑TC‑TAATAGTGTTCAAGTTCAGCTAAAAATGCGCGTTCATTCTCCAGTACTGTTTCTTTTGAATCACCATGCCAAATACAGCCa          <  1:5921890/91‑1 (MQ=255)
                                                           cgtgtgctcttCCGATC‑TC‑AAATAGTGTTCAAGTTCAGCTAAAAATGCGCGTTCATTCTCCAGTACTGTTTCTTTTGAATCACCATGCCAAATACAGCCa          <  1:5328761/91‑1 (MQ=255)
                                                            gtgtgctcttCCGATCTTC‑TAATAGTGTTCAAGTTCAGCTAAAAATGCGCGTTCATTCTCCAGTACTGTTTCTTTTGAATCACCATGCCAAATACAGCCa          <  1:6138872/92‑1 (MQ=255)
                                                            gtgtgctcttCCGATC‑TC‑TATTAGTGTTCAAGTTCAGCTAAAAATGCGCGTTCATTCTCCAGTACTGTTTCTTTTGAATCACCATGCCAAATACAGCCaa         <  1:1592622/92‑1 (MQ=255)
                                                            gtgtgctcttCCGATC‑TC‑TAATAGTGTTCAAGTTCAGCTAAAAATGCGCGTTCATTCTCCAGTACTGTTTCTTTTGAATCACCATGCCAAATACAGCCaa         <  1:6011610/92‑1 (MQ=255)
                                                             tgtgctcttCCGATC‑TCTTAATAGTGTTCAAGTTCAGCTAAAAATGCGCGTTCATTCTCCAGTACTGTTTCTTTTGAATCACCATGCCAAATACAGCCaa         <  1:4981433/93‑1 (MQ=255)
                                                              gtgctcttCCGATC‑TC‑TCTTAGTGTTCAAGTTCAGCTAAAAATGCGCGTTCATTCTCCAGTACTGTTTCTTTTGAATCACCATGCCAAATACAGCCAAAc       <  1:3578935/94‑1 (MQ=255)
                                                              gtgctcttCCGATC‑TC‑TAACAGTGTTCAAGTTCAGCTAAAAATGCGCGTTCATTCTCCAGTACTGTTTCTTTTGAATCACCATGCCAAATACAGCCAAAc       <  1:5726070/94‑1 (MQ=255)
                                                                gctcttCCGATC‑TC‑T‑ATAGTGTTCAAGTTCAGCTAAAAATGCGCGTTCATTCTCCAGTACTGTTTCTTTTGAATCACCATGCCAAATACAGCCAAACCtt    <  1:4217876/96‑1 (MQ=255)
                                                                 ctcttCCGATC‑TC‑TCTTAGTGTTCAAGTTCAGCTAAAAATGCGCGTTCATTCTCCAGTACTGTTTCTTTTGAATCACCATGCCAAATACAGCCAAACCtt    <  1:5597505/97‑1 (MQ=255)
                                                                  tcttCCGATC‑TC‑TAACAGTGTTCAAGTTCAGCTAAAAATGCGCGTTCATTCTCCAGTACTGTTTCTTTTGAATCACCATGCCAAATACAGCCAAACCTTg   <  1:2637421/98‑1 (MQ=255)
                                                                   gTTCCGAAC‑TC‑TAATAGTGTTCAAGTTCAGCTAAAAATGCGCGTTCATTCTCCAGTACTGTTTCTTTTGAATCACCATGCCAAATACAGCTAAAccagat  <  1:5520525/100‑5 (MQ=255)
                                                                          |                                                                                              
TCGAAAAAAGCTTGTTACAAACGCATATTCAGTCAACTACTGGCAACACAACATTGTAGCGCCTTATGTTCCGAAC‑TC‑TAATAGTGTTCAAGTTCAGCTAAAAATGCGCGTTCATTCTCCAGTACTGTTTCTTTTGAATCACCATGCCAAATACAGCCAAACCTTGA  >  NC_007793/2248935‑2249101

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