Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_007793 305,598 A→T D140V (GAT→GTT)  SAUSA300_RS01350 → iron‑sulfur cluster repair protein ScdA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_007793305,5980AT96.3% 71.4 / ‑3.7 27coding (419/675 nt)SAUSA300_RS01350iron‑sulfur cluster repair protein ScdA
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base T (25/1);  minor base . (1/0);  total (26/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 = 9.89e-01

GGACCTAATCATCCATATTTAGTTGAGTTAAAAGAAACATACGATACATTTAAAAATGGCATGTTAGAGCATATGCAAAAAGAAGACGATGTCGATTTTCCAAAACTCATTAAATATGAGCAAGGTGAGGTAGTAGACGATATTAATACTGTGATAGATGATTTAGTTTCAGACCACATTGC  >  NC_007793/305504‑305685
                                                                                              |                                                                                       
ggACCTAATCATCCATATTTAGTTGAGTTAAAAGAAACATACGATACATTTAAAAATGGCATGTTAGAGCATATGCAAAAAGAAGACGATGTCGTTTTTc                                                                                    >  1:824782/1‑100 (MQ=255)
ggACCTAATCATCCATATTTAGTTGAGTTAAAAGAAACATACGATACATTTAAAAATGGCATGTTAGAGCATATGCAAAAAGAAGACGATGTCGTTTTTc                                                                                    >  1:462638/1‑100 (MQ=255)
ggACCTAATCATCCATATTTAGTTGAGTTAAAAGAAACATACGATACATTTAAAAATGGCATGTTAGAGCATATGCAAAAAGAAGACGATGTCGTTTTTc                                                                                    >  1:412183/1‑100 (MQ=255)
 gACCTAATCATCCATATTTAGTTGAGTTAAAAGAAACATACGATACATTTAAAAATGGCATGTTAGAGCATATGCAAAAAGAAGACGATGTCGTTTTTcc                                                                                   >  1:1588830/1‑100 (MQ=255)
                   tAGTTGAGTTAAAAGAAACATACGATACATTTAAAAATGGCATGTTAGAGCATATGCAAAAAGAAGACGATGTCGTTTTTCCAAAACTCATTAAATATGa                                                                 >  1:1423180/1‑100 (MQ=255)
                        gAGTTAAAAGAAACATACGATACATTTAAAAATGGCATGTTAGAGCATATGCAAAAAGAAGACGATGTCGTTTTTCCAAAACTCATTAAATATGAGCAAg                                                            >  1:3879980/1‑100 (MQ=255)
                                                   aaaaaTGGCATGTTAGAGCATATGCAAAAAGAAGACGATGTCGTTTTTCCAAAACTCATTAAATATGAGCAAGGTGAGGTAGTAGACGATATTAATACtg                                 >  1:1403806/1‑100 (MQ=255)
                                                       aTGGCATGTTAGAGCATATGCAAAAAGAAGACGATGTCGTTTTTCCAAAACTCATTAAATATGAGCAAGGTGAGGTAGTAGACGATATTAATACTGTGAt                             <  1:2814975/100‑1 (MQ=255)
                                                        aggCATGTTAGAGCATATGCAAAAAGAAGACGATGTCGTTTTTCCAAAACTCATTAAATATGAGCAAGGTGAGGTAGTAGACGATATTAATACTGTGATa                            >  1:3569799/2‑100 (MQ=255)
                                                         ggCATGTTAGAGCATATGCAAAAAGAAGACGATGTCGTTTTTCCAAAACTCATTAAATATGAGCAAGGTGAGGTAGTAGACGATATTAATACTGTGATAg                           >  1:1133214/1‑100 (MQ=255)
                                                         ggCATGTTAGAGCATATGCAAAAAGAAGACGATGTCGTTTTTCCAAAACTCATTAAATATGAGCAAGGTGAGGTAGTAGACGATATTAATACTGTGATAg                           >  1:752479/1‑100 (MQ=255)
                                                         ggCATGTTAGAGCATATGCAAAAAGAAGACGATGTCGTTTTTCCAAAACTCATTAAATATGAGCAAGGTGAGGTAGTAGACGATATTAATACTGTGATAg                           >  1:4060220/1‑100 (MQ=255)
                                                         ggCATGTTAGAGCATATGCAAAAAGAAGACGATGTCGTTTTTCCAAAACTCATTAAATATGAGCAAGGTGAGGTAGTAGACGATATTAATACTGTGATAg                           >  1:3500733/1‑100 (MQ=255)
                                                         ggCATGTTAGAGCATATGCAAAAAGAAGACGATGTCGTTTTTCCAAAACTCATTAAATATGAGCAAGGTGAGGTAGTAGACGATATTAATACTGTGATAg                           >  1:2913518/1‑100 (MQ=255)
                                                         ggCATGTTAGAGCATATGCAAAAAGAAGACGATGTCGTTTTTCCAAAACTCATTAAATATGAGCAAGGTGAGGTAGTAGACGATATTAATACTGTGATAg                           >  1:2269202/1‑100 (MQ=255)
                                                         ggCATGTTAGAGCATATGCAAAAAGAAGACGATGTCGTTTTTCCAAAACTCATTAAATATGAGCAAGGTGAGGTAGTAGACGATATTAATACTGTGATAg                           >  1:1456513/1‑100 (MQ=255)
                                                         ggCATGTTAGAGCATATGCAAAAAGAAGACGATGTCGTTTTTCCAAAACTCATTAAATATGAGCAAGGTGAGGTAGTAGACGATATTAATACTGTGATAg                           >  1:1277957/1‑100 (MQ=255)
                                                         ggCATGTTAGAGCATATGCAAAAAGAAGACGATGTCGTTTTTCCAAAACTCATTAAATATGAGCAAGGTGAGGTAGTAGAAGATATTAATACTGTGATAg                           >  1:3721534/1‑100 (MQ=255)
                                                         ggCATGTTAGAGCATATGCAAAAAGAAGACGATGTCGTTTTTCCAAAACTCATTAAATATGAGCAAGGTGAGGCAGTAGACGATATTAATACCGTGATAg                           >  1:572827/1‑100 (MQ=255)
                                                                tAGAGCATATGCAAAAAGAAGACGATGTCGTTTTTCCAAAACTCATTAAATATGAGCAAGGTGAGGTAGTAGACGATATTAATACTGTGATAGATGAttt                    >  1:3182491/1‑100 (MQ=255)
                                                                   agCATATGCAAAAAGAAGACGATGTCGTTTTTCCAAAACTCATTAAATATGAGCAAGGTGAGGTAGTAGACGATATTAATACCGTGATAgatcggaagag                 >  1:1886366/1‑92 (MQ=255)
                                                                     cATATGCAAAAAGAAGACGATGTCGTTTTTCCAAAACTCATTAAATATGAGCAAGGTGAGGTAGTAGACGATATTAATACTGTGATAGATGATTTAGttt               >  1:2320187/1‑100 (MQ=255)
                                                                     cATATGCAAAAAGAAGACGATGTCGTTTTTCCAAAACTCATTAAATATGAGCAAGGTGAGGTAGTAGACGATATTAATACTGTGATAGATGATTTAGttt               >  1:1055204/1‑100 (MQ=255)
                                                                        atGCAAAAAGAAGACGATGTCGTTTTTCCAAAACTCATTAAATATGAGCAAGGTGAGGTAGTAGACGATATTAATACTGTGATAGATGATTTAGTTTCAg            >  1:325131/1‑100 (MQ=255)
                                                                        atGCAAAAAGAAGACGATGTCGTTTTTCCAAAACTCATTAAATATGAGCAAGGTGAGGTAGTAGACGATATTAATACTGTGATAGATGATTTAGTTTCAg            >  1:2887717/1‑100 (MQ=255)
                                                                        atGCAAAAAGAAGACGATGTCG‑TTTTCCAAAACTCATTAAATATGAGCAAGGTGAGGTAGTAGACGATATTAATACTGTGATAGATGATTTAGTTTCAGa           >  1:4473389/1‑100 (MQ=255)
                                                                                  aagaCGATGTCGTTTTTCCAAAACTCATTAAATATGAGCAAGGTGAGGTAGTAGACGATATTAATACTGTGATAGATGATTTAGTTTCAGACCACATTGc  >  1:4702274/1‑100 (MQ=255)
                                                                                              |                                                                                       
GGACCTAATCATCCATATTTAGTTGAGTTAAAAGAAACATACGATACATTTAAAAATGGCATGTTAGAGCATATGCAAAAAGAAGACGATGTCGATTTTCCAAAACTCATTAAATATGAGCAAGGTGAGGTAGTAGACGATATTAATACTGTGATAGATGATTTAGTTTCAGACCACATTGC  >  NC_007793/305504‑305685

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