Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,852,356 T→A L74Q (CTG→CAG)  hypD → hydrogenase maturation protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,852,3560TA97.0% 226.5 / ‑1.2 66L74Q (CTG→CAG) hypDhydrogenase maturation protein
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base A (38/26);  total (39/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.31e-01

GTGGCGGTCATACCCACGCTATCTTTAAATTCGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACTGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTACCTTTGGCGACGCCATGCGCGTGCCGGGGAAACA  >  NC_000913/2852257‑2852452
                                                                                                   |                                                                                                
gtgGCGGTCATACCCACGCTATCTTTAAATTCGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACt                                                                                                  >  1:1374210/1‑100 (MQ=255)
          tACCCACGCTATCTTTAAATTCGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGt                                                                                        <  2:712519/100‑1 (MQ=255)
          tACCCACGCTATCTTTAAATTCGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGt                                                                                        >  2:541237/1‑100 (MQ=255)
           aCCCACGCTATCTTTAAATTCGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTa                                                                                       >  1:1033640/1‑100 (MQ=255)
           aCCCACGCTATCTTTAAATTCGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTa                                                                                       <  1:789533/100‑1 (MQ=255)
            cccACGCTATCTTTAAATTCGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTa                                                                                       >  1:626854/1‑99 (MQ=255)
            cccACGCTATCTTTAAATTCGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAg                                                                                      >  2:1433233/1‑100 (MQ=255)
            cccACGCTATCTTTAAATTCGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAg                                                                                      <  1:1741272/100‑1 (MQ=255)
              cACGCTATCTTTAAATTCGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGaa                                                                                    >  1:563223/1‑100 (MQ=255)
                  cTATCTTTAAATTCGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGa                                                                                >  2:332902/1‑100 (MQ=255)
                   tATCTTTAAATTCGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCGACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGac                                                                               <  1:370740/100‑1 (MQ=255)
                     tCTTTAAATTCGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGacac                                                                             >  2:193596/1‑100 (MQ=255)
                     tCTTTAAATTCGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGacac                                                                             >  2:1139296/1‑100 (MQ=255)
                      cTTTAAATTCGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCAGATGGGTAGAATCGACAcc                                                                            >  1:1906449/1‑100 (MQ=255)
                       tttAAATTCGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCt                                                                           >  1:577749/1‑100 (MQ=255)
                       tttAAATTCGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCt                                                                           >  2:3088/1‑100 (MQ=255)
                       tttAAATTCGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCt                                                                           >  2:1131302/1‑100 (MQ=255)
                          aaaTTCGGCCTCGACCAGTTACTGCCGTAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCg                                                                        <  2:277481/100‑1 (MQ=255)
                             ttCGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTgg                                                                     <  2:1356912/100‑1 (MQ=255)
                              tCGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGga                                                                    <  1:358914/100‑1 (MQ=255)
                                ggCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGgaga                                                                  <  2:301412/100‑1 (MQ=255)
                                 gCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTTCAGCCGATGGGTAGAATCGACACCTGCGTGGAGAt                                                                 >  2:1212113/1‑100 (MQ=255)
                                  ccTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGAtt                                                                <  2:294944/100‑1 (MQ=255)
                                    tCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTgc                                                              <  2:1823229/100‑1 (MQ=255)
                                       aCCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACTGCCGATGGGTAGAATCGACACCTGCGTGGAGATTgccag                                                           <  2:1190082/100‑1 (MQ=255)
                                       aCCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTgccag                                                           <  2:1149970/100‑1 (MQ=255)
                                          aGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTgccagcca                                                        >  2:1011613/1‑100 (MQ=255)
                                          aGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCAGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTgccagcca                                                        >  1:791196/1‑100 (MQ=255)
                                             tACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATcc                                                     <  2:511601/100‑1 (MQ=255)
                                                 gCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGa                                                  >  1:1222810/1‑99 (MQ=255)
                                                  ccGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAg                                                >  1:699382/1‑100 (MQ=255)
                                                   cGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGt                                               >  1:616682/1‑100 (MQ=255)
                                                    ggAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTc                                              >  1:1016326/1‑100 (MQ=255)
                                                    ggAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTc                                              >  2:101029/1‑100 (MQ=255)
                                                    ggAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTc                                              >  1:683356/1‑100 (MQ=255)
                                                    ggAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGATAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTc                                              >  1:1242547/1‑100 (MQ=255)
                                                     gAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCa                                             <  1:8437/100‑1 (MQ=255)
                                                         aCGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCAtctt                                         <  2:563225/100‑1 (MQ=255)
                                                          cGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCAtcttc                                        >  1:1056427/1‑100 (MQ=255)
                                                          cGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCAtcttc                                        <  1:1440441/100‑1 (MQ=255)
                                                            ttGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTg                                      <  2:1056428/100‑1 (MQ=255)
                                                               aGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTa                                    >  2:1059815/1‑99 (MQ=255)
                                                               aGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTAc                                   >  2:1059923/1‑100 (MQ=255)
                                                               aGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTAc                                   >  2:1733277/1‑100 (MQ=255)
                                                                gTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGCTTGCCAGCCATCCGGAAGTCATCTTCTGTAcc                                  >  1:1700689/1‑100 (MQ=255)
                                                                   tATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTACCttt                               >  1:1009194/1‑100 (MQ=255)
                                                                     tCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTACCTTTgg                             >  1:419694/1‑100 (MQ=255)
                                                                     tCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTACCTTTgg                             <  1:193595/100‑1 (MQ=255)
                                                                      ccACGGTCCGGGGTGCCCGGTGTGCGTACTGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTACCTTTGGc                            >  1:1852699/1‑100 (MQ=255)
                                                                      ccACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGCAGTCATCTTCTGTACCTTTGGc                            <  1:1131301/100‑1 (MQ=255)
                                                                      ccACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTACCTTTGGc                            <  2:1338323/100‑1 (MQ=255)
                                                                      ccACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTACCTTTGGc                            >  1:1129647/1‑100 (MQ=255)
                                                                       cACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTACCTTTGGCg                           <  1:1746368/100‑1 (MQ=255)
                                                                          ggTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTACCTTTGGCGACg                        <  1:1930716/100‑1 (MQ=255)
                                                                               ggggTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTACCTTTGGCGACGCCATg                   <  2:341403/100‑1 (MQ=255)
                                                                               ggggTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTACCTTTGGCGACGCCATg                   >  1:1517701/1‑100 (MQ=255)
                                                                                   tGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTACCTTTGGCGACGCCATGcgcg               >  2:425862/1‑100 (MQ=255)
                                                                                   tGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTACCTTTGGCGACGCCATGcgcg               <  2:1016327/100‑1 (MQ=255)
                                                                                    gCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTACCTTTGGCGACGCCATGCGCGt              >  2:1085893/1‑100 (MQ=255)
                                                                                    gCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTACCTTTGGCGACGCCATGCGCGt              >  2:767989/1‑100 (MQ=255)
                                                                                     cccGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTACCTTTGGCGACGCCATGCGCGTg             <  1:697375/100‑1 (MQ=255)
                                                                                       cGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTACCTTTGGCGACGCCATGCGCGTGcc           >  2:1862985/1‑100 (MQ=255)
                                                                                       cGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTACCTTTGGCGACGCCATGCGCGTGcc           <  2:118798/100‑1 (MQ=255)
                                                                                          tgtgCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTACCTTTGGCGACGCCATGCGCGTGCCggg        <  2:881232/100‑1 (MQ=255)
                                                                                              cGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTACCTTTGGCGACGCCATGCGCGTGCCGGGGaaa    >  1:1626035/1‑100 (MQ=255)
                                                                                                tACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTACCTTTGGCGACGCCATGCGCGTGCCGGGGAAACa  >  1:1468138/1‑100 (MQ=255)
                                                                                                tACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTACCTTTGGCGACGCCATGCGCGTGCCGGGGAAACa  >  2:1748484/1‑100 (MQ=255)
                                                                                                   |                                                                                                
GTGGCGGTCATACCCACGCTATCTTTAAATTCGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACTGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTACCTTTGGCGACGCCATGCGCGTGCCGGGGAAACA  >  NC_000913/2852257‑2852452

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