Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,602,017 (C)7→8 coding (200/2019 nt) hyfB → hydrogenase 4 component B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,602,0101.C96.0% 168.5 / ‑0.6 50coding (193/2019 nt)hyfBhydrogenase 4 component B
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base C (19/29);  total (20/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

TCTGGCTTCGCTGTTTTTACTCGGTCTGGACAGACTGGCTATTAAGCTTTCCGGCATCACATCGCTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTG‑‑CCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCATTGACTTATATGCGCGAATACGA  >  NC_000913/2601871‑2602137
                                                                                                                                            |                                                                                                                                 
tCTGGCTTCGCTGTTTTTACTCGGTCTGGACAGACTGGCTATTAAGCTTTCCGGCATCACATCGCTGGTGGGCGGCGTGATTGGCATCATCGGCGGAATTACGCGATTACATGCTGGTGTAAATTTAGTCGCCCGTTTTg                                                                                                                                    >  2:272061/1‑140 (MQ=255)
    gCTTCGCTGTTTTTACTCGGTCTGGACAGACTGGCTATTAAGCTTTCCGGCATCACATCGCTGGTGGGCGGCGTGATTGGCATCATCAGCGGAAATACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTG‑‑cccc                                                                                                                              >  2:703078/1‑140 (MQ=255)
     cTTCGCTGTTTTTACTCGGTCTGGACAGACTGGCTATTAAGCTTTCCGGCATCACATCGCTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTG‑‑ccccc                                                                                                                             <  1:545488/140‑1 (MQ=255)
           tGTTTTTACTCGGTCTGGACAGACTGGCTATTAAGCTTTCCGGCATCACATCGCTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTG‑‑CCCCCCCcttt                                                                                                                      >  1:686234/1‑140 (MQ=255)
                    tCGGTCTGGACAGACTGGCTATTAAGCTTTCCGGCATCACATCGCTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGc                                                                                                              <  1:426987/140‑1 (MQ=255)
                                 aCTGGCTATTAAGCTTTCCGGCATCACATCGCTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTg                                                                                                 >  2:574607/1‑140 (MQ=255)
                                 aCTGGCTATTAAGCTTTCCGGCATCACATCGCTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTg                                                                                                 >  1:404954/1‑140 (MQ=255)
                                  cTGGCTATTAAGCTTTCCGGCATCACATCGCTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGc                                                                                                >  1:512056/1‑140 (MQ=255)
                                   tGGCTATTGAGCTTTCCGGCATTACATCGCTGGTGTGCGCCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCg                                                                                               <  2:352719/140‑1 (MQ=255)
                                      cTATTAAGCTTTCCGGCATCACATCGCTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAAt                                                                                            >  1:553878/1‑140 (MQ=255)
                                      cTATTAAGCTTTCCGGCATCACATCGCTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAAt                                                                                            >  1:480160/1‑140 (MQ=255)
                                            aGCTTTCCGGCATCACATCGCTTGTGGGCAGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAg                                                                                      <  1:272061/140‑1 (MQ=255)
                                                  ccGGCATCACATCGCTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCtctc                                                                                <  2:546734/140‑1 (MQ=255)
                                                     gCATCACATCGCTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGc                                                                             <  1:482223/140‑1 (MQ=255)
                                                          acaTCGCTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTa                                                                        <  1:185871/140‑1 (MQ=255)
                                                                cTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTg                                                                   >  2:596128/1‑139 (MQ=255)
                                                                cTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTg                                                                   <  1:596128/139‑1 (MQ=255)
                                                                cTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGc                                                                  <  1:501841/140‑1 (MQ=255)
                                                                 tggtggGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCt                                                                 >  2:384175/1‑140 (MQ=255)
                                                                 tggtggGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTAAACGCTGGGAATGGATAGCCTCTCGGCATTTATGGTGCt                                                                 >  2:287482/1‑140 (MQ=255)
                                                                   gtggGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTgg                                                               >  1:431928/1‑140 (MQ=255)
                                                                     gggCGGCGTGATTGGCATCATCCGCGGAATTCCGCAATTGCATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGtt                                                             <  2:303599/140‑1 (MQ=255)
                                                                       gcggcgTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTAt                                                           <  2:512056/140‑1 (MQ=255)
                                                                       gcggcgTGATTGGCAGCATCAGTGGTATTACGCAATTGCATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTAt                                                           <  1:344390/140‑1 (MQ=255)
                                                                           cgTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTcc                                                       >  2:113266/1‑140 (MQ=255)
                                                                                    catcatCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCtggtg                                              <  2:404954/140‑1 (MQ=255)
                                                                                       catcaGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCtggtggtg                                           >  2:430902/1‑140 (MQ=255)
                                                                                        atcaGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCAGCATTTATGGTGCTGGTTATCTCCTTGCtggtggtgg                                          >  1:250443/1‑140 (MQ=255)
                                                                                         tcaGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGATGGTTATCTACTTGCTggtggtggt                                         >  1:309520/1‑140 (MQ=255)
                                                                                            gCGGAATTACGCAATGACATGCTGGTGTATCTTTAGTTGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTg                                      <  1:240901/140‑1 (MQ=255)
                                                                                             cGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGt                                     <  2:686234/140‑1 (MQ=255)
                                                                                             cGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGCCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGt                                     >  1:46563/1‑140 (MQ=255)
                                                                                              ggAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGtt                                    <  1:633086/140‑1 (MQ=255)
                                                                                               gAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTtc                                   <  1:430902/140‑1 (MQ=255)
                                                                                                aaTTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTtct                                  <  1:113266/140‑1 (MQ=255)
                                                                                                 aTTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTtctc                                 <  1:574607/140‑1 (MQ=255)
                                                                                                  ttACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTtctct                                <  2:39675/140‑1 (MQ=255)
                                                                                                        caTTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTc                          <  1:488926/139‑1 (MQ=255)
                                                                                                        aaTTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTc                          <  1:384175/140‑1 (MQ=255)
                                                                                                        aaTTACATGCTGGTGCAACTTTTGTGGCCAGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTc                          <  1:287482/140‑1 (MQ=255)
                                                                                                          ttACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTggtggtggt                                         >  2:365586/1‑123 (MQ=255)
                                                                                                          ttACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTggtggtggt                                         <  1:365586/123‑1 (MQ=255)
                                                                                                          ttACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCAt                        <  2:480160/140‑1 (MQ=255)
                                                                                                          ttACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCAt                        <  2:553878/140‑1 (MQ=255)
                                                                                                          ttACATGCTGGTGTAACGTTAGACGCCCGTTATG‑‑CCCCCACCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCAt                        <  2:250443/140‑1 (MQ=255)
                                                                                                                    gtgtAACTTTAGTCGCCCGTTTTGCCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTACTGGTGGAGGTTTGTTCTCTCTATTCATTGACTtata               >  1:353505/1‑140 (MQ=255)
                                                                                                                    gtgtAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCATTGACTTatat              >  2:154142/1‑140 (MQ=255)
                                                                                                                    gtgtAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCATTGACTTatat              <  2:588633/140‑1 (MQ=255)
                                                                                                                    gtgtAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCATTGACTTatat              <  1:154142/140‑1 (MQ=255)
                                                                                                                            ttAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCATTGACTTATATGCGCGAAt      >  2:366268/1‑140 (MQ=255)
                                                                                                                              aGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCATTGACTTATATGCGCGAATAc    <  1:703078/140‑1 (MQ=255)
                                                                                                                               gTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCATTGACTTATATGCGCGAATACg   >  2:122679/1‑140 (MQ=255)
                                                                                                                                tCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCATTGACTTATATGCGCGAATACGa  <  1:122679/140‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                 
TCTGGCTTCGCTGTTTTTACTCGGTCTGGACAGACTGGCTATTAAGCTTTCCGGCATCACATCGCTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTG‑‑CCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCATTGACTTATATGCGCGAATACGA  >  NC_000913/2601871‑2602137

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