Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 198,135 T→A 39.0% S63R (AGT→AGA)  USA300HOU_0187 → ABC superfamily ATP binding cassette transporter, membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP000730198,1350TA39.0% 22.4 / 23.2 28S63R (AGT→AGA) USA300HOU_0187ABC superfamily ATP binding cassette transporter, membrane protein
Reads supporting (aligned to +/- strand):  ref base T (9/8);  new base A (6/5);  total (15/13)
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

AAATGGTCATTATTATTGGGCATTACCAACCTGTATTGTTACCGGGTCCTGCTCTTGTAGGAAAAAGTATATGGTCTTTCATTGTTACTGGAGAAATTTTCCAACATTTAGCAATTAGTTTATGGAGATTTGTAGCGGGCTTTGTTGTCGCATTG‑TTGGTTGCTATTCCATTGGGCTTCTTGCTTGGAAGGAATCGTTGG‑CTATACAACGCTATCGAACCGCTATTTCAATTGATTAGGCCGATATCTCCGAT  >  CP000730/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)
                    cATTACCAACCTGTATTGTTACCGGGTCCTGCTCTTGTAGGAAAAAGTATATGGTCTTTCATTGTTACTGGAGAAATTTTCCAACATTTAGCAATTAGTTTCTGGAGAGTTGTAGCGGGCCTTGTTGTCGCATTg‑ttggt                                                                                                >  2:330636/1‑140 (MQ=255)
                       tACCAACCTGTATTGTTACCGGGTCCTGCTCTTGTAGGAAAAAGTATATGGTCTTTCATTGTTACTGGAGAAATTTTCCAACATTTAGCAATTAGATTATGGAGATTTGTAGCGGGCTTTGTTGTCGCATTG‑TTGGTTGc                                                                                             <  1:240469/140‑1 (MQ=255)
                        aagatccgctATTGTTAGCGTGTCCTGTTCTTGTAAGAAAAAGTATATATTCCTGCATTGTTACTGGAGAAATTTTCCAACATTTAGCAATTAGATTATGGAGATTTGTAGCGGGCTTTGTTGTCGCATTG‑TTGGTTGCt                                                                                            <  1:47476/131‑1 (MQ=255)
                          cAACCTGTATTGTTACCGGGTCCTGCTCTTGTAGGAAAAAGTATATGGTCTTTCATTGTTACTGGAGAAATTTTCCAACATTTAGCAATTAGTTTATGGAGATTTGTAGCGGGCTTTGTTGTCGCATTG‑TTGGTTGCTAt                                                                                          >  1:63464/1‑140 (MQ=255)
                                      ttAACGGGTCGTGGGCTTGGAGGATAATGGAGATGGTCTTGCATTGTTTGTGGAGGAGTTTTCGAACATTTAGCAATTAGATTATGGAGATTTGTAGCGGGCTTTGTTGTAGCATTG‑TTTGTTGCTATTCCATTGGGttt                                                                              >  2:208450/1‑140 (MQ=255)
                                         ccGGGTCCTGCTCTTGTAGGAAAAAGTATATGGTCTTTCATTGTTACTGGAGAAATTTTCCAACATTTAGCAATTAGATTATGGAGATTTGTAGCGGGCTTTGTTGTCGCATTG‑TTGGTTGCTATTCCATTGGGcttctt                                                                           >  1:15439/1‑140 (MQ=255)
                                            ggTCCTGCTCTTGTAGTAAAAAGTATATGGTCTTTCATTGTTACTGGAGAAATTTTCCAACATTTAGCAATTAGATTATGGAGATTTGTAGCGGGCTTTGTTGTCGCATTG‑TTGGTTGCTATTCCATTGGGCTTcttgct                                                                        <  1:92439/140‑1 (MQ=255)
                                                     ctTGTAGGAAAAAGTATATGGTCCTTCATTGTTACTGGAGAAATTTTCCAACATTTAGCAATTAGTTTATGGAGATTTGTAGCGGGCTTTGTTGTCGCATTG‑TTGGTTGCTATTCCATTGGGCTTCTTGCTTggaaggaa                                                               <  1:260343/140‑1 (MQ=255)
                                                                aaGTATATGGTCTTTCATTGTTACTGCAGAAATTTTCCAACATTTAGCAATTAGTTTATGGAGATTTGTAGCGGGCTTTGTTGTCGCATTG‑TTGGTTGCTATTCCATTGGGCTTCTTGCTTGGAAGGAATCGTTGG‑Ctat                                                   <  1:96539/140‑1 (MQ=255)
                                                                      agggTTTTTCATTGTTACTGGAGAAATTTTCCAACATTTAGCAATTAGATTATGGAGATTTGTAGCGGGCTTTGTTGTCGCATTG‑TTGGTTGCTATTCCATTGGGCTTCTTGCTTGGAAGGAATCGTTGG‑CTATACAACg                                             <  1:66489/138‑1 (MQ=255)
                                                                        tgTCTTTCATTGTTACTGGAGAAATTTGCCAACATTTAGCAATTAGTTTATGGAGATTTGTAGCGGGCTTTGTTGTCGCATTG‑TTGGTTGCTATTCCATTGGGCTTCTTGCTTGGAAGGAATCGTTGG‑CTATACAACGCt                                           <  1:131807/139‑1 (MQ=255)
                                                                        ggTCTTTCATTGTTACTGGAGAAATTTTCCAACATTTAGCAATTAGTTTATGGAGATTTGTAGCGGGCTTTGTTGTCGCATTG‑TTGGTTGCTATTCCATTGGGCTTCTTGCTTGGAAGGAATCGTTTG‑CTATACAACGCt                                           >  1:214947/1‑140 (MQ=255)
                                                                                  tGTTACTGGAGAAATTTTCCAACATTTAGCAATTAGTTTATGGAGATTTGTAGCGGGCTTTGTTGTCGCATTG‑TTGGTTGCTATTCCATTGGGCTTCTTGCTTGGAAGGAATCGTTGG‑CTATACAACGCTATCGAACCGc                                 <  1:244269/140‑1 (MQ=255)
                                                                                   gTTACTGGAGAAATTTTCCAACATTTAGCAATTAGTTTATGGAGATTTGTAGCGGGCTTTGTTGTCGCATTG‑TTGGTTGCTATTCCATTGGGCTTCTTGCTTGGAAGGAATCGTTGG‑CTATACAACGCTATCGAACCGCt                                >  1:203307/1‑140 (MQ=255)
                                                                                   gTTACTGGAGAAATTTTCCAACATTTAGCAATTAGTTTATGGAGATTTGTAGCGGGCTTTGTTGTCGCATTG‑TTGGTTGCTATTCCATTGGGCTTCTTGCTTGGAAGGAATCGTTGG‑CTATACAACGCTATCGAACCGCt                                <  1:247907/140‑1 (MQ=255)
                                                                                      aCTGGAGAAATTTTCCAACATTTAGCAATTAGATTATGGAGATTTGTAGCGGGCTTTGTTGTCGCATTG‑TTGGTTGCTATTCCATTGGGCTTCTTGCTTGGAAGGAATCGTTGG‑CTATACAACGCTATCGAACCGCTAtt                             >  1:125923/1‑140 (MQ=255)
                                                                                            gaAATTTTCCAACATTTAGCAATTAGTTTACGGAGATTTGTAGCGGGCTTTGTCGTCGCATTG‑TTGGTTGCTATTCCATTGGGCTTCTTGCATGGAAGGAATCGTTGG‑CTATACAACGCTATCGAACCGCTATTTCAAtt                       >  2:217023/1‑140 (MQ=255)
                                                                                              aaTTTTCCAACATTTAGCAATTAGTTTATGGAGGTTTGTAGCGGGCTTTGTTGTAGCATTG‑TTGGTTGCTATTCCAGTGGGCTTCTTGTATGGAAGGAATCGTTGG‑CTATACAACGCTATCGAACCGCTATTTAAATTGa                     >  2:236076/1‑140 (MQ=255)
                                                                                                   tCCAACATTTAGCAATTAGATTATGGAGATTTGTAGCGGGCTTTGTTGTCGCATTG‑TTGGTTGCTATTCCATTGGGCTTCTTGCTTGGAAGGAATCGTTGG‑CTATACAACGCTATCGAACCGCTATTTCAATTGATTAgg                >  2:275058/1‑140 (MQ=255)
                                                                                                       aCATTTAGCAATTAGATTATGGAGATTTGTAGCGGGCTTTGTTGTCGCATTG‑TTGGTTGCTATTCCATTGGGCTTCTTGCTTGGAAGGAATCGTTGG‑CTATACAACGCTATCGAACCGCTATTTCAATTGATTAGGCCGa            <  1:278672/140‑1 (MQ=255)
                                                                                                           ttAGCAATTAGTTTATGGAGATTTGTAGCGGGCTTTGTTGTCGCATTG‑TTGGTTGCTATTCCATTGGGCTTCTTGCTTGGAAGGAATCGTTGG‑CTATACAACGCTATCGAACCGCTATTTCAATTGATTAGGCCGATAtc        <  1:333615/140‑1 (MQ=255)
                                                                                                             aGCAATTAGTTTATGGAGATTTGTAGCGGGCTTTGTTGTCGCATTG‑TTGGTCGATATTCCATTGGGCTTCTTGCTCGGCAGGAA‑AGTGGGCCGTTACAACTCTATCGAACCGCTATTTCATTTGAGTAGGCCGATAgctc      >  2:161961/1‑140 (MQ=255)
                                                                                                              gCAATTAGTTTATGGAGATTTGTGCCGGCCTTTGTT‑TCGCCTTGTTTGTTTGCTTTTCCTTGGGGCTCTTGGCTGGGAAGGACCCGTTGG‑CTTACCAACGCAACCGACCGGCTCTTTCAA‑TGGTTAGGGCGAGTTCTCCg    >  2:133980/1‑140 (MQ=255)
                                                                                                                 aTTAGATTATGGAGATTTGTAGCGGGCTTTGTTGTCGCATTG‑TTGGTTGCTATTCCATTGGGCTTCTTGCTTGGAAGGAATCGTTGG‑CTATACAACGCTATCGAACCGCTATTTCAATTGATTAGGCCGATATCTCCGAt  >  1:305565/1‑140 (MQ=255)
                                                                                                                      |                                                                                                                                        
AAATGGTCATTATTATTGGGCATTACCAACCTGTATTGTTACCGGGTCCTGCTCTTGTAGGAAAAAGTATATGGTCTTTCATTGTTACTGGAGAAATTTTCCAACATTTAGCAATTAGTTTATGGAGATTTGTAGCGGGCTTTGTTGTCGCATTG‑TTGGTTGCTATTCCATTGGGCTTCTTGCTTGGAAGGAATCGTTGG‑CTATACAACGCTATCGAACCGCTATTTCAATTGATTAGGCCGATATCTCCGAT  >  CP000730/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: ‑