Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,750,629 A→G 16.9% Y129H (TAT→CAT)  ydhX ← putative 4Fe‑4S ferredoxin‑type protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,750,6290AG16.9% 39.7 / 4.1 25Y129H (TAT→CAT) ydhXputative 4Fe‑4S ferredoxin‑type protein
Reads supporting (aligned to +/- strand):  ref base A (10/10);  new base G (2/2);  total (13/12)
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.49e-01

TTCCGGGCAGGCGCTGACGCAAATGGGCGGAAAGCCTTTTGCCAGTCGGGACTCAGCGCAGAAATCACATTTATCTGCCACTTTAGTGACAGGGTTAAGGTAGCGAACCTGGTACGGACAAGCGCCAATGCAATAACTGCAGCCGATGCACTGCGATTTTTCTACCCGCACGATCCCCTGCTCATCACGCCACGACGCCCCCGTCGGGCAAACGTCGATGCACGGAGCATCTTCGCAGTGCTGACATGACTGGCGAAAAAAGTGATACTGCGT  >  NC_000913/1750495‑1750767
                                                                                                                                      |                                                                                                                                          
ttCCGGGCAGGCGCTGACGCAAATGGGCGGAAAGCCTTTTGCCAGTCGGGACTCAGCGCAGAAATCACATTTATCTGCCACTTTAGTGACAGGGTTAAGGTAGCGAACCTGGTACGGACAAGCGCCAATGCAATAACTGCAGc                                                                                                                                    <  2:79940/143‑1 (MQ=255)
    gggCAGGCGCTGACGCAAATGGGCGGAAAGCCTTTTGACAGTCGAGACTCAGCGCAGAGATCACATTTATCTGCCACTTTAGTGACAGGGTTAAGGTAGCGAACCTGGTAGGGACAAGGGCCAAGGCAATAACTGACGCCgca                                                                                                                                >  2:291611/1‑141 (MQ=255)
             cggACGCAAATGGGCGGAAATCGTTTTGCCAGTCGGGACTCAGCGCAGAAATCACATTTATCTCCCACTTGAGTGACAGGGTTAAGGCAGCGAACCTGGTACGGACAAGCGCCAATCCAATAACTGCAGCCGATGCACTGCgg                                                                                                                       <  2:347733/141‑2 (MQ=255)
                  gCAAATGGGCGGAAAGCCTTTTGCCAGTCGGGACTCAGCGCAGAAATCACATTTATCTGCCACTTTAGTGACAGGGTTAAGGTAGCGAACCTGGTACGGGCAAGCGCCAATGCAATAACTGCAGCCGATGCACTGCGAttttt                                                                                                                  >  1:202417/1‑143 (MQ=255)
                          gCGGAAAGCCTTTTGCCAGTCGGGACTCAGCGCAGAAATCACATTTATCTGCCACTTTAGTGACAGGGTTAAGGTAGCGAACCTGGTACGGACAAGCGCCAATGCAATGACTGCAGCCGATGCACTGCGATTTTTCTACCCGc                                                                                                          >  2:212261/1‑143 (MQ=255)
                             gAAAGCCTTTTGCCAGTCGGGACTCAGCGCAGAAATCACATTTCTCTGCCACTTTAGTGACAGGGTTTAGGTCGCGGACCCGCCCCACAAATGGACTAAGTTTATAACTGCAGCCGATGCACTCCGACTTTTCTCCCCGCACg                                                                                                       >  2:75187/1‑143 (MQ=255)
                             gAAAGCCTTTTGCCAGTCGGGACTCAGCGCAGAAATCACATTTATCTGCCACTTTAGTGACAGGGTTAAGGTAGCGAACCTGGTACGGACAAGCGCCAATGCAATAACTGCAGCCGATGCACTGCGATTTTTCTACCCGCACg                                                                                                       >  1:293013/1‑143 (MQ=255)
                                   cTTTTGCCAGTCGGGACTCAGCGCAGAAATCACATTTATCTGCCACTTTAGTGACAGGGTTAAGGTAGCGAACCTGGTACGGACAAGCGCCAATGCAATAACTGCAGCCGATGCACTGCGCTTTTTCTACCCGCACGATcccc                                                                                                 >  1:52477/1‑143 (MQ=255)
                                   cTTTTACGAGTCGGGACACAGCGAAGAAAGCACATGTAGCGGCCACTGTAGTGACAGGGTTAAGGTAGCGAACCTGGTACGGACAAGCGCCAATGCAATAACTGCAGCCGATGCACTGCGATTTTTCTACCCGCACGATcccc                                                                                                 <  1:133986/143‑1 (MQ=255)
                                                           aGAAATCACATTTATCTGCCACTTTAGTGACAGGGTTAAGGTAGCGAACCTGGTACGGACAAGCGCCAATGCAATGACTGCAGCCGATGCACTGCGATTTTTCTACCCGCACGATCCCCTGACCGACACCCCACGCCGccccc                                                                         >  2:170942/1‑143 (MQ=255)
                                                                 cacaTTTATCTGCCACTTTAGTGACAGGGTTAAGGTAGCGAACCTGGTACGGACAAGCGCCAATGCAATAACTGCAGCCGATGCACTGCGATTTTTCTACCCGCACGATCCCCTGCTCATCACGCCACGACGACCCCGTCggg                                                                   >  1:288907/1‑143 (MQ=255)
                                                                      ttATCTGCCACTTTAGTGACAGGGTTAAGGTAGCGAACCTGGTACGGACAAGCGCCAATGCAATGACTGCAGCCGATGCACTGCGATTTTTCTACCCGCACGATCCCCTGCTCATCACGCCACGACGCCCCCGTCGGGCAAAc                                                              <  1:170942/143‑1 (MQ=255)
                                                                          cTGCCACTTTAGTGACAGGGTTAAGGTAGCGAACCTGGTACGGACAAGCGCCAATTCTTTATCTGCAGGCGACGCACTCCGACGTGTCTTCCCCCCCGATCCCCTGCCCATCCCGCCACGACGCCCCCGTCGGGCAAACGTCg                                                          <  2:164101/143‑1 (MQ=255)
                                                                               aCTTTAGTGACAGGGTTAAGGTAGCGAACCTGGTACGGACAAGCGCCAATGCAATAACTGCAGCCGATGCACTGCGATTTTTCTACCCGCACGATCCCCTGCTCATCACGCCACGACGCCCCCGTCGGGCAAACGTCGATGCa                                                     <  2:52477/143‑1 (MQ=255)
                                                                               aCTTTAGTGACAGAGTTAAGGTAGCGAACCTGGTACGGACAAGCGCCAATGCAATAACTGCAGCCGATGCACTGCGATTTTCCTACCCGCACGATCCCCTGCTCATCACGCAACGACGCCCCCGGCGGGCAAAAGTCGATGcc                                                     >  2:25037/1‑142 (MQ=255)
                                                                                 tgtcgtgtCAGGGTTATGGCATCGAACCTGGTACGGACAAGCGCCAATGCAATACCTGCCGCCGATGCACTGCGATTTTTCTACCCTCACGATCCCCTGCTCATCACGCCACGACGCCCCCGTCGGGCAAACGTCGATGCACg                                                   <  2:187506/139‑1 (MQ=255)
                                                                                    aGTGACAGGGTTAAGGTAGCGAACCTGGTACGGACAAGCGCCAATGCAATGACTGCAGCCGATGCACTGCGATTTTTCTACCCGCACGATCCCCTGCTCATCACGCCACGACGCCCCCGTCGGGCAAACGTCGATGCACGGAg                                                <  1:212261/143‑1 (MQ=255)
                                                                                      tGACAGGGTTAAGGTAGCGAACATGGTACGGACAAGCGCCAATGCAGTAACTGCAGCCGATGCACTGCGATTTTTCTACCCGCACGACCCCCTGCTCATCACGCCACGACGCCCCCGCAGAGCACACGTCGATGCACgcagca                                              >  2:56229/1‑143 (MQ=255)
                                                                                          aGGGTTAAGGTAGCGAACCTGGTACGGACAAGCGCCAATGCAATAACTGCAGCCGATGAACTACGATTTTTCTACCCGCACGATCCCCTGCTCATCATGCCACGACGCCCCCGTCGGGCAAACGTCGATGCACGGAGCATCtt                                          <  2:125814/143‑1 (MQ=255)
                                                                                              ttAAGGTAGCGAACCTGGTACGGACAAGCGCCAATGCAATAACTGCAGCCGATGCACTGCGATTTTTCTACCCGCACGATCCCCTGCTCATCACGCCACGACGCCCCCGTCGGGCAAACGTCGATGCACGGAGCATCTTCGCa                                      >  1:324981/1‑143 (MQ=255)
                                                                                                  ggTAGCGAACCTGGTACGGACAAGCGCCAATGCAATAACTGCAGCCGATGCACTGCGATTTTTCTACCCGCACGATCCCCTGCTCATCACGCCACGACGCCCCCGTCGGGCAAACGTCGACGCACGGACCACCTTCGCCGTGc                                  >  2:167812/1‑143 (MQ=255)
                                                                                                   gTAGCGAACCTGGTACGGACAAGAGCCAATCCAGTTACGGCATCCGATGCACTGCGATTTTTCTACCCGCACGATCCCCTGCTCATCACGCCACGACGCCCCCGTCGGGCAAACGTCGATGCACGGAGCATCTTCGCAGTGCt                                 >  1:249370/1‑143 (MQ=255)
                                                                                                                        aGCGCCAATGCAATAACTGCAGCCGATGCACTGCGATTTTTCTACCCGCACGATCCCCTGCTCATCACGCCACGACGCCCCCGTCGGGCAAACGTCGATGCACGGAGCATCTTCGCAGTGCTGACATGACTGGCGAAAAAAGt            <  1:241082/143‑1 (MQ=255)
                                                                                                                            ccAATGCAATAACTGCAGCCGTTGCACTACCAGTTTTGTGTCCCCACCACCCCCTGCTCATCACGTCACGACGCCCCCGGCGGGCAAACGTCGATGCACGGAGCATCTTCGCAGTGCTGACATGACTGGCGAAAAAAGTGATa        <  1:262793/143‑1 (MQ=255)
                                                                                                                                  ccaTAACTGCGCCCCATGCACCGCGATTTTTCTCGCCGGTTGTTCCCCTGCTCATCTCGCCACGACGCCCCCGCCTGGTAAACGTCGATGCACTGAGCATCTTCGCAGTGCTGACATGACTGGCGAAAAAAGTGATACTGCGt  <  2:97849/141‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                          
TTCCGGGCAGGCGCTGACGCAAATGGGCGGAAAGCCTTTTGCCAGTCGGGACTCAGCGCAGAAATCACATTTATCTGCCACTTTAGTGACAGGGTTAAGGTAGCGAACCTGGTACGGACAAGCGCCAATGCAATAACTGCAGCCGATGCACTGCGATTTTTCTACCCGCACGATCCCCTGCTCATCACGCCACGACGCCCCCGTCGGGCAAACGTCGATGCACGGAGCATCTTCGCAGTGCTGACATGACTGGCGAAAAAAGTGATACTGCGT  >  NC_000913/1750495‑1750767

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: