Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_007793 1,929,092 Δ1 bp intergenic (‑178/‑243) 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,0920T.80.0% 48.6 / 15.9 25intergenic (‑178/‑243)SAUSA300_RS09560/SAUSA300_RS09565transposase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (2/3);  new base . (13/7);  total (15/10)
Fisher's exact test for biased strand distribution p-value = 3.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.61e-02

GAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTA  >  NC_007793/1929020‑1929191
                                                                        |                                                                                                       
gAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCTGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑Tg                                                                          <  1:5920263/100‑1 (MQ=255)
                      tattcacagatggtATAATAATAGAGTCGCCTATCTCTCGGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCAcatcat                                                   <  1:4047085/92‑1 (MQ=255)
                       ctttGCCGATGATATAACAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCAcatcat                                                   <  1:4670545/99‑1 (MQ=255)
                          cacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATGGCa                                               >  1:6020183/5‑100 (MQ=255)
                            cagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATGGCAcc                                             >  1:2567913/3‑100 (MQ=255)
                              gatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAg                                           <  1:940964/100‑1 (MQ=255)
                              gatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAg                                           <  1:6029152/100‑1 (MQ=255)
                                tggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGtc                                         <  1:4559012/97‑1 (MQ=255)
                                 ggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCaagat                                        >  1:2345257/3‑96 (MQ=255)
                                 ggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGtca                                        >  1:5967744/3‑100 (MQ=255)
                                    ataataataGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGtcatct                                     <  1:2886011/100‑2 (MQ=255)
                                        ttctccggTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑C‑TTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTgg                                 <  1:4812918/93‑1 (MQ=255)
                                          gtagAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGc                                <  1:4747805/100‑1 (MQ=255)
                                              ggTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgtg                            >  1:5470383/2‑100 (MQ=255)
                                               cccGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGct                          >  1:1826417/3‑99 (MQ=255)
                                                 ggCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgt                        >  1:5383683/2‑97 (MQ=255)
                                                 cGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgt                        >  1:4392138/1‑97 (MQ=255)
                                                      aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGt                   >  1:2898774/1‑100 (MQ=255)
                                                      aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGt                   <  1:1245541/100‑1 (MQ=255)
                                                        gaaccAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGtt                  >  1:638519/5‑100 (MQ=255)
                                                          ctcAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtatta               >  1:4849075/1‑99 (MQ=255)
                                                          ctcAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtatta               >  1:537072/1‑99 (MQ=255)
                                                          ctcAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑CTT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtatta               >  1:3621474/1‑99 (MQ=255)
                                                          ctcAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑C‑TTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtatta               >  1:2786967/1‑99 (MQ=255)
                                                           tcAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtattat              >  1:2676618/1‑98 (MQ=255)
                                                                       ttAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATGGGCt   >  1:5366614/1‑100 (MQ=255)
                                                                       ttAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCG‑TTTTTTCTTATTggggt  <  1:4831372/100‑4 (MQ=255)
                                                                        tAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTa  >  1:4931874/1‑100 (MQ=255)
                                                                        |                                                                                                       
GAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTA  >  NC_007793/1929020‑1929191

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