Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 907,975 G→T 10.5% A432A (GCC→GCA 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
*CP009273907,9750GT10.5% 113.5 / 3.8 38A432A (GCC→GCAhcphybrid‑cluster [4Fe‑2S‑2O] protein in anaerobic terminal reductases
Reads supporting (aligned to +/- strand):  ref base G (18/16);  new base T (2/2);  total (20/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.70e-01

TCTGGTGTCAGGAAACCAGGTGCAGTCGGGCCGGTGACGATATTTTTCACGCCCAGAGAAAGCAGCGTCAGCAGAATGACGATCGCTTTCTGTTCAAACCAGGAGAGCACCAGCGACAGCGGCAGATCGTTGACACCGCAGCCCAGTTTCTCTGCCAGAGTGACAGCCAGAATAATCGCTGAGTAAGCATCGTTACATTGACCTGCATCTACCAGACGCGGCAGACCTTCGATATCGCCAAACTCAAGTTTGTTAAAGCGATATTTACCACAGGCGAGGGTCAGGATCAGGCAGTCATCCGGCACGCTGGTGGCGAAATCGGTGAAGTAGTGGCGCTCGCCGCGTGCGCCGTCACAGCCACCAAGCAGGAAGATATGACGCAGTTTTTCACGGCTCACCAGATCAATCAGCGTATCAGCAGCGCCAAGCAGCGTCTGGCGACCAAAACCCACGGTGATAAGGTGCGGAATTTCGCTGTACGGGAAGCCT  >  CP009273/907703‑908191
                                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                        
tCTGGTGTCAGGAAACCAGGTGCAGTCGGGCCGGTGACGATATTTTTCACGCCCAGAGAAAGCAGCGTCAGCAGAATGACGATCGCTTTCTGTTCAAACCAGGAGAGCACCAGCGACAGCGGCAGATCGTTGACACCGCAGCCCAGTTTCTCTGCCAGAGTGACAGCCAGAATAATCGCTGAGTAAGCATCGTTACATTGACCTGCATCTACCAGACGCGGCAGACCTTCGATATCGCCAAACTCAAGTTTGTTAAAGCGATATTTACCACAGGc                                                                                                                                                                                                                        <  1:193765/275‑1 (MQ=255)
  tGGTGCCAGGTACCCAGGGTCAGTCGGGCCGGTGACTATATTTTTCACGCCACGGGAAAGCAGCGTCAGCAGAATGACGATCGCTTGCTGTTCAAACCAGGAGAGCACCAGCGACAGCGGCAGATCGTTGCCACAGCAGCCCAGTTTCTCTGCCAGAGTGACAGCCAGAATAATCGCTGAGTAAGCATCGTTACATTGACCTGCATCTACCAGACGCGGCAGACCTTCGATATCGCCAAACTCAAGTTTGTTAAAGCGATATTTACCACAGGCg                                                                                                                                                                                                                       <  2:120969/274‑1 (MQ=255)
                                             acagcGCCCAGAGAAAGCAGCGTCAGCAGAATTATGATCGCTTTCTGTTCAAACCAGGAGAGCACCAGCGACAGCGGTAGATCGTTGACACCGCAGCCCAGTTTCTCTGCCAGAGTGACAGCCAGAATAATCGCTGAGTAAGCATCGTTACATTGACCTGCATCGACCAGACGCGGCAGACCTTCGATATCGCCAAACTCAAGTTTGTTAAAGCGATATTTACCACAGGCGAGGGTCAGGATCAGGCAGTCATCCGGCACGCTGGTGGCGAAAt                                                                                                                                                                            <  1:205098/270‑1 (MQ=255)
                                                            agcagcGTCAGCAGAATGACGATCGCTTTCTGTTCAAACCAGGAGAGCACCAGCGACAGCGGCAGATCGTTGACACCGCAGCCCAGTTTCTCTGCCAGAGTGACAGCCAGAATAATCGCTGAGTAAGCATCGTTACATTGACCTGCATCTACCAGACGCGGCAGACCTTCGATATCGCCAAACTCAAGTTTGTTAAAGCGATATTTACCACAGGCGa                                                                                                                                                                                                                      >  1:353831/1‑217 (MQ=255)
                                                                gcGTCAGCAGAATGACGATCGCTTTCTGTTCAAACCAGGAGAGCACCAGCGACAGCGGCAGATCGTTGACACCGCAGCCCAGTTTCTCTGCCAGAGTGACAGCCAGAATAATCGCTGAGTAAGCATCGTTACATTGACCTGCATCTACCAGACGCGGCAGACCTTCGATATCGCCAAACTCAAGTTTGTTAAAGCGATATTTACCACAGGCGAGGGt                                                                                                                                                                                                                  <  2:353826/217‑1 (MQ=255)
                                                                      gcagAATGACGATCGCTTTCTGTTCAAACCAGGAGAGCACCAGCGACAGCGGCAGATCGTTGACACCGCAGCCCAGTTTCTCTGCCAGAGTGACAGCCAGAATAATCGCTGAGTAAGCATCGTTACATTGACCTGCATCTACCAGACGCGGCAGACCTTCGATATCGCCAAACTCAAGTTTGTTAAAGCGATATTTACCACAGGCg                                                                                                                                                                                                                       <  2:172897/206‑1 (MQ=255)
                                                                                                                       cGGCAGATCGTTGACACCGCAGCCCAGTTTCTCTGCCAGAGTGACAGCCAGAATAATCGCTGAGTAAGCATCGTTACATTGACCTGCATCTACCAGACGCGGCAGACCTTCGATATCGCCAAACTCAAGTTTGTTAAAGCGATATTTACCACAGGCGAGGGTCAGGATCAGGCAGTCATcc                                                                                                                                                                                               >  1:264493/1‑181 (MQ=255)
                                                                                                                           aGATCGTTGACACCGCAGCCCAGTTTCTCTGCCAGAGTGACAGCCAGAATAATCGCTGAGTAAGCATCGTTACATTGACCTGCATCTACCAGACGCGGCAGACCTTCGATATCGCCAAACTCAAGTTTGTTAAAGCGATATTTACCACAGGCGAGGGTCAGGATCAGGCAGTCATCCGGCa                                                                                                                                                                                           <  2:264491/181‑1 (MQ=255)
                                                                                                                                     cacCGCCGCCCAGTTTCTCTGCCATAGTGACCGCCAGAATAATCGCTGAGTAAGCATCGTTACATTGACCTGCATCTACCAGACGCGGCAGACCTTCGATATCGCCAAACTCAAGTTTGTTAAAGCGATATTTACCACAGGCGAGGGTCAGGATCAGGCAGTCATCCGGCACGCTGGTGGCGAAATCGGTGAAGTAGTGGCGCTCGCCGCGTGCGCCGTCACAGCCACCAAGCAGGAAGATATGACGCAGTTTTTCACGGCTCACCAGatcaatc                                                                                   <  1:331231/275‑1 (MQ=255)
                                                                                                                                        cGCAGCCCAGTTTCTCTGCCAGAGTGACAGCCAGAATAATCGCTGAGTAAGCATCGTTACATTGACCTGCATCTACCAGACGCGGCAGACCTTCGATATCGCCAAACTCAAGTTTGTTAAAGCGATATTTACCACAGGCGAGGGTCAGGATCAGGCAGTCATCCGGCACGCTGGTGGCGAAATCGGTGAAGTAGTGGCGCTCGCCGCGTGCGCCGTCACAGCCACCAAGCAGGAAGATATGAAGCAGGTTTTCACGGCTCACCAGAGAAATCAGc                                                                                >  1:7238/1‑275 (MQ=255)
                                                                                                                                          cAGCCCAGTTTCTCTGCCAGAGTGACAGCCAGAATAATCGCTGAGTAAGCATCGTTACATTGACCTGCATCTACCAGACGCGGCAGACCTTCGATATCGCCAAACTCAAGTTTGTTAAAGCGATATTTACCACAGGCGAgg                                                                                                                                                                                                                    >  2:523677/1‑141 (MQ=255)
                                                                                                                                              ccAGTTTCTCTGCCAGAGTGACAGCCAGAATAATCGCTGAGTAAGCATCGTTACATTGACCTGCATCTACCAGACGCGGCAGACCTTCGATATCGCCAAACTCAAGTTTGTTAAAGCGATATTTACCACAGGCGAGGGTCa                                                                                                                                                                                                                <  1:523685/141‑1 (MQ=255)
                                                                                                                                                                                    gAGTAAGCATCGTTACATTGACCTGCATCTACCAGACGCGGCAGACCTTCGATATCGCCAAACTCAAGTTTGTTAAAGCGATATTTACCACATGCGAGGGTCAGGATCAGGCAGTCATCCGGCACGCTGGTGGCGAAATCGGTGAAGTAGTGgc                                                                                                                                                             >  1:133631/1‑154 (MQ=255)
                                                                                                                                                                                    gAGTAAGCATCGTTACATTGACCTGCATCTACCAGACGCGGCAGACCTTCGATATCGCCAAACTCAAGTTTGTTAAAGCGATATTTACCACAGGCGa                                                                                                                                                                                                                      >  1:352658/1‑97 (MQ=255)
                                                                                                                                                                                        aaGCATCGTTACATTGACCTGCATCTACCAGACGCGGCAGACCTTCGATATCGCCAAACTCAAGTTTGTTAAAGCGATATTTACCACATGCGAGGGTCAGGATCAGGCAGTCATCCGGCACGCTGGTGGCGAAATCGGTGAAGTAGTGGCGCTc                                                                                                                                                         <  2:133630/154‑1 (MQ=255)
                                                                                                                                                                                        aaGCATCGTTACATTGACCTGCATCTACCAGACGCGGCAGACCTTCGATATCGCCAAACTCAAGTTTGTTAAAGCGATATTTACCACAGGCGAGGGt                                                                                                                                                                                                                  <  2:352653/97‑1 (MQ=255)
                                                                                                                                                                                               gTTACATTGACCTGCATCTACCAGACCCGGCAGAACTTCGATATCACCAAACTCAAGTTTGTTAAAGCGATATTTACCACAGGCGAGGGTCAGGATCa                                                                                                                                                                                                          >  1:371386/1‑98 (MQ=255)
                                                                                                                                                                                                 tACATTGACCTGCATCTACCAGACGCGGCAGACCTTCGATATCGCCAAACTCAAGTTTGTTAAAGCGATATTTACCACAGGCGAGGGTCAGGATCAGGCAGTCATCCGGCACGCTGGTGGCGAAATCGGTGAAGTAGTGGCGCTCGCCGCGTGCGCCg                                                                                                                                            >  1:238989/1‑158 (MQ=255)
                                                                                                                                                                                                   cATTGACCTGCATCTACCAGACGCGGCAGACCTTCGATATCGCCAAACGCAAGTTTGTTAAAGCGATATTTACCACAGGCGAGGGTCAGGATCAGGCa                                                                                                                                                                                                      <  2:371381/98‑1 (MQ=255)
                                                                                                                                                                                                        aCCTGCATCTACCAGACGCGGCAGACCTTCGATATCGCCAAACTCAAGTTTGTTAAAGCGATATTTACCACAGGCGAGGGTCAGGATCAGGCAGTCATCCGGCACGCTGGTGGCGAAATCGGTGAAGTAGTGGCGCTCGCCGCGTGCGCCGTCACAGCCACCAAGCAGGAAGATATGACGCAGTTTTTCACGGCTCACCAGATCAATCAGCGTATCAGCAGCGCCAAGCAGCGTCTGGCGACCAAAACCCACGGTGATAAGGTGCGGAATTTCGc                >  1:306772/1‑275 (MQ=255)
                                                                                                                                                                                                        aCCTGCATCTACCAGACGCGGCAGACCTTCGATATCGCCAAACTCAAGTTTGTTAAAGCGATATTTACCACAGGCGAGGGTCAGGATCAGGCAGTCATCCGGCACGCTGGTGGCGAAATCGGTGAAGTAGTGGCGCTCGCCGCGTGCGCCGTCACAGCCACCAAGCAGGAAGATATGACGCAGTTTTTCACGGCTCACCAGATCAATCAGCGTATCAGCAGCGCCAAGCAGCGTCTGGCGACCAAAACACACGGTGATAAGGTGCGGAAtgt                   >  1:229790/1‑270 (MQ=255)
                                                                                                                                                                                                        aCCTGCATCTACCAGACGCGGCAGACCTTCGATATCGCCAAACTCAAGTTTGTTAAAGCGATATTTACCACAGGCGAGGGTCAGGATCAGGCAGTCATCCGGCACGCTGGTGGCGAAATCGGTGAAGTAGTGGCGCTCGCCGAGTGCGCCTTNACAGCCACCAAGCAGGAAGATATGACGCAGTTTTTCACGGCTAACCAGatc                                                                                       >  1:483629/1‑204 (MQ=255)
                                                                                                                                                                                                                 tcccAGACGCGGCAGACCTGCGATATCCCCAAACTCAAGTTGTGTAAGGCGTTATTTACCACAGGCGAGGGTCAGGATCAGGCAGTCATCCGGCACGCTGGTGGCGAAATCGGTGAAGTAGTGGCGCTCGCCGCGTGCGCCGTCACAGCCACCAAGCAGGAAGATATGACGCAGTTTTTCACGGCTCACCAGATCAATCAGCGTATCAGCAGCGCCAAGCAGCGTCTGGCGACCAAAACCCACGGTGATAAGGTGCGGAATTTCGCTGTACGGGa       <  2:229788/273‑1 (MQ=255)
                                                                                                                                                                                                                  aCCAGACGCGGCAGACCTTCGATATCGCCAAACTCAAGTTTGTTAAAGCGATATTTACCACATGCGAGGGTCAGGATCAGGCAGTCATCCGGCACGCTGGTGGCGAAATCGGTGAAGTAGTGGCGCTCGCCGCGTGCGCCGTCACAGCCACCAAGCAGGAAGatat                                                                                                                   >  2:43668/1‑166 (MQ=255)
                                                                                                                                                                                                                  aCCAGACGCGGCAGACCTTCGATATCGCCAAACTCAAGTTTGTTAAAGCGATATTTACCACAGGCGAGGGTCAg                                                                                                                                                                                                               >  1:28339/1‑74 (MQ=255)
                                                                                                                                                                                                                  aCCAGACGCGGCAGACCTTCGATATCGCCAAACTCAAGTTTGTTAAAGCGATATTTACCACAGGCGAGGGTCAGGATCAGGCAGTCATCCGGCACGCTGGTGGCGAAATCGGTGAAGTAGTGGCGCTCGCCGCGTGCGCCGTCACAGCCACCAAGCAGGAAGATATGACGCAGTTTTTCACGGCTCACCAGATCAATCAGCGTATCAGCAGCGCCAAGCAGCGTCTGGCACCCAAAACCAACGGAGCTAAAGGACGGAACTTCGCTGGACGGGaa      >  1:41020/1‑275 (MQ=255)
                                                                                                                                                                                                                    cAGACGCGGCAGACCTTCGATATCGCCAAACTCAAGTTTGTTAAAGCGATATTTACCACAGGCGAGGGTCAGGATCAGGCAGTCATCCGGCACGCTGGTGGCGAAATCGGTGAAGTAGTGGCGCTCGCcgcg                                                                                                                                                   >  1:375117/1‑132 (MQ=255)
                                                                                                                                                                                                                      gACGCGGCAGACCTTCGATATCGCCAAACTCAAGTTTGTTAAAGCGATATTTACCACATGCGAGGGTCAGGATCAGGCAGTCATCCGGCACGCTGGTGGCGAAATCGGTGAAGTAGTGGCGCTCGCCGCGTGCGCCGTCACAGCCACCAAGCAGGAAGATATGACg                                                                                                               <  1:43668/166‑1 (MQ=255)
                                                                                                                                                                                                                      gACGCGGCAGACCTTCGATATCGCCAAACTCAAGTTTGTTAAAGCGATATTTACCACAGGCGAGGGTCAGGATc                                                                                                                                                                                                           <  2:28339/74‑1 (MQ=255)
                                                                                                                                                                                                                       aCGCGGTAGACCTTCGATATCGCCAAACTCAAGTTTGTTAAAGCGATATTTACCACAGGCGAGGGTCAGGATCAGGCAGTCATCCGGCACGCTGGTGGCGAAATCGGTGAAGTAGTGGCGCTCGCCGCGTGCGCCGTCACAGCCACCAAGCAGGAAGATATGACACAGTTTTTCACGGCTCACCAGat                                                                                        >  2:341653/1‑188 (MQ=255)
                                                                                                                                                                                                                       aCGCGGCAGACCTTCGATATCGCCAAAGTCAAGTTTGTTAAAGCGATATTTACCACAGGCGAGGGTCAGGATCAGGCAGTCATCCGGCACGCTGGTGGCGAAATCGGTGAAGTAGTGGCGCTCGCCGCGTGCGCCGTCACAGCCACCAAGCAGGAAGATATGACGCAGTTTTTCACGGCTCACCAGATCAATCAGCGTATCAGCAGCGCCAAGCAGCGTCTGGCGACCAAAACCCACGGTGATAAAGTGCGGAAATTCGCTGTACGGGAAGCCt  >  1:340024/1‑274 (MQ=255)
                                                                                                                                                                                                                        cgcgGCAGACCTTCGATATCGCCAAACTCAAGTTTGTTAAAGCGATATTTACCACAGGCGAGGGTCAGGATCAGGCAGTCATCCGGCACGCTGGTGGCGAAATCGGTGAAGTAGTGGCGCTCGCCGCGTgcg                                                                                                                                               <  2:375112/132‑1 (MQ=255)
                                                                                                                                                                                                                           ggCAGACCTTCGATATCGCCAAACTCAAGTTTGTTAAAGCGATATTTACCACAGGCGAGGGTCAGGATCAGGCAGTCAt                                                                                                                                                                                                 >  1:366401/1‑79 (MQ=255)
                                                                                                                                                                                                                           ggCAGACCTTCGATATCGCCAAACTCAAGTTTGTTAAAGCGATATTTACCACAGGCGAGGGTCAGGATCAGGCAGTCATCCGGCACGCTGGTGGCGAAATCGGTGAAGTAGTGGCGCTCGCCGCGTGCGCCGTCACAGCCACCAAGCAGGAAGATATGACACAGTTTTTCACGGCTCACCAGatcaat                                                                                    <  1:341658/188‑1 (MQ=255)
                                                                                                                                                                                                                               gACCTTCGATATCGCCAAACTCAAGTTTGTTAAAGCGATATTTACCACAGGCGAGGGTCAGGATCAGGCAGTCATCCgg                                                                                                                                                                                             <  2:366396/79‑1 (MQ=255)
                                                                                                                                                                                                                                                          tGTTAAAGCGATATTTACCACAGGCGAGGGTCAGGATCAGGCAGTCATCCGGCACGCTGGTGGCGAAATCGGTGAAGTAGTGGCGCTCGCCGCGTGCGCCGTCACAGCCACCAAGCAGGAAGATATGACGCAGTTTTTc                                                                                                      >  2:471490/1‑139 (MQ=255)
                                                                                                                                                                                                                                                              aaaGCGATATTTACCACAGGCGAGGGTCAGGATCAGGCAGTCATCCGGCACGCTGGTGGCGAAATCGGTGAAGTAGTGGCGCTCGCCGCGTGCGCCGTCACAGCCACCAAGCAGGAAGATATGACGCAGTTTTTCACgg                                                                                                  <  1:471497/139‑1 (MQ=255)
                                                                                                                                                                                                                                                                               aGGCGAGGGTCAGGATCAGGCAGTCATCCGGCACGCTGGTGGCGaaa                                                                                                                                                                             >  1:485104/1‑47 (MQ=255)
                                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                        
TCTGGTGTCAGGAAACCAGGTGCAGTCGGGCCGGTGACGATATTTTTCACGCCCAGAGAAAGCAGCGTCAGCAGAATGACGATCGCTTTCTGTTCAAACCAGGAGAGCACCAGCGACAGCGGCAGATCGTTGACACCGCAGCCCAGTTTCTCTGCCAGAGTGACAGCCAGAATAATCGCTGAGTAAGCATCGTTACATTGACCTGCATCTACCAGACGCGGCAGACCTTCGATATCGCCAAACTCAAGTTTGTTAAAGCGATATTTACCACAGGCGAGGGTCAGGATCAGGCAGTCATCCGGCACGCTGGTGGCGAAATCGGTGAAGTAGTGGCGCTCGCCGCGTGCGCCGTCACAGCCACCAAGCAGGAAGATATGACGCAGTTTTTCACGGCTCACCAGATCAATCAGCGTATCAGCAGCGCCAAGCAGCGTCTGGCGACCAAAACCCACGGTGATAAGGTGCGGAATTTCGCTGTACGGGAAGCCT  >  CP009273/907703‑908191

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