Predicted mutation
evidence seq id position mutation annotation gene description
RA USA300TCH1516_ALE 2,373,698 +T coding (78/1335 nt) pbuG ← Guanine/hypoxanthine permease PbuG

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE2,373,6971.T97.3% 132.6 / NA 37S27T (TCT→ACT) pbuGGuanine/hypoxanthine permease PbuG
Reads supporting (aligned to +/- strand):  ref base . (0/0);  major base T (17/19);  minor base G (0/1);  total (17/20)
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.00e+00

GCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCATTTTCATATCTTCTGATACGCCTTTAACACCTGCTAAACTTAAAACTTGCGGGTTAACTGCTAAAATATAGGCCATAGA‑TAAGAAAGTTGTGATACCGCCTAAGATTTCTCTTTTAAAGTTTGTTCCATATTTATCGAACTGGAAATAGTTTTTCACGAATGGATGCTCCTTTTTTTAAAATACCTTAATAGTTTAATATCATTCT  >  USA300TCH1516_ALE/2373565‑2373824
                                                                                                                                     |                                                                                                                               
gCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCATTTTCATATCTTCTGATACGCCTTTAACACCTGCTAAACTTAAAACTTGCGGGTTAACTGCTAAAATATAGGCCATAGATTAAGaaa                                                                                                                          >  1:112560/1‑141 (MQ=255)
gCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCATTTTCATATCTTCTGATACGCCTTTAACACCTGCTAAACTTAAAACTTGCGGGTTAACTGCTAAAATATAGGCCATAGATTAAGaaa                                                                                                                          >  1:357909/1‑141 (MQ=255)
              gctgctAATGCAGTCGCTACAAAAATGGCACCTTGGTCCATTTTCATATCTTCTGATACGCCTTTAACACCTGCTAAACTTAAAACTTGCGGGTTAACTGCTAAAATATAGGCCATAGATTAAGAAAGTTGTGATACCGcc                                                                                                            >  1:50342/1‑141 (MQ=255)
                 gctAATGCAGTCGCTACAAAAATGGCACCTTGGTCCATTTTCATATCTTCAGATACGCCTTTAACACCTCCTAAACTTAAAACTTGCGGGTTAACTGCTAAAATATAGGCCATAGATTAAGAAAGTTGTGATACCGCCTaa                                                                                                         >  2:279760/1‑141 (MQ=255)
                                       tGGCACCTTGGTCCATTTTCATATCTTCTGATACGCCTTTAACACCTGCTAAACTTAAAACTTGCGGGTTAACTGCTAAAATATAGGCCATAGATTAAGAAAGTTGTGATACCGCCTAAGATTTCTCTTTTAAAGTTTGtt                                                                                   >  1:176898/1‑141 (MQ=255)
                                       tGGCACCTTGGTCCATTTTCATATCTTCTGATACGCCTTTAACACCTGCTAAACTTAAAACTTGCGGGTTAACTGCTAAAATATAGGCCATAGATTAAGAAAGTTGTGATACCGCCTAAGATTTCTCTTTTAAAGTTTGtt                                                                                   <  2:111902/141‑1 (MQ=255)
                                       tGGCACCTTGGTCCATTTTCATATCTTCTGATACGCCTTTAACACCTGCTAAACTTAAAACTTGCGGGTTAACTGCTAAAATATAGGCCATAGATTAAGAAAGTTGTGATACCGCCTAAGATTTCTCTTTTAAAGTTTGtt                                                                                   <  2:373329/141‑1 (MQ=255)
                                       tGGCACCTTGGTCCATTTTCATATCTTCTGATACGCCTTTAACACCTGCTAAACTTAAAACTTGCGGGTTAACTGCTAAAATATAGGCCATAGATTAAGAAAGTTGTGATACCGCCTAAGATTTCTCTTTTAAAGTTTGtt                                                                                   >  1:69770/1‑141 (MQ=255)
                                       tGGCACCTTGGTCCATTTTCATATCTTCTGATACGCCTTTAACACCTGCTAAAATTAAAACTTGCGGGTTAACTGCTAAAATATAGGCCATAGATTAAGAAAGTTGTGATACCGCCTAAGATTTCTCTTTTAAAGTTTGtt                                                                                   >  1:227881/1‑141 (MQ=255)
                                            ccTTGGTCCATTTTCATATCTTCTGATACGCCTTTAACACCTGCTAAACTTAAAACTTGCGGGTTAACTGCTAAAATATAGGCCATAGATTAAGAAAGTTGTGATACCGCCTAAGATTTCTCTTTTAAAGTTTGTTCCata                                                                              >  2:398950/1‑141 (MQ=255)
                                            ccTTGGTCCATTTTCATATCTTCTGATACGCCTTTAACACCTGCTAAACTTAAAACTTGCGGGTTAACTGCTAAAATATAGGCCATAGATTAAGAAAGTTGTGATACCGCCTAAGATTTCTCTTTTAAAGTTTGTTCCata                                                                              >  2:232393/1‑141 (MQ=255)
                                               tGGTCCATTTTCATATACGTTGATACGCCTTTAACACCTGCTAAACTTAAAACTTGCGGGTTAACTGCTAAAATATAGGCCATAGATTAAGAAAGTTGTGATACCGCCTAAGATTTCTCTTTTAAAGTTTGTTCCATAttt                                                                           <  1:3695/141‑1 (MQ=255)
                                                ggTCCATTTTCATATCTTCTGATACGCCTTTAACACCTGCTAAACTTAAAACTTGCGGGTTAACTGCTAAAATATAGGCCATAGATTAAGAAAGTTGTGATACCGCCTAAGATTTCTCTTTTAAAGTTTGTTCCATATTTa                                                                          <  2:157737/141‑1 (MQ=255)
                                                ggTCCATTTTCATATCTTCTGATACGCCTTTAACACCTGCTAAACTTAAAACTTGCGGGTTAACTGCTAAAATATAGGCCATAGATTAAGAAAGTTGTGATACCGCCTAAGATTTCTCTTTTAAAGTTTGTTCCATATTTa                                                                          <  2:294884/141‑1 (MQ=255)
                                                             atCTTCTGATACGCCTTTAACACCTGCTAAACTTAAAACTTGCGGGTTAACTGCTAAAATATAGGCCATAGATTAAGAAAGTTGTGATACCGCCTAAGATTTCTCTTTTAAAGTTTGTTCCATATTTATCGAACTGGAAAt                                                             <  1:277626/141‑1 (MQ=255)
                                                                          ccTTTAACACCTGCTAAACTTAAAACTTGCGGGTTAACTGCTAAAATATAGGCCATAGATTAAGAAAGTTGTGATACCGCCTAAGATTTCTCTTTTAAAGTTTGTTCCATATTTATCGAACTGGAAATAGTTTTTCACGaa                                                >  2:364028/1‑141 (MQ=255)
                                                                          ccTTTAACACCTGCTAAACTTAAAACTTGCGGGTTAACTGCTAAAATATAGGCCATAGATTAAGAAAGTTGTGATACCGCCTAAGATTTCTCTTTTAAAGTTTGTTCCATATTTATCGAACTGGAAATAGTTTTTCACGaa                                                >  2:10457/1‑141 (MQ=255)
                                                                           cTTTAACACCTGCTAAACTTAAGACTTGCGGGTTAACTGCTAAAATATAGGCCATAGATTAAGAAAGTTGTGATACCGCCTAAGATTTCTCTTTTAAAGTTTGTTCCATATTTATCGAACTGGAAATAGTTTTTCACGAAt                                               <  2:72588/141‑1 (MQ=255)
                                                                           cTTTAACACCTGCTAAACTTAAAACTTGCGGGTTAACTGCTAAAATATAGGCCATAGATTAAGAAAGTTGTGATACCGCCTAAGATTTCTCTTTTAAAGTTTGTTCCATATTTATCGAACTGGAAATAGTTTTTCACGAAt                                               <  2:58917/141‑1 (MQ=255)
                                                                                acacCTGCTAAACTTAAAACTTGCGGGTTAACTGCTAAAATATAGGCCATAGATTAAGAAAGTTGTGATACCGCCTAAGATTTCTCTTTTAAAGTTTGTTCCATATTTATCGAACTGGAAATAGTTTTTCACGAATGGATg                                          <  2:210440/141‑1 (MQ=255)
                                                                                acacCTGCTAAACTTAAAACTTGCGGGTTAACTGCTAAAATATAGGCCATAGATTAAGAAAGTTGTGATACCGCCTAAGATTTCTCTTTTAAAGTTTGTTCCATATTTATCGAACTGGAAATAGTTTTTCACGAATGGATg                                          <  2:348131/141‑1 (MQ=255)
                                                                                 cacCTGCTAAACTTAAAACTTGTGGGTTAACTGCTAAAATATAGGCCATAGATTAAGAAAGTTGTGATACCGCCTAAGATTTCTCTTTTAAAGTTTGTTCCATATTTATCGAACTGGAAATAGTTTTTCACGAATGGATGc                                         >  1:357765/1‑141 (MQ=255)
                                                                                 cacCTGCTAAACTTAAAACTTGCGGGTTAACTGCTAAAATATAGGCCATAGATTAAGAAAGTTGTGATACCGCCTAAGATTTCTCTTTTAAAGTTTGTTCCATATTTATCGAACTGGAAATAGTTTTTCACGAATGGATGc                                         <  1:422503/141‑1 (MQ=255)
                                                                                      gCTAAACTTAAGACTTGCGGGTTAACTGCTAAAATATAGGCCATAGAGTAAGAAAGTTGTGATACCGCCTAAGATTTCTCTGTTAAAGTTTGTTCCATATTTATCGAACTGGAAATAGTTTTTCACGAATGGATGCTCCtt                                    <  2:273819/141‑1 (MQ=255)
                                                                                         aaaCTTAAAACTTGCGGGTTAACTGCTAAAATATAGGCCATAGATTAAGAAAGTTGTGATACCGCCTAAGATTTCTCTTTTAAAGTTTGTTCCATATTTATCGAACTGGAAATAGTTTTTCACGAATGGATGCTCCttttt                                 <  2:279818/141‑1 (MQ=255)
                                                                                         aaaCTTAAAACTTGCGGGTTAACTGCTAAAATATAGGCCATAGATTAAGAAAGTTGTGATACCGCCTAAGATTTCTCTTTTAAAGTTTGTTCCATATTTATCGAACTGGAAATAGTTTTTCACGAATGGATGCTCCttttt                                 <  2:165236/141‑1 (MQ=255)
                                                                                           aCTTAAAACTTGCGGGTTAACTGCTAAAATATAGGCCATAGATTAAGAAAGTTGTGATACCGCCTAAGATTTCTCTTTTAAAGTTTGTTCCATATTTATCGAACTGGAAATAGTTTTTCACGAATGGATGCTCCttttttt                               <  1:66688/141‑1 (MQ=255)
                                                                                              tAAAACTTGCGGGTTAACTGCTAAAATATAGGCCATAGATTAAGAAAGTTGTGATACCGCCTAAGATTTCTCTTTTAAAGTTTGTTCCATATTTATCGAACTGGAAATAGTTTTTCACGAATGGATGCTCCTTTTTTTaaa                            <  1:398950/141‑1 (MQ=255)
                                                                                              tAAAACTTGCGGGTTAACTGCTAAAATATAGGCCATAGATTAAGAAAGTTGTGATACCGCCTAAGATTTCTCTTTTAAAGTTTGTTCCATATTTATCGAACTGGAAATAGTTTTTCACGAATGGATGCTCCTTTTTTTaaa                            <  1:232393/141‑1 (MQ=255)
                                                                                                     tGCGGGTTAACTGCTAAAATATAGGCCATAGATTAAGAAAGTTGTGATACCGCCTAAGATTTCTCTTTTAAAGTTTGTTCCATATTTATCGAACTGGAAATAGTTTTTCACGAATGGATGCTCCTTTTTTTAAAATACCtt                     <  2:364006/141‑1 (MQ=255)
                                                                                                      gCGGGTTAACTGCTAAAATATAGGCCATAGATTAAGAAAGTTGTGATACCGCCTAAGATTTCTCTTTTAAAGTTTGTTCCATATTTATCGAACTGGAAATAGTTTTTCACGAATGGATGCTCCTTTTTTTAAAATACCTTa                    >  2:434330/1‑141 (MQ=255)
                                                                                                      gCGGGTTAACTGCTAAAATATAGGCCATAGATTAAGAAAGTTGTGATACCGCCTAAGATTTCTCTTTTAAAGTTTGTTCCATATTTATCGAACTGGAAATAGTTTTTCACGAATGGATGCTCCTTTTTTTAAAATACCTTa                    >  2:100283/1‑141 (MQ=255)
                                                                                                      gCGGGTTAACTGCTAAAATATAGGCCATAGATTAAGAAAGTTGTGATACCGCCTAAGATTTCTCTTTTAAAGTTTGTTCCATATTTATCGAACTGGAAATAGTTTTTCACGAATGGATGCTCCTTTTTTTAAAATACCTTa                    >  1:208875/1‑141 (MQ=255)
                                                                                                         ggTTAACTGCTAAAATATAGGCCATAGATTAAGAAAGTTGTGATACCGCCTAAGATTTCTCTTTTAAAGTTTGTTCCATATTTATCGAACTGGAAATAGTTTTTCACGAATGGATGCTCCTTTTTTTAAAATACCTTAATa                 <  2:44653/141‑1 (MQ=255)
                                                                                                         ggTTAACTGCTAAAATATAGGCCATAGATTAAGAAAGTTGTGATACCGCCTAAGATTTCTCTTTTAAAGTTTGTTCCATATTTATCGAACTGGAAATAGTTTTTCACGAATGGATGCTCCTTTTTTTAAAATACCTTAATa                 <  2:44695/141‑1 (MQ=255)
                                                                                                                        tataGGCCATAGATTAAGAAAGTTGTGATACCGCCTAAGATTTCTCTTTTAAAGTTTGTTCCATATTTATCGAACTGGAAATAGTTTTTCACGAATGGATGCTCCTTTTTTTAAAATACCTTAATAGTTTAATATCATTCt  >  2:402385/1‑141 (MQ=255)
                                                                                                                        tataGGCCATAGATTAAGAAAGTTGTGATACCGCCTAAGATTTCTCTTTTAAAGTTTGTTCCATATTTATCGAACTGGAAATAGTTTTTCACGAATGGATGCTCCTTTTTTTAAAATACCTTAATAGTTTAATATCATTCt  >  2:44377/1‑141 (MQ=255)
                                                                                                                                     |                                                                                                                               
GCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCATTTTCATATCTTCTGATACGCCTTTAACACCTGCTAAACTTAAAACTTGCGGGTTAACTGCTAAAATATAGGCCATAGA‑TAAGAAAGTTGTGATACCGCCTAAGATTTCTCTTTTAAAGTTTGTTCCATATTTATCGAACTGGAAATAGTTTTTCACGAATGGATGCTCCTTTTTTTAAAATACCTTAATAGTTTAATATCATTCT  >  USA300TCH1516_ALE/2373565‑2373824

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