Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_007793 1,929,062 G→A intergenic (‑148/‑273) SAUSA300_RS09560 ← / → SAUSA300_RS09565 transposase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0077931,929,0620GA83.3% 18.9 / ‑0.9 13intergenic (‑148/‑273)SAUSA300_RS09560/SAUSA300_RS09565transposase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (1/0);  major base A (5/5);  minor base T (2/0);  total (8/5)
Fisher's exact test for biased strand distribution p-value = 4.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01

TTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATC  >  NC_007793/1928964‑1929155
                                                                                                  |                                                                                              
ttGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAAtagt                                                                                               <  1:4166739/100‑1 (MQ=255)
ttGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAAtagt                                                                                               <  1:628324/100‑1 (MQ=255)
                                                                              tattcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAA‑‑TCGACGAAGAGAGGAGGTGCATTTCTTTGCTAATCATTTTCGTTCAcatcat               <  1:1589450/92‑1 (MQ=255)
                                                                              tattcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAA‑‑TCGACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCAcatcat               <  1:4687136/92‑1 (MQ=255)
                                                                              tattcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAA‑‑TCGACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCAcatcat               <  1:4904910/92‑1 (MQ=255)
                                                                                cttGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAAGGTGATGCTTATTTTCGTTCACATCATa              >  1:373252/2‑100 (MQ=255)
                                                                                  cccagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAA‑‑TCGACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCa           >  1:4187742/5‑100 (MQ=255)
                                                                                  cccagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAA‑‑TCGACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCATTCACATCATAGCa           >  1:4996531/5‑100 (MQ=255)
                                                                                  cagagatgatATAATATACAAGTCGCCTATCTCTCAGTCGTCAACTTAGGCGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAg             >  1:5084312/5‑100 (MQ=255)
                                                                                  aacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAA‑‑TCGACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCa           >  1:2143166/5‑100 (MQ=255)
                                                                                       atgatATAATATACAAGTCGCCTATCTCTCAGTCGTCAACTTAGGCGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCa        >  1:4059432/1‑100 (MQ=255)
                                                                                        tggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAA‑‑TCGACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAgat     >  1:804867/4‑98 (MQ=255)
                                                                                        tggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAA‑‑TCGACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGtc     <  1:2803287/97‑1 (MQ=255)
                                                                                        tggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAA‑‑TCGACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTCCGTTCACATCATAGCACCAgat     >  1:109116/4‑98 (MQ=255)
                                                                                           tataATAATAGAGTCGCCTATCTCTCAGGCGTCAA‑‑TCGACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGtcatc  <  1:5032312/100‑1 (MQ=255)
                                                                                                  |                                                                                              
TTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATC  >  NC_007793/1928964‑1929155

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑