Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 313,191 C→T 10.8% C8Y (TGC→TAC)  fdhD ← protein sulfur transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947313,1910CT10.8% 106.1 / 3.7 37C8Y (TGC→TAC) fdhDprotein sulfur transferase
Reads supporting (aligned to +/- strand):  ref base C (16/17);  new base T (2/2);  total (18/19)
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.52e-01

TTGAGGCCGTTATAGACAATGGCCAGGGCCACTTCCTCGGCCAGCGGGGTGCTGTCGCGTTTGGCGTCGCTGAGCTGGACGTAGTCGTAGGTGTTGCTGGCGGCGGGCATGGACAAGGGCGATGAGGCCGCGCAGACCGGAGGCTTGCTGTTCATCGGGGCGCTGAATACCGGCGGCTTTCTGATAGCACAAGCCTAGGCGCGTGGGCCGGCGGCGTCTAATCGCTAGGGTCTATGCCTTGATAGACCGCGTCGATTGGTCGCGCAG  >  NC_002947/313059‑313325
                                                                                                                                    |                                                                                                                                      
ttGAGGCCGTTATAGACAATGGCCAGGGCCACTTCCTCGGCCAGCGGGGTGCTGTCGCGTTTGGCGTCGCTGAGCTGGACGTAGTCGTAGGTGTTGCTGGCGGCGGGCATGGACAAGGGCGATGAGGCCGCGCAGa                                                                                                                                     >  2:606743/1‑136 (MQ=255)
        gTTATAGACAATGGCCAGGGCCACTTCCTCGGCCAGCGGGGTGCTGTCGCGTTTGGCGTCGCTGAGCTGGACGTAGTCGTAGGTGTTGCTGGCGGCGGGCATGGACAAGGGCGATGAGGCCGCGCAGACCGGAGGc                                                                                                                             <  2:6591/136‑1 (MQ=255)
            taGACAATGGCCAGGGCCACTTCCTCGGCCAGCGGGGTGCTGTCGCGTTTGGCGTCGCTGAGCTGGACGTAGTCGTAGGTGTTGCTGGCGGCGGGCATGGACAAGGGCGATGAGGCCGCGTAGACCGGAGGCTTGc                                                                                                                         >  2:361755/1‑136 (MQ=255)
                      ccAGGGCCACTTCCTCGGCCAGCGGGGTGCTGTCGCGTTTGGCGTCGCTGAGCTGGACGTAGTCGTAGGTGTTGCTGGCGGCGGGCATGGACAAGGGCGATGAGGCCGCGCAGACCGGAGGCTTGCTGTTCATCgg                                                                                                               >  1:156202/1‑136 (MQ=255)
                      ccAGGGCCACTTCCTCGGCCAGCGGGGTGCTGTCGCGTTTGGCGTCGCTGAGCTGGACGTAGTCGTAGGTGTTGCTGGCGGCGGGCATGGACAAGGGCGATGAGGCCGCGCAGACCGGAGGCTTGCTGTTCATCgg                                                                                                               >  1:655031/1‑136 (MQ=255)
                       cAGGGCCACTTCCTCGGCCAGCGGGGTGCTGTCGCGTTTGGCGTCGCTGAGCTGGACGTAGTCGTAGGTGTTGCTGGCGGCGGGCATGGTCAAGGGCGATGAGGCCGCGCAGACCGGAGGCTTGCTGTTCATCggg                                                                                                              <  2:413813/136‑1 (MQ=255)
                       cAGGGCCACTTCCTCGGCCAGCGGGGTGCTGTCGCGTTTGGCGTCGCTGAGCTGGACGTAGTCGTAGGTGTTGCTGGCGGCGGGCATGGACAAGGGCGATGAGGCCGCGCAGACCGGAGGCTTGCTGTTCATCggg                                                                                                              <  2:310252/136‑1 (MQ=255)
                             cACTTCCTCGGCCAGCGGGGTGCTGTCGCGTTTGGCGTCGCTGAGCTGGACGTAGTCGTAGGTGTTGCTGGCGGCGGGCATGGACAAGGGCGATGAGGCCGCGCAGACCGGAGGCTTGCTGTTCATCGGGGCGCTg                                                                                                        >  1:443871/1‑136 (MQ=255)
                             cACTTCCTCGGCCAGCGGGGTGCTGTCGCGTTTGGCGTCGCGGAGCTGGACGTAGTCGTAGGTGTTGCTGGCGGCGGGCCTGGACAACGGCCATGAGGCCGCGCAGACCGGAGGCGTGCTGTTCATCGGGGCGCTg                                                                                                        >  1:533450/1‑136 (MQ=255)
                                  ccTCGGCCAGCGGGGTGCTGTCGCGTTTGGCGTCGCTGAGCTGGACGTAGTCGTAGGTGTTGCTGGCGGCGGGCATGGACAAGGGCGATGAGGCCGCGCAGACCGGAGGCTTGCTGTTCATCGGGGCGCTGAATAc                                                                                                   >  2:536520/1‑136 (MQ=255)
                                  ccTCGGCCAGCGGGGTGCTGTCGCGTTTGGCGTCGCTGAGCTGGACGTAGTCGTAGGTGTTGCTGGCGGCGGGCATGGACAAGGGCGATGAGGCCGCGCAGACCGGAGGCTTGCTGTTCATCGGGGCGCTGAATAc                                                                                                   >  1:566779/1‑136 (MQ=255)
                                   cTCGGCCAGCGGGGTGCTGTCGCGTTTGGCGTCGCTGAGCTGGACGTAGTCGTAGGTGTTGCTGGCGGCGGGCATGGACAAGGGCGATGAGGCCGCGCAGACCGGAGGCTTGCTGTTCATCGGGGCGCTGAATAcc                                                                                                  <  2:158204/136‑1 (MQ=255)
                                   cTCGGCCAGCGGGGTGCTGTCGCGTTTGGCGTCGCTGAGCTGGACGTAGTCGTAGGTGTTGCTGGCGGCGGGCATGGACAAGGGCGATGAGGCCGCGCAGACCGGAGGCTTGCTGTTCATCGGGGCGCTGAATAcc                                                                                                  <  1:606743/136‑1 (MQ=255)
                                   cTCGGCCAGCGGGGTGCTGTCGCGTTTGGCGTCGCTGAGCTGGACGTAGTCGTAGGTGTTGCTGGCGGCGGGCATGGACAAGGGCGATGAGGCCGCGCAGACCGGAGGCTTGCTGTTCATCGGGGCGCTGAATAcc                                                                                                  <  1:98969/136‑1 (MQ=255)
                                      ggCCAGCGGGGTGCTGTCGCGTTTGGCGTCGCTGAGCTGGACGTAGTCGTAGGTGTTGCTGGCGGCGGGCATGGACAAGGGCGATGAGGCCGCGCAGACCGGAGGCTTGCTGTTCATCGGGGCGCTGAATACcggc                                                                                               <  2:655031/136‑1 (MQ=255)
                                       gCCAGCGGGGTGCTGTCGCGTTTGGCGTCGCTGAGCTGGCCGTAGTCGTAGGGGTTGCTGGCGGCGGGCATGGACAAGGGCGATGAGGCCGCGCAGACCGGAGGGTTGGGGTTGATCCGGGCGCGGAAAACcggcg                                                                                              >  2:248786/1‑136 (MQ=255)
                                        ccAGCGGGGTGCTGTCGCGTTTGGCGTCGCTGAGCTGGACGTAGTCGTAGGTGTTGCTGGCGGCGGGCATGGACAAGGGCGATGAGGCCGCGCAGACCGGAGGCTTGCTGTTCATCGGGGCGCTGAATACcggcgg                                                                                             >  1:464487/1‑136 (MQ=255)
                                        ccAGCGGGGTGCTGTCGCGTTTGGCGTCGCTGAGCTGGACGTAGTCGTAGGTGTTGCTGGCGGCGGGCATGGACAAGGGCGATGAGGCCGCGCAGACCGGAGGCTTGCTGTTCATCGGGGCGCTGAATACcggcgg                                                                                             >  1:7578/1‑136 (MQ=255)
                                                   cTGTCGCGTTTGGCGTCGCTGAGCTGGACGTAGTCGTAGGTGTTGCTGGCGGCGGGCATGGACAAGGGCGATGAGGCCGCGCAGACCGGAGGCTTGCTGTTCATCGGGGCGCTGAATACCGGCGGCTTTCTGATAg                                                                                  <  1:471244/136‑1 (MQ=255)
                                                   cTGTCGCGCTTGGCGTCGCTGAGCTGGACGTAGTCGTAGGTGTTGCTGGCGGCGGGCATGGACAAGGGCGATGAGGCCGCGCAGACCGGAGGCTTGCTGTTCATCGGGGCGCTGAATACCGGCGGCTTTCTGATAg                                                                                  <  1:440400/136‑1 (MQ=255)
                                                     gTCGCGTTTGGCGTCGCTGAGCTGGACGTAGTCGTAGGTGTTGCTGGCGGCGGGCATGGACAAGGGCGATGAGGCCGCGTAGACCGGAGGCTTGCTGTTCATCGGGGCGCTGAATACCGGCGGCTTTCTGATAGca                                                                                >  2:56316/1‑136 (MQ=255)
                                                              ggCGTCGCTGAGCTGGACGTAGTAGTAGGTGTTGCTGGCGGCGGGCATGGACAAGGGCGATGAGGCCGCGCAGACCCGAGGCTTGCTGTTAATAGGGGCGCTGAATACCGGCGGCTTTATTATATCACAAAccaaa                                                                       >  2:229789/1‑133 (MQ=255)
                                                                cGTCGCTGAGCTGGACGTAGTCGTAGGTGTTGCTGGCGGCGGGCATGGACAAGGGCGATGAGGCCGCGCAGACCGGAGGCTTGCTGTTCATCGGGGCGCTGAATACCGGCGGCTTTCTGATAGCACAAGCCTAGgc                                                                     <  2:33186/136‑1 (MQ=255)
                                                                                 tAGTCGTAGGTGTTGCTGGCGGCGGGCATGGACAAGGGCGATGAGGCCGCGCAGACCGGAGGCTTGCTGTTCATCGGGGCGCTGAATACCGGCGGCTTTCTGATAGCACAAGCCTAGGCGCGTGGGCCGGCGGCGt                                                    <  2:566779/136‑1 (MQ=255)
                                                                                       tAGGTGTTGCTGGCGGCGGGCATGGACAAGGGCGATGAGGCCGCGCAGACCGGAGGCTTGCTGTTCATCGGGGCGCTGAATACCGGCGGCTTTCTGATAGCACAAGCCTAGGCGCGTGGGCCGGCGGCGTCTAATc                                              <  1:536520/136‑1 (MQ=255)
                                                                                        aGGTGTTGCTGGCGGCGGGCATGGACAAGGGCGATGAGGCCGCGTAGACCGGAGGCTTGCTGTTCATCGGGGCGCTGAATACCGGCGGCTTTCTGATAGCACAAGCCTAGGCGCGTGGGCCGGCGGCGTCTAATCg                                             <  1:56316/136‑1 (MQ=255)
                                                                                          gtgtTGCTGGCGGCGGGCATGGACAAGGGCGATGAGGCCGCGCAGACCGGAGGCTTGCTGTTCATCGGGGCGCTGAATACCGGCGGCTTTCTGATAGCACAAGCCTAGGCGCGTGGGCCGGCGGCGTCTAATCGCt                                           >  1:646378/1‑136 (MQ=255)
                                                                                          gtgtTGCTGGCGGCGGGCATGGACAAGGGCGATGAGGCCGCGCAGACCGGAGGCTTGCTGTTCATCGGGGCGCTGAATACCGGCGGCTTTCTGATAGCACAAGCCTAGGCGCGTGGGCCGGCGGCGTCTAATCGCt                                           >  1:585721/1‑136 (MQ=255)
                                                                                               gCTGGCGGCGGGCATGGACAAGGGCGATGAGGCCGCGCAGACCGGAGGCTTGCTGTTCATCGGGGCGCTGAATACCGGCGGCTTTCTGATAGCACAAGCCTAGGCGCGTGGGCCGGCGGCGTCTAATCGCTAGGGt                                      <  2:7578/136‑1 (MQ=255)
                                                                                               gCTGGCGGCGGGCATGGACAAGGGCGATGAGGCCGCGCAGACCGGAGGCTTGCTGTTCATCGGGGCGCTGAATACCGGCGGCTTTCTGATAGCACAAGCCTAGGCGCGTGGGCCGGCGGCGTCTAATCGCTAGGGt                                      <  2:464487/136‑1 (MQ=255)
                                                                                               gCTGGCGGCGGGCATGGACAAGGGCGATGAGGCCGCGCAGACCGGAGGCTTGCTGTTCATCGGGGCGCTGAATACCGGCGGCTTTCTGATAGCACAAGCCTAGGCGCGTGGGCCGGCGGCGTCTAATCGCTAGGGt                                      <  2:156202/136‑1 (MQ=255)
                                                                                                 tGGCGGCGGGCATGGACAAGGGCGATGAGGCCGCGCAGACCGGAGGCTTGCTGTTCATCGGGGCGCTGAATACCGGCGGCTTTCTGATAGCACAAGCCTAGGCGCGTGGGCCGgcggcg                                                     <  2:512945/119‑1 (MQ=255)
                                                                                                 tGGCGGCGGGCATGGACAAGGGCGATGAGGCCGCGCAGACCGGAGGCTTGCTGTTCATCGGGGCGCTGAATACCGGCGGCTTTCTGATAGCACAAGCCTAGGCGCGTGGGCCGgcggcg                                                     >  1:512945/1‑119 (MQ=255)
                                                                                                        gggCATGGACAAGGGCGATGAGGCCGCGCAGACCGGAGGCTTGCTGTTCATCGGGGCGCTGAATACCGGCGGCTTTCTGATAGCACAAGCCTAGGCGCGTGGGCCGGCGGCGTCTAATCGCTAGGGTCTATGCCtt                             >  2:355906/1‑136 (MQ=255)
                                                                                                        gggCATGGACAAGGGCGATGAGGCCGCGCAGACCGGAGGCTTGCTGTTCATCGGGGCGCTGAATACCGGCGGCTTTCTGATAGCACAAGCCTAGGCGCGTGGGCCGGCGGCGTCTAATCGCAAGGGTGTATGCCtt                             >  2:225401/1‑136 (MQ=255)
                                                                                                              ggACAAGGGCGATGAGGCCGCGTAGACCGGAGGCTTGCTGTTCATCGGGGCGCTGAATACCGGCGGCTTTCTGATAGCACAAGCCTAGGCGCGTGGGCCGGCGGCGTCTAATCGCTAGGGTCTATGCCTTGATAGa                       <  1:361755/136‑1 (MQ=255)
                                                                                                                             ggCCGCGCAGACCGGAGGCTTGCTGTTCATCGGGGCGCTGAATACCGGCGGCTTTCTGATAGCACAAGCCTAGGCGCGTGGGCCGGCGGCGTCTAATCGCTAGGGTCTATGCCTTGATAGACCGCGTCGATTGGTc        <  2:646378/136‑1 (MQ=255)
                                                                                                                                 gcgcAGACCGGAGGCTTGCTGTTCATCGGGGCGCTGAATACCGGCGGCTTTCTGATAGCACAAGCCTAGGCGCGTGGGCCGGCGGCGTCTAATAGCTAGGGTGTATTGCTTGATATACCGCGTCGATTGGTagcgg    >  1:405211/1‑135 (MQ=255)
                                                                                                                                   gcAGACCGGAGGCTTGCTGTTCATCGGGGCGCTGAATACCGGCGGCTTTCTGATAGCACAAGCCTAGGCGCGTGGGCCGGCGGCGTCTAATCGCTAGGGTCTATGCCTTGATAGACCGCGTCGATTGGTCGCGCAg  >  2:374238/1‑136 (MQ=255)
                                                                                                                                    |                                                                                                                                      
TTGAGGCCGTTATAGACAATGGCCAGGGCCACTTCCTCGGCCAGCGGGGTGCTGTCGCGTTTGGCGTCGCTGAGCTGGACGTAGTCGTAGGTGTTGCTGGCGGCGGGCATGGACAAGGGCGATGAGGCCGCGCAGACCGGAGGCTTGCTGTTCATCGGGGCGCTGAATACCGGCGGCTTTCTGATAGCACAAGCCTAGGCGCGTGGGCCGGCGGCGTCTAATCGCTAGGGTCTATGCCTTGATAGACCGCGTCGATTGGTCGCGCAG  >  NC_002947/313059‑313325

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