Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 638,335 A→G 8.3% V107A (GTC→GCC)  pdxK ← Pyridoxine kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE638,3350AG 8.3% 106.8 / 2.6 49V107A (GTC→GCC) pdxKPyridoxine kinase
Reads supporting (aligned to +/- strand):  ref base A (25/19);  new base G (2/2);  total (28/21)
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.35e-01

TTACCTAATCCTGATAATTGACCCGCTTCAAATAAGTTTGGTGTCACTACCGTTGCTTTAGGTAGTAAATATTTAATCATCGCCTCAGTATTTCCAGGATTAAGCACTTCATCTTCGCCTTTACAAACCATGACAGGATCTACTACAAAATATTGTGCATTAGATGCCTCATATACTTCTCCAGCACGTTTGATTATCTCCTCAGTACCTAACATACCTGTTTTAATAGCATCAGGTCCGATTGATAAAGCCGTTTCAAGTTGT  >  USA300TCH1516_ALE/638203‑638466
                                                                                                                                    |                                                                                                                                   
ttaCCTAATCCTGATAATTGACCCGCTTCAAATAAGTTTGGTGTCACTACCGTTGCTTTAGGTAGTAAATATTTAATCATCGCCTCAATATTTCCATGATTAAGCACTTCATCTTCGCCTTTACAAACCATGACCGGATct                                                                                                                             >  2:169686/1‑141 (MQ=255)
            gATAATTGACCCGCTTCAAATAAGTTTGGTGTCACTACCGTTGCTTTAGGTAGTAAATATTTAATCATCGCCTCAGTATTTCCAGGATTAAGCACTTCATCTTCGCCTTTACAAACCATGACAGGATCTACTACAAAatat                                                                                                                 <  2:387149/141‑1 (MQ=255)
              tAATTGACCCGCTTCAAATAAGTTTGGTGTCACTACCGTTGCTTTAGGTAGTAAATATTTAATCATCGCCTCAGTATTTCCAGGATTAAGCACTTCATCTTCGCCTTTACAAACCATGACAGGATCTACTACAAAATAttt                                                                                                               >  2:239574/1‑140 (MQ=255)
               aacttgcccGCTTCAAATAAGTTTGGTGTCACTACCGTTGCTTTAGGTAGTAAATATTTAATCATCGCCTCAGTATTTCCAGGATTAAGCACTTCATCTTCGCCTTTACAAACCATGACAGGATCTACTACAAAATATtgt                                                                                                              <  2:117331/135‑1 (MQ=255)
               aaTTGACCCGCTTCAAATAAGTTTGGTGTCACTACCGTTGCTTTAGGTAGTAAATATTTAATCATCGCCTCAGTATTTCCAGGATTAAGCACTTCATCTTCGCCTTTACAAACCATGACAGGATCTACTACAAAATATtgt                                                                                                              <  2:118859/141‑1 (MQ=255)
                   gtcccGTTTCAAATAAGTTTGGTGTCCCTACCGTTGCTTTAGGTAGTAACTATTTAATCATCGCCTCAGTATTTCCTGGATTAAGCACTTCATCTTCGCCTTTACAAACCATGACAGGATCTACTACTAAATATTGTGCAt                                                                                                          <  2:207214/139‑1 (MQ=255)
                   gACCCGCATCAAATAAGTTTGGTGTCACTACCGTTGCTTTAGGTAGTAAATATATGATGATGGTCGGAGTATTTCCAGTATTAAGCACTTCATCTTCGCCTTTACAAACCATGACAGGATATACTACAAAATATTGTGCAt                                                                                                          >  2:300083/1‑141 (MQ=255)
                    aCCCGCTTCAAATAAGTTTGGTGTCACTACCGTTGCTTTAGGTAGTAAATATTTAATCATCGCCTCAGTATTTCCAGGATTAAGCACTTCATCTTCGCCTTTACAAACCATGACAGGATCTACTACAAAATATTGTGCAtt                                                                                                         <  1:334265/141‑1 (MQ=255)
                        gCTTCAAATAAGTTTGGTGTCACTACCGTTGCTTTAGGTAGTAAATATTTAATCATCGCCTCAGTATTTCCAGGATTAAGCACTTCATCTTCGCCTTTACAAACCATGACAGGATCTACTACAAAATATTGTGCATTAGAt                                                                                                     <  2:93111/141‑1 (MQ=255)
                         cTTCAAATAAGTTTGGTGTCACTACCGTTGCTTTAGGTAGTAAATATTTAATCATCGCCTCAGTATTTCCAGGATTAAGCACTTCATCTTCGCCTTTACAAACCATGACAGGATCTACTACAAAATATTGTGCATTAGATg                                                                                                    >  1:356930/1‑141 (MQ=255)
                               aTAAGTTTGGTGTCACTACCGTTGCTTTAGGTAGTAAATATTTAATCATCGCCTCAGTATTTCCAGGATTAAGCACTTCATCTTCGCCTTTACAAACCATGACAGGATCTACTACAAAATATTGTGCATTAGATGCTTCat                                                                                              >  1:23774/1‑141 (MQ=255)
                               aTAAGTTTGGTGTCACTACCGTTGCTTTAGGTAGTAAATATTTAATCATCGCCTCAGTATTTCCAGGATTAAGCACTTCATCTTCGCCTTTACAAACCATGACAGGATCTACTACAAAATATTGTGCATTAGATGCCTCat                                                                                              >  1:205048/1‑141 (MQ=255)
                               aTAAGTTTGGTGTCACTACCGTTGCTTTAGGTAGTAAATATTTAATCATCGCCTCAGTATTTCCAGGATTAAGCACTTCATCTTCGCCTTTACAAACCATGACAGGATCTACTACAAAATATTGTGCATTAGATGCCTCat                                                                                              >  1:29448/1‑141 (MQ=255)
                                     ttGGTGTCACTACCGTTGCTTTAGGTAGTAAATATTTAATCATCGCCTCAGTATTTCCAGGATTAAGCACTTCATCTTCGCCTTTACAAACCATGACAGGATCTACTACAAAATATTGTGCATTAGATGCCTCATATACtt                                                                                        >  2:219315/1‑141 (MQ=255)
                                      tGGTGTCACTACCGTTGCTTTAGGTAGTAAATATTTAATCATCGCCTCAGTATTTCCAGGATTAAGCACTTCATCTTCGCCTTTACAAACCATGACAGGATCTACTACAAAATATTGTGCATTAGATGCCTCATATACTtc                                                                                       >  1:71065/1‑141 (MQ=255)
                                            cACTGCCGTTGCTTTAGGTAGTAAATATTTAATCATCGCCTCAGTATTTCCAGGATTAAGCACTTCATCTTCGCCTTTACAAACCATGACAGGATCTACTACAAAATATTGTGCATTAGATGCCTCATATACTTCTCCAGc                                                                                 <  2:337844/141‑1 (MQ=255)
                                               tACAGTGGCTTTAAGTAGTAAATATTTAATCATCGCCTCAGTATTTCCAGGATGAAGCACTTCATCTTCGCCTTTACAAACCATGACAGGATCTACTACAAAATATTGTGCATTAGATGCCTCATATACTTCTCCAGCACg                                                                              <  2:338102/141‑1 (MQ=255)
                                                aCCGTTGCTTTAGGTAGTAAATATTTAATCATCGCCTCAGTATTTCCAGGATTAAGCACTTCATCTTCGCCTTTAAAAACAATTACTTGATCTACTACAAAATATTGTGCAATAGATGCCTCATATACTTCGCCAGCACGt                                                                             >  2:45587/1‑141 (MQ=255)
                                                  cGTTGCTTTAGGTAGTAAATATTTAATCATCGCCTCAGTATTTCCAGGATTAAGCACTTCATCTTCGCCTTTACAAACCATGACAGGATCTACTACAAAATATTGTGCATTAGATGCCTCATATACTTCTCCAGCACGttt                                                                           <  1:299222/141‑1 (MQ=255)
                                                    ttGCTTTAGGTAGTAAATATTTAATCATCGCCTCAGTATTTCCAGGATTAAGCACTTCATCTTCGCCTTTACAAACCATGACAGGATCTACTACAAAATATTGTGCATTAGATGCCTCATATACTTCTCCAGCACGTTTGa                                                                         >  1:217572/1‑141 (MQ=255)
                                                      gCTTTAGGGAGTAAATATTTAATCATCGCCTCAGTATTTCCAGGATTAAGCACTTCATCTTCGCCTTTACAAACCATGACAGGATCTACTACAAAATATTGTGCATTAGATGCCTCATATACTTCTCCAGCACGTTTGAtt                                                                       >  1:245908/1‑141 (MQ=255)
                                                        tttAGGTAGTAAATATTTAATCATCGCCTCAGTATTTCCAGGATTAAGCACTTCATCTTCGCCTTTACAAACCATGACAGGATCTACTACAAAATATTGTGCATTACATGACTCATATACTTCTCCACCACGTTTGATTAt                                                                     >  2:153016/1‑141 (MQ=255)
                                                           aGGTAGTAAATATTTAATCATCGCCTCAGTATTTCCAGGATTAAGCACTTCATCTTCGCCTTTACAAACCATGACAGGATCTACTACAAAATATTGTGCATTAGATGCCTCATATACTTCTCCAGCACGTTTGATTAtctc                                                                  >  2:35555/1‑141 (MQ=255)
                                                            ggTATTAAATATTTAATCATCGCCTTAGTATTTCCAGGATTAAGCACTTCATCTTCGCCTTTACAAACCATGACAGAAGCTATTACCAAATGTATTCTGGTAGATGCCTCATATCCTCCTTGATTACTTTTGATTATctcc                                                                 <  2:386388/141‑1 (MQ=255)
                                                            ggTAGTAAATATTTAATCATCGCCTCAGTATTTCCAGGATTAAGCACTTCATCTTCGCCTTTACAAACCAGGACAGGATCTACTACAAAATATTGTGCATTTGATGCCTCATCTACTTTTCCGGCACCTTTGATTAtctct                                                                 >  1:67355/1‑140 (MQ=255)
                                                                   aaTATTCAATCATCGGCTCAGTCTTTCCATGATTAAGCACTTCATCTTCGCCTTTCCAATCCATGACAGCATCTACTACAAAATATTGTGCATTAGATGCCTCATATACTTCTCCAGCCCGTTTTATTATCTCCTCAGTAc                                                          <  2:359916/141‑1 (MQ=255)
                                                                      atTTAATCATCGCCTCAGTATTTCCAGGATTAAGCACTTCATCTTCGCCTTTACAAACCATGACAGGATCTACTACAAAATATTGTGCATTAGATGCCTCATATACTTCTCCAGCACGTTTGATTATCTCATCAGTACCTa                                                       >  1:162986/1‑141 (MQ=255)
                                                                              cgcGCCTCAGTGTATCCAGGATTAATCACTTCCTCTTCGCATTTACAAACCATGACAGGATCTACTACAAAATATTGTGCATTAGATGCCTCATATACTTCTCCAGCACGTTTGATTATCTCCTCAGTACCTAACATAACt                                               <  2:316153/139‑1 (MQ=255)
                                                                               tcGCCTCAGTATTTCCAGGATTAAGCACTTCATCTTCGCCTTTACAAACCATGGCAGGATCTACTACAAAATATTGTGCATTAGATGCCTCATATACTTCTCCAGCACGTTTGATTATCTCCTCAGTACCTAACATACCTg                                              >  1:301594/1‑141 (MQ=255)
                                                                               tcGCCTCAGTATTTCCAGGATTAAGCACTTCATCTTCGCCTTTACAAACCATGGCAGGATCTACTACAAAATATTGTGCATTAGATGCCTCATATACTTCTCCAGCACGTTTGATTATCTCCTCAGTACCTAACATACCTg                                              >  1:302101/1‑141 (MQ=255)
                                                                                cGCCTCAGTATTTCCAGGATTAAGCACTTCATCTTCGCCTTTACAAACCATGACAGGATCTACTACAAAATATTGTGCATTAGATGCCTCATATACTTCTCCAGCACGTTTGATTATCTCCTCAGTACCTAACATACCTGt                                             >  1:199810/1‑141 (MQ=255)
                                                                                 gCCTCAGTATTTCCAGGATTAAGCACTTCATCTTCGCCTTTACAAACCATGACAGGATCTACTACAAAATATTGTGCATTAGATGCCTCATATACTTCTCCAGCACGTTTGATTATCTCCTCAGTACCTAACATACCTGtt                                            <  2:90772/141‑1 (MQ=255)
                                                                                  ccTAAGTATTTCCAGGATTAAGCATTTCATCTTCGCCTTTACAAACCATGACAGGATCTACTACAAAATATTGTGCATTAGATGCCTCATATACTTCTCCAGCACGTTTGATTATCTCCTCAGTACCTAACATACCTGttt                                           <  1:359687/141‑1 (MQ=255)
                                                                                   cTCAGTATTTCCAGGATTAAGCACTTCATCTTCGCCTTTACAAACCATGACAGGATCTACTACAAAATATTGTCCATTAGATGCCTCATATACTTCTCCAGCACATTTGATTATGTCCTCAGTACCTAACATCCCTGtttt                                          >  2:151006/1‑141 (MQ=255)
                                                                                    tCAGTATTTCCAGGATTAAGCACTTCATCTTCGCCTTTACAAACAATGACAGGATCTACTACAAAATATTCTGCATTAGCTGCCGCATCTACTTCTCCAGCACGTTTGATTATCTCCTTCTTATCTTACCTACCTTTTTTa                                         >  2:362217/1‑141 (MQ=255)
                                                                                      aGTATTTACAGGATTAAGCACTTCATCTTCGCCTTTACAAACCATGACAGTATCTACTACAAAATATTGTGCATTAGATGCCTCATATACTTCTCCAGCACGTTTGATTATCTCCCAAGTACATAACATACCTGTTTTAAt                                       >  2:85734/1‑141 (MQ=255)
                                                                                        tATTTCCAGGATTAAGCACTTCATCTTCGCCTTTACAAACCATGACAGGATCTACTACAAAATATTGTGCATTAGATGCCTCATATACTTCTCCAGCACGTTTGATTATCTCCTCAGTACCTAACATACCTGTTTTAATAg                                     <  2:241523/141‑1 (MQ=255)
                                                                                                    tAAGCACTTCATCTTCGCCTTTACAAACCATGACAGGATCTACTACAAAATATTGTGCATTAGATGCCTCATATACTTCTCCAGCACGTTTGATTATCTCCTCAGTACCTAACATACCTGTTTTAATAGCATCAGGTCCGa                         <  2:349212/141‑1 (MQ=255)
                                                                                                     aaGCACTTCATCTTCGCCTTTACAAACCATGACAGGATCTACTACAAAATATTGTGCATTAGATGCCTCATATACTTCTCCAGCACGTTTGATTATCTCCTCAGTACCTAACATACCTGTTTTAATAGCATCAGGTCCGAt                        >  2:48826/1‑141 (MQ=255)
                                                                                                       gCACTTCATCTTCGCCTTTACAAACCATGACAGGATCTACTACAAAATATTGTGCATTAGATGCCTCATATACTTCTCCAGCACGTTTGATTATCTCCGCAGTACCTAACATACCTGTTTTAATGGCATCAGGTCCGATTg                      >  1:313903/1‑141 (MQ=255)
                                                                                                            tCATCTTCGCCTTTACAAACCATGACAGGATCTACTACAAAATATTGTGCATTAGATGCCTCATATACTTCTCCAGCACGTTTGATTATCTCCTCAGTGCCTAACATACCTGTTTTAATAGCATCAGGTCCGATTGATaaa                 >  1:275714/1‑141 (MQ=255)
                                                                                                               tCTTCGCCTTTACAAACCCTGCGAGGATCTACTACAAAATATTGTGCATTAGATGCCTCATATACTTCTCCAGCACGTTTGATTATCTCCGCACTACCTAACATACTTGTTTTAACAGCATCACACCCGTTTGATCAATcc              >  2:198729/1‑141 (MQ=255)
                                                                                                                cTTCGCCTTTACAAACCATGACAGGATCTACTACAAAATATTGTGCATTAGATGCCTCATATACTTCTCCAGCACGTTTGATTATCTCCTCACTACCTAACATACCTGTTTTACTAGCATCAGGCCCGATTGATAAACCCg             >  2:11416/1‑141 (MQ=255)
                                                                                                                 ttCGCCTTTACAAACCATGACAGGATCTACTACAAAATATTGTGCATTAGATGCCTCATATACTTCTCCAGCACGTTTGATTATCTCCTCAGTACATAACATACCTGTTTTAATAGCATAAGTTCCTACTTATAAAGccca            >  2:256106/1‑139 (MQ=255)
                                                                                                                  tCGCCTTTACAAACCATGACAGGATCTACTACAAAATATTGTGCATTAGATGCCTCATATACTTCTCCAGCACGTTTGATTATCTCCTCAGTACCTAACATACCTGTTTTAATAGCATCAGGTCCGATTGATAAAGCCGtt           <  2:80667/141‑1 (MQ=255)
                                                                                                                    gCCTTTACAAACCATGACAGGATCTACTACAAAATATTGTGCATTAGATGCCTCATATACTTCTCCAGCACGTTTGATTATCTCCTCAGTACCTAACATACCTGTTTTAATAGCATCAGGTCCGATTGATAAAGCCGTTTc         <  2:144115/141‑1 (MQ=255)
                                                                                                                     ccTTTACAAACCATGACAGGATCTACTACAAAATATTGTGCATTAGATGCCTCATATACTTCTCCAGCACGTTTGATTATCTCCTCAGTACCTAACATACCTGTTTTAATAGCATCAGGTCCGATTGATAAAGCCGTTTCa        <  1:175700/141‑1 (MQ=255)
                                                                                                                           cAAGCCATGGCAGGTTCTACTTCAAAATATTGTGCATTAGATGCTTCTTATACTTCTCCTGTTCGTTTGATTATCTCCTCAGTACCTAACATACCTGTTTTAATAGCTTCTGGTCCAATTGATAAAGCCGTTTCAAgttgt  <  2:302101/141‑1 (MQ=255)
                                                                                                                           cAAACCATGGCAGGATCTACTACAAAATATTGTGCATTAGATGCCTCATATACTTCTCCAGCACGTTTGATTATCTCCTCAGTACCTAACATACCTGTTTTAATAGCATCAGGTCCGATTGATAAAGCCGTTTCAAgttgt  <  2:301594/141‑1 (MQ=255)
                                                                                                                                    |                                                                                                                                   
TTACCTAATCCTGATAATTGACCCGCTTCAAATAAGTTTGGTGTCACTACCGTTGCTTTAGGTAGTAAATATTTAATCATCGCCTCAGTATTTCCAGGATTAAGCACTTCATCTTCGCCTTTACAAACCATGACAGGATCTACTACAAAATATTGTGCATTAGATGCCTCATATACTTCTCCAGCACGTTTGATTATCTCCTCAGTACCTAACATACCTGTTTTAATAGCATCAGGTCCGATTGATAAAGCCGTTTCAAGTTGT  >  USA300TCH1516_ALE/638203‑638466

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