Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0077931,929,0620GA72.2% 21.0 / 4.1 18intergenic (‑148/‑273)SAUSA300_RS09560/SAUSA300_RS09565transposase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (1/4);  new base A (6/7);  total (7/11)
Fisher's exact test for biased strand distribution p-value = 5.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.81e-01

TTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑GG‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATC  >  NC_007793/1928964‑1929155
                                                                                                  |                                                                                                
ttGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAAtagt                                                                                                 <  1:4238760/100‑1 (MQ=255)
ttGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAAtagt                                                                                                 <  1:2296237/100‑1 (MQ=255)
      aaTTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGta                                                                                           <  1:2581025/100‑2 (MQ=255)
                             tATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAAttt                                                                    <  1:3357040/100‑1 (MQ=255)
                                                                             atattcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TT‑TGCTAATCATTTTCGTTCAcatca                <  1:2066576/91‑1 (MQ=255)
                                                                              tattcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑CTTTGCTAATCATTTTCGTTCAcatcat               <  1:4730538/92‑1 (MQ=255)
                                                                                 ttttcGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑GG‑TGATGCTTATTTTCGTTcacagggctt             <  1:5032830/96‑7 (MQ=255)
                                                                                 ttGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑GG‑TGATGCTTATTTTCGTTCACATCATAg             >  1:3041997/1‑100 (MQ=255)
                                                                                  cacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TT‑TGCTAATCATTTTCGTTCACATCATAGcg           <  1:1922213/96‑2 (MQ=255)
                                                                                  cacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TT‑TGCTAATCATTTTCGTTCACATCATAGCa           >  1:235597/5‑100 (MQ=255)
                                                                                  cacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TT‑TGCTAATCATTTTCGTTCACATCATAGCa           >  1:202071/5‑100 (MQ=255)
                                                                                   ccagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TT‑TGCTAATCATTTTCGTTCACATCATGGCAc          >  1:2033958/4‑100 (MQ=255)
                                                                                   acagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TT‑TGCTAATCATTTTCGTTCAcatcagatcgg          >  1:1463683/4‑94 (MQ=255)
                                                                                   acagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACTT‑TGCTAATCATTTTCGTTCACATCATAGCAc          <  1:400538/97‑1 (MQ=255)
                                                                                    cagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACTT‑TGCTAATCATTTTCGTTCACATCATAGCAcc         <  1:6181703/98‑1 (MQ=255)
                                                                                     agatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TT‑TGCTAATCATTTTCGTTCACATCATAGCACCa        >  1:2248932/2‑100 (MQ=255)
                                                                                     agatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TT‑TGCTAATCATTTTCGTTCACATCATAGCACCa        <  1:4515294/99‑1 (MQ=255)
                                                                                       ctggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGt      >  1:613437/5‑100 (MQ=255)
                                                                                       atgatATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑GG‑TGATGCTTATTTTCGTTCACATCATAGCACCAg       <  1:4611488/100‑1 (MQ=255)
                                                                                           tataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAAGTGCA‑TTACTT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGtcatc  <  1:5856890/100‑1 (MQ=255)
                                                                                                  |                                                                                                
TTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑GG‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATC  >  NC_007793/1928964‑1929155

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