Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,105,994 T→C 85.3% T12T (ACT→ACC pdhD → Dihydrolipoyl dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,105,9940TC85.3% 140.8 / 12.0 60T12T (ACT→ACCpdhDDihydrolipoyl dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (6/3);  new base C (23/28);  total (29/31)
Fisher's exact test for biased strand distribution p-value = 2.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.95e-01

TTTGACCACAGACAAATTGATGGTGCAACTGGCCAAAATGCAATGAATCACATTAAACGTTTATTAAATAATCCAGAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACTATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCAATTAGGACAAAAAGTAACAATCGTTGAGAAAGGTAATCTTGGTGGTGTTTGCTTAAACGTAGGATGTATTCCTTCAAAA  >  USA300TCH1516_ALE/1105857‑1106126
                                                                                                                                         |                                                                                                                                    
tttGACCACAGACAAATTGATGGTGCAACTGGCCAAAATGCAATGAATCACATTAAACGTTTATTAAATAATCCAGAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACcata                                                                                                                                   <  1:307581/141‑1 (MQ=255)
 ttGACCACAGACAAATTGATGGTGCAACTGGCCAAAATGCAATGAATCACATTAAACGTTTATTAAATAATCCAGAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAg                                                                                                                                  <  1:517729/141‑1 (MQ=255)
         aGACAAATTGATGGTGCAACTGGCCAAAATGCAATGAATCACATTAAACGTTTATTAAATAATCCAGAATTATCATTAATGGAGGGGTAAAACATAGTGGTTGGAGATTTCCAAATTGAAACAAATTCCATTATAATCGGa                                                                                                                          >  2:420067/1‑141 (MQ=255)
               aTTGATGGTGCAACTGGCCAAAATGCAATGAATCACATTAAACGTTTATTAAATAATCCAGAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGt                                                                                                                    >  1:511250/1‑141 (MQ=255)
                ttGATGGTGCAACTGGCCAAAATGCAATGAATCACATTAAACGTTTATTAAATAATCCAGAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACTATAGTAATCGGAGCAGGTc                                                                                                                   >  2:475183/1‑141 (MQ=255)
                  gATGGTGCAACTGGCCAAAATGCAATGAATCACATTAAACGTTTATTAAATAATCCAGAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCt                                                                                                                 <  1:185246/141‑1 (MQ=255)
                     ggTGCAACTGGCCAAAATGCAATGAATCACATTAAACGTTTATTAAATAATCCAGAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCtggt                                                                                                              <  1:261346/141‑1 (MQ=255)
                      gTGCAACTGGCCAAAATGCAATGAATCACATTAAACGTTTATTAAATAATACAGAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGAAGGTCCtggtg                                                                                                             >  2:48993/1‑141 (MQ=255)
                        gCAACTGGCCAAAATGCAATGAATCACATTAAACGTTTATTAAATAATCCAGAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGa                                                                                                           >  1:158391/1‑141 (MQ=255)
                        gCAACTGGACAAAATGCAATGAATCACATTAAACGTTTATTAAATAATCCAGAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATGTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGa                                                                                                           <  2:44041/141‑1 (MQ=255)
                         caacgtgCCAAAATGCAATCAATCACATTAAACGTTTATTAAATAAGCCGGAATTATTATTAATGGAGGGGTAAAACAGTGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGGGCAGGTCCTGGTGGAt                                                                                                          <  2:44303/135‑1 (MQ=255)
                           aCTGGCCAAAATGCAATGAATCACATTAAACGTTTATTAAATAATCCAGAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATAc                                                                                                        >  2:104556/1‑141 (MQ=255)
                            cTGGCCAAAATGCAATGAATCACATTAAACGTTTATTAAATAATCCAGAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACg                                                                                                       <  1:567039/141‑1 (MQ=255)
                             ttggttaaaTGCAATGTATTACATTAAACCTTTATTAAATAATCCATAATTATGATTAATGTAGGGGTAAAACATGGGAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGt                                                                                                      <  2:288922/135‑1 (MQ=255)
                                ccAAAATGCAATGAATCACATTAAACGTTTATTAAATAATCCAGAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTgc                                                                                                   >  1:559879/1‑141 (MQ=255)
                                         aaTGAATCACATTAAACGTTAATTAAATAATCCAGAATTATTATCAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCg                                                                                          >  1:371784/1‑141 (MQ=255)
                                          aTGAATCACATTAAACGTTTATTAAATAATCCAGAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGt                                                                                         >  2:504540/1‑141 (MQ=255)
                                           tGAATCACATTAAACGTTTATTAAATAATCCAGAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTg                                                                                        >  2:589448/1‑141 (MQ=255)
                                                cacaTTAAACGTTTATTAAATAATCCAGAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACTATAGTAATCGGAGCAGGTCCAGGTGGATACGTTGCAGCAATTCGTGCAGcc                                                                                   >  2:564887/1‑140 (MQ=255)
                                                cacaTTAAACGTTTATTAAATAATCCAGAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCt                                                                                   <  2:43787/141‑1 (MQ=255)
                                                cacaTTAAACGTTTATTAAATAATCCAGAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCt                                                                                   <  1:499509/141‑1 (MQ=255)
                                                cacaTTAAACGTTTATTAAATAATCCAGAATCATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCt                                                                                   <  1:479007/141‑1 (MQ=255)
                                                 acaTTAAACGTTTATTAAATAATCCAGAATTATTATTAATGGCGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGTTACTATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGAATCAATTCGTGCCGCTc                                                                                  >  2:76219/1‑141 (MQ=255)
                                                 acaTTAAACGTTTATTAAATAATCCAGAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTGGTGCAGCta                                                                                  >  1:126673/1‑140 (MQ=255)
                                                   aTTAAACGTTTATTAAATAATCCAGAATTATTATTAATGGAGGGGTACAACATGGTAGTTGGAGATTTCACAATTGAAACCGATACCATAGTAATCGGAGCAGGACCTAGTGGATACGTTGCATCAAATGGTGCAGCGCaa                                                                                >  2:169425/1‑141 (MQ=255)
                                                   aTTAAACGTTTATTAAATAATCCAGAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCaa                                                                                <  1:503289/141‑1 (MQ=255)
                                                   aTTAAACGTTTATTAAATAATCCAGAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGACCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCaa                                                                                >  1:353921/1‑141 (MQ=255)
                                                    ttAAACGCTTATTAAATAATCCAGAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCAAt                                                                               <  2:382109/141‑1 (MQ=255)
                                                        aCGTTTATTAAATAATCCAGAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTAGTGCAGCTCAATTAgg                                                                           >  2:146737/1‑141 (MQ=255)
                                                         cGTTTATTAAATAATCCAGAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACTATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCAGTTAGGa                                                                          >  2:194111/1‑141 (MQ=255)
                                                          gTTTATTAAATAATCCAGAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCAATTAGGAc                                                                         >  2:546967/1‑141 (MQ=255)
                                                          gTTTATTAAATAATCCAGAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCAATTAGGAc                                                                         <  1:384069/141‑1 (MQ=255)
                                                              atgaaaTAATCCAGAATTATTATTAATGGAGGGGTAAAACAGGGTAGTTGGAGATTGCCCAATTGAAACAGATACCATAGTAATCGGGGGAGGTCCTGGTGGATGCGTTGCAGAAATTCGTGCAGCTCAAATAGAAGaaaa                                                                     >  2:160110/4‑141 (MQ=255)
                                                               gtaAATAATCCAGAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGACTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGTTACGTTGCAGCAATTCGTGCAGCTCAATTAGGACaaaaa                                                                    <  1:279233/140‑1 (MQ=255)
                                                                taAATAATCCAGAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCAATTAGGACAAAAAg                                                                   >  2:445774/1‑141 (MQ=255)
                                                                taAATAATCCAGAATTATTATTAATGGAGGGGTAAAAAATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGAAGGTCCTGGTGAATACGTTGCAGCAATTCGTGCAGCTCAATTAGGACAAAAAg                                                                   >  2:114549/1‑141 (MQ=255)
                                                                   ataatCCAGAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCAATTAGGACAAAAAGtca                                                                >  2:495014/1‑139 (MQ=255)
                                                                         cAGAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCAATTAGGACAAAAAGTAACAATCg                                                          <  1:589448/141‑1 (MQ=255)
                                                                           gAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCAATTAGGACAAAAAGTAACAATCGtt                                                        <  1:445774/141‑1 (MQ=255)
                                                                               tattattaATGGAGGGGTAAAACATGGTAGTTGGGGATTTCCCAATTGAAACAGATACTATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCAATTAGGACAAAAAGTAACAATCGTTgaga                                                    <  2:70278/141‑1 (MQ=255)
                                                                               tattattaATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCAATTAGGACAAAAAGTAACAATCGTTgaga                                                    <  2:270774/141‑1 (MQ=255)
                                                                               tattattaATGAAGTGTTAATACATTGTAGTTGGATATTTCCCAATTGAAACAGATACTATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCAATTAGGACAAAAAGTAACAATCGTTgaga                                                    <  2:475536/141‑1 (MQ=255)
                                                                                 ttgttaATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACTATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCAATTAGGACAAAAAGTAACAATCGTTGAGaaa                                                  <  1:259626/138‑1 (MQ=255)
                                                                                  tattaAGGGAGGGGTACAACATGGTAGTTGGGGATTTACACATTGAACCAGATACCAAAGTAATCGGAGCAGGTCCTGGTGGAAACGTTGCAGCAATTCGTGCAGCTCAATTTGGACAAAAATTAACAATCGTTGAGAAAg                                                 <  2:110484/141‑1 (MQ=255)
                                                                                    ttaATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCAATTAGGACAAAAAGTAACAATCGTTGAGAAAGGt                                               <  2:182538/141‑1 (MQ=255)
                                                                                    ttaATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCAATTAGGACAAAAAGTAACAATCGTTGAGAAAGGt                                               <  2:468711/141‑1 (MQ=255)
                                                                                          gAGGGGTAAAACATGGTAGTTGGAGATATCCCAATTGAAACAGATACCATAGTCATCGGAGCAGGTCAAGGTGGATACGTTGCAGCAATTCGTGCAGCTCAATTAGGACAAAAAGTAACAATCGTTGAGAAAGGTAATCtt                                         >  1:571383/1‑141 (MQ=255)
                                                                                              ggTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCAATTAGGACAAAAAGTAACAATCGTTGAGAAAGGTAATCTtggtg                                     <  2:371784/141‑1 (MQ=255)
                                                                                                tAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCAATTAGGACAAAAAGTAACAATCGTTGAGAAAGGTAATATtggtgag                                   <  2:158391/141‑3 (MQ=255)
                                                                                                  aaaCATGGTAGTTGGAGATTTCCCAATTGAAACAGATACTATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCAATTAGGACAAAAAGTAACAATCGTTGAGAAAGGTAATCTTGGTGgtgt                                 >  2:164319/1‑141 (MQ=255)
                                                                                                  aaaCATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCCATTCGTGCAGCTCAATTAGGACAAAAAGTAACAATCGTTGAGAAAGGTAATCTTGGTGgtgt                                 <  1:571465/141‑1 (MQ=255)
                                                                                                           aGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCAATTAGGACAAAAAGTAACAATCGTTGAGAAAGGTAATCTTGGTGGTGTTTGCTTaaa                        >  2:271178/1‑141 (MQ=255)
                                                                                                            gTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCAATTAGGACAAAAAGTAACAATCGTTGAGAAAGGTAATCTTGGTGGTTGTTTCTTaaaa                       >  1:354087/1‑140 (MQ=255)
                                                                                                              tGGAGATTTCCCAATTTAAACAGATACCATAGTAATGTGATCAGTTGCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCAATTAGGACAAAAAGTAACAATCGTTGAGAAAGGTAATCTTGGTGGTGTTTGCTTAAACGt                     <  2:126673/141‑1 (MQ=255)
                                                                                                               ggAGATTTCCCAATTGAAACAGATACTATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCAATTAGGACAAAAAGTAACAATCGTTGAGAAAGGTAATCTTGGTGGTGTTTGCTTAAACGTa                    >  2:353782/1‑141 (MQ=255)
                                                                                                                 agaTTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCAATTAGGACAAAAAGTAACAATCGTTGAGAAAGGTAATCTTGGTGGTGTTTGCTTAAACGTAgg                  >  1:326284/1‑141 (MQ=255)
                                                                                                                   aTTTCCCAATTGAAACACATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCAATTAGGACAAAAAGTAACAATCGTTGAGAAAGGTAATCTTGGTGGTGTTTGCTTAAACGTAGGAt                <  1:264885/141‑1 (MQ=255)
                                                                                                                    tttCCCAATTGAAACAGATACCATAGTAATCGGAGCATGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCAATTAGGACAAAAAGTAACAATCGTTGAGAAAGGTAATCTTGGTGGTGTTTGCTTAAACGTAGGATg               <  1:262412/141‑1 (MQ=255)
                                                                                                                    tttCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCAATTAGGACAAAAAGTAACAATCGTTGAGAAAGGTAATCTTGGTGGTGTTTGCTTAAACGTAGGATg               <  2:353921/141‑1 (MQ=255)
                                                                                                                              gAAACAGATACCATAGTAATCGGAGCAGGTCCTGGGGGATACGTTGCAGCAATTCGTGCAGCTCAATTAGGACAAAAAGTAACAATCGTTGAGAAAGGTAATCTTGGTGGTGTTTGCTTAAACGTAGGATGTATTCCTTCa     <  1:91139/141‑1 (MQ=255)
                                                                                                                                 aCAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCAATTAGGACGAAAGGTAACAATCGTTGAGAAAGGTAATCTTGGTGGTGTTTGCTTAAACGAAGGAGGTATTCCTTCaaag  >  2:281/1‑140 (MQ=255)
                                                                                                                                         |                                                                                                                                    
TTTGACCACAGACAAATTGATGGTGCAACTGGCCAAAATGCAATGAATCACATTAAACGTTTATTAAATAATCCAGAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACTATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCAATTAGGACAAAAAGTAACAATCGTTGAGAAAGGTAATCTTGGTGGTGTTTGCTTAAACGTAGGATGTATTCCTTCAAAA  >  USA300TCH1516_ALE/1105857‑1106126

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