Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,425,673 A→G 11.6% F837L (TTT→CTT)  USA300HOU_RS12500 ← Putative formate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE2,425,6730AG11.6% 137.5 / 11.2 52F837L (TTT→CTT) USA300HOU_RS12500Putative formate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (25/21);  new base G (3/3);  total (28/24)
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.23e-01

CAACAAATGCACGTGGCACTTTGTATTTAATCATTGGTGTTTGATAAGTCATATTACCTTCATGGAAGTGCTCTAACAATCGACCATTGTTTACATGAATATCATAAATTTCATCTTGCTTAAAGTAATTATCAAATGATAATGGGAATAATTTTGCTTTACCATTATCAAAATTGAATCCTTCTAAGTATAGAATAGGCTCATCAGTACCATCAGGTTGTACTGGCCATTGTAAACTATTGAATCCTTCTAAACGATCATAACTTACCCCAGC  >  USA300TCH1516_ALE/2425538‑2425811
                                                                                                                                       |                                                                                                                                          
caacaaATGCACGTGGCACTTTGTATTTAATCATTGGTGTTTGATAAGTCATATTACCTTCATGGAAGTGCTCTAACAATCGACCATTGTTTACATGAATATCATAAATTTCATCTTGCTTAAAGTAATTATCAAGTGATa                                                                                                                                       >  2:189443/1‑141 (MQ=255)
caacaaATGCACGTGGCACTTTGTATTTAATCATTGGTGTTTGATAAGTCATATTACCTTCATGGAAGTGCTCTAACAATCGACCATTGTTTACATGAATATCATAAATTTCATCTTGCTTAAAGTAATTATCAAGTGATa                                                                                                                                       >  2:30449/1‑141 (MQ=255)
caacaaATGCACGTGGCACTTTGTATTTAATCATTGGTGTTTGATAAGTCATATTACCTTCATGGAAGTGCTCTAACAATCGACCATTGTTTACATGAATATCATAAATTTCATCTTGCTTAAAGTAATTATCAAGTGATa                                                                                                                                       >  2:66725/1‑141 (MQ=255)
       tGCACGTGGCACTTTGTATTTAATCATTGGTGTTTGATAAGTCATATTACCTTCATGGAAGTGCTCTAACAATCGACCATTGTTTACATGAATATCATAAATTTCATCTTGCTTAAAGTAATTATCAAATGATAATGGGaa                                                                                                                                >  2:369216/1‑141 (MQ=255)
       tGCACGTGGCACTTTGTATTTAATCATTGGTGTTTGATAAGTCATATTACCTTCATGGAAGTGCTCTAACAATCGACCATTGTTTACATGAATATCATAAATTTCATCTTGCTTAAAGTAATTATCAAATGATAATGGGaa                                                                                                                                <  2:372527/141‑1 (MQ=255)
        gCACGTGGCACTTTGTATTTAATCATTGGTGTTTGATAAGTCATATTACCTTCATGGAAGTGCTCTAACAATCGACCATTGTTTACATGAATATCATAAATTTCATCTTGCTTAAAGTAATTATCAAATGATAATGGGaat                                                                                                                               >  1:412786/1‑141 (MQ=255)
        gCACGTGGCACTTTGTATTTAATCATTGGTGTTTGATAAGTCATATTACCTTCATGGAAGTGCTCTAACAATCGACCATTGTTTACATGAATATCATAAATTTCATCTTGCTTAAAGTAATTATCAAATGATAATGGGaat                                                                                                                               >  1:110463/1‑141 (MQ=255)
        gCACGTGGCACTTTGTATTTAATCATTGGTGTTTGATAAGTCATATTACCTTCATGGAAGTGCTCTAACAATCGACCATTGTTTACATGAATATCATAAATTTCATCTTGCTTAAAGTAATTATCAAATGATAATGGGaat                                                                                                                               >  1:287762/1‑141 (MQ=255)
          aCGTGGCACTTTGTATTTAATCATTGGTGTTTGATAAGTCATATTACCTTCATGGAAGTGCTCTAACAATCGACCATTGTTTACATGAATATCATAAATTTCATCTTGCTTAAAGTAATTATCAAATGATAATGGGaataa                                                                                                                             <  1:331581/141‑1 (MQ=255)
                                 ttGGTGTTTGATAAGTCATATTACCTTCATGGAAGTGCTCTAACAATCGACCATTGTTTACATGAATATCATAAATTTCATCTTGCTTAAAGTAATTATCAAATGATAATGGGAATAATTTTGCTTTACCATTATCAAAAt                                                                                                      <  2:281494/141‑1 (MQ=255)
                                   ggTGTTTGATAAGTCATATTACCTTCATGGAAGTGCTCTAACAATCGACCATTGTTTACATGAATATCATAAATTTCATCTTGCTTAAAGTAATTATCAAATGATAATGGGAATAATTTTGCTTTACCATTATCAAAATTg                                                                                                    <  1:127433/141‑1 (MQ=255)
                                   ggTGTTTGATAAGTCATATTACCTTCATGGAAGTGCTCTAACAATCGACCATTGTTTACATGAATATCATAAATTTCATCTTGCTTAAAGTAATTATCAAATGATAATGGGAATAATTTTGCTTTACCATTATCAAAATTg                                                                                                    <  1:375740/141‑1 (MQ=255)
                                      gtTTGATAAGTCATATTACCTTCATGGAAGTGCTCTAACAATCGACCATTGTTTACGTGAGTAGCGTAAATAGCATCTTGCTTAAATGAAGTAGCAAATGATAATGTGAATAATTTTGCTTTACCATTATCAAAATTGAAt                                                                                                 >  1:275482/1‑141 (MQ=255)
                                      gtTTGATAAGTCATATTACCTTCATGGAAGTGCTCTAACAATCGACCATTGTTTACATGAATATCATAAATTTCATCTTGCTTAAAGTAATTATCAAATGATAATGGGAATAATTTTGCTTTACCATTATCAAAATTGAAt                                                                                                 >  2:66137/1‑141 (MQ=255)
                                       tttGATAAGTCATATTACCTTCATGGAAGTGCTCTAACAATCGACCATTGTTTACATGAATATCATAAATTTCATCTTGCTTAAAGTAATTATCAAATGATAATGGGAATAATTTTGCTTTACCATTATCAAAATTGAATc                                                                                                <  2:147321/141‑1 (MQ=255)
                                       tttGATAAGTCATATTACCTTCATGGAAGTGCTCTAACAATCGACCATTGTTTACATGAATATCATAAATTTCATCTTGCTTAAAGTAATTATCAAATGATAATGGGAATAATTTTGCTTTACCATTATCAAAATTGAATc                                                                                                <  2:147288/141‑1 (MQ=255)
                                           aTAAGTCATATTACCTTCATGGAAGTGCTCTAACAATCGACCATTGTTTACATGAATATCATAAATTTCATCTTGCTTAAAGTAATTATCAAATGATAATGGGAATAATTTTGCTTTACCATTATCAAAATTGAATCCTTc                                                                                            <  2:63533/141‑1 (MQ=255)
                                           aTAAGTCATATTACCTTCATGGAAGTGCTCTAACAATCGACCATTGTTTACATGAATATCATAAATTTCATCTTGCTTAAAGTAATTATCAAATGATAATGGGAATAATTTTGCTTTACCATTATCAAAATTGAATCCTTc                                                                                            <  2:70112/141‑1 (MQ=255)
                                            tAAGTCATATTACCTTCATGGACGTGCTCTAACAATCGACCATTGTTTACATGAATATCATAAATTTCATCTTGCTTAAAGTAATTATCAAATGATAATGGGAATAATTTTGCTTTACCATTATCAAAATTGAATCCTTCt                                                                                           >  2:283324/1‑141 (MQ=255)
                                            tAAGTCATATTACCTTCATGGAAGTGCTCTAACAATCGACCATTGTTTACATGAATATCATAAATTTCATCTTGCTTAAAGTAATTATCAAATGATAATGGGAATAATTTTGCTTTACCATTATCAAAATTGAATCCTTCt                                                                                           >  2:73591/1‑141 (MQ=255)
                                            tAAGTCATATTACCTTCATGGAAGTGCTCTAACAATCGACCATTGTTTACATGAATATCATAAATTTCATCTTGCTTAAAGTAATTATCAAATGATAATGGGAATAATTTTGCTTTACCATTATCAAAATTGAATCCTTCt                                                                                           >  2:214145/1‑141 (MQ=255)
                                               gTCATATTACCTTCATGGAAGTGCTCTAACAATCGACCATTGTTTACATGAATATCATAAATTTCATCTTGCTTAAAGTAATTATCAAATGATAATGGGAATAATTTTGCTTTACCATTATCAAAATTGAATCCTTCTAAg                                                                                        <  1:177166/141‑1 (MQ=255)
                                                      tACCTTCATGGAAGTGCTCTAACAATCGACCATTGTTTACATGAATATCATAAATTTCATCTTGCTTAAAGTAATTATCAAATGATAATGGGAATAATTTTGCTTTACCATTATCACAATTGAAACCTTCTAAGTATAGaa                                                                                 >  2:97849/1‑141 (MQ=255)
                                                      tACCTTCATGGAAGTGCTCTAACAATCGACCATTGTTTACATGAATATCATAAATTTCATCTTGCTTAAAGTAATTATCAAATGATAATGGGAATAATTTTGCTTTACCATTATCAAAATTGAATCCTTCTAAGTATAGaa                                                                                 >  2:371660/1‑141 (MQ=255)
                                                                   gTGCTCTAACAATCGACCATTGTTTACATGAATATCATAAATTTCATCTTGCTTAAAGTAATTATCAAATGATAATGGGAATAATTTTGCTTTACCATTATCAAAATTGAATCCTTCTAAGTATAGAATAGGCTCATCAGt                                                                    <  2:100337/141‑1 (MQ=255)
                                                                        ctAACAATCGACCATTGTTTACGTGAATATCATAAATTTCATCTTGCTTAAAGTAATTATCAAATGATAATGGGAATAATTTTGCTTTACCATTATCAAAATTGAATCCTTCTAAGTATAGAATAGGCTCATCAGTACCAt                                                               >  2:61896/1‑141 (MQ=255)
                                                                        ctAACAATCGACCATTGTTTACATGAATATCATAAATTTCATCTTGCTTAAAGTAATTATCAAATGATAATGGGAATAATTTTGCTTTACCATTATCAAAATTGAATCCTTCTAAGTATAGAATAGGCTCATCAGTACCAt                                                               >  2:58254/1‑141 (MQ=255)
                                                                        ctAACAATCGACCATTGTTTACATGAATATCATAAATTTCATCTTGCTTAAAGTAATTATCAAATGATAATGGGAATAATTTTGCTTTACCATTATCAAAATTGAATCCTTCTAAGTATAGAATAGGCTCATCAGTACCAt                                                               >  1:200256/1‑141 (MQ=255)
                                                                        ctAACAATCGACCATTGTTTACATGAATATCATAAATTTCATCTTGCTTAAAGTAATTATCAAATGATAATGGGAATAATTTTGCTTTACCATTATCAAAATTGAATCCTTCTAAGTATAGAATAGGCTCATCAGTACCAt                                                               >  2:199675/1‑141 (MQ=255)
                                                                        ctAACAATCGACCATTGTTTACATGAATATCATAAATTTCATCTTGCTTAAAGTAATTATCAAATGATAATGGGAATAATTTTGCTTTACCATTATCAAAATTGAATCCTTCTAAGTATAGAATAGGCTCATCAGTACCAt                                                               >  2:85249/1‑141 (MQ=255)
                                                                             aaTCGACCATTGTTTACATGAATATCATAAATTTCATCTTGCTTAAAGTAATTATCAAATGATAATGGGAATAATTTTGCTTTACCATTATCAAAATTGAATCCTTCTAAGTAGAGCATAGGCTCATCAGTACCATCAGGt                                                          <  2:157315/141‑1 (MQ=255)
                                                                                  aCCATTGTTTACATGAATATCATAAATTTCATCTTGCTTAAAGTAATTATCAAATGATAATGGGAATAATTTTGCTTTACCATTATCAAAATTGAATCCTTCTAAGTATAGAATAGGCTCATCAGTACCATCAGGTTGTAc                                                     >  1:410028/1‑141 (MQ=255)
                                                                                       tGTTTACATGAATATCATAAATTTCATCTTGCTTAAAGTAATTATCAAATGATAATGGGAATAATTTTGCTTTACCATTCTCAAAATTGAATCCTTCTAAGTATAGAATAGGCTCACCAGTACCATCAGGTTGTACTGGcc                                                >  2:341947/1‑141 (MQ=255)
                                                                                       tGTTTACATGAACATCATAAATTTCACCTTGCTTAAAGTAAGTATCAAATGATAATGGGAATAATTTTGCTTTACCATTATCAAAATTGAATCCTTCTAAGTATAGAATACGCTCCTCACTACCCTCAGGTTGTACTggtc                                                >  2:181328/1‑139 (MQ=255)
                                                                                        gTTTACATGAATATCATAAATTTCATCTTGCTTAAAGTAATTATCAAATGATAATGGGAATAATTTTGCTTTACCATTATCAAAATTGAATCCTTCTAAGTATAGAATAGGCTCATCAGTACCATCAGGTTGTACTGGCCa                                               >  1:427798/1‑141 (MQ=255)
                                                                                           tACATGAATATCATAAATTTCATCTTGCTTAAAGTAATTATCAAATGATAATGGGAATAATTTTGCTTTACCATTATCAAAATTGAATCCTTCTAAGTATAGAATAGGCTCATCAGTACCATCAGGTTGTACTGGCCATTg                                            <  1:390543/141‑1 (MQ=255)
                                                                                           tACATGAATATCATAAATTTCATCTTGCTTAAAGTAATTATCAAATGATAATGGGAATAATTTTGCTTTACCATTATCAAAATTGAATCCTTCTAAGTATAGAATAGGCTCATCAGTACCATCAGGTTGTACTGGCCATTg                                            <  2:4590/141‑1 (MQ=255)
                                                                                                gAATATCATAAATTTCATCTTGCTTAAAGTAATTATCAAATGATAATGGGAATAATTTTGCTTTACCATTATCAAAATTGAATCCTTCTAAGTATAGAATAGGCTCATCAGTACCATCAGGTTGTACTGGCCATTGTAAAc                                       <  1:237401/141‑1 (MQ=255)
                                                                                                 aaTATCATAAATTTCATCTTGCTTAAAGTAATTATCAAATGATAATGGGAATAATTTTGCTTTACCATTATCAAAATTGAATCCTTCTAAGTATAGAATAGGCTCATCAGTACCATCAGGTTGTACTGGCCATTGTAAACt                                      <  1:81678/141‑1 (MQ=255)
                                                                                                 aaTATCATAAATTTCATCTTGCTTAAAGTAATTATCAAATGATAATGGGAATAATTTTGCTTTACCATTATCAAAATTGAATCCTTCTAAGTATAGAATAGGCTCATCAGTACCATCAGGTTGTACTGGCCATTGTAAACt                                      <  1:36425/141‑1 (MQ=255)
                                                                                                         aaaTTTCATCTTGCTTAAAGTAATTATCAAATGATAATGGGAATAATTTTGCTTTACCATTATCAAAATTGAATCCTTCTAAGTATAGAATAGGCTCATCAGTACCATCAGGTTGTACTGGCCATTGTAAACTATTGAATc                              >  2:257357/1‑141 (MQ=255)
                                                                                                         aaaTTTCATCTTGCTTAAAGTAATTATCAAATGATAATGGGAATAATTTTGCTTTACCATTATCAAAATTGAATCCTTCTAAGTATAGAATAGGCTCATCAGTACCATCAGGTTGTACTGGCCATTGTAAACTATTGAATc                              >  2:79292/1‑141 (MQ=255)
                                                                                                           aTTTCATCTTGCTTAAAGTAATTATCAAGTGATAATGGGAATAATTTTGTTTTACCATTATCAAAATTGAATCCTTCTAAGTATAGAATAGGCTCATCAGTACCATCAGGTTGTACTGGCCATTGTAAACTATTGAATCCt                            <  1:66725/141‑1 (MQ=255)
                                                                                                           aTTTCATCTTGCTTAAAGTAATTATCAAGTGATAATGGGAATAATTTTGCTTTACCATTATCAAAATTGAATCCTTCTAAGTATAGAATAGGCTCATCAGTACCATCAGGTTGTACTGGCCATTGTAAACTATTGAATCCt                            <  1:30449/141‑1 (MQ=255)
                                                                                                           aTTTCATCTTGCTTAAAGTAATTATCAAGTGATAATGGGAATAATTTTGCTTTACCATTATCAAAATTGAATCCTTCTAAGTATAGAATAGGCTCATCAGTACCATCAGGTTGTACTGGCCATTGTAAACTATTGAATCCt                            <  1:189443/141‑1 (MQ=255)
                                                                                                              tCATCTTGCTTAAAGTAATTATCAAATGATAATGGGAATAATTTTGCTTTACCATTATCAAAATTGAATCCTTCTAAGTATAGAATAGGCTCATCAGTACCATCAGGTTGTACTGGCCATTGTAAACTATTGAATCCTTCt                         >  1:29402/1‑141 (MQ=255)
                                                                                                                   ttGCTTAAAGTAATTATCAAATGATAATGGGAATAATTTTGCTTTACCATTATCAAAATTGAATCCTTCTAAGTATAGAATAGGCTCATCAGTACCATCAGGTTGTACTGGCCATTGTAAACTATTGAATCCTTCTAAACg                    >  2:147251/1‑141 (MQ=255)
                                                                                                                            gTAATTATCAAATGATAATGGGAATAATTTTGCTTTACCATTATCAAAATTGAATCCTTCTAAGTATAGAATAGGCTCATCAGTACCATCAGGTTGTACTGGCCATTGTAAACTATTGAATCCTTCTAAACGATCATAACt           <  1:369216/141‑1 (MQ=255)
                                                                                                                              aaTTATCAAATGATAATGGGAATAATTTTGCTTTACCATTATCAAAATTGAATCCTTCTAAGTATAGAATAGGCTCATCAGTACCATCAGGTTGTACTGGCCATTGTAAACTATTGAATCCTTCTAAACGATCATAACTTa         <  2:46346/141‑1 (MQ=255)
                                                                                                                              aaTTATCAAATGATAATGGGAATAATTTTGCTTTACCATTATCAAAATTGAATCCTTCTAAGTATAGAATAGGCTCATCAGTACCATCAGGTTGTACTGGCCATTGTAAACTATTGAATCCTTCTAAACGATCATAACTTa         <  1:303344/141‑1 (MQ=255)
                                                                                                                              aaTTATCAAATGATAATGGGAATAATTTTGCTTTACCATTATCAAAATTGAATCCTTCTAAGTATAGAATAGGCTCATCAGTACCATCAGGTTGTACTGGCCATTGTAAACTATTGAATCCTTCTAAACGATCATAACTTa         <  2:42812/141‑1 (MQ=255)
                                                                                                                                  aTCAAATGATAATGGGAATAATTTTGCTTTACCATTATCAAAATTGAATCCTTCTAAGTATAGAATAGGCTCATCAGTACCATCAGGTTGTACTGGCCATTGTAAACTATTGAATCCTTCTAAACGATCATAACTTAcccc     >  1:3285/1‑141 (MQ=255)
                                                                                                                                     aaaTGATAATGGGAATAATTTTGCTTTACCATTATCAAAATTGAATCCTTCTAAGTATAGAATAGGCTCATCAGTACCATCAGGTTGTACTGGCCATTGTAAACTATTGAATCCTTCTAAACGATCATAACTTACCCCAGc  >  2:367961/1‑141 (MQ=255)
                                                                                                                                       |                                                                                                                                          
CAACAAATGCACGTGGCACTTTGTATTTAATCATTGGTGTTTGATAAGTCATATTACCTTCATGGAAGTGCTCTAACAATCGACCATTGTTTACATGAATATCATAAATTTCATCTTGCTTAAAGTAATTATCAAATGATAATGGGAATAATTTTGCTTTACCATTATCAAAATTGAATCCTTCTAAGTATAGAATAGGCTCATCAGTACCATCAGGTTGTACTGGCCATTGTAAACTATTGAATCCTTCTAAACGATCATAACTTACCCCAGC  >  USA300TCH1516_ALE/2425538‑2425811

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