Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 129,360 G→T intergenic (+24/+47) lpd → / ← yacH dihydrolipoyl dehydrogenase; E3 component of pyruvate and 2‑oxoglutarate dehydrogenases complexes; glycine cleavage system L protein; dihydrolipoamide dehydrogenase/DUF3300 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913129,3600GT93.7% 205.7 / 3.7 63intergenic (+24/+47)lpd/yacHdihydrolipoyl dehydrogenase; E3 component of pyruvate and 2‑oxoglutarate dehydrogenases complexes; glycine cleavage system L protein; dihydrolipoamide dehydrogenase/DUF3300 family protein
Reads supporting (aligned to +/- strand):  ref base G (0/4);  new base T (41/18);  total (41/22)
Fisher's exact test for biased strand distribution p-value = 1.23e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

TCTGTGGGCCTGGCGGCAGAAGTGTTCGAAGGTAGCATTACCGACCTGCCGAACCCGAAAGCGAAGAAGAAGTAATTTTTCGTTTGCCGGAACATCCGGCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCTGCAATTTACCTTTCCAGTCTTCTTGCTCCACGTTCAGAGAGACGTTCGCATAC  >  NC_000913/129262‑129454
                                                                                                  |                                                                                              
tCTGTGGGCCTGGCGGCAGAAGTGTTCGAAGGTAGCATTACCGACCTGCCGAACCCGAAAGCGAAGAAGAAGTAATTTTTCGTTTGCCGGAACATCCGGc                                                                                               <  1:174581/100‑1 (MQ=255)
   gtgGGCCTGGCGGCAGAAGTGTTCGAAGGTAGCATTACCGACCTGCCGAACCCGAAAGCGAAGAAGAAGTAATTTTTCGTTTGCCGGAACATCCGTCAAt                                                                                            >  2:516302/1‑100 (MQ=255)
      ggCCTGGCGGCAGAAGTGTTCGAAGGTAGCATTACCGACCTGCCGAACCCGAAAGCGAAGAAGAAGTAATTTTTCGTTTGCCGGAACATCCGTCAATTaa                                                                                         >  2:1293822/1‑100 (MQ=255)
      ggCCTGGCGGCAGAAGTGTTCGAAGGTAGCATTACCGACCTGCCGAACCCGAAAGCGAAGAAGAAGTAATTTTTCGTTTGCCGGAACATCCGTCAATTaa                                                                                         >  1:402283/1‑100 (MQ=255)
               gcAGAAGTGTTCGAAGGTAGCATTACCGACCTGCCGAACCCGAAAGCGAAGAAGAAGTAATTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGc                                                                                >  1:969320/1‑100 (MQ=255)
                cAGAAGTGTTCGAAGGTAGCATTACCGACCTGCCGAACCCGTAAGCGAAGAAGAAGTAATTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCt                                                                               >  2:371181/1‑100 (MQ=255)
                cAGAAGTGTTCGAAGGTAGCATTACCGACCTGCCGAACCCGAAAGCGAAGAAGAAGTAATTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCt                                                                               >  2:1513730/1‑100 (MQ=255)
                       gtTCGAAGGTAGCATTACCGACCTGCCGAACCCGAAAGCGAAGAAGAAGTAATTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCAcg                                                                        >  1:1467752/1‑100 (MQ=255)
                         tCGAAGGTAGCATTACCGACCTGCCGAACCCGAAAGCGAAGAAGAAGTAATTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCAcgcc                                                                      <  2:691824/100‑1 (MQ=255)
                          cGAAGGTAGCATTACCGACCTGCCGAACCCGAAAGCGAAGAAGAAGTAATTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCAcgccg                                                                     >  2:568769/1‑100 (MQ=255)
                             aGGTAGCATTACCGACCTGCCGAACCCGAAAGCGAAGAAGAAGTAATTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCtt                                                                  >  1:1245357/1‑100 (MQ=255)
                             aGGTAGCATTACCGACCTGCCGAACCCGAAAGCGAAGAAGAAGTAATTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCtt                                                                  >  1:1295709/1‑100 (MQ=255)
                              ggTAGCATTACCGACCTACCGAACCCGAAAGCGAAGAAGAAGTAATTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCttt                                                                 >  2:195591/1‑100 (MQ=255)
                                 aGCATTACCGACCTGCCGAACCCGAAAGCGAAGAAGCAGTAATTTTTCGTTTGCCGGACCATCCGTCAATTACCAAAGCGGCTAACCAcgccgc                                                                    >  2:1530997/1‑94 (MQ=255)
                                  gCATTACCGACCTGCCGAACCCGAAAGCGAAGAAGAAGTAATTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCttttttt                                                             <  1:921730/100‑1 (MQ=255)
                                  gCATTACCGACCTGCCGAACCCGAAAGCGAAGAAGAAGTAATTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCttttttt                                                             >  2:783102/1‑100 (MQ=255)
                                    aTTACCGACCTGCCGAACCCGAAAGCGAAGAAGAAGTAATTTTTCGTTTGCCGGAACATCCGGCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTAc                                                           <  2:201738/100‑1 (MQ=255)
                                       aCCGACCTGCCGAACCCGAAAGCGAAGAAGAAGTAATTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGt                                                         >  2:1033959/1‑99 (MQ=255)
                                        ccGACCTGCCGAACCCGAAAGCGAAGAAGAAGTAATTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCt                                                       >  1:1173880/1‑100 (MQ=255)
                                         cGACCTGCCGAACCCGAAAGCGAAGAAGAAGTAATTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCTg                                                      <  2:637520/100‑1 (MQ=255)
                                          gACCTGCCGAACCCGAAAGCGAAGAAGAAGTAATTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCTGc                                                     >  1:783857/1‑100 (MQ=255)
                                           aCCTGCCGAACCCGAAAGCGAAGAAGAAGTAATTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCTGCa                                                    >  2:679361/1‑100 (MQ=255)
                                            ccTGCCGAACCCGAAAGCGAAGAAGAAGTAATTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCTGCaa                                                   >  2:583091/1‑100 (MQ=255)
                                               gCCGAACCCGAAAGCGAAGAAGAAGTAATTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCTGCAAttt                                                >  1:922431/1‑100 (MQ=255)
                                               gCCGAACCCGAAAGCGAAGAAGAAGTAATTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCTGCAAttt                                                <  1:1328613/100‑1 (MQ=255)
                                                ccGAACCCGAAAGCGAAGAAGAAGTAATTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCTGCAATTTa                                               >  1:1407846/1‑100 (MQ=255)
                                                ccGAACCCGAAAGCGAAGAAGAAGTAATTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCATTTTTTACGTCTGCAATTTa                                               <  2:969315/100‑1 (MQ=255)
                                                 cGAACCCGAAAGCGAAGAAGAAGTAATTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCCAGCCGCTTTTTTTACGTCTGCAAttt                                                >  2:592042/1‑98 (MQ=255)
                                                 cGAACCCGAAAGCGAAGAAGAAGTAATTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCTGCAATTTAc                                              >  1:465531/1‑100 (MQ=255)
                                                 cGAACCCGAAAGCGAAGAAGAAGTAATTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCTGCAATTTAc                                              <  1:1513739/100‑1 (MQ=255)
                                                     cccGAAAGCGAAGAAGAAGTAATTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCTGCAATTTACCttt                                          >  1:56185/1‑100 (MQ=255)
                                                     cccGAAAGCGAAGAAGAAGTAATTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCTGCAATTTACCttt                                          >  1:177213/1‑100 (MQ=255)
                                                      ccGAAAGCGAAGAAGAAGTAATTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCAGCTAACCACGCCGCTTTTTTTACGTCTGCAATTTACCttt                                          >  1:1228149/1‑99 (MQ=255)
                                                       cGAAAGCGAAGAAGAAGTAATTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCTGCAATTTACCTTTcc                                        >  2:817003/1‑100 (MQ=255)
                                                       cGAAAGCGAAGAAGAAGTAATTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCTGCAATTTACCTTTcc                                        <  2:1247902/100‑1 (MQ=255)
                                                             cGAAGAAGAAGTAATTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCTGCAATTTACCTTTCCAGtctt                                  >  1:927117/1‑100 (MQ=255)
                                                              gaagaagaaGTAATTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCTGCAATTTACCTTTCCAGtcttc                                 <  2:1334378/100‑1 (MQ=255)
                                                              gaagaagaaGTAATTTTTCGTTTGCCGGAACATCCGGCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCTGCAATTTACCTTTCCAGtcttc                                 <  2:431416/100‑1 (MQ=255)
                                                               aagaagaagTAATTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCTGCAATTTACCTTTCCAGtcttct                                <  1:936317/100‑1 (MQ=255)
                                                                 gaagaagTAATTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCTGCTAACCACGCCGCTTTTTTTACGTCTGCAATTTACCTTTCCAGTcttctt                               >  2:561512/1‑99 (MQ=255)
                                                                   agaagTAATTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAGCCACGCCGCTTTTTTTACGTCTGCAATTTACCTTTCCAGTCTTCTTGCt                            <  2:1521179/100‑1 (MQ=255)
                                                                      agTAATTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCTGCAATTTACCTTTCCAGTCTTCTTGCTCCa                         <  2:1442367/100‑1 (MQ=255)
                                                                      agTAATTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCTGCAATTTACCTTTCCAGTCTTCTTGCTCCa                         >  1:647857/1‑100 (MQ=255)
                                                                       gTAATTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCTGCAATTTACCTTTCCAGTCTTCTTGCTCCAc                        >  1:322513/1‑100 (MQ=255)
                                                                        tAATTTTTCGTTTGCCGGAACATCCGGCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCTGCAATTTACCTTTCCAGTCTTCTTGCTCCACg                       <  1:1014272/100‑1 (MQ=255)
                                                                         aaTTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCTGCAATTTACCTTTCCAGTCTTCTTGCTCCACGt                      >  1:981614/1‑100 (MQ=255)
                                                                         aaTTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCTGCAATTTACCTTTCCAGTCTTCTTGCTCCACGt                      >  1:300793/1‑100 (MQ=255)
                                                                         aaTTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCTGCAATTTACCTTTCCAGTCTTCTTGCTCCACGt                      >  1:1136425/1‑100 (MQ=255)
                                                                          aTTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCTGCAATTTACCTTTCCAGTCTTCTTGCTCCACGtt                     >  1:372243/1‑100 (MQ=255)
                                                                           tttttCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCTGCAATTTACCTTTCCAGTCTTCTTGCTCCACGTTc                    <  2:986018/100‑1 (MQ=255)
                                                                             tttCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCTGCAATTTACCTTTCCAGTCTTCTTGCTCCACGTTCag                  <  1:975436/100‑1 (MQ=255)
                                                                             tttCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCTGCAATTTACCTTTCCAGTCTTCTTGCTCCACGTTCag                  <  1:195591/100‑1 (MQ=255)
                                                                                 gTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCTGCAATTTACCTTTCCAGTCTTCTTGCTCCACGTTCagagag              >  1:396206/1‑100 (MQ=255)
                                                                                 gTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCTGCAATTTACCTTTCCAGTCTTCTTGCTCCACGTTCagagag              >  1:1531160/1‑100 (MQ=255)
                                                                                 gTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCTGCAATTTACCTTTCCAGTCTTCTTGCTCCACGTTCagagag              >  1:1042663/1‑100 (MQ=255)
                                                                                    tGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCTGCAATTTACCTTTCCAGTCTTCTTGCTCCACGTTCAGAGAGACg           >  2:168880/1‑100 (MQ=255)
                                                                                    tGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCTGCAATTTACCTTTCCAGTCTTCTTGCTCCACGTTCAGAGAGACg           >  2:876932/1‑100 (MQ=255)
                                                                                      ccGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCTGCAATTTACCTTTCCAGTCTTCTTGCTCCACGTTCAGAGAGACGt          >  2:970774/1‑99 (MQ=255)
                                                                                        ggAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCTGCAATTTACCTTTCCAGTCTTCTTGCTCCACGTTCAGAGAGACGTTCg       <  2:234754/100‑1 (MQ=255)
                                                                                        ggAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTGACGTCTGCAATTTACCTTTCCAGTCTTCTTGCTCCACGTTCAGAGAGACGTTCg       >  2:1173869/1‑100 (MQ=255)
                                                                                          aaCATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCTGCAATTTACCTTTCCAGTCTTCTTGCTCCACGTTCAGAGAGACGTTCGCa     <  2:322513/100‑1 (MQ=255)
                                                                                             aTCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCTGCAATTTACCTTTCCAGTCTTCTTGCTCCACGTTCAGAGAGACGTTCGCATAc  <  1:516303/100‑1 (MQ=255)
                                                                                             aTCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCTGCAATTTACCTTTCCAGTCTTCTTGCTCCACGTTCAGAGAGACGTTCGCATAc  <  2:775829/100‑1 (MQ=255)
                                                                                                  |                                                                                              
TCTGTGGGCCTGGCGGCAGAAGTGTTCGAAGGTAGCATTACCGACCTGCCGAACCCGAAAGCGAAGAAGAAGTAATTTTTCGTTTGCCGGAACATCCGGCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCTGCAATTTACCTTTCCAGTCTTCTTGCTCCACGTTCAGAGAGACGTTCGCATAC  >  NC_000913/129262‑129454

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