Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,646,189 G→C 100% A155G (GCA→GGA)  ypdF ← Aminopeptidase YpdF

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,646,1890GC99.5% 70.3 / ‑7.0 26A155G (GCA→GGA) ypdFAminopeptidase YpdF
Reads supporting (aligned to +/- strand):  ref base G (0/1);  major base C (13/11);  minor base T (1/0);  total (14/12)
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.

TCACTTGCAACACCATGTGGTAATGCACCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTGCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTTATCTACAGTATTAGAAATGCTTATTAATGATATAC  >  USA300TCH1516_ALE/1646057‑1646324
                                                                                                                                    |                                                                                                                                       
tCACTTGCAACACCATGTGGTAATGCACCTCTATGACCAGATGCTACAATGGTATGACATGATGGGCCAGCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTTATTCTTTGTCAGTACTGCCTTCTTTTaaa                                                                                                                                 >  1:132228/1‑139 (MQ=255)
tCACTTGCAACACCATGTGGTAATGCACCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTaca                                                                                                                                 <  2:15361/141‑1 (MQ=255)
        aaCACCATGTGGTAATGCACCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATTCCTCCTTTTACAACCTTTaa                                                                                                                         >  2:100142/1‑141 (MQ=255)
                      aTGCACCGCTATGTCCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATg                                                                                                           <  2:151201/141‑1 (MQ=255)
                        gCACCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCCGATCCTAATGCTCGCATTTTTCTGTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGtt                                                                                                         >  1:147630/1‑141 (MQ=255)
                                          gCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAgct                                                                                       <  2:221254/141‑1 (MQ=255)
                                                      tCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAAtt                                                                           >  2:48344/1‑141 (MQ=255)
                                                        gAATGATGGTCCATCTGCTGCTAATTCTAGCATGTTGCTTTCTAATTTTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTaa                                                                         <  2:27367/141‑1 (MQ=255)
                                                        gAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACTGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTaa                                                                         <  2:46706/141‑1 (MQ=255)
                                                              tgGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTTTCTAATGCCCCCTTTTCCAGGTTCAAAAGGAAAGGCATATGTTTCATTAACACTATTATCTTATTTTTGAATTAAAACAAt                                                                   >  2:74124/1‑141 (MQ=255)
                                                                gTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGATACTTTTTTAATTAAAGCAAttt                                                                 >  2:102066/1‑141 (MQ=255)
                                                                  ccATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTTAATTAAATCAATTTCg                                                               >  2:77604/1‑141 (MQ=255)
                                                                          tCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATc                                                       <  1:260377/141‑1 (MQ=255)
                                                                            cTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCtt                                                     <  1:100142/141‑1 (MQ=255)
                                                                                    aGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGtct                                             <  2:146183/141‑1 (MQ=255)
                                                                                    aGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGtct                                             <  1:121108/141‑1 (MQ=255)
                                                                                     tctttttGCTTTCTAATAGTTGCTTTAATTCTTTTTCAGTCATGCCTCCTTTTTCAACAGTTAGTATATTTTCTTATGGTTCAAACTCTCGACATCCGGCTTTTTGAATTAAAGCAATTTCGTGTGCATCTTTGACGtctc                                            <  2:333768/138‑1 (MQ=255)
                                                                                                 cTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTTAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACTTCTCTAATTTTATCTa                                >  2:164668/1‑141 (MQ=255)
                                                                                                       ttGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTTATCTACAGTAt                          <  2:27559/141‑1 (MQ=255)
                                                                                                             ttAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTTATCTACAGTATTAGaaa                    >  1:221196/1‑141 (MQ=255)
                                                                                                               aaTTCTTTTTCAGTCATGCTTGCTTTTACAACAGTTAAAATATATTCCTATGTCTCGTCAAGTATATGAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGAGGTCTCTAATTTTATCTACAGTATTAGAAATg                  <  1:112931/141‑1 (MQ=255)
                                                                                                                     ttttCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTTATCTACAGTATTAGAAATGCttatt            >  2:76713/1‑141 (MQ=255)
                                                                                                                        tCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTTATTTACAGTATTAGAAATGCTTATTAAt         >  1:73707/1‑141 (MQ=255)
                                                                                                                        tCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTTATCTACAGTATTAGAAATTCTTATTAAt         >  2:22337/1‑141 (MQ=255)
                                                                                                                             cATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACCATATTAGCTGCTTTTTGAATTAAAGCACTTTCGTCAGCATCTTTGACGTCTCTAATTTCATCTACAGTATTAGAACTGCTTATTAATCatat    >  2:155095/1‑141 (MQ=255)
                                                                                                                               tGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTGATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTTATCTACAGTATTAGAAGTGCTTATTAGTGATCTAc  >  1:277863/1‑141 (MQ=255)
                                                                                                                                    |                                                                                                                                       
TCACTTGCAACACCATGTGGTAATGCACCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTGCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTTATCTACAGTATTAGAAATGCTTATTAATGATATAC  >  USA300TCH1516_ALE/1646057‑1646324

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: