Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,103,408 A→G 10.8% L310L (TTA→TTG pdhA → Pyruvate dehydrogenase E1 component subunit alpha

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,103,4080AG10.8% 102.7 / 5.8 37L310L (TTA→TTGpdhAPyruvate dehydrogenase E1 component subunit alpha
Reads supporting (aligned to +/- strand):  ref base A (18/15);  new base G (2/2);  total (20/17)
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.74e-01

ACTATGACATATCGTTATGGTCCTCATACAATGGCTGGTGACGATCCAACTCGTTACAGAACTTCAGACGAAGATGCTGAATGGGAGAAAAAAGACCCATTAGTACGTTTCCGTAAATTCCTTGAAAACAAAGGTTTATGGAATGAAGACAAAGAAAATGAAGTTATTGAACGTGCAAAAGCTGATATTAAAGCAGCAATTAAAGAGGCTGATAACACTGAAAAACAAACTGTTACTTCTCTAATGGAAATTATGTATGAAGATATGCCTCAAAACTT  >  USA300TCH1516_ALE/1103271‑1103548
                                                                                                                                         |                                                                                                                                            
aCTATGACATATCGTTATGGTCCTCATACAATGGCTGGTGACGATCCAACTCGTTACAGAACTTCAGACGAAGATGCTGAATGGGAGAAAAAAGACCCATTAGTACGTTTCCGTAAATTCCTTGAAAACAAAGGTTTATgg                                                                                                                                           >  1:27421/1‑141 (MQ=255)
aCTATGACATATCGTTATGGTCCTCAGACAATGGCTGGTGACGATCCAACTCGTTACAGAACTTCAGACGAAGATGCTGAATGGGAGAAAAAAGACCCATTAGTACGTTTCCGTAAATTCCTTGAAAACAAAGGTTTATgg                                                                                                                                           <  1:100363/141‑1 (MQ=255)
                      cTCATACAATGGCTGGTGACGATCCAACTCGTTACAGAACTTCAGACGAAGATGCTGAATGGGAGAAAAAAGACCCATTAGTACGTTTCCGTAAATTCCTTGAAAACAAAGGTTTATGGAATGAAGACAAAGAAAATGAAg                                                                                                                     <  2:245759/141‑1 (MQ=255)
                       tCATACAATGGCTGGTGACGATCCAACTCGTTACAGAACTTCAGACGAAGATGCTGAATGGGAGAAAAAAGACCCATTAGTACGTTTCCGTAAATTCCTTGAAAACAAAGGTTTATGGAATGAAGACAAAGAAAATGAAGt                                                                                                                    >  1:342075/1‑141 (MQ=255)
                               tGGCTGGTGACGATCCAACTCGTTACAGAACTTCAGACGAAGATGCTGAATGGGAGAAAAAAGACCCATTAGTACGTTTCCGTAAATTCCTTGAAAACAAAGGTTTATGGAATGAAGACAAAGAAAATGAAGTTATTGAAc                                                                                                            >  1:268951/1‑141 (MQ=255)
                                 gCTGGTGACGATCCAACTCGTTACAGAACTTCAGACGAAGATGCTGAATGGGAGAAAAAAGACCCATTAGTACGTTTCCGTAAATTCCTTGAAAACAAAGGTTTATGGAATGAAGACAAAGAAAATGAAGTTATTGAACGt                                                                                                          >  2:325357/1‑141 (MQ=255)
                                     gTGACGATCCAACTCGTTACAGAACTTCAGACGAAGATGCTGAATGGGAGAAAAAAGACCCATTAGTACGTTTCCGTAAATTCCTTGAAAACAAAGGTTTGTGGAATGAAGACAAAGAAAATGAAGTTATTGAACGTGCaa                                                                                                      >  2:36247/1‑141 (MQ=255)
                                     gTGACGATCCAACTCGTTACAGAACTTCAGACGAAGATGCTGAATGGGAGAAAAAAGACCCATTAGTACGTTTCCGTAAATTCCTTGAAAACAAAGGTTTGTGGAATGAAGACAAAGAAAATGAAGTTATTGAACGTGCaa                                                                                                      >  2:36261/1‑141 (MQ=255)
                                      agACGATCCAACTCGTTACAGAACTTCAGACGAAGATGCTGAAAGGAAGAAAAAAGACCCATTAGTACGTTTCCGTAAATTCCTTGAAAACAAAGGTTTATGGAATGAAGACAAAGAAAATGAAGTTATTGAACGTGCaaa                                                                                                     <  2:347030/140‑1 (MQ=255)
                                       gACGATCCAACTCGTTACAGAACTTCAGACGAAGATGCTGAATGGGAGAAAAAAGACCCATTAGTACGTTTCCGTAAATTCCTTGAAAACAAAGGTTTATGGAATGAAGACAAAGAAAATGAAGTTATTGAACGTGCaaaa                                                                                                    >  1:301312/1‑141 (MQ=255)
                                        aCGATCCAACTCGTTACAGAACTTCAGACGAAGATGCTGAATGGGAGAAAAAAGACCCATTAGTACGTTTCCGTAAATTCCTTGAAAACAAAGGTTTATGGAATGAAGACAAAGAAAATGAAGTTATTGAACGTGCAAAAg                                                                                                   <  1:226502/141‑1 (MQ=255)
                                        aCGATCCAACTCGTTACAGAACTACAGACCAAGATGCTGAATAGGAGAAAAAAGACCCATTAGTACGTTCCCGTAAATTCCTTGAAAACAAAGGTTTATGGAATGAAGACAAAGAAAATGAAGTTATTGAACGTGCAAAAg                                                                                                   <  1:208158/141‑1 (MQ=255)
                                         cGATCCAACTCGTTACAGAACTTCAGACGAAGATGCTGAATGGGAGAAAAAAGACCCATTAGTACGTTTCCGTAAATTCCTTGAAAACAAAGGTTTATGGAATGAAGACAAAGAAAATGAAGTTATTGAACGTGCAAAAGc                                                                                                  >  2:161124/1‑141 (MQ=255)
                                         cGATCCAACTCGTTACAGAACTTCAGACGAAGATGCTGAATGGGAGAAAAAAGAACAGTTAGTACGTTTCCGTAAATTCCTTGAAAACAAAGGTTTATGGAATGAAGACAAAGAAAATGAAGTTATTGAACGTGCAAAAGc                                                                                                  >  1:291008/1‑141 (MQ=255)
                                          gATCCAACTCGTTACAGAACTTCAGACGAAGATGCTGAATGGGAGAAAAAAGACCCATTAGTACGTTTCCGTAAATTCCTTGAAAACAAAGGTTTATGGAATGAAGACAAAGAAAATGAAGTTATTGAACGTGCAAAAGCt                                                                                                 >  1:184930/1‑141 (MQ=255)
                                                   cGTTACAGAACTTCAGACGAAGATGCTGAATGGGAGAAAAAAGACCCATTAGTACGTTTCCGTAAATTCCTTGAAAACAAAGGTTTATTGAATGAAGACAAAGAAAATGAAGTTATTGAACGTGCAAAAGCTGATATTaaa                                                                                        <  2:55661/141‑1 (MQ=255)
                                                   cGTTACAGAACTTCAGACGAAGATGCTGAATGGGAGAAAAAAGACCCATTAGTACGTTTCCGTAAATTCCTTGAAAACAAAGGTTTATGGAATGAAGACAAAGAAAATGAAGTTATTGAACGTGCAAAAGCTGATATTaaa                                                                                        <  2:27421/141‑1 (MQ=255)
                                                                 aGACGAAGATGCTGAATGGGAGAAAAAAGACCCATTAGTACGTTTCCGTAAATTCCTTGAAAACAAAGGTTTATGGAATGAAGACAAAGAAAATGAAGTTATTGAACGTGCAAAAGCTGATATTAAAGCAGCAATTAAaga                                                                          >  2:106875/1‑141 (MQ=255)
                                                                   aCGAAGATGCTGAATGGGAGAAAAAAGACCCATTAGTACGTTTCCGTAAATTCCTTGAAAACAAAGGTTTATGGAATGAAGACAAAGAAAATGAAGTTATTGAACGTGCAAAAGCTGATATTAAAGCAGCAATTAAAGAgg                                                                        <  1:67504/141‑1 (MQ=255)
                                                                            cTGAATGGGAGAAAAAAGACCCATTAGTACGTTTCCGTAAATTCCTTGAAAACAAAGGTTTATGGAATGAAGACAAAGAAAATGAAGTTATTGAACGTGCAAAAGCTGATATTAAAGCAGCAATTAAAGAGGCTGATAaca                                                               >  1:160130/1‑141 (MQ=255)
                                                                                  gggAGAAAAAAGACCCATTAGTACGTTTCCGTAAATTCCTTGAAAACAAAGGTTTATGGAATGAAGACAAAGAAAATGAAGTTATTGAACGTGCAAAAGCTGATATTAAAGCAGCAATTAAAGAGGCTGATAACACTGaaa                                                         <  2:276766/141‑1 (MQ=255)
                                                                                             gACCCATTAGTACGTTTCCGTAAATTCCTTGAAAACAAAGGTTTATGGAATGAAGACAAAGAAAATGAAGTTATTGAACGTGCAAAAGCTGATATTAAAGCAGCAATTAAAGAGGCTGATAACACTGAAAAACAAACTGtt                                              >  2:246711/1‑141 (MQ=255)
                                                                                              aCCCATTAGTACGTTTCCGTAAATTCCTTGAAAACAAAGGTTTATGGAATGAAGACAAAGAAAATGAAGTTATTGAACGTGCAAAAGCAGATATTAAAGCAGCAATTAAAGAGGCTGATAACACTGAAAAACAAACTGTTa                                             <  1:79118/141‑1 (MQ=255)
                                                                                                     aGTACGTTTCCGTAAATTCCTTGAAAACAAAGGTTTGTGGAATGAAGACAAAGAAAATGAAGTTATTGAACGTGCAAAAGCTGATATTAAAGCAGCAATTAAAGAGGCTGATAACACTGAAAAACAAACTGTTACTTctct                                      <  1:36261/141‑1 (MQ=255)
                                                                                                     aGTACGTTTCCGTAAATTCCTTGAAAACAAAGGTTTGTGGAATGAAGACAAAGAAAATGAAGTTATTGAACGTGCAAAAGCTGATATTAAAGCAGCAATTAAAGAGGCTGATAACACTGAAAAACAAACTGTTACTTctct                                      <  1:36247/141‑1 (MQ=255)
                                                                                                     aGTACGTTTCCGTAAATTCCTTGAAAACAAAGGTTTATGGAATGAAGACAAAGAAAATGAAGTTATTGAACGTGCAAAAGCTGATATTAAAGCAGCAATTAAAGAGGCTGATAACACTGAAAAACAAACTGTTACTTctct                                      <  1:111067/141‑1 (MQ=255)
                                                                                                      gTACGTTTCCGTAAATTCCTTGAAAACAAAGGTTTATGGAATGAAGACAAAGAAAATGAAGTTATTGAACGTGCAAAAGCTGATATTAAAGCAGCAATTAAAGAGGCTGATAACACTGAAAAACAAACTGTTACTTCTCTa                                     <  1:325357/141‑1 (MQ=255)
                                                                                                         cGTTTCCGTAAATTCCTTGAAAACAAAGGTTTATGGAATGAAGACAAAGAAAATGAAGTTATTGAACGTGCAAAAGCTGATATTAAAGCAGCAATTAAAGAGGCTGATAACACTGAAAAACAAACTGTTACTTCTCTAATg                                  >  1:209471/1‑141 (MQ=255)
                                                                                                           tttCCGTAAATTCCTTGAAAACAAAGGTTTATGGAATGAAGACAAAGAAAATGAAGTTATTGAACGTGCAAAAGCTGATATTAAAGCAGCAATTAAAGAGGCTGATAACACTGAAAAACAAACTGTTACTTCTCTAATGGa                                >  1:207742/1‑141 (MQ=255)
                                                                                                                gTAAATTCCTTGAAAACAAAGGTTTATGGAATGAAGACAAAGAAAATGAAGTTATTGAACGTGCAAAAGCTGATATTAAAGCAGCAATTAAAGAGGCTGATAACACTGAAAAACAAACTGTTACTTCTCTAATGGAAATta                           >  1:270517/1‑141 (MQ=255)
                                                                                                                  aaaTTCCTTGAAAACAAAGGTTTATGGAATGAAGACAAAGAAAATGAAGTTATTGAACGTGCAAAAGCTGATATTAAAGCAGCAATTAAAGAGGCTGATAACACTGAAAAACAAACTGTTACTTCTCTAATGGAAATtatg                         <  2:342075/141‑1 (MQ=255)
                                                                                                                  aaaTTCCGTGAAAACAAAGGTTTATGGAATGAAGACAAAGAAAATGAAGTTATTGAACGTGCAAAAGCTGATATTAAAGCAGCAATTAAAGAGGCTGATAACACTGAAAAACAAACTGTTACTTCTCTAATGGAAATtatg                         <  1:149169/141‑1 (MQ=255)
                                                                                                                        cTTGAAAACAAAGGTTTATGGAATGAAGACAAAGAAAATGAAGTTATTGAACGTGCAAAAGCTGATATTAAAGCAGCAATTAAAGAGGCTGATAACACTGAAAAACAAACTGTTACTTCTCTAATGGAAATTATGTATGaa                   >  2:10946/1‑141 (MQ=255)
                                                                                                                        cTTGAAAACAAAGGTTTATGGAATGAAGACAAAGAAAATGAAGTTATTGAACGTGCAAAAGCTGATATTAAAGCAGCAATTAAAGAGGCTGATAACACTGAAAAACAAACTGTTACTTCTCTAATGGAAATTATGTATGaa                   >  2:289634/1‑141 (MQ=255)
                                                                                                                                cAAAGGTTTATGGAATGAAGACAAAGAAAATGAAGTTATTGAACGTGCAAAAGCTGATATTAAAGCAGCAATTAAAGAGGCTGATAACACTGAAAAACAAACTGTTACTTCTCTAATGGAAATTATGTATGAAGATATGcc           <  1:53979/141‑1 (MQ=255)
                                                                                                                                  aaGGTTTATGGAATGAAGACAAAGAAAATGAAGTTATTGAACGTGCAAAAGCTGATATTAAAGCAGCAATTAAAGAGGCTGATAACACTGAAAAACAAACTGTTACTTCTCTAATGGAAATTATGTATGAAGATATGCCTc         >  1:269250/1‑141 (MQ=255)
                                                                                                                                   aGGTTTATGGAATGAAGACAAAGAAAATGAAGTTATTGAACGTGCAAAAGCTGATATTAAAGCAGCAATTAAAGAGGCTGATAACACTGAAAAACAAACTGTTACTTCTCTAATGGAAATTATGTATGAAGATATGCCTCa        >  2:288046/1‑141 (MQ=255)
                                                                                                                                         aTGGAATGAAGACAAAGAAAATGAAGTTATTGAACGTGCAAAAGCTGATATTAAAGCAGCAATTAAAGAGGCTGATAACACTGAAAAACAAACTGTTACTTCTCTAATGGAAATTATGTATGAAGATATGCCTCAAAACtt  <  1:10946/141‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                            
ACTATGACATATCGTTATGGTCCTCATACAATGGCTGGTGACGATCCAACTCGTTACAGAACTTCAGACGAAGATGCTGAATGGGAGAAAAAAGACCCATTAGTACGTTTCCGTAAATTCCTTGAAAACAAAGGTTTATGGAATGAAGACAAAGAAAATGAAGTTATTGAACGTGCAAAAGCTGATATTAAAGCAGCAATTAAAGAGGCTGATAACACTGAAAAACAAACTGTTACTTCTCTAATGGAAATTATGTATGAAGATATGCCTCAAAACTT  >  USA300TCH1516_ALE/1103271‑1103548

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: