Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 328,487 C→A 19.4% I210I (ATC→ATA)  USA300HOU_0294 → hypothetical membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP000730328,4870CA19.4% 68.5 / 12.7 31I210I (ATC→ATA) USA300HOU_0294hypothetical membrane protein
Reads supporting (aligned to +/- strand):  ref base C (13/12);  new base A (4/2);  total (17/14)
Fisher's exact test for biased strand distribution p-value = 6.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.30e-02

ATATAAGTCTAATAGTTCCTGTAACACAGCAAGAAGAGATTTTTATAGGGTTAAGTATATTCATCAATATAATTGTTGATGCTCTTCTTGTTAGGTCATATTTTAATTACGCGTTATATAAGAGTATTAAGAAAGATATCGAAAATGAAGGTAAGACAATGTAAAGGAAAGACATAAAAAATAATAGGAGTCTAACCATGATAATTTTTATTTTAATAACAATATTTGCTATTTATTATATAGCTATGATTGCTAGTTTGTTTAAAAGTGAAGGTTTTT  >  CP000730/328348‑328626
                                                                                                                                           |                                                                                                                                           
atataAGTCTAATAGTTCCTGTAACACAGCAAGAAGAGATTTTTATAGGGTTAAGTATATTCATCAATATAATTGTTGATGCTCTTCTTGTTAGGTCATATTTTAATTACGCGTTATATAAGAGTATTAAGAAAGATATCg                                                                                                                                            <  2:161130/141‑1 (MQ=255)
   taAGTCTAATAGTTCCTGTAACACAGCAAGAAGAGATTTTTATAGGGTTAAGTATATTCATCAATATAATTGTTGATGCTCTTCTTGTTAGGTCATATTTTAATTACGCGTTATATAAGAGTATTAAGAAAGATATCGaaa                                                                                                                                         >  1:472625/1‑141 (MQ=255)
                tCCTGTAACACAGCAAGAAGAGATTTTTATAGGGTTAAGTATATTCATCAATATAATTGTTGATGCTCTTCTTGTTAGGTCATATTTTAATTACGCGTTATATAAGAGTATTAAGAAAGATATCGAAAATGAAGGTAAGAc                                                                                                                            <  1:66068/141‑1 (MQ=255)
                tCCTGTAACACAGCAAGAAGAGATTTTTATAGGGTTAAGTATATTCATCAATATAATTGTTGATGCTCTTCTTGTTAGGTCATATTTTAATTACGCGTTATATAAGAGTATTAAGAAAGATATCGAAAATGAAGGTAAGAc                                                                                                                            >  2:331775/1‑141 (MQ=255)
                    gTAACACAGCAAGAAGAGATTTTTATAGGGTTAAGTATATTCATCAATATAATTGTTGATGCTCTTCTTGTTAGGTCATATTTTAATTACGCGTTATATAAGAGTATTAAGAAAGATATCGAAAATGAAGGTAAGACAATg                                                                                                                        <  1:201587/141‑1 (MQ=255)
                       acacAGCAAGAAGAGATTTTTATAGGGTTAAGTATATTCATCAATATAATTGTTGATGCTCTTCTTGTTAGGTCATATTTTAATTACGCGTTATATAAGAGTATTAAGAAAGATATCGAAAATGAAGGTAAGACAATGTaa                                                                                                                     >  1:363864/1‑141 (MQ=255)
                                          ttATAGGGTTAAGTATATTCATCAATATAATTGTTGATGCTCTTCTTGTTAGGTCATATTTTAATTACGCGTTATATAAGAGTATTAAGAAAGATATCGAAAATGAAGGTAAGACAATGTAAAGGAAAGACATAAAAaata                                                                                                  >  2:178856/1‑141 (MQ=255)
                                             taGGGTTAAGTATATTTATCAATATAATTGTTGATGCTCTTCTTGTTAGGTCATATTTCAATTACGCCTTATATAAGAGTAGTAATAAAGATATCGAAAATGAAGGTAAGACAATGTAAAGGAAAGACATAAAAAataata                                                                                               <  1:260438/141‑1 (MQ=255)
                                             taGGGTTAAGTATATTCATCAATATAATTGTTGATGCTCTTCTTGTTAGGTCATATTTTAATTACGCGTTATATAAGAGTATTAAGAAAGATATCGAAAATGAAGGTAAGACAATGTAAAGGAAAGACATAAAAAataata                                                                                               <  1:227907/141‑1 (MQ=255)
                                               gggTTAAGTATATTCATCAATATAATTGTTGATGCTCTTCTTGTTAGGTCATATTTTAATTACGCGTTATATAAGAGTATTAAGAAAGATATCGAAAATGAAGGTAAGACAATGTAAAGGAAAGACATAAAAAATAATAgg                                                                                             >  1:505405/1‑141 (MQ=255)
                                                     aGTATATTCATCAATATAATTGTTGATGCTCTTCTTGTTAGGTCATATTTTAATTACGCGTTATATAAGAGTATTAAGAAAGATATCGAAAATGAAGGTAAGACAATGTAAAGGAAAGACATAAAAAATAATAGGAGTCTa                                                                                       <  1:164021/141‑1 (MQ=255)
                                                            tcatcaATATAATTGTTGATGCTCTTCTTGTTAGGTCATATTTTAATTACGCGTTATATAAGAGTATTAAGAAAGATATCGAAAATGAAGGTAAGACAATGTAAAGGAAAGACATAAAAAATAATAGGAGTCTAACCATGa                                                                                <  2:470842/141‑1 (MQ=255)
                                                             catcaATATAATTGTTGATGCTCTTCTTGTTAGGTCATATTTTAATTACGCGTTATATAAGAGTATTAAGAAAGATATCGAAAATGAAGGGAAGACAATGTAAAGGAAAGACATAAAAAATAATAGGAGTCTAACCATGAt                                                                               >  2:432936/1‑141 (MQ=255)
                                                               tcaATATAATTGTTGATGCTCTTCTTGTTAGGTCATATTTTAATTACGCGTTATATAAGAGTATTAAGAAAGATATCGAAAATGAAGGTAAGACAATGTAAAGGAAAGACATAAAAAATAATAGGAGTCTAACCATGATaa                                                                             >  1:75441/1‑141 (MQ=255)
                                                               tcaATATAATTGTTGATGCTCTTCTTGTTAGGTCATATTTTAATTACGCGTTATATAAGAGTATTAAGAAAGATATCGAAAATGAAGGTAAGACAATGTAAAGGAAAGACATAAAAAATAATAGGAGTATAACCATGATaa                                                                             >  2:375499/1‑141 (MQ=255)
                                                                     taATTGTTGATGCTCTTCTTGTTAGGTCATATTTTAATTACGCGTTATATAAGAGTATTAAGAAAGATATCGAAAATGAAGGTAAGACAATGTAAAGGAAAGACATAAAAAATAATAGGAGTCTAACCATGATAATtttta                                                                       >  2:49206/1‑141 (MQ=255)
                                                                     taATTGTTGATGCTCTTCTTGTTAGGTCATATTTTAATTACGCGTTATATAAGAGTATTAAGAAAGATATCGAAAATGAAGGTAAGACAATGTAAAGGAAAGACATAAAAAATAATAGGAGTCTAACCATGATAATtttta                                                                       <  1:55076/141‑1 (MQ=255)
                                                                         tgttgATGCTCTTCTTGTTAGGTCATATTTTAATTACGCGTTATATAAGAGTATTAAGAAAGATATAGAAAATGAAGGTAAGACAATGTAAAGGaa                                                                                                                >  2:250283/1‑96 (MQ=255)
                                                                         tgttgATGCTCTTCTTGTTAGGTCATATTTTAATTACGCGTTATATAAGAGTATTAAGAAAGATATAGAAAATGAAGGTAAGACAATGTAAAGGaa                                                                                                                <  1:250283/96‑1 (MQ=255)
                                                                             gATGCTCTTCTTGTTAGGTCATATTTTAATTACGCGTTATATAAGAGTATTAAGAAAGATATCGAAAATGAAGGTAAGACAATGTAAAGGAAAGACATAAAAAATAATAGGAGTCTAACCATGATAATTTTTATTTtaata                                                               <  2:363864/141‑1 (MQ=255)
                                                                                                 acatgttAATTACGCGTGATATAAGAGTATTAAGACAGATATCGAATATGATGGCAAGTTGATGTAAAGGAAAAACATAACCAATAATAGATGGTTTACCATGATAATTTTTATTTTAATATCAACGTTTGCTATttatta                                           <  2:317917/136‑1 (MQ=255)
                                                                                                    ttttAATTACGCGTTATATAAGAGTATTAAGAAAGATATAGAAAATGAAGGTAAGACAATGTAAAGGAAAGACATAAAAAATAATAGGAGTCTAACCATGATAATTTTTATTTTAATAACAATATTTGCTATTTATtatat                                        <  1:442203/141‑1 (MQ=255)
                                                                                                       tAATTACGCGTTATATAAGAGTATTAAGAAAGATATAGAAAATGAAGGTAAGACAATGTAAAGGAAAGACATAAAAAATAATAGGAGTCTAACCATGATAATTTTTATTTTAATAACAATATTTGCTATTTATTATATAGc                                     >  1:203935/1‑141 (MQ=255)
                                                                                                          ttACGCGTTATATAAGAGTATTAAGAAAGATATAGAAAATGAAGGTAAGACAATGTAAAGGAAAGACATAAAAAATAATAGGAGTCTAACCATGATAATTTTTATTTTAATAACAATATTTGCTATTTATTATATAGCTAt                                  >  2:423262/1‑141 (MQ=255)
                                                                                                             cgcgTTATATAAGAGTATTAAGAAAGATATCGAAAATGAAGGTAAGACAATGTAAAGGAAAGACATAAAAAATAATAGGAGTCTAACCATGATAATTTTTATTTTAATAACAATATTTGCTATTTATTATATAGCTATGAt                               >  1:495752/1‑141 (MQ=255)
                                                                                                              gcgTTATATAAGAGTATTAAGAAAGATATAGAAAATGAAGGTAAGACAATGTAAAGGAAAGACATAAAAAATAATAGGAGTCTAACCATGATAATTTTTATTTTAATAACAATATTTGCTATTTATTATATAGCTATGAtt                              >  1:489250/1‑141 (MQ=255)
                                                                                                               cgTTATATAAGAGTATTAAGAAAGATATCGAAAATGAAGGTAAGACAATGTAAAGGAAAGACATAAAAAATAATAGGAGTCTAACCATGATAATTTTTATTTTAATAACAATATTTGCTATTTATTATATAGCTATGATTg                             <  2:356359/141‑1 (MQ=255)
                                                                                                               aattctaaaaGTGTATTAAGTAAGATATCGAACATGAAGTGAAGACAGTATAAAGCAAAGACATATTAAATAAAGGGAGTCTGACCATGATAATTTTTATTTTAATATCATTTTTTGCTATTTATTCTTTATCTCTAATTg                             <  2:117538/133‑1 (MQ=255)
                                                                                                                         gagTATTAAGAAAGATATCGAAAATGAAGGTAAGACAATGTAAAGGAAAGACATGAAAAATAATTGGAGTCTAACCATGATAATTTTTATTTTAATAACAATATTTGCTATTTATTATATAGCTATGATTGCTAgtttgtt                   >  2:379873/1‑141 (MQ=255)
                                                                                                                                          tCGAAAATGAAGGTAAGACAATGTAAAGGAAAGACATAAAAAATAATAGGAGTCTAACCATGATAATTTTTATTTTAATAACAATATTTGCTATTTATTATATAGCTATGATTGCTAGTTTGTTTAAAAGTGAAGGttttt  >  2:428632/1‑141 (MQ=255)
                                                                                                                                          tCGAAAATGAAGGTAAGACAATGTAAAGGAAAGACATAAAAAATAATAGGAGTCTAACCATGATAATTTTTATTTTAATAACAATATTTGCTATTTATTATATAGCTATGATTGCTAGTTTGTTTAAAAGTGAAGGttttt  >  1:231234/1‑141 (MQ=255)
                                                                                                                                           |                                                                                                                                           
ATATAAGTCTAATAGTTCCTGTAACACAGCAAGAAGAGATTTTTATAGGGTTAAGTATATTCATCAATATAATTGTTGATGCTCTTCTTGTTAGGTCATATTTTAATTACGCGTTATATAAGAGTATTAAGAAAGATATCGAAAATGAAGGTAAGACAATGTAAAGGAAAGACATAAAAAATAATAGGAGTCTAACCATGATAATTTTTATTTTAATAACAATATTTGCTATTTATTATATAGCTATGATTGCTAGTTTGTTTAAAAGTGAAGGTTTTT  >  CP000730/328348‑328626

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