Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 913,779 A→G 68.1% 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,7790AG68.1% 59.7 / 38.7 47L153P (CTG→CCG) hcphybrid‑cluster [4Fe‑2S‑2O] protein in anaerobic terminal reductases
Reads supporting (aligned to +/- strand):  ref base A (10/5);  new base G (18/14);  total (28/19)
Fisher's exact test for biased strand distribution p-value = 5.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.49e-01

ACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTAAATTC  >  W3110S.gb/913718‑913837
                                                             |                                                          
aCCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCAGACGCAg                                                     >  1:481889/1‑69 (MQ=255)
aCCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCAGACGCAg                                                     >  1:485843/1‑69 (MQ=255)
aCCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCAGACGCAg                                                     >  1:1867361/1‑69 (MQ=255)
aCCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCAGACGCAg                                                     >  1:1909547/1‑69 (MQ=255)
aCCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCAGACGCAg                                                     >  1:3395531/1‑69 (MQ=255)
      aaCATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCGGACGCa                                                      >  1:2927200/1‑62 (MQ=255)
      aaCATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCAGACGCAGGCCgag                                               <  1:3403422/69‑1 (MQ=255)
         aTGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCAGACGCAGGCCGAGAAt                                            >  1:3297211/1‑69 (MQ=255)
         aTGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCAGACGCAGGCCGAGAAt                                            >  1:2492769/1‑69 (MQ=255)
            cgcgTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCGGACGCAGGCCGAGAATGtt                                         >  1:707345/1‑69 (MQ=255)
            cgcgTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCGGACGCAGGCCGAGAATGtt                                         >  1:3022500/1‑69 (MQ=255)
            cgcgTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCGGACGCAGGCCGAGAATGtt                                         >  1:2976308/1‑69 (MQ=255)
                  ttCCATATAGGCCGCCTCACCTTTCAGGCCATACAGGCACAGCGGATGCAGGCCGAGAATGTTTTCGcc                                   >  1:22534/1‑69 (MQ=255)
                  ttCCATATAGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCGGACGCAGGCCGAGAATGTTTTCGcc                                   >  1:1924822/1‑69 (MQ=255)
                   tCCATATAGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCGGACGCAGGCCGAGAATGTTTTCGcc                                   <  1:1942751/68‑1 (MQ=255)
                         taGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCGGACGCAGGCCGAGAATGTTTTCGCCAATcgcc                            >  1:2064943/1‑69 (MQ=255)
                         taGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATcgcc                            >  1:459107/1‑69 (MQ=255)
                              cgccgcACCTTTCAGGCCATACAGGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCttt                       <  1:1942977/69‑1 (MQ=255)
                              cgccgcACCTTTCAGGCCATACAGGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCttt                       <  1:2273641/69‑1 (MQ=255)
                                ccgcACCTTTCAGGCCATACAGGCACAGCGGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTa                      <  1:754397/68‑1 (MQ=255)
                                 cgcACCTTTCAGGCCATACAGGCACAGCGGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATc                    <  1:3153293/69‑1 (MQ=255)
                                  gcACCTTTCAGGCCATACAGGCACAGCAGACGCAGGCCGAGAATGtt                                         <  1:856009/47‑1 (MQ=255)
                                      cTTTCAGGCCATACAGGCACAGCGGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGt               <  1:392398/69‑1 (MQ=255)
                                      cTTTCAGGCCATACAGGCACAGCGGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGt               <  1:608591/69‑1 (MQ=255)
                                      cTTTCAGGCCATACAGGCACAGCGGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGt               <  1:700193/69‑1 (MQ=255)
                                      cTTTCAGGCCATACAGGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGt               <  1:68594/69‑1 (MQ=255)
                                        ttCAGGCCATACAGGCACAGCGGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTa             >  1:688626/1‑69 (MQ=255)
                                        ttCAGGCCATACAGGCACAGCGGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTa             >  1:903143/1‑69 (MQ=255)
                                         tCAGGCCATACAGGCACAGCGGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAg            >  1:926669/1‑69 (MQ=255)
                                         tCAGGCCATACAGGCACAGCGGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAg            >  1:1750975/1‑69 (MQ=255)
                                          cgggCCATACAGGCACAGCGGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGt               <  1:2552826/63‑1 (MQ=255)
                                          cAGGCCATACAGGCACAGCGGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGt               >  1:213941/1‑65 (MQ=255)
                                          cAGGCCATACAGGCACAGCGGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGt               <  1:1145596/65‑1 (MQ=255)
                                          cAGGCCATACAGGCACAGCGGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGt               >  1:80149/1‑65 (MQ=255)
                                          cAGGCCATACAGGCACAGCGGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGt               <  1:2216148/65‑1 (MQ=255)
                                                  aCAGGCACAGCGGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTAt                     >  1:913188/1‑51 (MQ=255)
                                                   cAGGCACAGCGGACGCAGGCCGAGAATGTTTTCGCCAAt                                >  1:642242/1‑39 (MQ=255)
                                                   cAGGCACAGCGGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTa       <  1:2477253/64‑1 (MQ=255)
                                                   cAGGCACAGCGGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTa       <  1:915163/64‑1 (MQ=255)
                                                   cAGGCACAGCGGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTAAAtt   >  1:2056393/1‑68 (MQ=255)
                                                   cAGGCACAGCGGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTAAAt    <  1:3488343/67‑1 (MQ=255)
                                                   cAGGCACAGCGGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTAAATTc  <  1:2435653/69‑1 (MQ=255)
                                                   cAGGCACAGCGGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTAAATTc  <  1:1829976/69‑1 (MQ=255)
                                                   cAGGCACAGCGGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTAAATTc  >  1:2789228/1‑69 (MQ=255)
                                                   cAGGCACAGCGGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTAAATTc  >  1:1624123/1‑69 (MQ=255)
                                                   cAGGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTAAAtt   >  1:469378/1‑68 (MQ=255)
                                                   cAGGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTAAAtt   >  1:2898301/1‑68 (MQ=255)
                                                             |                                                          
ACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTAAATTC  >  W3110S.gb/913718‑913837

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: