Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 599,226 G→C 11.1% P228R (CCA→CGA)  USA300HOU_0542 ← M20D subfamily peptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP000730599,2260GC11.1% 101.4 / 5.5 36P228R (CCA→CGA) USA300HOU_0542M20D subfamily peptidase
Reads supporting (aligned to +/- strand):  ref base G (17/15);  new base C (2/2);  total (19/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.70e-01

TATCCATTTTGTCCATTACATGATTCTGTATATCTGAATCGAAAGTTCTAACTGTACCTTTACAAAATGCTTGATCAGGAATAACGCTATCTGTGGTGCCTGCTTGAATCATTCCAAATGAAAGTACAGCTTGTTTAACTGGATCGATCGTACGTGAAATTATTTTTTGTGCACTTAAAATGAACTCTGCCATGATTACTATTGGGTCAATGGTTTCATGAGGTTTGGCACCATGACCACCACGACCTTTAAATGTGACGCTAAATTCA  >  CP000730/599086‑599354
                                                                                                                                            |                                                                                                                                
tATCCATTTTGTCCATTACATGATTCTGTATATCTGAATCGAAAGTTCTAACTGTACCTTTACAAAATGCTTGATCAGGAATAACGCTATCTGTGGTGCCTGCTTGAATCATTCCAAATGAAAGTACAGCTTGTTTAACTg                                                                                                                                  >  1:88956/1‑141 (MQ=255)
     aTTTTGTCCATTACATGATTCTGTATATCTGAATCGAAAGTTCTAACTGTACCTTTACAAAATGCTTGATCAGGAATAACGCTATCTGTGGTGCCTGCTTGAATCATTCCAAATGAAAGTACAGCTTGTTTAACTGgatcg                                                                                                                             <  1:305430/141‑1 (MQ=255)
        ttGTCCATTACATGATTCTGTATATCTGAATCGAAACTTCTAACTGTACCTTTACAAAATGCTTGATCAGGAATAACGCTATCTGTGGTGCCTGCTTGAATCATTCCAAATGAAAGTACAGCTTGTTTAACTGgatcgatc                                                                                                                          <  1:184707/141‑1 (MQ=255)
         tGTCCATTACATGATTCTGTATATCTGAATCGAAAGTTCTAACTTTACCTTTACAAAATGCTTGATCAGGAATAACGCTATCTGTGGTGCCTGCTTGAATCATTCCAAATGAAAGTACAGCTTGTTTAACTGGatcgatcg                                                                                                                         <  1:172800/141‑1 (MQ=255)
            ccATTACATGATTCTGTATATCTGAATCGAAAGTTCTAACTGTACCTTTACAAAATGCTTGATCAGGAATAACGCTATCTGTGGTGCCTGCTTGAATCATTCCAAATGAAAGTACAGCTTGTTTAACTGGATCGATCGTag                                                                                                                      >  1:109677/1‑140 (MQ=255)
             cATTACATGATTCTGTATATCTGAATCGAAAGTTCTAACTGTACCTTTACAAAATGCTTGATCAGGAATAACGCTATCTGTGGTGCCTGCTTGAATCATTCCAAATGAAAGTACAGCTTGTTTAACTGGATCGATCGTACg                                                                                                                     >  1:267294/1‑141 (MQ=255)
                   aTGATTCTGTATATCTGAATCGAAAGTTCTAACTGTACCTTTACAAAATGCTTGATCAGGAATAACGCTATCTGTGGTGCCTGCTTGAATCATTCCAAATGAAAGTACAGCTTGTTTAACTGGATCGATCGTACGTGAAat                                                                                                               >  1:166805/1‑141 (MQ=255)
                   aTGATTCTGTATATCTGAATCCAAAGTTCTAACTGTACCTTTACAAAATGCTTGATCAGGAATAACGCTAGCTGTGGTGCCTGCTTGAATCATTCCAAATGAAAGTACAGCTTGTTTAACTGGATCGATCGTACGTGAAat                                                                                                               <  2:65147/141‑1 (MQ=255)
                         cTGTATATCTGAATCGAAAGTTCTAACTGTACCTTTACAAAATGCTTGATCAGGAATAACGCTATCTGTGGTGCCTGCTTGAATCATTCCAAATGAAAGTACAGCTTGTTTAACTGGATCGATCGTACGTGAAATTAtttt                                                                                                         >  2:291983/1‑141 (MQ=255)
                              tatCTGAATCGAAAGTTCTAACTGTACCTTTACAAAATGCTTGATCAGGAATAACGCTATCTGTGGTGCCTGCTTGAATCATTCCAAATGAAAGTACAGCTTGTTTAACTCGATCGATCGTACGTGAAATTATTTTTtgtg                                                                                                    >  1:320715/1‑141 (MQ=255)
                              tatCTGAATCGAAAGTTCTAACTGTACCTTTACAAAATGCTTGATCAGGAATAACGCTATCTGTGGTGCCTGCTTGAATCATTCCAAATGAAAGTACAGCTTGTTTAACTCGATCGATCGTACGTGAAATTATTTTTtgtg                                                                                                    >  1:29820/1‑141 (MQ=255)
                              tatCTGAATCGAAAGTTCTAACTGTACCTTTACAAAATGCTTGATCAGGAATAACGCTAGCTGTGGTGCCTGCTTGAATCATTCCAAATGAAAGTACAGCTGGTTTAACTGGATCGATCGTACGTGAAATTATGTTTtgtg                                                                                                    >  1:36597/1‑141 (MQ=255)
                                 cTGAATCGAAAGTTCTAACTGTACATTTACAAAATGCTTGATCAGGAATAACGCTATCTGTGGTGCCTGCTTGAATCATTCCAAATGAAAGTAAAGCTTGTTTAACTGGATCGATCGTACGTGAAATTATTCTTTGTGCAc                                                                                                 >  1:309635/1‑141 (MQ=255)
                                       cGAAAGTTCTAACTGTACCTTTACAAAATGCTTGATCAGGAATAACGCTATCTGTGGTGCCTGCTTGAATCATTCCAAATGAAAGTACAGCTTGTTTAACTGGATCGATCGTACGTGAAATTATTTTTTGTGCACTTaaaa                                                                                           <  2:158154/141‑1 (MQ=255)
                                       cGAAAGTTCTAACTGTACCTTTACAAAATGCTTGATCAGGAATAACGCTATCTGTGGTGCCCGCTTGAATCATTCCAAATGAAAGTACAGCTTGTTTAACTGGATCGATCGTACGTGAAATTATTTTTTGTGCACTTaaaa                                                                                           <  1:74668/141‑1 (MQ=255)
                                         aaaGTTCTAACTGTACCTTTACAAAATGCTTGATCAGGAATAACGCTATCTGTGGTGCCTGCTTGAATCATTCCAAATGAAAGTACAGCTTGTTTAACTGGATCGATCGTACGTGAAATTTTTTTTTGTGCACTTAAAATg                                                                                         >  1:27597/1‑141 (MQ=255)
                                                 aaCTGTACCTTTACAAAATGCTTGATCAGGAATAACGCTATCTGTGGTGCCTGCTTGAATCATTCCAAATGAAAGTACAGCTTGTTTAACTGGATCGATCGTACGTGAAATTATTTTTTGTGCACTTAAAATGAACTCTGc                                                                                 >  2:114645/1‑141 (MQ=255)
                                                     gTACCTTTACAAAATGCTTGATCAGGAATAACGCTATCTGTGGTGCCTGCTTGAATCATTCCAAATGAAAGTACAGCTTGTTTAACTGGATCGATCGTACGTGAAATTATTTTTTGTGCACTTAAAATGAACTCTGCCATg                                                                             <  2:259243/141‑1 (MQ=255)
                                                       aCCTTTAAAAAATGGTTGATCAGGAATAACGCTATCTGTGGTGCCTGCTTGAATCATTCCAAATGAAAGTACAGCTTGTTTAACTGGATCGATCGTACGTGAAATTATTTTTTGTGCACTTAAAATGAACTCTGCCATGAt                                                                           <  2:24036/141‑1 (MQ=255)
                                                              cAAAATGCTTGATCAGGAATAACGCTATCTGTGGTGCCTGCTTGAATCATTCCAAATGAAAGTACAGCTTGTTTAACTGGATCGATCGTACGTGAAATTATTTTTTGTGCACTTAAAATGAACTCTGCCATGATTACTAtt                                                                    <  1:184504/141‑1 (MQ=255)
                                                                         aTCAGGAATAACGCTATCTGTGGTCCCTGCTTGAATCATTCCAAATGAAAGTACAGCTTGTTTAACTGGATCGATCGTACGTGAAATTATTTTCTGTGCACTTAAAATGAACTCTGCCATGATTACTATTGGGTCAATGGt                                                         <  2:208339/141‑1 (MQ=255)
                                                                           cAGGAATAACGCTATCTGTGGTGCCTGCTTGAATCATTCCAAATGAAAGTACAGCTTGTTTAACTGGATCGATCGTACGTGAAATTATTTTTTGTGCCCTTAAAATGAACTCTGCCATGATTACTATTGGGTCACTGGttt                                                       >  1:153522/1‑141 (MQ=255)
                                                                             ggAATAACGCTATCTGTGGTGCCTGCTTGAATCATTCCAAATGAAAGTACAGCTTGTTTAACTGGATCGATCGTACGTGAAATTATTTTTTGTGCACTTAAAATGAACTCTGCCATGATTATTTTTGGGTCATTGGtttgt                                                     >  2:200158/1‑139 (MQ=255)
                                                                               aaTAACGCTATCTGTGGTGCCTGCTTGAATCATTCCAAATGAAAGTACAGCTTGTTTAACTGGATCGATCGTACGTGAAATTATTTTTTGTGCACTTAAAATGAACTCTGCCATGATTACTATTGGGTCAATGGTTTCATg                                                   >  1:327432/1‑141 (MQ=255)
                                                                                             tggtgCCTGCTTGAATCATTCCAAATGAAAGTACAGCTTGTTTAACTGGATCGATCGTACGTGAAATTATTTTTTGTGCACTTAAAATGAACTCTGCCATGATTACTATTGGGTCAATGGTTTCATGAGGTTTGGCACCAt                                     >  2:177727/1‑141 (MQ=255)
                                                                                                   cTGCTTGAATCATTCCAAATGAAAGTACAGCTTGTTTAGCTGGATCGATCGTACGTGAAATTATTTTTTGTGCACTTAAAATGAACTCTGCCATGATTACTATTGGGTCAATGGTTTCATGAGGTTTGGCACCATGaccac                               <  1:45391/141‑1 (MQ=255)
                                                                                                   cTGCTTGAATCATTCCAAATGAAAGTACAGCTTGTTTAACTGGATCGATCGTACGTGAAATTATTTTTTGTGCACTTAAAATGAACTCTGCCATGATTACTATTGGGTCAATGGTTTCATGAGGTTTGGCACCATGaccac                               <  1:308574/141‑1 (MQ=255)
                                                                                                     gCTTGAATCATTCCAAATGAAAGTACAGCTTGTTTAACTGGATCGATCGTACGTGAAATTATTTTTTGTGCACTTAAAATGAACTCTGCCATGATTACTATTGGGTCAATGGTTTCATGAGGTTTGTCACCCTGAccacca                             >  1:29055/1‑141 (MQ=255)
                                                                                                      cTTGAATCATTCCAAATGAAAGTACAGCTTGTTTAACTGGATCGATCGTACGTGAAATTATTTTTTGTGCACTTAAAATGAACTCTGCCATGATTACTATTGGGTCAATGGTTTCATGAGGTTTGGCACCATGACcaccac                            >  1:284456/1‑141 (MQ=255)
                                                                                                                    aaTGAAAGTACAGCTTGTTTAACTGGATCGATCGTACGTGAAATTATTTTTTGTGCACTTAAAATGAACTCTGCCATGATTACTATTGGGTCAATGGTTTCATGAGGTTTGGCACCATGACCACGACGACGTTTAAAtgtg              >  2:1164/1‑141 (MQ=255)
                                                                                                                        aaaGTACAGCTTGTTTAACTGGATCGATCGTACGTGAAATTATTTTTTGTGCACTTAAAATGAACTCTGCCATGATTACTATTGGGTCAATGGTTTCATGAGGTTTGGCACCATGACCACCACGACCTTTAAATGTGACGc          <  1:114645/141‑1 (MQ=255)
                                                                                                                        aaaGTACAGCTTGTTTAACTGGATCGATCGTACGTGAAATTATTTTTTGTGCACTTAAAATGAACTCTGCCATGATTACTATTGGGTCAATGGTTTCATGAGGTTTGGCACCATGACCACCACGACCTTTAAATGTGACGc          <  1:112692/141‑1 (MQ=255)
                                                                                                                        aaaGTACAGCTTGTTTAACTGGATCGATCGTACGTGAAATTATTTTTTGTGCACTTAAAATGAACTCTGCCATGATTACTATTGGGTCAATGGTTTCATGAGGTTTGGCACCATGACCACCACGACCTTTAAATGTGACGc          >  2:299057/1‑141 (MQ=255)
                                                                                                                        aaaGTACAGCTTGTTGAACTGGATCGATCGTACGTGAAATTATTTTTTGTGCACTTAAAATGAACTCTGCCATGATTACTATTGGGTCAATGGTTTCATGAGGTTTGGCACCATCACCACCACGACCTTTAAATGTGACGc          <  2:102617/141‑1 (MQ=255)
                                                                                                                          aGTACAGCTTGTTTAACTCGATCGATCGTACGTGCAATTAATTTTTGCGCCCACAAAAGGAACTCTGCCATGATTACTATTGGGTCAATGGTTTCATGAGGTTTGGCACCATGACCACCACGACCTTTAAATGTGACGCTa        <  2:320715/141‑1 (MQ=255)
                                                                                                                          aGTACAGCTTGTTTAACTCGATCGATCGTACGTGAAATTATTTTTTGTGCACTTAAAATGAACTCTGCCATGATTACTATTGGGTCAATGGTTTCATGAGGTTTGGCACCATGACCACCACGACCTTTAAATGTGACGCTa        <  2:29820/141‑1 (MQ=255)
                                                                                                                                gCTTGTTTAACTGGATCGATCGTACGTGAAATTATTTTTTGTGCACTTAAAATGAACTCTGCCATGATTACTATTGGGTCAATGGTTTCATGAGGTTTTGCACCATGACCACCACGACCTTTAAATGTGACGCTAAATTCa  >  2:110076/1‑141 (MQ=255)
                                                                                                                                            |                                                                                                                                
TATCCATTTTGTCCATTACATGATTCTGTATATCTGAATCGAAAGTTCTAACTGTACCTTTACAAAATGCTTGATCAGGAATAACGCTATCTGTGGTGCCTGCTTGAATCATTCCAAATGAAAGTACAGCTTGTTTAACTGGATCGATCGTACGTGAAATTATTTTTTGTGCACTTAAAATGAACTCTGCCATGATTACTATTGGGTCAATGGTTTCATGAGGTTTGGCACCATGACCACCACGACCTTTAAATGTGACGCTAAATTCA  >  CP000730/599086‑599354

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: