Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 1,105,994 T→C 16.9% T12T (ACT→ACC)  pdhD → dihydrolipoyl dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007301,105,9940TC16.9% 61.6 / 5.7 29T12T (ACT→ACC) pdhDdihydrolipoyl dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (11/13);  new base C (2/3);  total (13/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTGACCACAGACAAATTGATGGTGCAACTGGCCAAAATGCAATGAATCACATTAAACGTTTATTAAATAATCCAGAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACTATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCAATTAGGACAAAAAGTAACAATCGTTGAGAAAGGTAATCTTGGTGGTGTTTGCTTAAACGTAGGATGTATTCCTTCAAAAGCAT  >  CP000730/1105857‑1106130
                                                                                                                                         |                                                                                                                                        
tttGACCACAGACAAATTGATGGTGCAACTGGCCAAAATGCAATGAATCACATTAAACGTTTATTAAATAATCCAGAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACtata                                                                                                                                       >  1:108820/1‑141 (MQ=255)
   gACCACAGACAAATTGATGGTGCAACTGGCCAAAATGCAATGAATCACATTAAACGTTTATTAAATAATCCAGAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACTATAGTa                                                                                                                                    >  2:206787/1‑141 (MQ=255)
    aCCACAGACAAATTGATGGTGCAACTGGCCAAAATGCAATGAATCACATTAAACGTTTATTAAATAATCCAGAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACTATAGTaa                                                                                                                                   >  2:286261/1‑141 (MQ=255)
             gttcaattattGCAACTGGCCAAAATGCAATGAATCACATTAAACGTTTATTAAATAATCCAGAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACTATAGTAATCGGAGCAg                                                                                                                          <  2:40031/131‑1 (MQ=255)
               aTTGATGGTGCAACTGGCCAAAATGCAATGAATCACATTAAACGTTTATTAAATAATCCAGAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACTATAGTAATCGGAGCAGGt                                                                                                                        <  1:24519/141‑1 (MQ=255)
                    tGGTGCAACTGGCCAAAATGCAATGAATCACATTAAACNTTTATTAAATAATCCAGAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGCTACTATAGTAATCGGAGCAGGTCCtgg                                                                                                                   >  1:152871/1‑141 (MQ=255)
                    tGGTGCAACTGGCCAAAATGCAATGAATCACATTAAACGTTTATTAAATAATCCAGAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACTATAGTAATCGGAGCGGGTCCtgg                                                                                                                   >  2:158232/1‑141 (MQ=255)
                    tGGTGCAACTGGCCAAAATGCAATGAATCACATTAAACGTTTATTAAATAATCCAGAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACTATAGTAATCGGAGCAGGTCCtgg                                                                                                                   >  2:91097/1‑141 (MQ=255)
                      gTGCAACTGGCCAAAATGCAATGAATCACATTAAACGTTTATTAAATAATCCAGAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCtggtg                                                                                                                 >  2:142913/1‑141 (MQ=255)
                             tGGCCAAAATGCAATGAATCACATTAAACGTTTATTAAATAATCCAGAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGt                                                                                                          <  1:205263/141‑1 (MQ=255)
                                  aaaaTGCAATGAATCACATTAAACGTTTATTAAATAATCCAGAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACTATAGTAATCGGAGCAGGTCCTGGTGGATACGTTgcag                                                                                                     <  1:138378/141‑1 (MQ=255)
                                         aaTGAATCACATTAAACGTTTATTAAATAATCCAGAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACTATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCg                                                                                              >  2:189999/1‑141 (MQ=255)
                                          aTGAATCACATTAAACGTTTATTAAATAATCCAGAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACTATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGt                                                                                             >  1:31391/1‑141 (MQ=255)
                                              aTCACATTAAACGTTTATTAAATAATCCAGAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACTATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAg                                                                                         <  1:35615/141‑1 (MQ=255)
                                                cacaTTAAACGTTTATTAAATAATCCAGAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAACCAGATACTATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCt                                                                                       <  1:257078/141‑1 (MQ=255)
                                                      aaaCGTTTATTAAATAATCCAGAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTGCCAAGTGAAACAGCGACTGTAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCAATTa                                                                                 >  1:171898/1‑141 (MQ=255)
                                                             tattaAATAATCCAGAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCAATTAGGACaaa                                                                          >  2:186700/1‑141 (MQ=255)
                                                               ttaAAGAATCCAGAGTTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCAATTAGGACaaaaa                                                                        <  2:330609/141‑1 (MQ=255)
                                                                               tattattaATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACTATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCAATTAGGACAAAAAGTAACAATCGTTgaga                                                        <  1:116481/141‑1 (MQ=255)
                                                                                       aTGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACTATAGTAATCGGAGCAGGTCATGGTGGATACGTTGCAGCAATTCGTGCAGCGCAATTAGGACAAAAAGTAACAATCGTTGAGAAAGGTAAt                                                <  2:108820/141‑1 (MQ=255)
                                                                                                   aaCATGGTAGTTGGAGATTTCCCAATTGAAACAGATACTATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCAATTAGGACAAAAAGTAACAATCGTTGAGAAAGGTAATCTTGGTGGTGtt                                    <  1:185577/141‑1 (MQ=255)
                                                                                                       tGGTAGTTGGAGATTTCCCAATTGAAACAGATACTATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCAATTAGGACAAAAAGTAACAATCGTTGAGAAAGGTAATCTTGGTGGTGTTTGCt                                >  1:33352/1‑141 (MQ=255)
                                                                                                              tGGAGATTTCCCAATTGAAACAGATACTATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCAATTAGGACAAAAAGTAACAATCGTTGAGAAAGGTAATCTTGGTGGTGTTTGCTTAAACGt                         <  1:318774/141‑1 (MQ=255)
                                                                                                                gagaTTTCCCAATTGAAACAGATACTATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCAATTAGGACAAAAAGTAACAATCGTTGAGAAAGGTAATCTTGGTGGTGTTTGCTTAAACGTAg                       <  1:305779/141‑1 (MQ=255)
                                                                                                                 agaTTTCCCAATTGAAACAGATACTATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCAATTAGGACAAAAAGTAACAATCGTTGAGAAAGGTAATCTTGGTGGTGTTTGCTTAAACGTAgg                      <  1:5748/141‑1 (MQ=255)
                                                                                                                   aTTTCCCAATTGAAACAGATACTATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCAATTAGGACAAAAAGTAACAATCGTTGAGAAAGGTAATCTTGGTGGTGTTTGCTTAAACGTAGGAt                    <  1:206787/141‑1 (MQ=255)
                                                                                                                           aTTGAAACAGATACTATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCAATTAGGACAAAAAGTAACAATCGTTGAGAAAGGTAATCTTGGTGGTGTTTGCTTAAACGTAGGATGTATTCCt            >  2:92768/1‑141 (MQ=255)
                                                                                                                                aaCAGATACTATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCAATTAGGACAAAAAGTAACAATCGTTGAGAAAGGTAATCTTGGTGGTGTTTGCTTAAACGTAGGATGTATTCCTTCaaa       >  2:317228/1‑141 (MQ=255)
                                                                                                                                  cAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCAATTAGGACAAAAAGTAACAATCGTTGAGAAAGGTAATCTTGGTGGTGTTTGCTTAAACGTAGGATGTATTCCTTCAAAAg     <  1:142913/141‑1 (MQ=255)
                                                                                                                                     aTACTATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCAATTAGGACAAAAAGTAACAATCGTTGAGAAAGGTAATCTTGGTGGTGTTTGCTTAAACGTAGGATGTATTCCTTCAAAAGCAt  <  1:214205/141‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                        
TTTGACCACAGACAAATTGATGGTGCAACTGGCCAAAATGCAATGAATCACATTAAACGTTTATTAAATAATCCAGAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACTATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCAATTAGGACAAAAAGTAACAATCGTTGAGAAAGGTAATCTTGGTGGTGTTTGCTTAAACGTAGGATGTATTCCTTCAAAAGCAT  >  CP000730/1105857‑1106130

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