Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_007793 1,193,217 (T)7→8 intergenic (+68/‑331) SAUSA300_RS05910 → / → SAUSA300_RS05915 hypothetical protein/bifunctional protein PyrR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0077931,193,2101.T88.9% 138.4 / 10.7 55intergenic (+61/‑338)SAUSA300_RS05910/SAUSA300_RS05915hypothetical protein/bifunctional protein PyrR
Reads supporting (aligned to +/- strand):  ref base . (6/0);  new base T (47/1);  total (54/1)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.97e-01

AGATACAATTCGAAAAAGAGATGCATAATTGTGTCGTGCTATAATTACGATAACATTATTATGTAAACTATAAGTGTTTATCGTTTTCG‑‑T‑TTTTTTAAATGAAAACACATTATAAAATAGATTAAGGTTAACATTTTTACTGAGTATGGTTAATTTTTCAATAT  >  NC_007793/1193122‑1193285
                                                                                         |                                                                             
aGATACAATTCGAAAAAGAGATGCATAATTGTGTCGTGCTATAATTACGATAACATTATTATGTAAACTATAAGTGTTTATCGTTTTCG‑‑T‑TTTTTTtaaa                                                                  >  1:5137781/1‑99 (MQ=255)
aGATACAATTCGAAAAAGAGATGCATAATTGTGTCGTGCTATAATTACGATAACATTATTATGTAAACTATAAGTGTTTATCGTTTTCG‑‑T‑TTTTTTtaaa                                                                  >  1:436563/1‑99 (MQ=255)
aGATACAATTCGAAAAAGAGATGCATAATTGTGTCGTGCTATAATTACGATAACATTATTATGTAAACTATAAGTGTTTATCGTTTTCG‑‑T‑TTTTTTtaaa                                                                  >  1:2113638/1‑99 (MQ=255)
  aTACAATTCGAAAAAGAGATGCATAATTGTGTCGTGCTATAATTACGATAACATTATTATGTAAACTATAAGTGTTTATCGTTTTCGT‑T‑TTTTTTAAATg                                                                 >  1:2341241/1‑100 (MQ=255)
    aCAATTCGAAAAAGAGATGCATAATTGTGTCGTGCTATAATTACGATAACATTATTATGTAAACTATAAGTGTTTATCGTTTTCGT‑T‑TTTTTTAAATGaa                                                               >  1:1832849/1‑100 (MQ=255)
    aCAATTCGAAAAAGAGATGCATAATTGTGTCGTGCTATAATTACGATAACATTATTATGTAAACTATAAGTGTTTATCGTTTTCGT‑T‑TTTTTTAAATGaa                                                               >  1:3340835/1‑100 (MQ=255)
    aCAATTCGAAAAAGAGATGCATAATTGTGTCGTGCTATAATTACGATAACATTATTATGTAAACTATAAGTGTTTATCGTTTTCG‑‑t‑ttttttttaaatga                                                              >  1:6247355/1‑92 (MQ=255)
    aCAATTCGAAAAAGAGATGCATAATTGTGTCGTGCTATAATTACGATAACATTATTATGTAAACTATAAGTGTCTATCGTTTTCGT‑T‑TTTTTTAAATGaa                                                               >  1:6058599/1‑100 (MQ=255)
                tgagaTGCATAATTGTGTCGTGCTATAATTACGATAACATTATTATGTAAACTATAAGTGTTTATCGTTTTCGT‑T‑TTTTTTAAATGAAAACACATTATaa                                                   >  1:1359664/2‑100 (MQ=255)
                 ttgaTGCATAATTGTGTCGTGCTATAATTACGATAACATTATTATGTAAACTATAAGTGTTTATCGTTTTCGT‑T‑TTTTTTAAATGAAAACACATTATaaa                                                  >  1:84986/3‑100 (MQ=255)
                 atgaTGCATAATTGTGTCGTGCTATAATTACGATAACCTTATTATGTAAACTATAAGTGTTTATCGTTTTCG‑‑T‑TTTTTTTAATGAAAACACATTATaaaa                                                 >  1:3575359/3‑100 (MQ=255)
                   tgtGCATAATTGTGTCGTGCTATAATTACGATAACATTATTATGTAAACTATAAGTGTTTATCGTTTTCG‑‑T‑TTTTTTTAATGAAAACACATTATAAAATa                                               >  1:4231389/3‑100 (MQ=255)
                   gaTGCATAATTGTGTCGTGCTATAATTACGATAACATTATTACGTAAACTATAAGCGTTTATCGTTTTCGT‑T‑TTTTTTAAATGAAAACACATTATAAAAt                                                >  1:499049/1‑100 (MQ=255)
                   ctcgCATAATTGTGTCGTGCTATAATTACGATAACATTATTATGTAAACTATAAGTGTTTATCGTTTTCGT‑T‑TTTTTTAAATGAAAACACATTATAAAAt                                                >  1:1883424/4‑100 (MQ=255)
                   aaTGCATAATTGTGTCGTGCTATAATTACGATAACATTATTATGTAAACTATAAGTGTTTATCGTTTTCGT‑T‑TTTTTTAAATGAAAACACATTATAAAAt                                                >  1:5434665/2‑100 (MQ=255)
                    ttGCATAATTGTGTCGTGCTATAATTACGATAACATTATTATGTAAACTATAAGTGTTTATCGTTTTCGT‑T‑TTTTTTAAATGAAAACACATTATAAAATa                                               >  1:448819/2‑100 (MQ=255)
                    aTGCATAATTGTGTCGTGCTATAATTACGATAACATTATTATGTAAACTATAAGTGTTTATCGTTTTCGT‑T‑TTTTTTAAATGAAAACACATTATAAAATa                                               >  1:2201614/1‑100 (MQ=255)
                    aTGCATAATTGTGTCGTGCTATAATTACGATAACATTATTATGTAAACTATAAGTGTTTATCGTTTTCGT‑T‑TTTTTTAAATGAAAACACATTATAAAATa                                               >  1:3921742/1‑100 (MQ=255)
                    aTGCATAATTGTGTCGTGCTATAATTACGATAACATTATTATGTAAACTATAAGTGTTTATCGTTTTCGT‑T‑TTTTTTAAATGAAAACACATTATAAAATa                                               >  1:3987794/1‑100 (MQ=255)
                    aTGCATAATTGTGTCGTGCTATAATTACGATAACATTATTATGTAAACTATAAGTGTTTATCGTTTTCGT‑T‑TTTTTTAAATGAAAACACATTATAAAATa                                               >  1:5596977/1‑100 (MQ=255)
                    aTGCATAATTGTGTCGTGCTATAATTACGATAACATTATTATGTAAACTATAAGTGTTTATCGTTTTCGT‑T‑TTTTTTAAATGAAAACACATTATAAAATa                                               >  1:2884944/1‑100 (MQ=255)
                    aTGCATAATTGTGTCGTGCTATAATTACGATAACATTATTATGTAAACTATAAGTGTTTATCGTTTTCGT‑T‑TTTTTTAAATGAAAACACATTATAAAATa                                               >  1:1402314/1‑100 (MQ=255)
                    aTGCATAATTGTGTCGTGCTATAATTACGATAACATTATTATGTAAACTATAAGTGTTTATCGTTTTCG‑‑TCTTTTTTAAATGAAACCACATTATAAAATa                                               >  1:4918731/1‑100 (MQ=255)
                      gCATAATTGTGTCGTGCTATAATTACGATAACATTATTATGTAAACTATAAGTGTTTATCGTTTTCGT‑T‑TTTTTTAAATGAAAACACATTATAATATAGa                                             >  1:3792786/1‑100 (MQ=255)
                      gCATAATTGTGTCGTGCTATAATTACGATAACATTATTATGTAAACTATAAGTGTTTATCGTTTTCGT‑T‑TTTTTTAAATGAAAACACATTATAAAATAgg                                             >  1:4569324/1‑99 (MQ=255)
                      gCATAATTGTGTCGTGCTATAATTACGATAACATTATTATGTAAACTATAAGTGTTTATCGTTTTCGT‑T‑TTTTTTAAATGAAAACACATTATAAAATAGa                                             >  1:4301560/1‑100 (MQ=255)
                      gCATAATTGTGTCGTGCTATAATTACGATAACATTATTATGTAAACTATAAGTGTTTATCGTTTTCGT‑T‑TTTTTTAAATGAAAACACATTATAAAATAGa                                             >  1:2215524/1‑100 (MQ=255)
                      gCATAATTGTGTCGTGCTATAATTACGATAACATTATTATGTAAACTATAAGTGTTTATCGTTTTCGT‑T‑TTTTTTAAATGAAAACACATTATAAAATAGa                                             >  1:1603627/1‑100 (MQ=255)
                      gCATAATTGTGTCGTGCTATAATTACGATAACATTATTATGTAAACTATAAGTGTTTATCGTTTTCGT‑T‑TTTTTTAAATGAAAACACATTATAAAATAGa                                             >  1:4589393/1‑100 (MQ=255)
                      gCATAATTGTGTCGTGCTATAATTACGATAACATTATTATGTAAACTATAAGTGTTTATCGTTTTCGT‑T‑TTTTTTAAATGAAAACACATTATAAAATAGa                                             >  1:3656461/1‑100 (MQ=255)
                      gCATAATTGTGTCGTGCTATAATTACGATAACATTATTATGTAAACTATAAGTGTTTATCGTTTTCGT‑T‑TTTTTTAAATGAAAACACATTATAAAATAGa                                             >  1:3553160/1‑100 (MQ=255)
                      gCATAATTGTGTCGTGCTATAATTACGATAACATTATTATGTAAACTATAAGTGTTTATCGTTTTCGT‑T‑TTTTTTAAATGAAAACACATTATAAAATAGa                                             >  1:5975996/1‑100 (MQ=255)
                      gCATAATTGTGTCGTGCTATAATTACGATAACATTATTATGTAAACTATAAGTGTTTATCGTTTTCGT‑T‑TTTTTTAAATGAAAACACATTATAAAATAGa                                             >  1:3322941/1‑100 (MQ=255)
                      gCATAATTGTGTCGTGCTATAATTACGATAACATTATTATGTAAACTATAAGTGTTTATCGTTTTCGT‑T‑TTTTTTAAATGAAAACACATTATAAAATAGa                                             >  1:1469641/1‑100 (MQ=255)
                      gCATAATTGTGTCGTGCTATAATTACGATAACATTATTATGTAAACTATAAGTGTTTATCGTTTTCGT‑T‑TTTTTTAAATGAAAACACATTATAAAATAGa                                             >  1:3221597/1‑100 (MQ=255)
                      gCATAATTGTGTCGTGCTATAATTACGATAACATTATTATGTAAACTATAAGTGTTTATCGTTTTCGT‑T‑TTTTTTAAATGAAAACACATTATAAAATAGa                                             >  1:2794819/1‑100 (MQ=255)
                      gCATAATTGAGTCGTGCTATAATTACGATAACATTATTATGTAAACTATAAGTGTTTATCGTTTTCGT‑T‑TTTTTTAAATGAAAACACATTATAAAATAGa                                             >  1:3630723/1‑100 (MQ=255)
                        agccttGTGTCGTGCTATAATTACGATAACATTATTATGTAAACTATAAGTGTTTATCGTTTTCGT‑T‑TTTTTTAAATGAAAACACATTATAAAATAGAtt                                           >  1:4655422/5‑100 (MQ=255)
                          aaTTGTGTCGTGCTATAATTACGATAACATTATTATGTAAACTATAAGTGTTTATCGTTTTCGT‑T‑TTTTTTAAATGAAAACACATTATAAAATAGATTaa                                         >  1:631042/1‑100 (MQ=255)
                           gatgtgTCGTGCTATAATTACGATAACATTATTATGTAAACTATAAGTGTTTATCGTTTTCGT‑T‑TTTTTTAAATGAAAACACATTATAAAATAGATTAAg                                        >  1:1747928/3‑100 (MQ=255)
                                aacGTGCTATAATTACGATAACATTATTATGTAAACTATAAGTGTTTATCGTTTTCGT‑T‑TTTTTTAAATGAAAACACATTATAAAATAGATTAAGGTTaa                                   >  1:3251017/3‑100 (MQ=255)
                                 ccGTGCTATAATTACGATAACATTATTATGTAAACTATAAGTGTTTATCGTTTTCGT‑T‑TTTTTTAAATGAAAACACATTATAAAATAGATTAAGGTTAAc                                  >  1:188092/2‑100 (MQ=255)
                                 acGTGCTATAATTACGATAACATTATTATGTAAACTATAAGTGTTTATCGTTTTCGT‑T‑TTTTTTAAATGAAAACACATTATAAAATAGATTAAGGTTAAc                                  >  1:4608610/2‑100 (MQ=255)
                                 acGTGCTATAATTACGATAACATTATTATGTAAACTATAAGTGTTTATCGTTTTCGT‑T‑TTTTTTAAATGAAAACACATTATAAAATAGATTAAGGTTAAc                                  >  1:4400068/2‑100 (MQ=255)
                                  cGTGCTATAATTACGATAACATTATTATGTAAACTATAAGTGTTTATCGTTTTCGT‑T‑TTTTTTCAATGAAAACACATTATAAAATAGATTAAGGTTAACa                                 >  1:591989/1‑100 (MQ=255)
                                  cGTGCTATAATTACGATAACATTATTATGTAAACTATAAGTGTTTATCGTTTTCGT‑T‑TTTTTTAAATGAAAACACATTATAAAATAGATTAAGGTTAACa                                 >  1:5009630/1‑100 (MQ=255)
                                  cGTGCTATAATTACGATAACATTATTATGTAAACTATAAGTGTTTATCGTTTTCGT‑T‑TTTTTTAAATGAAAACACATTATAAAATAGATTAAGGTTAACa                                 >  1:4220188/1‑100 (MQ=255)
                                  cGTGCTATAATTACGATAACATTATTATGTAAACTATAAGTGTTTATCGTTTTCGT‑T‑TTTTTTAAATGAAAACACATTATAAAATAGATTAAGGTTAACa                                 >  1:3656282/1‑100 (MQ=255)
                                  cGTGCTATAATTACGATAACATTATTATGTAAACTATAAGTGTTTATCGTTTTCGT‑T‑TTTTTTAAATGAAAACACATTATAAAATAGATTAAGGTTAACa                                 >  1:1108415/1‑100 (MQ=255)
                                    tGCTATAATTACGATAACATTATTATGTACACTATAAGTGTTTATCGTTTTCGC‑T‑TTTTTTAAATGAAAACACATTATAAAATAGATTAAGGTTAACAtt                               >  1:3493840/1‑100 (MQ=255)
                                    ggCTATAATTACGATAACATTATTATGTAAACTATAAGTGTTTATCGTTTTCGT‑T‑TTTTTTAAATGAAAACACATTATAAAATAGATTAAGGTTAACAtt                               >  1:776646/2‑100 (MQ=255)
                                                       atattatGTAAACTATAAGTGTTTATCGTTTTCGT‑T‑TTTTTTAAATGAAAACACATTATAAAATAGATTAAGGTTAACATTTTTACTGAGTATGGTTAAt            >  1:4602672/2‑100 (MQ=255)
                                                         tttatGTAAACTATAAGTGTTTATCGTTTTCGT‑T‑TTTTTTAAATGAAAACACATTATAAAATAGATTAAGGTTAACATTTTTACTGAGTATGGTTAAttt          >  1:4574465/2‑100 (MQ=255)
                                                         attatGTAAACTATAAGTGTTTATCGTTTTCGTTT‑TTTTTTAAATGAAAACACATTATAAAATAGATTAAGGTTAACATTTTTACTGAGTATGGTTAAtt           >  1:2503059/1‑100 (MQ=255)
                                                                aaaCTATAAGTGTTTATCGTTTTCGTTT‑TTTTTTAAATGAAAACACATTATAAAATAGATTAAGGTTAACATTTTTACTGAGTATGGTTAATTTTTCAat    >  1:2452375/1‑100 (MQ=255)
                                                                 aaCTATAAGTGTTTATCGTTTTCGT‑T‑TTTTTTAAATGAAAACACATTATAAAATAGATTAAGGTTAACATTTTTACTGAGTATGGTTAATTTTTCAatat  <  1:5678432/100‑1 (MQ=255)
                                                                                         |                                                                             
AGATACAATTCGAAAAAGAGATGCATAATTGTGTCGTGCTATAATTACGATAACATTATTATGTAAACTATAAGTGTTTATCGTTTTCG‑‑T‑TTTTTTAAATGAAAACACATTATAAAATAGATTAAGGTTAACATTTTTACTGAGTATGGTTAATTTTTCAATAT  >  NC_007793/1193122‑1193285

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