Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 913,779 A→G 79.0% 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,7790AG79.0% 79.7 / 16.4 38L153P (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 (16/14);  total (19/19)
Fisher's exact test for biased strand distribution p-value = 6.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.04e-01

TGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTAAATTCTG  >  W3110S.gb/913728‑913839
                                                   |                                                            
tGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCAGACGCAGGCCGATAATGtt                                           <  1:3707/71‑1 (MQ=255)
tGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCAGACGCAGGCCGAGAATGtt                                           <  1:3521076/71‑1 (MQ=255)
tGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCAGACGCAGGCCGAGAATGtt                                           <  1:2935803/71‑1 (MQ=255)
    cgTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCAGACGCAGGCCGAGAATGTTTTCg                                       >  1:3186108/1‑71 (MQ=255)
    cgTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCAGACGCAGGCCGAGAATGTTTTCg                                       >  1:2482115/1‑71 (MQ=255)
    cgTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCAGACGCAGGCCGAGAATGTTTTCg                                       >  1:2247648/1‑71 (MQ=255)
               taGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCGGa                                                            >  1:857585/1‑39 (MQ=255)
               taGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCGGa                                                            >  1:588984/1‑39 (MQ=255)
               taGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCGGa                                                            >  1:1378630/1‑39 (MQ=255)
               taGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCGGa                                                            >  1:1656896/1‑39 (MQ=255)
               taGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCGGACGCAGGCCGAGAATGTTTTCGCCAATcgccg                             >  1:3467078/1‑70 (MQ=255)
                          aCCTTTCAGGCCATACAGGCACATCGGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGt                 <  1:2916516/71‑1 (MQ=255)
                             tttCAGGCCATACAGGCACAGCGGAGGCAGGCCGAGAATGTTTTCGCCaa                                   >  1:2937396/1‑50 (MQ=255)
                             tttCAGGCCATACAGGCACAGCGGACGCAGGCCGAGAAt                                              >  1:1433600/1‑39 (MQ=255)
                             tttCAGGCCATACAGGCACAGCGGACGCAGGCCGAGAATGTTTTCGCCaa                                   >  1:1102801/1‑50 (MQ=255)
                             tttCAGGCCATACAGGCACAGCGGACGCAGGCCGAGAATGTTTTCGCCaa                                   >  1:741960/1‑50 (MQ=255)
                             tttCAGGCCATACAGGCACAGCGGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAg              >  1:478507/1‑71 (MQ=255)
                             tttCAGGCCATACAGGCACAGCGGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAg              >  1:3678201/1‑71 (MQ=255)
                               tCAGGCCATACAGGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAg              <  1:3363921/69‑1 (MQ=255)
                                cAGGCCATATAGGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAg              <  1:2059332/68‑1 (MQ=255)
                                     cATACAGGCACAGCGGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTaaa       <  1:3322666/70‑1 (MQ=255)
                                       tACAGGCACAGCGGACGCAGGCCGAGAATGTTTTCGcc                                     <  1:3718281/38‑1 (MQ=255)
                                       tACAGGCACAGCGGACGCAGGCCGAGAATGTTTTCGcc                                     <  1:708857/38‑1 (MQ=255)
                                       tACAGGCACAGCGGACGCAGGCCGAGAATGTTTTCGcc                                     <  1:1411229/38‑1 (MQ=255)
                                       tACAGGCACAGCGGACGCAGGCCGAGAATGTTTTCGcc                                     <  1:3191143/38‑1 (MQ=255)
                                       tACAGGCACAGCGGACGCAGGCCGAGAATGTTTTCGcc                                     <  1:1546258/38‑1 (MQ=255)
                                       tACAGGCACAGCGGACGCAGGCCGAGAATGTTTTCGCCa                                    <  1:2042085/39‑1 (MQ=255)
                                       tACAGGCACAGCGGACGCAGGCCGAGAATGTTTTCGCCa                                    <  1:3255296/39‑1 (MQ=255)
                                       tACAGGCACAGCGGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGt                 <  1:173273/58‑1 (MQ=255)
                                        aCAGGCACAGCGGACGCAGGCCGAGAATGTTTTCGcc                                     <  1:456297/37‑1 (MQ=255)
                                        aCAGGCACAGCGGACGCAGGCCGAGAATGTTTTCGCCAATc                                 >  1:2656023/1‑41 (MQ=255)
                                        aCAGGCACAGCGGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTAAATTct   >  1:77921/1‑71 (MQ=255)
                                         cAGGCACAGCGGACGCAGGCCGAGAATGTTTTCGCCaa                                   >  1:2386527/1‑38 (MQ=255)
                                         cAGGCACAGCGGACGCAGGCCGAGAATGTTTTCGCCaa                                   >  1:412761/1‑38 (MQ=255)
                                         cAGGCACAGCGGACGCAGGCCGAGAATGTTTTCGCCaa                                   >  1:2125418/1‑38 (MQ=255)
                                         cAGGCACAGCGGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTAAATTctg  <  1:2240848/71‑1 (MQ=255)
                                         cAGGCACAGCGGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTAAATTctg  <  1:826729/71‑1 (MQ=255)
                                         cAGGCACAGCGGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTAAATTctg  <  1:844738/71‑1 (MQ=255)
                                               cagcggACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTa         <  1:2860377/58‑1 (MQ=255)
                                               cagcggACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTa         <  1:2004204/58‑1 (MQ=255)
                                                   |                                                            
TGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTAAATTCTG  >  W3110S.gb/913728‑913839

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