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,4540CT99.0% 55.2 / ‑6.7 23G73D (GGC→GAC) aptAdenine phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (8/14);  minor base G (1/0);  total (9/14)
Fisher's exact test for biased strand distribution p-value = 3.91e-01
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.

AATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCAT  >  USA300TCH1516_ALE/1744324‑1744586
                                                                                                                                  |                                                                                                                                    
aaTTAACACACGTTAACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTGCCCTTCTTTTCTAACAGGTGCAAAGTCAATCCCCa                                                                                                                             >  1:161285/1‑140 (MQ=255)
      cacacGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGTCAATCCCCATTGAAt                                                                                                                       <  2:65426/140‑1 (MQ=255)
         acGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACAATTGTACCATACTCTAGGTCATACTCATAACGAATGACTTCACGCGGTAATTTCCCTTCTTTTCTAAGAGGCACAAAGTCAATCCCCATTGAATaaa                                                                                                                    >  2:93440/1‑139 (MQ=255)
                    gTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGTCAATCCCCATTGAATAAGCTACAGGACAg                                                                                                         <  2:215197/140‑1 (MQ=255)
                                    gtgCTTTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGTCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCAc                                                                                         <  2:211178/140‑1 (MQ=255)
                                     tgCATTGTTAAAACATTTGTACCATAATCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGTCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCAcg                                                                                        >  2:174334/1‑140 (MQ=255)
                                              aaaaCATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGTCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAgg                                                                               <  1:279062/140‑1 (MQ=255)
                                              aaaaCATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGTCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAgg                                                                               <  2:193806/140‑1 (MQ=255)
                                                             tatTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGTCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATc                                                                <  1:66374/140‑1 (MQ=255)
                                                                tCTAGGTCATTCTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGTCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAAc                                                             <  1:263628/140‑1 (MQ=255)
                                                                   aGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGTCAATCCCCATTGAATAAGCTACAGGCCAGCCAATGATAAAGCCACGCGCGTCAGGTCCTACAACGCTCTCAACAtc                                                          >  2:261149/1‑140 (MQ=255)
                                                                     gTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGTCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACAtctc                                                        <  1:139249/140‑1 (MQ=255)
                                                                     gTCATACTCATAACGAATGACTTCACGAGGTAATTTACCTTCTTTTCTACCAGGTGCAAAGTCAATCTCCATTGAATAAGCTACAGGACAGCAAATAATAAAGCAACGCGATTCAGGGCTTATAACGATATCAACAtttt                                                        >  2:131293/1‑137 (MQ=255)
                                                                      tCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGTCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATctct                                                       <  1:83809/140‑1 (MQ=255)
                                                                          aCTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGTCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCt                                                   <  1:174334/140‑1 (MQ=255)
                                                                             cATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGTCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACAtct                                                         <  1:250856/131‑1 (MQ=255)
                                                                             cATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGTCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACAtct                                                         >  2:250856/1‑131 (MQ=255)
                                                                               tAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGTCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTCCAACGATATCAACATCTCTGTCTTTTGc                                              >  2:289733/1‑140 (MQ=255)
                                                                                   gAATGACGTCACTAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGTCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTAt                                          <  2:201087/140‑1 (MQ=255)
                                                                                      tGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGTCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCt                                       >  1:257050/1‑140 (MQ=255)
                                                                                          ttaaCGAGGTAATTTCTCTACTTTTCTAAAAGGTGCAAAGGCAATGCCCATTGAATAAGCTACAGGACAGCCAATGGTAAAGCCACGCGCTTCAAAACATAAAACAATATCAACATCTATGTCTTTTTCGTATTATACaa                                   >  2:166923/4‑140 (MQ=255)
                                                                                                                 ttCTAACAGGTGCAAAGTCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGTTGCATAGCCatat            <  2:161978/140‑1 (MQ=255)
                                                                                                                           tGCAAAGTCAATCCCCATGGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCAt  <  2:96084/140‑1 (MQ=255)
                                                                                                                                  |                                                                                                                                    
AATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCAT  >  USA300TCH1516_ALE/1744324‑1744586

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: