Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 913,349 G→A L156L (CTG→TTG)  hcp ← hybrid‑cluster [4Fe‑2S‑2O] subunit of anaerobic terminal reductases

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913913,3490GA98.3% 186.1 / NA 59L156L (CTG→TTG) hcphybrid‑cluster [4Fe‑2S‑2O] subunit of anaerobic terminal reductases
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (31/27);  total (31/28)
Fisher's exact test for biased strand distribution p-value = 4.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAAGCGCGTTCATATCGGCAGGCCAGGTCCCCAGCCACGCCATGATTTTATGGTACTGGGCATAAATATCGTTGTCGTATTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTAAATTCTGCTGCCTGACGTTGCAGCTCGCCGAGATCGTCGCTCACCAGTTGCAGGTCAGCCATCGGGTTATCG  >  NC_000913/913214‑913482
                                                                                                                                       |                                                                                                                                     
gaagCGCGTTCATATCGGCAGGCCAGGTCCCCAGCCACGCCATGATTTTATGGTACTGGGCATAAATATCGTTGTCGTATTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGca                                                                                                                                    <  1:434468/139‑1 (MQ=255)
   gcgcGTTCATATCGGCAGGCCAGGTCCCCAGCCACGCCATGATTTTATGGTACTGGGCATAAATATCGTTGTCGTATTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCAcag                                                                                                                                 <  2:633924/139‑1 (MQ=255)
    cgcgTTCATATCGGCAGGCCAGGTCCCCAGCCACGCCATGATTTTATGGTACTGGGCATAAATATCGTTGTCGTATTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCAcagc                                                                                                                                >  1:653373/1‑139 (MQ=255)
    cgcgTTCATATCGGCAGGCCAGGTCCCCAGCCACGCCATGATTTTATGGTACTGGGCATAAATATCGTTGTCGTATTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCAcagc                                                                                                                                >  1:354751/1‑139 (MQ=255)
                               cAGCCACGCCATGATTTTATGGTACTGGGCATAAATATCGTTGTCGTATTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCa                                                                                                     >  2:1147059/1‑139 (MQ=255)
                                      gCCATGATTTTATGGTACTGGGCATAAATATCGTTGTCGTATTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATcgccg                                                                                              <  1:1147059/139‑1 (MQ=255)
                                        cATGATTTTATGGTACTGGGCATAAATATCGTTGTCGTATTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCt                                                                                            <  1:312808/139‑1 (MQ=255)
                                            aTTTTATGGTACTGGGCATAAATATCGTTGTCGTATTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGc                                                                                                                      <  1:1432144/109‑1 (MQ=255)
                                            aTTTTATGGTACTGGGCATAAATATCGTTGTCGTATTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGc                                                                                                                      >  2:1432144/1‑109 (MQ=255)
                                                      aCTGGGCATAAATATCGTTGTCGTATTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCgag                                                                                                                  <  2:124234/103‑1 (MQ=255)
                                                      aCTGGGCATAAATATCGTTGTCGTATTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCgag                                                                                                                  >  1:124234/1‑103 (MQ=255)
                                                      aCTGGGCATAAATATCGTTGTCGTATTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGtt                                                                                 >  1:842828/1‑136 (MQ=255)
                                                      aCTGGGCATAAATATCGTTGTCGTATTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGtt                                                                                 <  2:842828/136‑1 (MQ=255)
                                                       cTGGGCATAAATATCGTTGTCGTATTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGa                                                                                                                             >  2:92711/1‑91 (MQ=255)
                                                       cTGGGCATAAATATCGTTGTCGTATTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGa                                                                                                                             <  1:92711/91‑1 (MQ=255)
                                                         gggCATAAATATCGTTGTCGTATTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATc                                                                                                  <  2:1385998/116‑1 (MQ=255)
                                                         gggCATAAATATCGTTGTCGTATTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATc                                                                                                  >  1:1385998/1‑116 (MQ=255)
                                                          ggCATAAATATCGTTGTCGTATTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCtt                                                                                           >  2:2071/1‑122 (MQ=255)
                                                          ggCATAAATATCGTTGTCGTATTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCtt                                                                                           <  1:2071/122‑1 (MQ=255)
                                                          ggCATAAATATCGTTGTCGTATTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTa                                                                          >  2:479071/1‑139 (MQ=255)
                                                          ggCATAAATATCGTTGTCGTATTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTa                                                                          >  2:160898/1‑139 (MQ=255)
                                                               aaaTATCGTTGTCGTATTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCt                                                                                            >  2:1370096/1‑116 (MQ=255)
                                                               aaaTATCGTTGTCGTATTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCt                                                                                            <  1:1370096/116‑1 (MQ=255)
                                                                    tCGTTGTCGTATTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATcgccg                                                                                              <  1:32659/109‑1 (MQ=255)
                                                                    tCGTTGTCGTATTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATcgccg                                                                                              >  2:32659/1‑109 (MQ=255)
                                                                      gTTGTCGTATTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTAAATTCTGCTGcc                                                              <  2:403207/139‑1 (MQ=255)
                                                                         gTCGTATTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCg                                                                                                                    >  2:265103/1‑82 (MQ=255)
                                                                         gTCGTATTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCg                                                                                                                    <  1:265103/82‑1 (MQ=255)
                                                                         gTCGTATTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTa                                                                                >  1:114942/1‑118 (MQ=255)
                                                                         gTCGTATTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTa                                                                                <  2:114942/118‑1 (MQ=255)
                                                                                 gACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTAAATTCTGCTGCCTGACGTTGCAg                                                   <  1:479071/139‑1 (MQ=255)
                                                                                 gACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTAAATTCTGCTGCCTGACGTTGCAg                                                   >  1:1282884/1‑139 (MQ=255)
                                                                                  aCCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGt                                                                                  <  1:835177/107‑1 (MQ=255)
                                                                                  aCCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGt                                                                                  >  2:835177/1‑107 (MQ=255)
                                                                                  aCCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTAAATTCTGCTGCCTGACGTTGCAGc                                                  >  1:188895/1‑139 (MQ=255)
                                                                                    cGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTAAATTCTGCTGCCTg                                                            <  1:1532697/127‑1 (MQ=255)
                                                                                    cGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTAAATTCTGCTGCCTg                                                            >  2:1532697/1‑127 (MQ=255)
                                                                                       gaACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTAAATTCTGCTGCCTGACGTTGCAGCTCGcc                                             >  2:39491/1‑139 (MQ=255)
                                                                                            tGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATATAAGCACAGCGGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTAAATTCTGCTGCCTGACGTTGCAGCTCGCCGAGAt                                        >  2:945673/1‑139 (MQ=255)
                                                                                             gcgcGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCtt                                                                                           >  2:379099/1‑87 (MQ=255)
                                                                                             gcgcGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCtt                                                                                           <  1:379099/87‑1 (MQ=255)
                                                                                                  tgtTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCa                                                                                                                         <  2:1298461/52‑1 (MQ=255)
                                                                                                  tgtTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCa                                                                                                                         >  1:1298461/1‑52 (MQ=255)
                                                                                                             ggCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTAAATTCTGCTGCCTGACGTTGCAGCTCGCCGAGATCGTCGCTCACCAGTTGc                       <  1:160898/139‑1 (MQ=255)
                                                                                                              gccgccGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATc                                                                                                  >  1:623239/1‑63 (MQ=255)
                                                                                                              gccgccGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATc                                                                                                  <  2:623239/63‑1 (MQ=255)
                                                                                                              gccgccGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTAAATTCTGCTGCCTGACGTTGCAGCTCGCCGAGATCGTCGCTCACCAGTTGCa                      >  2:415465/1‑139 (MQ=255)
                                                                                                                 gcggcACCTTTCAGGCCATACAGGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTAAATTCTGCTGCCTGACATTGCAGCTCGCCGAGATCGTCGCTCACCAGTTGCAGGt                   <  2:354751/136‑1 (MQ=255)
                                                                                                                 gccgcACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATcgc                                                                                                >  1:604205/1‑62 (MQ=255)
                                                                                                                 gccgcACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATcgc                                                                                                <  2:604205/62‑1 (MQ=255)
                                                                                                                 gccgcACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTAAATTctg                                                                   >  1:513357/1‑91 (MQ=255)
                                                                                                                 gccgcACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTAAATTctg                                                                   <  2:513357/91‑1 (MQ=255)
                                                                                                                 gccgcACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTAAATTCTGCTGCCTGACGTTGCAGCTCGCCGAGATCGTCGCTCACCAGTTGCAGGt                   <  1:415465/139‑1 (MQ=255)
                                                                                                                   cgcACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGt                                                                           >  1:1068954/1‑81 (MQ=255)
                                                                                                                   cgcACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGt                                                                           <  2:1068954/81‑1 (MQ=255)
                                                                                                                   cgcACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTAAATTCTGCTGCCTGACGTTGCAGCTCGCCGAGATCGTCGCTCACCAGTTGCAGGTCa                 >  2:789133/1‑139 (MQ=255)
                                                                                                                                  aTACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTAAATTCTGCTGCCTGACGTTGCAGCTCGCCGAGATCGTCGCTCACCAGTTGCAGGTCAGCCATCGGGTTATCg  >  2:1412868/1‑139 (MQ=255)
                                                                                                                                   tACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTAAATTCTGCTGCCTGACGTTGCAGCTCGCCGAGATCGTCGCTCAc                              <  2:496441/110‑1 (MQ=255)
                                                                                                                                   tACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTAAATTCTGCTGCCTGACGTTGCAGCTCGCCGAGATCGTCGCTCAc                              >  1:496441/1‑110 (MQ=255)
                                                                                                                                       |                                                                                                                                     
GAAGCGCGTTCATATCGGCAGGCCAGGTCCCCAGCCACGCCATGATTTTATGGTACTGGGCATAAATATCGTTGTCGTATTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTAAATTCTGCTGCCTGACGTTGCAGCTCGCCGAGATCGTCGCTCACCAGTTGCAGGTCAGCCATCGGGTTATCG  >  NC_000913/913214‑913482

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