Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,744,454 C→T 100% G73D (GGC→GAC)  apt ← Adenine phosphoribosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,744,4540CT98.4% 44.3 / ‑6.6 16G73D (GGC→GAC) aptAdenine phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (7/8);  minor base G (0/1);  total (7/9)
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
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGTTGCATAGCCAT  >  USA300TCH1516_ALE/1744317‑1744574
                                                                                                                                         |                                                                                                                        
catcTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACGGGTGCAAAGTCaa                                                                                                                       <  2:284232/141‑1 (MQ=255)
      tAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGTCAATCCCCa                                                                                                                 <  1:134712/141‑1 (MQ=255)
             cacacGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCTAAGTCATTCCCCATTGAATa                                                                                                          >  1:89553/1‑141 (MQ=255)
               cacGTTGCACTTGTAGAATTGCATCTTTGTGCATTGTTAAAACATTTATACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCATTCTTTTCTAACAGGTGCAAAGTCAATCCCCATTGAATAAg                                                                                                        >  2:335828/1‑141 (MQ=255)
                 cggtGATCTGGTTTAAGTGCATCTTTGTGCATTTTTAAAACATTTGTACCATATTCCAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGTCAATCCCCATTGAATAAGCt                                                                                                      <  2:190446/138‑1 (MQ=255)
                 cGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGTCAATCCCCATTGAATAAGCt                                                                                                      <  2:258181/141‑1 (MQ=255)
                     gACCTCGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATGACGAATGACTTCACGAGGTACTTTCCTGTCTTTGCTAACAGGTGCAAAGTCAATCCCCATTGAATAAGCTACAg                                                                                                  <  1:66551/141‑1 (MQ=255)
                        cTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGTCAATCCCCATTGAATAAGCTACAGGAc                                                                                               >  1:182702/1‑141 (MQ=255)
                        cTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTGTCTCACAGGTGCAAAGTCAATCCCCAGTGGATATGCTACAGGAc                                                                                               >  1:34575/1‑141 (MQ=255)
                         tGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGTCAATCCCCATTGAATAAGCTACAGGACa                                                                                              <  1:228805/141‑1 (MQ=255)
                                         ttGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGTCAATCCCCATTGAATAAGCTACAGGACAGCCAATGGTAAAGCCa                                                                              >  1:263096/1‑141 (MQ=255)
                                              cATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGTCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCAagcat                                                                         >  2:256634/1‑139 (MQ=255)
                                                            ttGTACCATATTCTAGGTCATACTCATAACGAGTGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGTCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACaa                                                           <  1:197461/141‑1 (MQ=255)
                                                                                               aCTTCACGAGTGAATTTCCCTTCTCTTCTAACAGGTGCAAAGTCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCATCATCTCTGTCTTTTGCGTATTCTACa                        <  2:206930/141‑1 (MQ=255)
                                                                                                              ttCCCTTCTTTTCTAACAGGTGCAAAGGCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGTTGCa         <  1:3945/141‑1 (MQ=255)
                                                                                                                     cttTTCTAACAGGTGCAAAGTCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGTTGCATAGCCat  >  1:17803/1‑141 (MQ=255)
                                                                                                                                         |                                                                                                                        
CATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGTTGCATAGCCAT  >  USA300TCH1516_ALE/1744317‑1744574

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