Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,744,428 T→C 37.0% K82E (AAA→GAA)  apt ← Adenine phosphoribosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,744,4280TC37.0% 6.3 / 15.3 19K82E (AAA→GAA) aptAdenine phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base T (6/6);  new base C (4/3);  total (10/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.49e-01

CAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGT  >  USA300TCH1516_ALE/1744302‑1744563
                                                                                                                              |                                                                                                                                       
caGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGACGTCATTTCCCTTCTTTACtca                                                                                                                           >  2:33470/1‑139 (MQ=255)
       taataaATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTCCCCTTCTTTTCTAACAGGTGc                                                                                                                    <  1:9134/141‑1 (MQ=255)
       taataaATCATCTGTAATTAACACACGTTGACCGGGTGTAAGTGCATGTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTCCCCTTCTTTTCTAACAGGTGc                                                                                                                    <  1:232294/141‑1 (MQ=255)
                      aaTTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTTCAAAGCCAATCCCCAt                                                                                                     >  1:184095/1‑141 (MQ=255)
                          aaCACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTCCCCTTCTTTTCTAACAGGTGCAAAGCCAATACCCATTGaa                                                                                                 >  1:66407/1‑141 (MQ=255)
                                      ccTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTCCAGGa                                                                                     >  1:9753/1‑141 (MQ=255)
                                                   cATCTTTGTGCATTGTTAAAACATTTGTATCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCACAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATaa                                                                        >  2:182492/1‑141 (MQ=255)
                                                                       aCATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATACAGCCACGCGCTTCAGGTCCt                                                    >  1:87641/1‑141 (MQ=255)
                                                                       aCATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCt                                                    <  2:186782/141‑1 (MQ=255)
                                                                                       cTAGGTCATACTCATAACGAATGACTTCACGAGGTACTTCCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACACGACAGCCAACGATAAAGCCACGCCCTTCACATCCTACAACAATACCAAcct                                    >  2:58590/1‑139 (MQ=255)
                                                                                                 cTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCttt                          <  1:147417/141‑1 (MQ=255)
                                                                                                   cATAACGAATGACTTCACGAGGTAATTTCCCTACTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGTCAGCCAATGATAAAGCCACGCGCTGCAGGTCCTACAACGATATCAACATCTCTGTCTTTTg                        <  1:30918/141‑1 (MQ=255)
                                                                                                     tAACGAATGACTTCACGAGGTAATTCCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCg                      >  1:242204/1‑141 (MQ=255)
                                                                                                             taCTTCATGAGGTAATTCTCTTTCTTTTCGAACATGTTCAAAGCCAATCCCAATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTAc              <  2:154534/140‑1 (MQ=255)
                                                                                                                ttCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACAATATCAACATCTCTGTCTTTTGCGTCTTCTACAAt           >  2:144535/1‑141 (MQ=255)
                                                                                                                    cGAGGTAATTTACCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACCGGACAGACAATGATAAAGCCACGCGCTTCAGGTCCTACAACTATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTa       <  1:95647/141‑1 (MQ=255)
                                                                                                                         tAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGt  <  1:52293/141‑1 (MQ=255)
                                                                                                                         tAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGt  <  1:52286/141‑1 (MQ=255)
                                                                                                                         tAATTCCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGt  >  1:192930/1‑141 (MQ=255)
                                                                                                                              |                                                                                                                                       
CAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGT  >  USA300TCH1516_ALE/1744302‑1744563

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: