Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,744,472 G→A 100% A67V (GCT→GTT)  apt ← Adenine phosphoribosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,744,4720GA90.7% 80.5 / 0.4 32A67V (GCT→GTT) aptAdenine phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base G (2/1);  new base A (13/16);  total (15/17)
Fisher's exact test for biased strand distribution p-value = 5.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.17e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATCCATAATTGT  >  USA300TCH1516_ALE/1744338‑1744599
                                                                                                                                      |                                                                                                                               
caccggctttCATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAACTACa                                                                                                                            <  1:115360/133‑1 (MQ=255)
   cTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAACTACAGGa                                                                                                                         <  2:305666/140‑1 (MQ=255)
    tGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGTATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAACTACAGGAc                                                                                                                        <  2:66107/140‑1 (MQ=255)
         tAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCa                                                                                                                   <  2:239497/140‑1 (MQ=255)
         tAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAACTACAGGACAGCCa                                                                                                                   <  1:265025/140‑1 (MQ=255)
         tAATTGCATCTTTGTGCATTGTTAAAACACTAACACACTTTTCTAGGGCTTGCTCATAACGAATTACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAACTACAGGACAGCCa                                                                                                                   <  1:240243/140‑1 (MQ=255)
            ttGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTAATTAAACTCAAGGACAGCCAATg                                                                                                                >  1:557723/1‑140 (MQ=255)
                 tCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATcaa                                                                                                           >  2:332787/1‑140 (MQ=255)
                 tCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAACTAAAGTACAGCCAATGATaaa                                                                                                           <  2:70900/140‑1 (MQ=255)
                   tttGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAACTACAGGACAGCCCATGATAAAGc                                                                                                         >  1:361376/1‑140 (MQ=255)
                               tAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAACTACAGGACAGCCAATGATAAAGCCACGCGCTTCAg                                                                                             >  1:380033/1‑140 (MQ=255)
                                 aaaCATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAACTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGt                                                                                           >  1:435492/1‑140 (MQ=255)
                                    cATTTGTACCATATTGTAGGTCATACTCAGAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAACTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCt                                                                                        <  2:210536/140‑1 (MQ=255)
                                       ttGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAACTACAGGACAGCCAATGATAAAGCCCCGCGCTTCAGGTCCTACa                                                                                     >  1:491526/1‑140 (MQ=255)
                                             cATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAACTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGata                                                                               <  1:247188/140‑1 (MQ=255)
                                                   cTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAAAAGGTGCAAAGCCAATCCACATTGAATAAACTACACGACAGACAATGATAACGCCAAGCGCTTCAGGTCCTACAACGATACAAACa                                                                         >  2:509815/1‑140 (MQ=255)
                                                     aGGTCATACTCATTACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAACTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACAtc                                                                       <  2:301958/140‑1 (MQ=255)
                                                        tCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAACTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATctct                                                                    <  1:299886/140‑1 (MQ=255)
                                                            aCTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAACTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTcc                                                                >  2:148562/1‑139 (MQ=255)
                                                             cTCATACCGAATGTCTTGACGCGCTAATTACCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCCTTGAATAAACTACAGGCCTGCCAATGATAAAGGCACGCGCTTCAGGTCCCACAACGATATCAACATCTCTGTCtt                                                               <  2:182726/140‑1 (MQ=255)
                                                                      aaTGACTTCACGAGGTAATTTCCCTTCTTTTCTAGCAGGTGCAAAGCCAATCCCCATTGAATAAACTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTAtt                                                      <  1:118212/140‑1 (MQ=255)
                                                                       aTGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAACTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTc                                                     >  1:330526/1‑140 (MQ=255)
                                                                       aTGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAACTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTACTACAACGATATCAACATCTCTGTCTTTTGCGTATTc                                                     >  2:23847/1‑140 (MQ=255)
                                                                       aTGACTGCACGTGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGCCTTTTGCGTATTc                                                     >  2:469311/1‑140 (MQ=255)
                                                                                   ggTAATTTCCCTTCCTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAACTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATc                                         <  1:120367/140‑1 (MQ=255)
                                                                                                 tttgcaaaCGGGTGCAAAGCCAACCCCCGTGGAATAAACTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGTTGCATAGCCat                           <  2:51296/134‑1 (MQ=255)
                                                                                                   ttCTAACAGGTGCAAAGCCAATCTCCATTGAATAAACTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGTTGCATAGCCatat                         <  2:421377/140‑1 (MQ=255)
                                                                                                        aCAGGTGCAAAGCCAATCCCCATTGAATAAACTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGTTGCATAGCCATATGCTTc                    >  1:554568/1‑140 (MQ=255)
                                                                                                           ggTGCAAAGCCAATCCCCATTGAATAAACTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGTTGCATCGCCATATGCTTCAcc                 >  2:15594/1‑140 (MQ=255)
                                                                                                             tGCAAAGCCAATCCCCATTGAATAAACTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCAt               <  2:190058/140‑1 (MQ=255)
                                                                                                                        tCCCCATTGAATAAACTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTATACAATTTTATCTGTTGCATAGCCATATGCTTAACCATTTTCCTTGAtt    >  1:548628/1‑140 (MQ=255)
                                                                                                                          cccATTGAATAAAATACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTATGTATTTTGCGTATTCTACAATTTTACATGTTGAATCGACCTATGTTTCATCATCATCCATAATTgt  >  2:358808/1‑140 (MQ=255)
                                                                                                                                      |                                                                                                                               
GACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATCCATAATTGT  >  USA300TCH1516_ALE/1744338‑1744599

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: