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,3560TA93.4% 197.6 / 5.0 61L74Q (CTG→CAG) hypDhydrogenase maturation protein
Reads supporting (aligned to +/- strand):  ref base T (3/1);  new base A (24/33);  total (27/34)
Fisher's exact test for biased strand distribution p-value = 3.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.44e-01

GGCGGTCATACCCACGCTATCTTTAAATTCGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACTGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTACCTTTGGCGACGCCATGCGCGTGCCGGGGAAACAGG  >  NC_000913/2852259‑2852454
                                                                                                 |                                                                                                  
ggCGGTCATACCCACGCTATCTTTAAATTCGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGc                                                                                                  >  1:570488/1‑100 (MQ=255)
  cGGTCATACCCACGCTATCTTTAAATTCGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCg                                                                                                <  1:1324564/100‑1 (MQ=255)
     tCATACCCACGCTATCTTTAAATTCGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATg                                                                                             <  1:266170/100‑1 (MQ=255)
       aTACCCACGCTATCTTTAAATTCGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATggg                                                                                           <  2:401681/100‑1 (MQ=255)
         aCCCACGCTATCTTTAAATTCGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCACGGTGTGCGTACAACCGATGGGTa                                                                                         >  2:1229336/1‑100 (MQ=255)
          cccACGCTATCTTTAAATTCGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAg                                                                                        <  1:560452/100‑1 (MQ=255)
           ccACGCTATCTTTAAATTCGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGa                                                                                       >  2:373402/1‑100 (MQ=255)
            cACGCTATCTTTAAATTCGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGaa                                                                                      >  1:208648/1‑100 (MQ=255)
            cACGCTATCGTTAAATTCGGCCTCGACCAGTTCCTGCCGGAAAACCTTGAGTTTACCCACGGTCCAGGGTGCCCGCTGTGCGTACACCCGATGGGGAg                                                                                        >  2:28874/1‑98 (MQ=255)
                 tATCTTTAAATTCGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGac                                                                                 <  1:948053/100‑1 (MQ=255)
                  aTCTTTAAATTCGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGaca                                                                                <  1:1168102/100‑1 (MQ=255)
                   tCTTTAAATTCGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGacac                                                                               >  1:563475/1‑100 (MQ=255)
                   tCTTTAAATTCGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGacac                                                                               <  1:1176146/100‑1 (MQ=255)
                            tCGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGga                                                                      <  2:321733/100‑1 (MQ=255)
                            tCGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGga                                                                      >  1:344841/1‑100 (MQ=255)
                             cGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGgag                                                                     >  2:931426/1‑100 (MQ=255)
                             cGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGgag                                                                     <  2:180365/100‑1 (MQ=255)
                             cGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGgag                                                                     <  2:190538/100‑1 (MQ=255)
                              ggCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGgaga                                                                    <  1:1394460/100‑1 (MQ=255)
                                   cGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACTGCCGATGGGTAGAATCGACACCTGCGTGGAGATTgcc                                                               >  2:75647/1‑100 (MQ=255)
                                     aCCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTgccag                                                             <  1:1293650/100‑1 (MQ=255)
                                      ccAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTGCAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTgccagc                                                            >  2:1110040/1‑100 (MQ=255)
                                      ccAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTgccagc                                                            <  2:1390978/100‑1 (MQ=255)
                                       cAGTTACTGCCGGAAAACGTTGAGTCTATCCACGGTCCGGGGCGCCCGGTGTGCGTACAGCCGAGGGGTAGAATCGACCCGTGCGTGGAGATTg                                                                 <  2:240825/94‑1 (MQ=255)
                                        aGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTgccagcca                                                          >  1:109393/1‑100 (MQ=255)
                                         gTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCAt                                                         >  1:902871/1‑100 (MQ=255)
                                           tACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATcc                                                       >  2:1229963/1‑100 (MQ=255)
                                            aCTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCg                                                      <  1:1119351/100‑1 (MQ=255)
                                            aCTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCg                                                      >  1:1302767/1‑100 (MQ=255)
                                             cTGCCGGAAATCGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGACTCGACACCTGCGTGGAGATTGCCAGCCATCCgg                                                     <  2:1035123/100‑1 (MQ=255)
                                              tGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGa                                                    >  2:1501265/1‑100 (MQ=255)
                                              tGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGa                                                    <  2:859829/100‑1 (MQ=255)
                                                 cGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGt                                                 >  1:97441/1‑100 (MQ=255)
                                                   gAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCa                                               <  2:667086/100‑1 (MQ=255)
                                                     aaaCGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCAtc                                             <  1:1185619/100‑1 (MQ=255)
                                                        cGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCAtcttc                                          <  1:1055165/100‑1 (MQ=255)
                                                             aGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTAc                                     >  1:1126617/1‑100 (MQ=255)
                                                             aGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTAc                                     >  1:1150695/1‑100 (MQ=255)
                                                                 tATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTACCttt                                 >  1:1078259/1‑100 (MQ=255)
                                                                   tCCACGGTCCGGGGTGCCCGGTGTGCGTACTGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTACCTTTgg                               <  1:75647/100‑1 (MQ=255)
                                                                   tCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTACCTTTgg                               <  1:1448889/100‑1 (MQ=255)
                                                                    ccACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGTAAGTCATCTTCTGTACCTTTgg                               <  2:920096/99‑1 (MQ=255)
                                                                    ccACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTACCTTTGGc                              <  2:1569472/100‑1 (MQ=255)
                                                                    ccACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTACCTTTGGc                              >  1:1369081/1‑100 (MQ=255)
                                                                     cACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTACCTTTGGCg                             >  2:767905/1‑100 (MQ=255)
                                                                     cACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTACCTTTGGCg                             <  1:1229968/100‑1 (MQ=255)
                                                                        ggTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTACCTTTGGCGACg                          >  2:74388/1‑100 (MQ=255)
                                                                            cGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTACCTTTGGCGACGCCAt                      >  1:260629/1‑100 (MQ=255)
                                                                                 tGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTACCTTTGGCGACGCCATGcgcg                 >  1:1569437/1‑100 (MQ=255)
                                                                                    ccGGTGTGCGTACTGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTACCTTTGGCGACGCCATGCGCGTGc              >  1:1194545/1‑100 (MQ=255)
                                                                                       gtgtGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTACCTTTGGCGACGCCATGCGCGTGCCgg           <  1:610297/100‑1 (MQ=255)
                                                                                       gtgtGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTACCTTTGGCGACGCCATGCGCGTGCCgg           <  1:212452/100‑1 (MQ=255)
                                                                                        tgtgCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTACCTTTGGCGACGCCATGCGCGTGCCggg          <  1:858512/100‑1 (MQ=255)
                                                                                        tgtgCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTACCTTTGGCGACGCCATGCGCGTGCCggg          <  2:244483/100‑1 (MQ=255)
                                                                                          tgCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTACCTTTGGCGACGCCATGCGCGTGCCGGGGa        <  2:1284934/100‑1 (MQ=255)
                                                                                          tgCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTACCTTTGGCGACGCCATGCGCGTGCCGGGGa        <  2:326928/100‑1 (MQ=255)
                                                                                          tgCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTACCTTTGGCGACGCCATGCGCGTGCCGGGGa        <  1:934819/100‑1 (MQ=255)
                                                                                           gCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTACCTTTGGCGACGCCATGCGCGTGCCGGGGaa       <  2:78419/100‑1 (MQ=255)
                                                                                           gCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTACCTTTGGCGACGCCATGCGCGTGCCGGGGaa       <  2:1203504/100‑1 (MQ=255)
                                                                                           gCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTACCTTTGGCGACGCCATGCGCGTGCCGGGGaa       >  1:1172010/1‑100 (MQ=255)
                                                                                                cTGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTACCTTTGGCGACGCCATGCGCGTGCCGGGGAAACAgg  >  1:203076/1‑100 (MQ=255)
                                                                                                 |                                                                                                  
GGCGGTCATACCCACGCTATCTTTAAATTCGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACTGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTACCTTTGGCGACGCCATGCGCGTGCCGGGGAAACAGG  >  NC_000913/2852259‑2852454

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: ‑