Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 638,614 G→T 100% T14K (ACA→AAA)  pdxK ← Pyridoxine kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE638,6140GT100.0% 118.4 / NA 35T14K (ACA→AAA) pdxKPyridoxine kinase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (20/15);  total (20/15)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CATCCATTGGTAATGGTGTAACATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTGTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCAATATTTGATGATTGTTGAAA  >  USA300TCH1516_ALE/638477‑638749
                                                                                                                                         |                                                                                                                                       
cATCCATTGGTAATTGTGTAACATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTttgt                                                                                                                                      <  1:464363/141‑1 (MQ=255)
       tGGTAATGGTGTAACATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAAcc                                                                                                                               <  1:276707/141‑1 (MQ=255)
       tGGTAATGGTGTAACATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAAcc                                                                                                                               <  1:554841/141‑1 (MQ=255)
         gTAATGGTGTAACATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCGCTTTTGTCGGAGCCgg                                                                                                                             >  1:9032/1‑141 (MQ=255)
               gtgtAACATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTg                                                                                                                       >  1:329261/1‑141 (MQ=255)
                tgtAACATCGTGTGACCATGTATCTTTATCCATATCCACGACGGCAGTTAAAGCGACCATGCCATACGTATCCAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGt                                                                                                                      <  1:517461/141‑1 (MQ=255)
                 gtAACATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCTCTTTTGTCAGAACCGGCAATTGtt                                                                                                                     <  1:490791/141‑1 (MQ=255)
                 gtAACATCGTGTGACCATGTATCTCTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGtt                                                                                                                     <  1:249801/141‑1 (MQ=255)
                      aTCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAAc                                                                                                                >  1:410821/1‑141 (MQ=255)
                       tCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGGAATTGTTAAAAct                                                                                                               >  1:321925/1‑141 (MQ=255)
                        cGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAActt                                                                                                              >  1:167959/1‑141 (MQ=255)
                           gtgACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAActttct                                                                                                           >  1:529395/1‑141 (MQ=255)
                           gtgACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAActttct                                                                                                           >  1:529275/1‑141 (MQ=255)
                                   gTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTGTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGcc                                                                                                   <  1:245731/141‑1 (MQ=255)
                                      tCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTACTTGGGTACTGGCTCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCAtt                                                                                                >  1:659166/1‑141 (MQ=255)
                                         ttATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAg                                                                                             <  1:215780/141‑1 (MQ=255)
                                          tATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCGTACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGc                                                                                            <  1:342847/141‑1 (MQ=255)
                                                      cGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACa                                                                                >  1:469263/1‑141 (MQ=255)
                                                                tAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACAGAATAATAtt                                                                      >  1:418490/1‑141 (MQ=255)
                                                                          aTGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATACTATTGTATtcatca                                                            >  1:398635/1‑141 (MQ=255)
                                                                           tGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTcatcat                                                           >  1:159956/1‑141 (MQ=255)
                                                                             ccATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTATATAATAATATTGTATTCatcatat                                                         <  1:430505/141‑1 (MQ=255)
                                                                             ccATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCatcatat                                                         >  1:567650/1‑141 (MQ=255)
                                                                              cATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCATCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCatcatatc                                                        <  1:197517/141‑1 (MQ=255)
                                                                                      tCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTaa                                                <  1:453593/141‑1 (MQ=255)
                                                                                                 gAACGTTTTCAAAGCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAAAAATATTGTATTCATCATATCATTTTTAACCTAATTGaaa                                     <  1:159111/141‑1 (MQ=255)
                                                                                                  aaCGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGaaaa                                    >  1:514259/1‑141 (MQ=255)
                                                                                                    cGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTAAACCCTGACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAat                                  >  1:458323/1‑141 (MQ=255)
                                                                                                        ttCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTa                              <  1:360501/141‑1 (MQ=255)
                                                                                                          cAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAg                            >  1:330165/1‑141 (MQ=255)
                                                                                                               cTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCCTTTTTAACCTAATTGAAAAATATTAAGCATTc                       >  1:590277/1‑141 (MQ=255)
                                                                                                                tGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCa                      >  1:608969/1‑141 (MQ=255)
                                                                                                                           cTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCAATATTtgatga           <  1:294208/141‑1 (MQ=255)
                                                                                                                               gcCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCAATATTTGATGAttgt       <  1:148936/141‑1 (MQ=255)
                                                                                                                                   gCACTTTTGTCAGAACGGGCAATTGTTAAAACTTTCTTAAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGCAAAATATTAAGCATTCAATATTTGATGATTGTTGaa   >  1:374420/1‑141 (MQ=255)
                                                                                                                                    cACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCAATATTTGATGATTGTTGaaa  >  1:569689/1‑141 (MQ=255)
                                                                                                                                         |                                                                                                                                       
CATCCATTGGTAATGGTGTAACATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTGTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCAATATTTGATGATTGTTGAAA  >  USA300TCH1516_ALE/638477‑638749

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