Predicted mutation
evidence seq id position mutation annotation gene description
RA W3110S.gb 2,698,480 C→T C54C (TGC→TGT yfhL → predicted 4Fe‑4S cluster‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,698,4800CT76.5% 90.0 / 26.8 51C54C (TGC→TGTyfhLpredicted 4Fe‑4S cluster‑containing protein
Reads supporting (aligned to +/- strand):  ref base C (0/12);  new base T (39/0);  total (39/12)
Fisher's exact test for biased strand distribution p-value = 6.30e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.53e-01

AGCGATAAGTGTACCGAATGCGTAGGGCACTACGAGACACCAACCTGCCAGAAGGTGTGCCCGATCCCCAATACTATTGTGAAAGATCCGGCGCATGTCGAGACAG  >  W3110S.gb/2698421‑2698526
                                                           |                                              
aGCGATAAGTGTACCGAATGCGTAGGGCACTACGAGACACCAACCTGCCAGAAGGTGTGccc                                              <  1:348644/62‑1 (MQ=255)
aGCGATAAGTGTACCGAATGCGTAGGGCACTACGAGACACCAACCTGCCAGAAGGTGTGccc                                              <  1:1701488/62‑1 (MQ=255)
                 aTGCGTAGGGCACTACGAGACACCAACCTGCCAGAAGGTGTGCCCGATCCCCAATACTATtg                             <  1:2509271/62‑1 (MQ=255)
                 aTGCGTAGGGCACTACGAGACACCAACCTGCCAGAAGGTGTGCCCGATCCCCAATACTATtg                             <  1:2679822/62‑1 (MQ=255)
                                            cTGCCAGAAGGTGTGCCCGATCCCCAATACTTTTGTGAAAGATCCGGCGCATGTCGATACag  <  1:2411279/62‑1 (MQ=255)
                                            cTGCCAGAAGGTGTGCCCGATCCCCAATACTATTGTGAAAGATCCGGCGCATGTCGAGACag  <  1:197124/62‑1 (MQ=255)
                                            cTGCCAGAAGGTGTGCCCGATCCCCAATACTATTGTGAAAGATCCGGCGCATGTCGAGACag  <  1:219946/62‑1 (MQ=255)
                                            cTGCCAGAAGGTGTGCCCGATCCCCAATACTATTGTGAAAGATCCGGCGCATGTCGAGACag  <  1:3022708/62‑1 (MQ=255)
                                            cTGCCAGAAGGTGTGCCCGATCCCCAATACTATTGTGAAAGATCCGGCGCATGTCGAGACag  <  1:2048980/62‑1 (MQ=255)
                                            cTGCCAGAAGGTGTGCCCGATCCCCAATACTATTGTGAAAGATCCGGCGCATGTCGAGACag  <  1:34262/62‑1 (MQ=255)
                                            cTGCCAGAAGGTGTGCCCGATCCCCAATACTATTGTGAAAGATCCGGCGCATGTCGAGACag  <  1:597034/62‑1 (MQ=255)
                                            cTGCCAGAAGGTGTGCCCGATCCCCAATACTATTGTGAAAGATCCGGCGCATGTCGAGACag  <  1:600722/62‑1 (MQ=255)
                                            cTGCCAGAAGGTGTGCCCGATCCCCAATACTATTGTGAAAGATCCGGCGCATGTCGAGACag  <  1:670278/62‑1 (MQ=255)
                                            cTGCCAGAAGGTGTGCCCGATCCCCAATACTATTGTGAAAGATCCGGCGCATGTCGAGACag  <  1:1371545/62‑1 (MQ=255)
                                                cAGAAGGTGTGTCCGATCCCCAATACTATTGTGaaa                        >  1:822378/1‑36 (MQ=255)
                                                cAGAAGGTGTGTCCGATCCCCAATACTATTGTGAAAg                       >  1:2369862/1‑37 (MQ=255)
                                                cAGAAGGTGTGTCCGATCCCCAATACTATTGTGAAAg                       >  1:68307/1‑37 (MQ=255)
                                                cAGAAGGTGTGTCCGATCCCCAATACTATTGTGAAAg                       >  1:700198/1‑37 (MQ=255)
                                                cAGAAGGTGTGTCCGATCCCCAATACTATTGTGAAAg                       >  1:791078/1‑37 (MQ=255)
                                                cAGAAGGTGTGTCCGATCCCCAATACTATTGTGAAAg                       >  1:3062284/1‑37 (MQ=255)
                                                cAGAAGGTGTGTCCGATCCCCAATACTATTGTGAAAg                       >  1:608325/1‑37 (MQ=255)
                                                cAGAAGGTGTGTCCGATCCCCAATACTATTGTGAAAg                       >  1:2706930/1‑37 (MQ=255)
                                                cAGAAGGTGTGTCCGATCCCCAATACTATTGTGAAAg                       >  1:830705/1‑37 (MQ=255)
                                                cAGAAGGTGTGTCCGATCCCCAATACTATTGTGAAAg                       >  1:2586527/1‑37 (MQ=255)
                                                cAGAAGGTGTGTCCGATCCCCAATACTATTGTGAAAg                       >  1:863909/1‑37 (MQ=255)
                                                cAGAAGGTGTGTCCGATCCCCAATACTATTGTGAAAg                       >  1:2560921/1‑37 (MQ=255)
                                                cAGAAGGTGTGTCCGATCCCCAATACTATTGTGAAAg                       >  1:880660/1‑37 (MQ=255)
                                                cAGAAGGTGTGTCCGATCCCCAATACTATTGTGAAAg                       >  1:2429225/1‑37 (MQ=255)
                                                cAGAAGGTGTGTCCGATCCCCAATACTATTGTGAAAg                       >  1:96001/1‑37 (MQ=255)
                                                cAGAAGGTGTGTCCGATCCCCAATACTATTGTGAAAg                       >  1:2011040/1‑37 (MQ=255)
                                                cAGAAGGTGTGTCCGATCCCCAATACTATTGTGAAAg                       >  1:1101441/1‑37 (MQ=255)
                                                cAGAAGGTGTGTCCGATCCCCAATACTATTGTGAAAg                       >  1:1104939/1‑37 (MQ=255)
                                                cAGAAGGTGTGTCCGATCCCCAATACTATTGTGAAAg                       >  1:1195300/1‑37 (MQ=255)
                                                cAGAAGGTGTGTCCGATCCCCAATACTATTGTGAAAg                       >  1:1239633/1‑37 (MQ=255)
                                                cAGAAGGTGTGTCCGATCCCCAATACTATTGTGAAAg                       >  1:1269728/1‑37 (MQ=255)
                                                cAGAAGGTGTGTCCGATCCCCAATACTATTGTGAAAg                       >  1:1298922/1‑37 (MQ=255)
                                                cAGAAGGTGTGTCCGATCCCCAATACTATTGTGAAAg                       >  1:135691/1‑37 (MQ=255)
                                                cAGAAGGTGTGTCCGATCCCCAATACTATTGTGAAAg                       >  1:1366821/1‑37 (MQ=255)
                                                cAGAAGGTGTGTCCGATCCCCAATACTATTGTGAAAg                       >  1:1532568/1‑37 (MQ=255)
                                                cAGAAGGTGTGTCCGATCCCCAATACTATTGTGAAAg                       >  1:1855622/1‑37 (MQ=255)
                                                cAGAAGGTGTGTCCGATCCCCAATACTATTGTGAAAg                       >  1:1931532/1‑37 (MQ=255)
                                                cAGAAGGTGTGTCCGATCCCCAATACTATTGTGAAAg                       >  1:2384206/1‑37 (MQ=255)
                                                cAGAAGGTGTGTCCGATCCCCAATACTATTGTGAAAg                       >  1:2035755/1‑37 (MQ=255)
                                                cAGAAGGTGTGTCCGATCCCCAATACTATTGTGAAAg                       >  1:2046788/1‑37 (MQ=255)
                                                cAGAAGGTGTGTCCGATCCCCAATACTATTGTGAAAg                       >  1:2049185/1‑37 (MQ=255)
                                                cAGAAGGTGTGTCCGATCCCCAATACTATTGTGAAAg                       >  1:2108928/1‑37 (MQ=255)
                                                cAGAAGGTGTGTCCGATCCCCAATACTATTGTGAAAg                       >  1:215235/1‑37 (MQ=255)
                                                cAGAAGGTGTGTCCGATCCCCAATACTATTGTGAAAg                       >  1:2187651/1‑37 (MQ=255)
                                                cAGAAGGTGTGTCCGATCCCCAATACTATTGTGAAAg                       >  1:2191363/1‑37 (MQ=255)
                                                cAGAAGGTGTGTCCGATCCCCAATACTATTGTGAAAg                       >  1:2283148/1‑37 (MQ=255)
                                                cAGAAGGTGTGTCCGATCCCCAATACTATTGTGAAAg                       >  1:103824/1‑37 (MQ=255)
                                                cAGAAGGTGTGTCCGATCCCCAATAATATTGTGAAAg                       >  1:2952724/1‑37 (MQ=37)
                                                cAGAAGGAGTGTCCGATCCCCAATACTATTGTGaaa                        >  1:2689998/1‑36 (MQ=38)
                                                           |                                              
AGCGATAAGTGTACCGAATGCGTAGGGCACTACGAGACACCAACCTGCCAGAAGGTGTGCCCGATCCCCAATACTATTGTGAAAGATCCGGCGCATGTCGAGACAG  >  W3110S.gb/2698421‑2698526

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