Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 198,135 T→A 36.7% S63R (AGT→AGA ssuC → Putative aliphatic sulfonates transport permease protein SsuC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE198,1350TA36.7% 27.2 / 19.0 27S63R (AGT→AGAssuCPutative aliphatic sulfonates transport permease protein SsuC
Reads supporting (aligned to +/- strand):  ref base T (9/8);  new base A (6/4);  total (15/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

AAATGGTCATTATTATTGGGCATTACCAACCTGTATTGTTACCGGGTCCTGCTCTTGTAGGAAAAAGTATATGGTCTTTCATTGTTACTGGAGAAATTTTCCAACATTTAGCAATTAGTTTATGGAGATTTGTAGCGGGCTTTGTTGTCGCATTGTTGGTTGCTATTCCATTGGGCTTCTTGCTTGGAAGGAATCGTTGGCTATACAACGCTATCGAACCGCTATTTCAATTGATTAGGCCGATATCTCCGAT  >  USA300TCH1516_ALE/198017‑198269
                                                                                                                      |                                                                                                                                        
aaaTGGTCATTATTATTGGGCATTACCAACCTGTATTGTTACCGGGTCCTGCTCTTGTAGGAAAAAGTATATGGTCTTTCATTGTTACTGGAGAAATTTTCCAACATTTAGCAATTAGTTTATGGAGATTTGTAGCGGGc                                                                                                                     >  2:333615/1‑140 (MQ=255)
 aaTGGGCATTATTATTGGGCATTACCAACCTGTATTGTTACCGGGTCCTGCTCTTGTAGGAAAAAGTATATGGTCTTTCATTGTTACTGGAGAAATTTTCCAACATTTAGCAATTAGTTTATGGAGATTTGTAGCGGGCt                                                                                                                    <  2:278793/140‑1 (MQ=255)
    ggTCATTATTATTGGGCATTACCAACCTGTATTGTTACCGGGTCCTGCTCTTGTAGGAAAAAGTATATGGTCTTTCATTGTTACTGGAGAAATTTTCCAACATTTAGCAATTAGTTTATGGAGATTTGTAGCGGGCTttg                                                                                                                 <  1:76327/140‑1 (MQ=255)
              attGGGCATTACCAACCTGTATTGTTACCGGGTCCTGCTCTTGTAGGAAAAAGTATATGGTCTTTCATTGTTACTGGAGAAATTTTCCAACATTTAGCAATTAGATTAGGTAGATTTTTAGCGATCTTTGTTGTCGCAtt                                                                                                       >  1:148085/1‑140 (MQ=255)
                    cATTACCAACCTGTATTGTTACCGGGTCCTGCTCTTGTAGGAAAAAGTATATGGTCTTTCATTGTTACTGGAGAAATTTTCCAACATTTAGCAATTAGTTTCTGGAGAGTTGTAGCGGGCCTTGTTGTCGCATTgttggt                                                                                                >  2:330636/1‑140 (MQ=255)
                       tACCAACCTGTATTGTTACCGGGTCCTGCTCTTGTAGGAAAAAGTATATGGTCTTTCATTGTTACTGGAGAAATTTTCCAACATTTAGCAATTAGATTATGGAGATTTGTAGCGGGCTTTGTTGTCGCATTGTTGGTTGc                                                                                             <  1:240469/140‑1 (MQ=255)
                          cAACCTGTATTGTTACCGGGTCCTGCTCTTGTAGGAAAAAGTATATGGTCTTTCATTGTTACTGGAGAAATTTTCCAACATTTAGCAATTAGTTTATGGAGATTTGTAGCGGGCTTTGTTGTCGCATTGTTGGTTGCTAt                                                                                          >  1:63464/1‑140 (MQ=255)
                                      ttAACGGGTCGTGGGCTTGGAGGATAATGGAGATGGTCTTGCATTGTTTGTGGAGGAGTTTTCGAACATTTAGCAATTAGATTATGGAGATTTGTAGCGGGCTTTGTTGTAGCATTGTTTGTTGCTATTCCATTGGGttt                                                                              >  2:208450/1‑140 (MQ=255)
                                         ccGGGTCCTGCTCTTGTAGGAAAAAGTATATGGTCTTTCATTGTTACTGGAGAAATTTTCCAACATTTAGCAATTAGATTATGGAGATTTGTAGCGGGCTTTGTTGTCGCATTGTTGGTTGCTATTCCATTGGGcttctt                                                                           >  1:15439/1‑140 (MQ=255)
                                            ggTCCTGCTCTTGTAGTAAAAAGTATATGGTCTTTCATTGTTACTGGAGAAATTTTCCAACATTTAGCAATTAGATTATGGAGATTTGTAGCGGGCTTTGTTGTCGCATTGTTGGTTGCTATTCCATTGGGCTTcttgct                                                                        <  1:92439/140‑1 (MQ=255)
                                                     ctTGTAGGAAAAAGTATATGGTCCTTCATTGTTACTGGAGAAATTTTCCAACATTTAGCAATTAGTTTATGGAGATTTGTAGCGGGCTTTGTTGTCGCATTGTTGGTTGCTATTCCATTGGGCTTCTTGCTTggaaggaa                                                               <  1:260343/140‑1 (MQ=255)
                                                                aaGTATATGGTCTTTCATTGTTACTGCAGAAATTTTCCAACATTTAGCAATTAGTTTATGGAGATTTGTAGCGGGCTTTGTTGTCGCATTGTTGGTTGCTATTCCATTGGGCTTCTTGCTTGGAAGGAATCGTTGGCtat                                                   <  1:96539/140‑1 (MQ=255)
                                                                      agggTTTTTCATTGTTACTGGAGAAATTTTCCAACATTTAGCAATTAGATTATGGAGATTTGTAGCGGGCTTTGTTGTCGCATTGTTGGTTGCTATTCCATTGGGCTTCTTGCTTGGAAGGAATCGTTGGCTATACAACg                                             <  1:66489/138‑1 (MQ=255)
                                                                        tgTCTTTCATTGTTACTGGAGAAATTTGCCAACATTTAGCAATTAGTTTATGGAGATTTGTAGCGGGCTTTGTTGTCGCATTGTTGGTTGCTATTCCATTGGGCTTCTTGCTTGGAAGGAATCGTTGGCTATACAACGCt                                           <  1:131807/139‑1 (MQ=255)
                                                                        ggTCTTTCATTGTTACTGGAGAAATTTTCCAACATTTAGCAATTAGTTTATGGAGATTTGTAGCGGGCTTTGTTGTCGCATTGTTGGTTGCTATTCCATTGGGCTTCTTGCTTGGAAGGAATCGTTTGCTATACAACGCt                                           >  1:214947/1‑140 (MQ=255)
                                                                                  tGTTACTGGAGAAATTTTCCAACATTTAGCAATTAGTTTATGGAGATTTGTAGCGGGCTTTGTTGTCGCATTGTTGGTTGCTATTCCATTGGGCTTCTTGCTTGGAAGGAATCGTTGGCTATACAACGCTATCGAACCGc                                 <  1:244269/140‑1 (MQ=255)
                                                                                   gTTACTGGAGAAATTTTCCAACATTTAGCAATTAGTTTATGGAGATTTGTAGCGGGCTTTGTTGTCGCATTGTTGGTTGCTATTCCATTGGGCTTCTTGCTTGGAAGGAATCGTTGGCTATACAACGCTATCGAACCGCt                                >  1:203307/1‑140 (MQ=255)
                                                                                   gTTACTGGAGAAATTTTCCAACATTTAGCAATTAGTTTATGGAGATTTGTAGCGGGCTTTGTTGTCGCATTGTTGGTTGCTATTCCATTGGGCTTCTTGCTTGGAAGGAATCGTTGGCTATACAACGCTATCGAACCGCt                                <  1:247907/140‑1 (MQ=255)
                                                                                      aCTGGAGAAATTTTCCAACATTTAGCAATTAGATTATGGAGATTTGTAGCGGGCTTTGTTGTCGCATTGTTGGTTGCTATTCCATTGGGCTTCTTGCTTGGAAGGAATCGTTGGCTATACAACGCTATCGAACCGCTAtt                             >  1:125923/1‑140 (MQ=255)
                                                                                            gaAATTTTCCAACATTTAGCAATTAGTTTACGGAGATTTGTAGCGGGCTTTGTCGTCGCATTGTTGGTTGCTATTCCATTGGGCTTCTTGCATGGAAGGAATCGTTGGCTATACAACGCTATCGAACCGCTATTTCAAtt                       >  2:217023/1‑140 (MQ=255)
                                                                                              aaTTTTCCAACATTTAGCAATTAGTTTATGGAGGTTTGTAGCGGGCTTTGTTGTAGCATTGTTGGTTGCTATTCCAGTGGGCTTCTTGTATGGAAGGAATCGTTGGCTATACAACGCTATCGAACCGCTATTTAAATTGa                     >  2:236076/1‑140 (MQ=255)
                                                                                                   tCCAACATTTAGCAATTAGATTATGGAGATTTGTAGCGGGCTTTGTTGTCGCATTGTTGGTTGCTATTCCATTGGGCTTCTTGCTTGGAAGGAATCGTTGGCTATACAACGCTATCGAACCGCTATTTCAATTGATTAgg                >  2:275058/1‑140 (MQ=255)
                                                                                                       aCATTTAGCAATTAGATTATGGAGATTTGTAGCGGGCTTTGTTGTCGCATTGTTGGTTGCTATTCCATTGGGCTTCTTGCTTGGAAGGAATCGTTGGCTATACAACGCTATCGAACCGCTATTTCAATTGATTAGGCCGa            <  1:278672/140‑1 (MQ=255)
                                                                                                           ttAGCAATTAGTTTATGGAGATTTGTAGCGGGCTTTGTTGTCGCATTGTTGGTTGCTATTCCATTGGGCTTCTTGCTTGGAAGGAATCGTTGGCTATACAACGCTATCGAACCGCTATTTCAATTGATTAGGCCGATAtc        <  1:333615/140‑1 (MQ=255)
                                                                                                             aGCAATTAGTTTATGGAGATTTGTAGCGGGCTTTGTTGTCGCATTGTTGGTCGATATTCCATTGGGCTTCTTGCTCGGCAGGAAAGTGGGCCGTTACAACTCTATCGAACCGCTATTTCATTTGAGTAGGCCGATAgctc      >  2:161961/1‑140 (MQ=255)
                                                                                                              gCAATTAGTTTATGGAGATTTGTGCCGGCCTTTGTTTCGCCTTGTTTGTTTGCTTTTCCTTGGGGCTCTTGGCTGGGAAGGACCCGTTGGCTTACCAACGCAACCGACCGGCTCTTTCAATGGTTAGGGCGAGTTCTCCg    >  2:133980/1‑140 (MQ=255)
                                                                                                                 aTTAGATTATGGAGATTTGTAGCGGGCTTTGTTGTCGCATTGTTGGTTGCTATTCCATTGGGCTTCTTGCTTGGAAGGAATCGTTGGCTATACAACGCTATCGAACCGCTATTTCAATTGATTAGGCCGATATCTCCGAt  >  1:305565/1‑140 (MQ=255)
                                                                                                                      |                                                                                                                                        
AAATGGTCATTATTATTGGGCATTACCAACCTGTATTGTTACCGGGTCCTGCTCTTGTAGGAAAAAGTATATGGTCTTTCATTGTTACTGGAGAAATTTTCCAACATTTAGCAATTAGTTTATGGAGATTTGTAGCGGGCTTTGTTGTCGCATTGTTGGTTGCTATTCCATTGGGCTTCTTGCTTGGAAGGAATCGTTGGCTATACAACGCTATCGAACCGCTATTTCAATTGATTAGGCCGATATCTCCGAT  >  USA300TCH1516_ALE/198017‑198269

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