Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913_3_pae_tpiA 2,355,843 T→C 9.4% F108S (TTC→TCC)  glpC → anaerobic sn‑glycerol‑3‑phosphate dehydrogenase, C subunit, 4Fe‑4S iron‑sulfur cluster

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_pae_tpiA2,355,8430TC 9.4% 100.1 / 2.6 42F108S (TTC→TCC) glpCanaerobic sn‑glycerol‑3‑phosphate dehydrogenase, C subunit, 4Fe‑4S iron‑sulfur cluster
Reads supporting (aligned to +/- strand):  ref base T (21/17);  new base C (2/2);  total (23/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 = 3.49e-01

GCAAACGTTGTGAAGTCGCCTGCCCGTCCGATGTGAAGATTGGCGATATTATCCAGCGCGCGCGGGCGAAATATGACACCACGCGCCCGTCGCTGCGTAATTTTGTGTTGAGTCATACCGACCTGATGGGTAGCGTTTCCACGCCGTTCGCACCAATCGTCAACACCGCTACCTCGCTGAAACCGGTGCGGCAGCTGCTTGATGCGGCGTTAAAAATCGATCATCGCCGCACGCTACCGAAATACTCCTTCGGCACGTTCCGTCGCTGGTATCGCAGCGTGGCGGCTCAGCAAGC  >  NC_000913_3_pae_tpiA/2355696‑2355990
                                                                                                                                                   |                                                                                                                                                   
gCAAACGTTGTGAGTTCGCCGGCCCGGCCGATGTGAAGATTGGCGATATTATCCAGCGCGCGCGGGAGAAATATGAAACCACGCGCCAGTCGCTGTGTAATTTTGTGTTGAGTCATACCGATCTGATGGGTAGCGTTTCCACGCCGTTc                                                                                                                                                    <  2:468171/149‑1 (MQ=255)
  aaaCGTTGTGAAGTCGCCTGCCCGTCCGATGTGAAGATTGGCGATATTATCCAGCGCGCGCGGGCGAAATATGACACCACGCGCCCGTCGCTGCGTAATTTTGTGTTGCGTAATAGCGACGTGATGGGTAGCGTTTCCACGCCGTTCGc                                                                                                                                                  >  1:635916/1‑149 (MQ=255)
  aaaCGTTGTGAAGTCGCCTGCCCGTCCGATGTGAAGATTGGCGATATTATCCAGCGCGCGCGGGCGAAATATGACACCACGCGCCCGTCGCTGCGTAATTTTGTGTTGAGTCATACCGACCTGATGGGTAGCGTTTCCACGCCGTCCGc                                                                                                                                                  >  2:353320/1‑149 (MQ=255)
  aaaCGTTGTGAAGTCGACTGCCCGTCCGATGTGAAGATTGGCGATATTATCCAGCGCGCGCGGGCGAAATATGACACCACGCGCCCGTCGCTGCGTAATTTTGTGTTGAGTCATACCGACCTGATGGGTAGCGTTTCCACGCCGTCCGc                                                                                                                                                  >  2:352877/1‑149 (MQ=255)
   aaCGTTGTGAATTCGCATGCCCGTCCGATGGGAAGATTGGCGATATTATCCAGCGCGCGCGGGCGAAATATGACCCCACGCGCCCGTCGCTGCGTAATTTTGTGTTGAGTCATACCGACCTGATGGGTAGCGTTTCCACGCCGTTCGCa                                                                                                                                                 <  2:394765/149‑1 (MQ=255)
   aaCGTTGTGAAGTCGCCTGCCCGTCCGATGTGAAGATTGGCGATATTATCCAGCGCGCGCGGGCGAAATATGACACCACGCGCCCGTCGCTGCGTAATTTTGTGTTGAGTCATACCGACCTGAAGGGTAGCGTTTCCACGCCGTTCGCa                                                                                                                                                 >  1:140028/1‑149 (MQ=255)
      gTTGTGAAGTCGCCTGCACGTCCGATGTGAAGATTGGCGATATTATCCAGCGCGCGCGGGCGAAATATGACACCACGCGCCCGTCGCTGCGTAATTTTGTGTTGAGTCATACCGACCTGATGGGTGGCGTTTCCACGCCGTTCGCAACa                                                                                                                                              >  2:321916/1‑149 (MQ=255)
       ttGTGAAGTCGCCTGCCCGTCCGATGTGAAGATTGGCGATATTATCCAGCGCGCGCGGGCGAAATATGACACCACGCGCCCGTCGCTGCGTAATTTTGTGTTGAGTCATACCGACCTGATGGGTAGCGTTTCCACGCCGTTCGCACCaa                                                                                                                                             >  1:56856/1‑149 (MQ=255)
         gtgGAGTCGCCTGCCCGTCCGATGTGACGATTGGCGATATTATCCAGCGCGCGCGGGCGAAATATGACACCACGCGCCCGTCGCTGCGTAATTTTGTGTTGAGTCATACCGACCTGATGGGTAGCGTTTCCACGCCGTTCGCACCAATc                                                                                                                                           <  1:99598/149‑1 (MQ=255)
         gtgAAGTCGCCTGCCCGTCCGATGTGAAGATTGGCGATATTATCCAGCGCGCGCGGGCGAAATATGACACCACGCGCCCGTCGCTGCGTAATTTTGTGTTGAGTTATACCGACCTGATGGGTAGCGTTTCCACGCCGTTCGCACCAATc                                                                                                                                           <  1:99793/149‑1 (MQ=255)
                    tGCCCGTCCGATGTGAAGATTGGCGATATTATCCAGCGCGCGCGGGCGAAATATGACACCACGCGCCCGTCGCTGCGTAATTTTGTGTTGAGTCATACCGACCTGATGGGTAGCGTTTCCACGCCGTTCGCACCAATCGTCAACACGGc                                                                                                                                >  1:341418/1‑149 (MQ=255)
                            cGATGTGAAGATTGGCGATATTATCCAGCGCGCGCGGGCGAAATATGACACCACGCGCCCGTCGCTGCGTAATTTTGTGTTGAGTCATACCGACCTGATGGGTAGCGTTTCCACGCCGTTCGCACCAATCATCAACACCGCTACCCCGc                                                                                                                        >  2:700446/1‑149 (MQ=255)
                             gATGTGAAGATTGGCGATATTATCCAGCGCGCGCGGGCGAAATATGACACCACGCGCCCGTCGCTGCGTAATTTTGTGTTGAGTCATACCGACCTGATGGGTAGCGTTTCCACGCCGTTCGCACCAATCGTCAACACCGCTACCTCGCt                                                                                                                       <  2:120378/149‑1 (MQ=255)
                                        tGGCGATATTATCCAGCGCGCGCGGGCGAAATATGACACCACGCGCCCGTCGCTGCGTAATTTTGTGTTGAGTCATACCGACCTGATGGGTAGCGTTTCCACGCCGTTCGCACCAATCGTCAACACCGCTACCTCGCTGAAACCGGTGc                                                                                                            <  1:584619/149‑1 (MQ=255)
                                          gCGATATTATCCAGCGCGCGCGGGCGAAATATGACACCACGCGCCCGTCGCTGCGTAATTTTGTGTTGAGTCATACCGACCTGATGGGTAGCGTTTCCACGCCGTTCGCACCAATCGTCAACACCGCTACCTCGCCGAAACAGCTGCgg                                                                                                          >  2:32862/1‑149 (MQ=255)
                                             atatTATCCTGCGCGCGCGGGCGAAATATGACACCACGCGCCCGTCGCTGCGTAATTTTGTGTTGAGTCATACCGACCTGATGGGTAGCGTTTCCACGCCGTTCGCACCAATCGTCAACACCGCTACCTCGCTGAAACCGGTGCGGCAg                                                                                                       <  1:630206/149‑1 (MQ=255)
                                                 tgtCCAGCGCGCCCGGGCGAAATATGACACCACGCGCCCGTCGCTGCGTAATTTTGTGTTGAGTCATACCGACCTGATGGGTAGCGTTTCCACGCCGTTCGCACCAATCGTCAACACCGCTACCTCGCTGAAACCGGTGCGGCAgctgc                                                                                                   <  1:476137/147‑1 (MQ=255)
                                                      aGCGCGCGCGGGCGAAATATGACACCACGCGCCCGTCGCTGCGTAATTTTGTGTTGAGTCATACCGACCTGATGGGTAGCGTTTCCACGCCGTTCGCACCAATCGTCAACACCGCTACCTCGCTGAAACCGGTGCGGCAGCTGCTTGAt                                                                                              >  1:550569/1‑149 (MQ=255)
                                                                   gAAATATGACACCACGCGCCCGTCGCTGCGTAATTTTGTGTTGAGTCATACCGACCTGATGGGTAGCGTTTCCACGCCGTTCGCACCAATCGTCAACACCGCTACCTCGCTGAAACCGGTGCGGCAGCTGCTTGATGCGGCGTTaaaaa                                                                                 >  1:711713/1‑149 (MQ=255)
                                                                   gAAATATGACACCACGCGCCCGACGCTGCGTAATTTTGTGATAAATCATACCTACCTGATGGGTAGCGTTGCCACGCCGTTCGCACCAACCGTCAACACCGCTACCTCGCTGAACCCCGTCCGGCAGCTGCTTGATGCGGCGTAgaaaa                                                                                 >  2:157971/1‑149 (MQ=255)
                                                                       tatGACACCACGCGCCCGTCGCTGCGTAATTTTGTGTTGAGTCATACCGCCCTGATGGGTAGCGTTTCCACGCCGTTCGCACCAATCGTCAACACCGCTACCTCGCTGAAACCGGTGCGGCAGCTGCTTGATGCGGCGTTAAAAATCGa                                                                             >  1:695306/1‑149 (MQ=255)
                                                                        catgcgccagGCGCCCGTCGTCGCTTAATTTTGTGTTGAGTCCTACCGACCTGATGGGTAGCGTTTCCACGCCGTTCGCACCAATCGTCAACACCGCTACCTCGCTGAAACCGGTGCGGCAGCTGCTTGATGCGGCGTTAAAAATCGat                                                                            <  2:416440/143‑1 (MQ=255)
                                                                            cacctcGCTCGCGGCGCTGCGTTATTTTGTGGTGTGTCATACCGGCCTGCTGGGTAGCGTTTCCCCGCCGTTCGCGCCAATCGTCATCACCGCTGCCTCGCTGGAACCGGTGCGGCAGCTGCTTGATGCGGCGTTAAAAATCGatcatc                                                                        <  2:77620/149‑1 (MQ=255)
                                                                              ccacGCGCCCGTCGCTGCGTAATTTTGTGTTGAGTCATACCGACCTGATGGGTAGCGTTTCCACGCCGTTCGCACCAATCGTCAACACCGCTACCTCGCTGAAACCGGTGCGGCAGCTGCTTGATGCGGCGTTAAAAATCGATCATcgc                                                                      <  2:757241/149‑1 (MQ=255)
                                                                                acGCGCCCGTCGCTGCGTAATTTTGTGTTGAGTCATACCGACCTGATGGGTAGCGTTTCCACGCCGTTCGCACCAATCGTCAACACCGCTACCTCGCTGAAACCTGTGCGGCAGCTGCTTGATGCGGCGTTAAAAATCGATCATcgccg                                                                    >  1:605248/1‑149 (MQ=255)
                                                                                acGCGCCCGTCGCTGCGTAATTTTGTGTTGAGTCATACCGACCTGATGGGTAGCGTTTCCACGCCGTTCGCACCAATCGTCAACACCGCTACCTCGCTGAAACCGGTGCGGCAGCTGCTTGATGCGGCGTTAAAAAACGATAATcgccg                                                                    >  2:10661/1‑149 (MQ=255)
                                                                                         tCGCTGCGTAATTTTGTGTTGAGTCATACCGACCTGATGGGTAGCGTTTCCACGCCGTTCGCACCAATCGTCAACACCGCTACCTCGCTGAAACCGGTGCGGCAGCTGCTTGATGCGGCGTTATAAATCGATCATCGCCGCACGCTAcc                                                           >  1:831595/1‑149 (MQ=255)
                                                                                         tCGCTGCGTAATTTTGTGTTGAGTCATACCGACCTGATGGGTAGCGTTTCCACGCCGTTCGCACCAATCGTCAACACCGCTACCTCGCTGAAACCGGTGCGGCAGCTGCTTGATGCGGCGTTAAAAATCGATCATCGCCGCGCGCTAcc                                                           >  1:551196/1‑149 (MQ=255)
                                                                                             tGCGTAATTTTGTGTTGAGTCATACCGACCTGATGGGTAGCGTTTCCACGCCGTTCGCACCAATCGTCAACACCGCTACCTCGCTGAAACCGGTGCGGCAGCTGCTTGATGCGGCGTTAAAAATCGATCATCGCCGCACGCTACCGaaa                                                       <  1:238645/149‑1 (MQ=255)
                                                                                                   aTTTTGTGTTGAGTCATACCGACCTGATGGGTAGCGTTTCCACGCCGTTCGCACCAATCGTCAACACCGCTACCTCGCTGAAACCGGTGCGGCAGCTGCTTGATGCGGCGTTAAAAATCGATCATTGCCGCACGCTACCGAAATACTcc                                                 <  1:806470/149‑1 (MQ=255)
                                                                                                                tCATACCGACCTGATGGGTAGCGTTTCCACGCCGTTCGCACCAATCGTCAAAACCGCTACCTCGCTGAAACCGGTGCGGCAGCTGCTTGATGCGGCGTTAAAAATCGATCATCGCCGCACGCTACCGAAATACTCATTCGGCCGGTTcc                                    >  2:215351/1‑149 (MQ=255)
                                                                                                                   tACCGACCTGATGGGTAGCGTTTCCACGCCGTCCGCACCAATCGTCAACACCGCTACCTCGCTGAAACCGGTGCGGCAGCTGCTTGATGCGGCGTTAAAAATCGATCATCGCCGCACGCTACCGAAATACTCCTTCGGCACGTTCCGTc                                 <  1:353320/149‑1 (MQ=255)
                                                                                                                   tACCGACCTGATGGGTAGCGTTTCCACGCCGTCCGCACCAATCGTCAACACCGCTACCTCGCTGAAACCGGTGCGGCAGCTGCTTGATGCGGCGTTAAAAATCGATCATCGCCGCACGCTACCGAAATACTCCTTCGGCACGTTCCGTc                                 <  1:352877/149‑1 (MQ=255)
                                                                                                                      cGACCTGATGGGTAGCGTTTCCACGCCGTTCGCACCAATCGTCAACACCGCTACCTCGCTGAAACCGGTGCGGCAGCTGCTTGATGCGGCGTTAAAAATCGATCATCGCCGCACGCTACCGAAATACTCCTTCGGCACGTTCCGTCGCt                              >  2:58882/1‑149 (MQ=255)
                                                                                                                         ccTGATGGGTAGCGTTTCCACGCCGTTCGCACCAATCGTCAACACCGCTACCTCGCTGAAACCGGTGCGGCAGCTGCTTGATGCGGCGTTAAAAATCGATCATCGCCGCACGCTACCGAAATACTCCTTCGGCACGTTCCGTCGCTGGt                           >  2:785829/1‑149 (MQ=255)
                                                                                                                           tGATGGGTAGCGTTTCCACGCCGTTCGCACCAATCGTCAACACCGCTACCTCGCTGAAACCGGTGCGGCAGCTGCTTGATGCGGCGTTAAAAATCGATCATCGCCGCACGCTACCGAAATACTCCTTCGGCACGTTCCGTCGCTGGTAt                         >  1:285627/1‑149 (MQ=255)
                                                                                                                                   aGCGTTTCCACGCCGTTCGCACCAATCGTCAACACCGCTACCTCGCTGAAACCGGTGCGGCAGCTGCTTGATGCGGCGTTAAAAATCGATCATCGCCGCACGCTACCGAAATACTCCTTCGGCACGTTCCGTCGCTGGTATCGCAGCGt                 >  1:259765/1‑149 (MQ=255)
                                                                                                                                    gCGTTTCCACGCCGTTCGCACCAATNGTCAACACCGCTACCTCGCTGAAACCGGTGCGGCAGCTGCTTGATGCGGCGTTAAAAATCGATCATCGCCGCACGCTACCGAAATACTCCTTCGGCACGTTCCGTCGCTGGTATCGCAGCGTg                <  1:617110/149‑1 (MQ=255)
                                                                                                                                       tttCCACGCCGTTCGCACCAATCGTCAACACCGCTACCTCGCTGAAACCGGTGCGGCAGCTGCTTGATGCGGCGTTAAAAATCGATCATCGCCGCACGCTACCGAAATACTCCTTCGGCACGTTCCGTCGCTGGTATCGCAGCGTggcg             <  1:755824/149‑1 (MQ=255)
                                                                                                                                         tCCACGCCGTTCGCACCAATCGTCAACACCGCTACCTCGCTGAAACCGGTGCGGCAGCTGCTTGATGCGGCGTTAAAAATCGATCATCGCCGCACGCTACCGAAATACTCCTTCGGCACGTTCCGTCGCTGGTATCGCAGCGTggcggc           >  2:503043/1‑149 (MQ=255)
                                                                                                                                                cGTTCGTACCAATCGTCAACACCGCGTCCTCGCTGAAGCCGGTGCGGCAGCTGCTTGATGCGGCGTTAAAAATCGATCATCGCCGCACGCTACCGAAATACTCCTTCGGCACGTTCCGTCGCTGGTATCGCAGCGTGGCGGCTCagcaa    <  2:456830/149‑1 (MQ=255)
                                                                                                                                                cGTTCGCACCAATCGTCAACACCGCTACCTCGCTGAAACCGGTGCGGCAGCTGCTTGATGCGGCGTTAAAAATCGATCATCGCCGCACGCTACCGAAATACTCCTTCGGCACGTTCCGTCGCTGGTATCGCAGCGTGGCGGCTCagcaa    <  2:56856/149‑1 (MQ=255)
                                                                                                                                                  ttCGCACCAATCGTCAACACCGCTACCTCGCTGAAACCGGTGCGGCAGCTGCTTGATGCGGCGTTAAAAATCGATCATCGCCGCACGCTACCGAAATACTCCTTCGGCACGTTCCGTCGCTGGTATCGCAGCGTGGCGGCTCagcaagc  >  1:317997/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GCAAACGTTGTGAAGTCGCCTGCCCGTCCGATGTGAAGATTGGCGATATTATCCAGCGCGCGCGGGCGAAATATGACACCACGCGCCCGTCGCTGCGTAATTTTGTGTTGAGTCATACCGACCTGATGGGTAGCGTTTCCACGCCGTTCGCACCAATCGTCAACACCGCTACCTCGCTGAAACCGGTGCGGCAGCTGCTTGATGCGGCGTTAAAAATCGATCATCGCCGCACGCTACCGAAATACTCCTTCGGCACGTTCCGTCGCTGGTATCGCAGCGTGGCGGCTCAGCAAGC  >  NC_000913_3_pae_tpiA/2355696‑2355990

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