Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_007793 1,589,875 2 bp→CG coding (5‑6/126 nt) SAUSA300_RS07845 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0077931,589,8750TC100.0% 97.6 / NA 40V2A (GTA→GCA) SAUSA300_RS07845hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (37/3);  total (37/3)
*NC_0077931,589,8760AG97.5% 95.7 / ‑4.7 40V2V (GTA→GTG) SAUSA300_RS07845hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (37/2);  total (37/3)
Fisher's exact test for biased strand distribution p-value = 7.50e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.42e-01

AAAAAAATCAAAAAACGCCTACTAGTGTAGACGTTGAATGGTGGTGAGAATTTTATGGTAGATAAAAACAAAAAACAAGAAACTACTCGTAGTAACCCATTAAACAAAAGTTTTGAAAAGTCAGGCGCCAGCGAAAAATTAAAAAGCACTTTATC  >  NC_007793/1589817‑1589971
                                                          ||                                                                                               
ttactaaTCAAAAAACGCCTACAAGTGTAGACGTTGAATGGTGGTGAGAATTTTATGGCGGATAAAAACAAAAAACAAGAAGCTACTCGTAGTAACCCaa                                                         >  1:2422635/6‑99 (MQ=21)
     aaTCAAAAAACGCCTACAAGTGTAGACGTTGAATGGTGGTGAGAATTTTATGGCGGATAAAAACAAAAAACAAGAAGCTACTCGTAGTAACCCAATAAAc                                                    >  1:3584357/1‑100 (MQ=25)
                   tacaaGTGTAGACGTTGAATGGTGGTGAGAATTTTATGGCGGATAAAAACAAAAAACAAGAAGCTACTCGTAGTAACCCAATAAACAAAAGTTTTGaaaa                                      >  1:2046978/1‑100 (MQ=38)
                    acaaGTGTAGACGTTGAATGGTGGTGAGAATTTTATGGCGGATAAAAACAAAAAACAAGAAGCTACTCGTAGTAACCCAATAAACCAAAGTTTTGAAAAg                                     >  1:4948690/4‑100 (MQ=38)
                     caaGTGTAGACGTTGAATGGTGGTGAGAATTTTATGGCGGATAAAAACAAAAAACAAGAAGCTACTCGTAGTAACCCAATAAACAAAAGTTTTGAAAAgc                                    >  1:3996684/3‑99 (MQ=38)
                     aaaGTGTAGACGTTGAATGGTGGTGAGAATTTTATGGCGGATAAAAAAAAAAAACAAGAAGCTACTCGTAGGAACCCAATAAACAAAAGTTTTGAAAAgc                                    >  1:298778/3‑99 (MQ=38)
                         tgtAGACGTTGAATGGTGGTGAGAATTTTATGGCGGATAAAAACAAAAAACAAGAAGCTACTCGTAGTAACCCAATAAAAAAAAGTTTTGCAAAgccggg                                >  1:3375461/1‑95 (MQ=255)
                         cgtAGACGTTGAATGGTGGTGAGAATTTTATGGCGGATAAAAACAAAAAACAAGAAGCTACTCGTAGTAACCCAATAAACAAAAGTTTTGAAAAgccggg                                >  1:4670977/2‑95 (MQ=255)
                               tccgatctggtggTGAGAATTTTATGGCGGATAAAAACAAAAAACAAGAAGCTACTCGTAGTAACCCAATAAACAAAAGTTTTGAAAAGCCGGGTGCCAg                          <  1:4220980/93‑1 (MQ=255)
                                 ctGAATGGTGGTGAGAATTTTATGGCGGATAAAAACAAAAAACAAGAAGCTACTCGTAGTAACCCAATAAACAAAAGTTTTGAAAAGCCGGGTGCCAGCg                        >  1:3347942/2‑100 (MQ=255)
                                 ccgatctgtggtgAGAATTTTATGGCGGATAAAAACAAAAAACAAGAAGCTACTCGTAGTAACCCAATAAACAAAAGTTTGGAAAAGCCGGGTGCCAGCg                        <  1:3142452/93‑1 (MQ=255)
                                  tGAATGGTGGTGAGAATTTTATGGCGGATAAAAACAAAAAACAAGAAGCTACTCGTAGTAACCCAATAAACAAAAGTTTTGAAAAGCCGGGTGCCAGCGa                       >  1:1680375/1‑100 (MQ=255)
                                  tGAATGGTGGTGAGAATTTTATGGCGGATAAAAACAAAAAACAAGAAGCTACTCGTAGTAAACCAATAAACAAAAGTTTTGAAAAGCCGGGCGCCAGCGa                       >  1:1104808/1‑100 (MQ=255)
                                   gAATGGTGGTGAGAATTTTATGGCGGATAAAAACAAAAAACAAGAAGCTACTCGTAGTAACCCAATAAACAAAAGTTTTGAAAAGCCGGGTGCCAGCGaa                      >  1:5777351/1‑100 (MQ=255)
                                    tctggtggTGAGAATTTTATGGCGGATAAAAACAAAAAACAAGAAGCTACTCGTAGTAACCCAATAAACAAAAGTTTTGAAAAGCCGGGTGCCAGCGaaa                     >  1:3452206/3‑100 (MQ=255)
                                    taTGGTGGTGAGAATTTTATGGCGGATAAAAACAAAAAACAAGAAGCTACTCGTAGTAACCCAATAAACAAAAGTTTTGAAAAGCCGGGTGCCAGCGaaa                     >  1:5444817/2‑100 (MQ=255)
                                    gctggtggTGAGAATTTTATGGCGGATAAAAACAAAAAACAAGAAGCTACTCGTAGTAACCCAATAAACAAAAGTTTTGAAAAGCCGGGGGCCAGCGaaa                     >  1:3312694/3‑100 (MQ=255)
                                    gaTGGTGGTGAGAATTTTATGGCGGATAAAAACAAAAAACAAGAAGCTATTTGTAGCAACCCAATAAACAAAAGTTTTGAAAAGCCGGGTGCCAGCGaag                     >  1:953029/2‑99 (MQ=255)
                                    aaTGGTGGTGAGAATTTTATGGCGGATAAAAACAAAAAACAAGAAGCTACTCGTAGTAACCCAATAAACAAAAGTTTTGAAAAGCCGGGGGCCAGCGaaa                     >  1:828730/1‑100 (MQ=255)
                                    aaTGGTGGTGAGAATTTTATGGCGGATAAAAACAAAAAACAAGAAGCTACTCGTAGGAACCCAATAAACAAAAGTTTTGAAAAGCCGGGTGCCAGCGaaa                     >  1:3914006/1‑100 (MQ=255)
                                     ttggtggTGAGAATTTTATGGCGGATAAAAACAAAAAACAAGAAGCTACTCGTAGTAACCCAATAAACAAAAGTTTTGAAGAGCCGGGTGCCAGCGaaaa                    >  1:2617275/2‑100 (MQ=255)
                                     aTGGTGGTGAGAATTTTATGGCGGATAAAAACAAAAAACAAGAAGCTACTCGTAGTAACCCAATAAACAAAAGTTTTGAAAAGCCGGGTGCCAGCGaaaa                    >  1:4588470/1‑100 (MQ=255)
                                     aTGGTGGTGAGAATTTTATGGCGGATAAAAACAAAAAACAAGAAGCTACTCGTAGTAACCCAATAAACAAAAGTTTTGAAAAGCCGGGGGCCAGCGaaaa                    >  1:4740688/1‑100 (MQ=255)
                                     aTGGTGGTGAGAATTTTATGGCGGATAAAAACAAAAAACAAAAAGCTACTCGTAGTAACCCAATAAACAAAAGTTTTGAAAAGCCGGGTGCCAGCGaaaa                    >  1:626814/1‑100 (MQ=255)
                                      gtgtggtgAGAATTTTATGGCGGATAAAAACAAAAAACAAGAAGCTACTCGTAGTAACCCAATAAACAAAAGTTTTGAAAAGCCGGGGGCCAGCGaaaac                   >  1:2295226/3‑99 (MQ=255)
                                      gcctggtgAGAATTTTATGGCGGATAAAAACAAAAAACAAGAAGCTACTCGTAGTAACCCAATAAACAAAAGTTTTGAAAAGCCGGGGGCCAGCGaaaac                   >  1:2287713/4‑99 (MQ=255)
                                      ctgtggtgAGAATTTTATGGCGGATAAAAACAAAAAACAAGAAACTACTCGTAATAACCCAATAAACAAAAGCTTTGAAAAACCGGGGGGCAGCGaagat                   >  1:445744/3‑97 (MQ=255)
                                      acatggtgAGAATTTTATGGCGGATAAAAACAAAAAACAAGGAACTACTTGTAGTAAACCAATAAACAAAAGGGTTGAAAAGCCGGGGGCCAGCGgaaac                   >  1:6321703/4‑99 (MQ=255)
                                       ggtggtGAGAATTTTATGGCGGATAAAAACAAAAAACAAGAAGCTACTCGTAGTAACCCAATAAACAAAAGTTTTTAAAAGCCGGGTGCCAGCGaaaact                  >  1:5815377/1‑98 (MQ=255)
                                       ggtggtGAGAATTTTATGGCGGATAAAAACAAAAAACAAGAAGCTACTCGCAGTAACCCAATAAACAAAAGCTTTGAAAAGCCGGGTGCCAGCGaaaact                  >  1:5568292/1‑98 (MQ=255)
                                        gtggtgAGAATTTTATGGCGGATAAAAACAAAAAACAAGGAGCCACTCGTAGTAATCCAATAAACAAAAGCTTTGAAAAGCCCGGCGCCAgatcggaaga                 >  1:496662/1‑91 (MQ=255)
                                        gtggtgAGAATTTTATGGCGGATAAAAACAAAAAACAAGAAGCTACTCGTAGTAAACCAATAAACAAAAGCTTTGAAAAGCCGGGTGCCAGCGAAAACtt                 >  1:2350372/1‑100 (MQ=255)
                                        gtggtgAGAATTTTATGGCGGATAAAAACAAAAAACAAGAAGATACTCGTAGTAACCCAATAAACAAAAGTTTTGAAAAGCCGGGAGCCAGCGAAAACtt                 >  1:4786197/1‑100 (MQ=255)
                                          ggtgAGAATTTTATGGCGGATAAAAACAAAAAACAAGAAGCTACTCGTAGAAACCCAATAAACAAAAGTTTTGAAAAGCCGGGGGCCAGCGAAAACTTaa               >  1:5107228/1‑100 (MQ=255)
                                          ggtgAGAATTTTATGGCGGATAAAAACAAAAAACAAGAACCTCCTCGTTGTAACCCCATTAACAAAAGTTTTGAAAAGCCGGGAGCCAGCGAAAACTTaa               >  1:3111504/1‑100 (MQ=255)
                                                    ttATGGCGGATAAAAACAAAAAACAAGAAGCTACTCGTAGTAACCCAATAAACAAAAGTTTTGAAAAGCCGGGGGCCAGCGAAAACTTAAAAAGCACttt     >  1:3428327/1‑100 (MQ=255)
                                                    ctATGGCGGATAAAAACAAAAAACAAGAAGCTACTCGTAGTAACCCAATAAACAAAAGTTTTGAAAAGCCGGGTGCCAGCGAAAACTTAAAAAGCACttt     >  1:3021768/2‑100 (MQ=255)
                                                     tATGGCGGATAAAAACAAAAAACAAGAAGCTACTCGTAGTAACCCAATAAACAAAAGTTTTGAAAAGCCGGGTGCCAGCGAAAACTTAAAAAGCACTTTa    >  1:4147905/1‑100 (MQ=255)
                                                     tATGGCGGATAAAAACAAAAAACAAGAAGCTACTCGCAGTAACCCAATAAACAAAAGTTTTGAAAAGCCGGGTGCCAGCGAAAACTTAAAAAGCACTTTa    >  1:4847659/1‑100 (MQ=255)
                                                       tGGCAGATAAAAACAAAAAACAAGAAGCTACTCGTAGTAACCCAATAAACAAAAGTTTTGAAAAGCCGGGTGCCAGCGAAAACTTAAAAAGCACTTTATc  <  1:1395034/100‑1 (MQ=255)
                                                          ||                                                                                               
AAAAAAATCAAAAAACGCCTACTAGTGTAGACGTTGAATGGTGGTGAGAATTTTATGGTAGATAAAAACAAAAAACAAGAAACTACTCGTAGTAACCCATTAAACAAAAGTTTTGAAAAGTCAGGCGCCAGCGAAAAATTAAAAAGCACTTTATC  >  NC_007793/1589817‑1589971

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