Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,605,627 C→T 8.6% G273G (GGC→GGT)  hyfD → hydrogenase 4, membrane subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,605,6270CT 8.6% 118.0 / 2.9 49G273G (GGC→GGT) hyfDhydrogenase 4, membrane subunit
Reads supporting (aligned to +/- strand):  ref base C (26/18);  new base T (2/2);  total (29/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.80e-01

CGATGCAAGCGTGGCTACCGGACGCAATGGAAGCGCCAACACCGATCAGCGCCTATCTCCACGCCGCATCGATGGTGAAAGTGGGCGTTTACATTTTTGCCCGCGCTATTATCGACGGCGGCAATATCCCGCATGTGATTGGCGGCGTTGGCATGGTCATGGCACTGGTCACCATTCTTTATGGCTTTCTGATGTATTTGCCACAGCAGGATATGAAGCGGTTGCTAGCCTGGTCGACCATCACTCAACTTGGCTGGATGTTCTTCGGCTTGTCGCTCTCCATCTTCGGCTCGC  >  NC_000913/2605482‑2605775
                                                                                                                                                 |                                                                                                                                                    
cGATGCAAGCGGGGCTATCGGACGCAGTGGAAGCGCCAACACCGAGCGGGGCCTCTCTCCAGGCCGCATCGATAGTGAAAGTGGGCGTTTACATTTTTGCCCGCGCGCTTATCGACGGCGGCACTATCCCGCATGTGATTAGCGGCGtt                                                                                                                                                   >  2:663576/1‑149 (MQ=255)
 ggtGCAAGCGTGGCTACCGGACGCAATGGAAGCGCCAACACCGATCAGCGCCTATCTCCACGCCGCATCGATGGTGAAAGTGGGCGTTTACATTTTTGCCCGCGCTATTATCGACGGCGGCAATATCCCGCATGTGATTGGCGGCGTTg                                                                                                                                                  <  2:96109/147‑1 (MQ=255)
  aTGCAAGCGTGGCTACCGGACGCAATGGAAGCGCCGACACCGATCAGCGCCTATCTCCACGCCGCATCGATGGTGAAAGTGGGCGTTTACATTTTTGCCCGCGCTATTATCGACGGCGGCAATATCCCGCATGTGATTGGCGGCGTTgg                                                                                                                                                 >  2:453824/1‑149 (MQ=255)
  aTGCAAGCGTGGCTACCGGACGCAATGGAAGCGCCAACACCGATCAGCGCCTATCTCCACGCCGCATCGATGGTGAAAGTGGGCGTTTACATTTTTGCCCGCGCTATTATCGACGGCGGCAATATCCCGCATGTGATTGGCGGGGTTgg                                                                                                                                                 >  1:216712/1‑149 (MQ=255)
  aTGCAAGCGTGGCTACCGGACGCAATGGAAGCGCCAACACCGATCAGCGCCTATCTCCACGCCGCATCGATGGTGAAAGTGGGCGTTTACATTTTTGCCCGCGCTATTATCGACGGCGGCAATATCCCGCATGTGATTGGCGGCGTTgg                                                                                                                                                 >  1:216532/1‑149 (MQ=255)
  aTGCAAGCGCGGCCACCGGACGCAATGGAAGCGCCAACACCGATCAGCGCCTATCTCCACGCCGCATCGATGGTGAAAGAGGGCGTTTACATTTTTGCCCGCGCGATTATCGACGGCGGCAATATCCCGCATGTGATTGGCGGCGTTgg                                                                                                                                                 <  2:230761/149‑1 (MQ=255)
        gCGTGGCTACCGGACGCAATGGAAGCGCCAACACCGATCAGCGCCTATCTCCACGCCGCATCGATGGTGAAAGTGGGCGTTTACATTTTTGCCCGCGCTATTATCGACGGCGGCAATATCCCGCATGTGATCGGCGGCGTTGGCATGGt                                                                                                                                           >  2:749974/1‑149 (MQ=255)
               tACCGGACGCAATGGAAGCGCCAACACCGATCAGCGCCTATCTCCACGCCGCATCGATGGTGAAAGTGGGCGTTTACATTTTTGCCCGCGCTATTATCGACGGCGGCAATATCCCGCATGTGATTGGCGGCGTTGGCATGGTCATGGCa                                                                                                                                    >  1:277831/1‑149 (MQ=255)
                 gcGGAAGCAATGGAAGCGCCAACAGCGACCAGCGCCTATCTCAACGTTGCCTCGATGGTGAAAGTGGCGGTTTGCGTCTTTGCCCGCGCTACTCACGAAGGCAGCAAAAACACCAATGTGATTGGCGGCGTTGGCATGGTCATggaaat                                                                                                                                  >  2:280410/2‑145 (MQ=255)
                 ccGGACGCAATGGAAGCGCCAACACCGATCAGCGCCTATCTCCACGCCGCATCGATGGTGAAAGTGGGCGTTTACATTTTTGCCCGCGCTATTATCGACGGCGGCAATATCCCGCATGTGATTGGCGGCGTTGGCATGGTCATGGCACt                                                                                                                                  >  2:355501/1‑149 (MQ=255)
                       gCAATGGAAGCGCCAACACCGATCCGCGCCTATCTCCACGCCGCATCGATGGTGAAAGTGGGCGTTTACATTTTTGCCCGCGCTATTATCGACGGCGGCAATATCCCGCATGTGATTGGCGGCGTTGGCATGGTCATGGCACTGGTCAc                                                                                                                            <  1:48044/149‑1 (MQ=255)
                       gCAATGGAAGCGCCAACACCGATCAGCGCCTATCTCCACGCCGCATCGATGGTGAAAGTGGGCGTTTACATTTTTGCCCGCGCTATTATCGACGGCGGCAATATCCCGCATGTGATTGGCGGCGTTGGCATGGTCATGGCGCTGGTCAc                                                                                                                            >  2:297296/1‑149 (MQ=255)
                         aaTGGAAGCGCCAACACCGATCAGCGCCTATCTCTTCTTTGCATCGATGGTGAAAGTGGGCGTTTACATTTTTGCCCGCGCTATTATCGACGGCGGCAATATCCCGCATGTGATTGGCGGCGTTGGCATGGTCATGGCACTGGTCACCa                                                                                                                          <  2:601207/149‑1 (MQ=255)
                           tGGAAGCGCCAACACCGATCAGCGCCTATCTCCACGCCGCATCGATGGTGAAAGTGGGCGTTTACATTTTTTCCCGCGCTATTATCGACGGCGTCGATATCCCGCATGTGATTGTCGGCGTTGACATGGTCATGGCACTGGTCACCAtt                                                                                                                        >  1:723158/1‑149 (MQ=255)
                                       cacCGATCAGCGCCTATCTCCACGCCGCATCGATGGTGAAAGTGGGCGTTTACATTTTTGCCCGCGCTATTATCGACGGCGGCAATATCCCGCATGTGATTGGCGGCGTTGGCATGGTCATGGCACTGGTCACCATTCTTTATGGCttt                                                                                                            <  1:265907/149‑1 (MQ=255)
                                            aTCAGCGCCTATCTCCACGCCGCATCGATGGTGAAAGTGGGCGTTTACATTTTTGCCCGCGCTATTATCGACGGCGGCAATATCCCGCATGTGATTGACGGCGTTGGCATGGGCATGGCCCTGGTCACCACTATTTATGGCTTTcggcg                                                                                                       >  2:280867/1‑145 (MQ=255)
                                                 ggcCTATCCCCACGCCGCATCGATGGTGAAAGTGGGCGTTTACATTGTTGCCCGCGCTATTATCGACGGCGGCAATATCCCGCATGTGATTGGCGGCGTTGGCATGGTCATGGCACTGGTCACCATTCTTTATGGCTTTCTGATGTAtt                                                                                                  <  2:266431/148‑1 (MQ=255)
                                                    cTATCTCCCCGCCGCATCGATGGTGAAAGTGGGCGTTTACATTTTTGCCCGCGCTATTATCGACGGCGGCAATATCCCGCATGTGATTGGCGGCGTTGGCATGGTCATGGCACTGGTCACCATTCTTTATGGCTTTCTGATGTATTTGc                                                                                               <  1:189157/149‑1 (MQ=255)
                                                     tATCTCCACGCCGCATCGATGGTGAAAGTGGGCGTTTACATTTTTGCCCGCGCCATCATCAACGCCGGCAATATCCCGCAGGTGATTGGCGGCGTTGGCATGGTCATGGCCCTGGTCCCCATTCTTTTTGGCTGTCTGTTGTTTTTgac                                                                                              >  1:255232/1‑147 (MQ=255)
                                                      aTCTCCACGCCGCATCGATGGTGAAAGTGGGCGTTTACATTTTTGCCCGCGCTATTATCGACGGCGGCAATATCCCGCATGTGATTGGCGGCGTTGGCATGGTCATGGCACTGGTCACCATTCTTTATGGCTTTCTGATGGATTTGCca                                                                                             >  1:395078/1‑149 (MQ=255)
                                                      aTCTCCACGCAGCAGCGATGGTGAAAGTGGGCGTTTACATTTTTGCCCGCGCTATTATCGACGGCGGCAATATCCCGCATGTGATTGGCGGCGTTGGCATGGTCATGGCACTGGTCACCATTCTTTATGGCTTTCTGATGTATTTGCca                                                                                             <  1:568056/149‑1 (MQ=255)
                                                           cACGCCGCATCGATGGTGAAAGTGGGCGTTTACATTTTTGCCCGCGCTATTATCGACGGCGGCAATATCCCGCATGTGATTGGCGGCGTTGGCATGGTCATGGCACTGGTCACCATTCTTTATGGCTTTGTGATGTATTTGCCAcagca                                                                                        >  1:442008/1‑149 (MQ=255)
                                                                cgcATCGATGGTGAAAGTGGGCGTTTACATTTTTGCCCGCGCTATTATCGACGGCGGCAATATCCCGCATGTGATTGGCGGCGTTGGCATGGTCATGGCACTGGTCACCATTCTTTATGGCTTTCTGATGTATTTGCCACAGCAGGata                                                                                   <  1:606851/149‑1 (MQ=255)
                                                                       aTGGTGAAAGTGGGCGTTTACATTTTTGCCCGCGCTATTATCGACGGCGGCAATATCCCGCATGTGATTGGCGGCGTTGGCCTGGTCATGGCACTGGTCACCATTCTTTATGGCTTTCTGATGTATTTGCCACAGCAGGATATTAAGCg                                                                            >  1:114425/1‑149 (MQ=255)
                                                                       aTGGTGAAAGTGGGCGTTTACATTTTTGCCCGCGCTATTATCGACGGCGGCAATATCCCGCATGTGATTGGCGGCGTTGGCATGGTCATGGCACTGGTCACCATTCTTTATGGCTTTCTGATGTATTTGCCACAGCAGGATATGAAGCg                                                                            <  2:10073/149‑1 (MQ=255)
                                                                              aaGTGGGCGTTTACATTTTTTTTTTTTTTTTTATCGACGGCGGCAATATCCCGCATGTGATTGGCGGTGTTGGCATGGTCATGGCACTGGTCACCATTCTTTATGGCTTTCTGATGTATTTGCCACAGCAGGATATGAAGCGGTTGCTa                                                                     >  2:749098/1‑149 (MQ=255)
                                                                              aaGTGGGCGTTTACATTTTTGCCCGCGCTATTATCGACGGCGGCAATATCCCGCATGTGATTGGCGGTGTTGGCATGGTCATGGCACTGGTCACCATTCTTTATGGCTTTCTGATGTATTTGCCACAGCAGGATATGAAGCGGTTGCTa                                                                     >  2:749167/1‑149 (MQ=255)
                                                                                 tGGGCGTTTACATTTTTGCCCGCGCTATTATCGACGGCGGCAATATCCCGCATGTGATTGGCGGCGTTGGCATGGTCATGGCACTGGTCACCATTCTTTATGGCTTTCTGATGTATTTGCCACAGCAGGATATGAAGCGGTTGCTAGcc                                                                  >  1:426449/1‑149 (MQ=255)
                                                                                         tACATTTTTGCCCGCGCTATTATCGACGGCGGCAATATCCCGCATGTGATTGGCGGCGTTGGCATGGTCATGGCACTGGTCACCATTCTTTATGGCTTTCTGATGTATTTGCCACAGCAGGATATGAAGCGGTTGCTAGCCTGGTCGAc                                                          <  1:715994/149‑1 (MQ=255)
                                                                                          aCATTTTGGCCCGCGCTATTATCGCCGGCGGCAATATGCCTCATGTTATTGGCGGCGTTGGCATGGTCATGGCACTGGGCACCATTCTTTAGGGCTTTCTCATGTATTTGCCACAGCAGGATATGAACCGGTTGCTAGCCTGGTCGAcc                                                         <  2:255232/149‑1 (MQ=255)
                                                                                             tttttGCCCGCGCTATTATCGACGGCGGCAATATCCCGCATGTGATTGGCGGCGTTGGCATGGTCATGGCACTGGTCACCATTCTTTATGGCTTTCTGATGTATTTGCCACAGCAGGATATGAAGCGGTTGCTAGCCTGGTCGACCATc                                                      <  1:503847/149‑1 (MQ=255)
                                                                                               tttGCCCGCGCTATTATCGACGGCGGCAATATCCCGCATGTGATTGGCGGCGTTGGCATGGTCATGGCACTGGTCACCATTCTTTATGGCTTTCTGATGTATTTGCCACAGCAGGATATGAAGCGGTTGCTAGCCTGGTCGACCATCAc                                                    >  2:582740/1‑149 (MQ=255)
                                                                                                   cccGCGCTATTATCGACGGCGGCAATATCCCGCATGTGATTGGCGGCGTTGGCATGGTCATGGCACTGGTCACCATTCTTTATGGCTTTCTGATGTATTTGCCACAGCAGGATATGAAGCGGTTGCTAGCCTGGTCGACCATCACTCaa                                                >  1:123176/1‑149 (MQ=255)
                                                                                                    ccGCGCTATTAGCGACGGCGGCAATATCCCGCATGTGATTGGCGGCGTTGGCATGGTCATGGCACTGGTCACCATTCTTTATGGCTTTCTGATGTATTTGCCACAGCAGGATATGAAGCGGTTGCTAGCCTGGTCGACCATCACTCAAc                                               <  1:370844/149‑1 (MQ=255)
                                                                                                     cgcgCTATTATCGACGGCGGCAATATCCCGCATGTGATTGGCGGCGTTGGCATGGTCATGGCACTGGTCACCATTCTTTATGGCTTTCTGATGTATTTGCCACAGCAGGATATGAAGCGGTTGCTAGCCTGGTCGACCATCACTaaatc                                              >  2:560321/1‑147 (MQ=255)
                                                                                                     cgcgCTATTATCGACGGCGGCAATATCCCGCATGTGATTGGCGGCGTTGGCATGGTCATGGCACTGGTCACCATTCTTTATGGCTTTCTGATGTATTTGCCACAGCAGGATATGAAGCGGTTGCTAGCCTGGTCGACCATCACTCAACt                                              >  2:283466/1‑149 (MQ=255)
                                                                                                                    ggcggcAATATCCCGCATGTGATTGGCGGCGTTGGCATGGTCATGGCACTGGTCACCATTCTTTATGGCTTTCTGATGTATTTGCCACAGCAGGATATGAAGCGGTTGCTAGCCTGGTCGACCATCAATCAACTTGGCTGGATGttctt                               >  2:462939/1‑149 (MQ=255)
                                                                                                                         cAATATCCCGCCTGTGATTGGCGGCGTTGGCATGGTCATGGCACTGGTCACCATTCTTTATGGCTTTCTGCTGTATTTGCCACCGCAGGATCTGACGCGGTTGCTCGCCTGGTCGACCATCACTCAACTTGGCTGGATGTTCTTCGGCt                          <  2:441113/149‑1 (MQ=255)
                                                                                                                          gatagCCCGCATGTGATTGGCGGCGTTGTCATGGTCATGTCACTGGTCACCATTCTTTATGGCTTTCTGATGTATTTGCCACATCAGGATATGAAGCGGTTGCTAGCCTGGTCGACCATCACTCAACTTGGCTGGATGTTCTTCGGCtt                         <  2:395078/148‑1 (MQ=255)
                                                                                                                            tatCCCGCATGTGATTGGCGGCGTTGGAAGGGTCATGGCACTGGTCACCATTCTTTATGGCTTTCTGATGTATTTGCCACAGCAGGATATGAACCGGGTGCTAGCCTGATCGACCATCACTCAACTTGACTGGATGTTCTTCGGCGTgg                       >  2:273455/1‑148 (MQ=255)
                                                                                                                              tCCCGCATGTGATTGGCGGTGTTGGCATGGTCATGGCACTGGTCACCATTCTTTATGGCTTTCTGATGTATTTGCCACAGCAGGATATGAAGCGGTTGCTAGCCTGGTCGACCATCACTCAACTTGGCTGGATGTTCTTCGGCTTGTCg                     <  1:749167/149‑1 (MQ=255)
                                                                                                                              tCCCGCATGTGATTGGCGGTGTTGGCATGGTCATGGCACTGGTCACCATTCTTTATGGCTTTCTGATGTATTTGCCACAGCAGGATATGAAGCGGTTGCTAGCCTGGTCGACCATCACTCAACTTGGCTGGATGTTCTTCGGCTTGTCg                     <  1:749098/149‑1 (MQ=255)
                                                                                                                                ccGCATGTGATTGGCGGCGTTGGCATGGTCATGGCACTGGTCACCATTCTTTATGGCTTTCTGATGTATTTGCCACAGCAGGATATGAAGCGGTTGCTAGCCTGGTCGACCATCACTCAACTTGGCTGGATGTTCTTCGGCTTGTCGct                   >  1:53429/1‑149 (MQ=255)
                                                                                                                                     tgtgATTGGCGGCGTTGGCATGGTCATGGCACTGGTCACCATTCTTTATGGCTTTCTGATGTATTTGCCACAGCAGGATGTGAAGCGGTTGCTGGCCTGGTCGACCGGCACTGAAGTTGGCTGGAGGTTCTTCGGCTTGTCGCtctcag              >  2:415221/1‑147 (MQ=255)
                                                                                                                                     tgtgATTGGCGGCGTTGGCATGGTCATGGCACTGGTCACCATTCTTTATGGCTTTCTGATGTATTTGCCACAGCAGGATATGAAGCGGTTGCTAGCCTGGTCGACCATCACTCAACTTGGCTGGATGTTCTTCGGCTTGTCGCTCTCCa              >  1:602148/1‑149 (MQ=255)
                                                                                                                                         aTTGGCGGCGTTGGCATGGTCATGGCACTGGTCACCATTCTTTATGGCTTTCTGATGTATTTGCCACAGCAGGATATGAAGCGGTTGCTAGCCTGGTCGACCATCACTCAACTTGGCTGGATGTTCTTCGGCTTGTCGCTCTCCATCtt          >  2:397252/1‑149 (MQ=255)
                                                                                                                                         aTTGGCCGCGTTGGCATGGTCATGGCACTGGTCACCATTCTTTATGGCTTTCTGATGTATTTGCCACAGCAGGATATGAAGCGGTTGCTAGCCTGGTCGACCATCACTCAACTTGGCTGGATGTTCTTCGGCTTGTCGCTCTCCATCtt          <  1:749974/149‑1 (MQ=255)
                                                                                                                                          ttGGCGGCTTTGGCATGGTCATGGCACTGGTCACCATTCTTTATGGCTTTCTGATGTATTTGCCACAGCAGGATATGAAGCGGTTGCTAGCCTGGTCGACCATCACTCAACTTGGCTGGATGTTCTTCGGCTTGTTGCTCCCCATCTTc         <  2:277062/149‑1 (MQ=255)
                                                                                                                                          ntGGCGGCGTTGGCATGGTCATGGCACTGGTCACCATTCTTTATGGCTTTCTGATGTATTTGCCACAGCAGGATATGAAGCGGTTGCTAGCCTGGTCGACCATCACTCAACTTGGCNGGATGTTCTTCGGCTTGTCGCTCTCCATCTTc         >  1:50365/2‑149 (MQ=255)
                                                                                                                                                 cgTTGGCATGGTCATGGCACTGGTCACCATTCTTTATGGCTTTCTGATGTATTTGCCACAGCAGGATATGAAGCGGTTGCTAGCCTGGTCGACCATCACTCAACTTGGCTGGATGTTCTTCGGCTTGTCGCTCTCCATCTTCGGCTcgc  <  1:453824/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
CGATGCAAGCGTGGCTACCGGACGCAATGGAAGCGCCAACACCGATCAGCGCCTATCTCCACGCCGCATCGATGGTGAAAGTGGGCGTTTACATTTTTGCCCGCGCTATTATCGACGGCGGCAATATCCCGCATGTGATTGGCGGCGTTGGCATGGTCATGGCACTGGTCACCATTCTTTATGGCTTTCTGATGTATTTGCCACAGCAGGATATGAAGCGGTTGCTAGCCTGGTCGACCATCACTCAACTTGGCTGGATGTTCTTCGGCTTGTCGCTCTCCATCTTCGGCTCGC  >  NC_000913/2605482‑2605775

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑