Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 913,779 A→G 72.4% L153P (CTG→CCG)  hcp ← hybrid‑cluster [4Fe‑2S‑2O] protein in anaerobic terminal reductases

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb913,7790AG72.4% 46.1 / 17.9 29L153P (CTG→CCG) hcphybrid‑cluster [4Fe‑2S‑2O] protein in anaerobic terminal reductases
Reads supporting (aligned to +/- strand):  ref base A (3/5);  new base G (2/19);  total (5/24)
Fisher's exact test for biased strand distribution p-value = 1.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

AGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTAAATTCTG  >  W3110S.gb/913722‑913839
                                                         |                                                            
agaACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCGGACGCAGGCCgag                                                 <  1:788137/71‑1 (MQ=255)
agaACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCGGACGCAGGCCgag                                                 <  1:189332/71‑1 (MQ=255)
agaACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCGGACGCAGGCCgag                                                 <  1:1597993/71‑1 (MQ=255)
agaACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCGGACGCAGGCCgag                                                 <  1:1297050/71‑1 (MQ=255)
agaACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCGGACGCAGGCCgag                                                 <  1:1485155/71‑1 (MQ=255)
 gaACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCGGACGCAGGCCgag                                                 <  1:1343686/70‑1 (MQ=255)
              ttCCATATAGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCGGACGCAGGCCGAGAATGTTTTCGCCa                                    <  1:1129203/70‑1 (MQ=255)
              ttCCATATAGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCGGACGCAGGCCGAGAATGTTTTCGCCa                                    <  1:585987/70‑1 (MQ=255)
              ttCCATATAGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCGGACGCAGGCCGAGAATGTTTTCGCCa                                    <  1:1540338/70‑1 (MQ=255)
              ttCCATATAGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCGGACGCAGGCCGAGAATGTTTTCGCCa                                    <  1:1329742/70‑1 (MQ=255)
              ttCCATATAGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCGGACGCAGGCCGAGAATGTTTTCGCCa                                    <  1:1327910/70‑1 (MQ=255)
                     taGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATcgccgc                            <  1:1237130/71‑1 (MQ=255)
                     taGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATcgccgc                            <  1:443578/71‑1 (MQ=255)
                     taGGCCGCCGCACCTTTCAGGCCATACAGGCACACCAGACGCAGGCCGAGAATGTTTTCGCCAATcgccgc                            <  1:876373/71‑1 (MQ=255)
                             cgcACCTTTCAGGCCATACAGGCACAGCGGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCtt                    <  1:889314/71‑1 (MQ=255)
                             cgcACCTTTCAGGCCATACAGGCACAGCGGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCtt                    <  1:470669/71‑1 (MQ=255)
                                aCCTTTCAGGCCATACAGGCACAGCGGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGt                 <  1:1951560/71‑1 (MQ=255)
                                             tACAGGCACAGCGGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGt                 <  1:1347224/58‑1 (MQ=255)
                                             tACAGGCACAGCGGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGt                 <  1:1845424/58‑1 (MQ=255)
                                             tACAGGCACAGCGGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGt                 <  1:816953/58‑1 (MQ=255)
                                             tACAGGCACAGCGGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGt                 <  1:863627/58‑1 (MQ=255)
                                             tACAGGCACAGCGGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGt                 <  1:1520162/58‑1 (MQ=255)
                                             tACAGGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGt                 <  1:90015/58‑1 (MQ=255)
                                              aCAGGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGt                 <  1:179575/57‑1 (MQ=255)
                                               cAGGCACAGCGGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTAAATTctg  >  1:128188/1‑71 (MQ=255)
                                               cAGGCACAGCGGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTAAATTctg  >  1:1088730/1‑71 (MQ=255)
                                               cAGGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTAAATTctg  >  1:1932158/1‑71 (MQ=255)
                                               cAGGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTAAATTctg  >  1:1295174/1‑71 (MQ=255)
                                               cAGGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTAAATTctg  >  1:1530800/1‑71 (MQ=255)
                                                         |                                                            
AGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTAAATTCTG  >  W3110S.gb/913722‑913839

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