Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 730,535 T→C 11.0% intergenic (+142/‑187) lipA → / → USA300HOU_0677 possible S9C family peptidase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP000730730,5350TC11.0% 88.0 / 2.7 38intergenic (+142/‑187)lipA/USA300HOU_0677possible S9C family peptidase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (20/13);  new base C (2/2);  total (23/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

AGTAATTTATTGCTTAAGCAACATTTAAGATTTTCAAATTAAAAACGAAGAATTTAAAACATGTGGTGCTAATGTGTAAGAGTCTGAGATATAATAAATTGTTAGCGGTTCTTTATCATTTCTATCTCACTCTATTATGATGTGACATTTATTTTACAAACTAATTTGTTTTGAACCTGAAAATAACTTTTTACAACTAAATTGTCGAAAAACAGTGTGTACATGATATAATAAATTTATAAATTGAAAAGAATTCAAAAGA  >  CP000730/730405‑730666
                                                                                                                                  |                                                                                                                                   
agtaATTTATTACTTAAGCAACATTTAAGATCTTCAAATTGAAAACGAAGAATTTAAAACATGTGGGGCTAATGTGTAAGAGTCTGAGATATAATCAATTGTTAGCGGTTCTTTATCATTTCTATCTCACCCTATgatgat                                                                                                                           >  1:240913/1‑141 (MQ=255)
     tttATTGCTTAAGCAACATTTAAGATTTTCAAATTAAAAACGAAGAATTTAAAACATGTGGTGCTAATGTGTAAGAGTCTGAGATATAATAAATTGTTAGCGGTTCTTTATCATTTCTATCTCACTCTATTATGATGTGAc                                                                                                                      >  2:4698/1‑141 (MQ=255)
                    aCATTTAAGATTTTCAAATTAAAAACGAAGAATTTAAAACATGTGGTGCTAATGTGTAAGAGTCTGAGATATAATAAATTGTTAGCGGTTCTTTATCATTTCTATCTCACTCTATTATGATGTGACATTTATTTTACAAAc                                                                                                       >  1:173665/1‑141 (MQ=255)
                     cATTTAAGATCTTCAAATTGAAAACGAAGAATTTAAAACATGTGGGGCTAATGTGTAAGAGTCTGAGATATAATAAATTGTTAGCGGTTCTTTATCATTTCTATCTCACCCTATTATGATGTGACATTTATTTTACAAACt                                                                                                      >  1:164403/1‑141 (MQ=255)
                         tAAGATTTTCAAATTAAAAACGAAGAATTTAAAACATGTGGTGCTAATGTGTAAGAGTCTGAGATATAATAAATTGTTAGCGGTTCTTTATCATTTCTATCTCACTCTATTATGATGTGACATTTATTTTACAAACTAAtt                                                                                                  >  2:30492/1‑141 (MQ=255)
                         tAAGATTTTCAAATTAAAAACGAAGAATTTAAAACATGTGGTGCTAATGTGTAAGAGTCTGAGATATAATAAATTGTTAGCGGTTCTTTATCATTTCTATCTCACTCTATTATGATGTGACATTTATTTTACAAACTAAtt                                                                                                  <  1:117193/141‑1 (MQ=255)
                            gATTTTCAAATTAAAAACGAAGAATTTAAAACATGTGGTGCTAATGTGTAAGAGTCTGAGATATAATAAATTGTTAGCGGTTCTTTATCATTTCTATCTCACTCTATTATGATGTGACATTTATTTTACAAACTAATTTGt                                                                                               >  1:456663/1‑141 (MQ=255)
                                ttCAAATTAAAAACGAAGAATTTAAAACATGTGGTGCTAATGTGTAAGAGTCTGAGATATAATAAATTGTTAGCGGTCCTTTATCATTTCTATCTCACTCTATTATGATGTGACATTTATTTTACAAACTAATTTGTTTTg                                                                                           <  2:131360/141‑1 (MQ=255)
                                  cAAATTAAAAACGAAGAATTTAAAACATGTGGTGCTAATGTGTAAGAGTCTGAGATATAATAAATTGTTAGCGGTTCTCTATCATTTCTATCTCACTCTATTATGATGTGACATTTATTTTACAAACTAATTTGTTTTGaa                                                                                         >  2:385052/1‑141 (MQ=255)
                                    aaTTAAAAACGAAGAATTTAAAACATGTGGTGCTAATGTGTAAGAGTCTGAGATATAATAACTACTTAGCGGTACATTAACCATACCATCACGAACTTTAATACTAAGGAATTTATTTTACAAACTAATTTGTTTTGaaac                                                                                       >  1:293369/1‑139 (MQ=255)
                                    aaTTAAAAACGAAGAATTTAAAACATGTGGTGCTAATGTGTAAGAGACTGAGATATAATAAATTGTTAGCGGTTCTTTATCATTTCTATCTCACTCTATTATGATGTGACATTTATTTTACAAACTAATTTGTTTTTAAcc                                                                                       >  2:140339/1‑141 (MQ=255)
                                          aaaCGAAGAATTTAAAACAGGTGGTGCTAATGTGTAAGAGTCTGGGATATAATAAATTGTGAGCGGTTCTTTATCATTTCTATCTCACTCTATTATGATGTGACATTTATTTTACAAACTAATTTGTTTTGAACCTGaaaa                                                                                 >  2:11419/1‑141 (MQ=255)
                                                     ttAAAACATGTGGTGCTAATGTGTAAGAGTCTGAGATATAATAAATTGTTAGCGGTTCTTTATCATTTCTATCTCACTCTATTATGATGTGACATTTATTTTACAAACTAATTTGTTTTGAACCTGAAAATAACTTTTTAc                                                                      >  1:158894/1‑141 (MQ=255)
                                                      tAAAACATGTGGTGCTAATGTGTAAGAGTCTGAGATATAATAAATTGTTAGCGGTTCTTTATCATTTCTATCTCACTCTATTATGATGTGACATTTATTTTACAAACTAATTTGTTTTGAACCTGAAAATAACTTTTTACa                                                                     <  1:39736/141‑1 (MQ=255)
                                                         aaCATGTGGTGCTAATGTGTAAGAGTCTGAGATATAATAAATTGTTAGCGGTTCTTTATCATTTCTATCTCACTCTATTATGATGTGACATTTATTTTACAAACTAATTTGTTTTGAACCTGAAAATAACTTTTTACAACt                                                                  >  2:324098/1‑141 (MQ=255)
                                                         aaCATGTGGTGCTAATGTGTAAGAGTCTGAGATATAATAAATTGTTAGCGGTTCTTTATCATTTCTATCTCACTCTATTATGATGTGACATTTATTTTACAAACTAATTTGTTTTGAACCTGAAAATAACTTTTTACAACt                                                                  >  2:242351/1‑141 (MQ=255)
                                                             tgtgGTGCTAATGTGTAAGAGTCTGAGATATAATAAATTGTTAGCGGTTCTTTATCATTTCTATCTCACTCTATTATGATGTGACATTTATTTTACAAACTAATTTGTTTTGAACCTGAAAATAACTTTTTTCTATTAAAt                                                              >  2:189301/1‑141 (MQ=255)
                                                             tgtgGTGCTAATGTGTAAGAGTCTGAGATATAATAAATTGTTAGCGGTTCTTTATCATTTCTATCTCACTCTATTATGATGTGACATTTATTTTACAAACTAATTTGTTTTGAACCTGAAAATAACTTTTTACAACTAAAt                                                              >  2:288075/1‑141 (MQ=255)
                                                               tggggCTAATGTGTAAGAGTCTGAGATATAATAAATTGTTAGCGGTTCTTTATCATTTCTATCTCACCCTATTATGATGTGACATTTATTTTACAAACTAATTTGTTTTGAACCTGAAAATAACTTTTTACAACTAAATTg                                                            <  1:353057/137‑1 (MQ=255)
                                                                         gtgtAAGAGTCTGAGGTATAATAAATTGTTAGCGGTTCTTTATCATTTCTATCTCACTCTATTATGATGTGACATTTATTTTACAAACTAATTTGTTTTGAACCTGAAAATAACTTTTTACAACTAAATTGTCGAAAAACa                                                  <  1:191137/141‑1 (MQ=255)
                                                                          tgtAAGAGTCTGAGATATAATAAATTGTTAGCGGTTCTTTATCATTTCTATCTCACTCTATTATGATGTGACATTTATTTTACAAACTAATTTGTTTTGAACCTGAAAATAACTTTTTACAACTAAATTGTCGAAAAACAg                                                 <  1:135014/141‑1 (MQ=255)
                                                                                    tGAGATATAATAAATTGTTAGCGGTTCTTTATCATTTCTATCTCACTCTATTATGATGTGACATTTATTTTACAAACTAATTTGTTTTGAACCTGAAAATAACTTTTTACAACTAAATTGTCGAAAAACAGTGTGTACATg                                       >  2:47394/1‑141 (MQ=255)
                                                                                    tGAGATATAATAAATTGTTAGCGGTTCTTTATCATTTCTATCTCACTCTATTATGATGTGACATTTATTTTACAAACTAATTTGTTTTGAACCTGAAAATAACTTTTTACAACTAAATTGTCGAAAAACAGTGTGTACATg                                       <  2:260235/141‑1 (MQ=255)
                                                                                    tGAGATATAATAAATTGTTAGCGGTTCTTTATCATTTCTATCTCACTCTATTATGATGTGACATTTATTTTACAAACTAATTTGTTTTGAACCTGAAAATAACTTTTTACAACTAAATTGTCGAAAAACAGTGTGTACATg                                       <  1:31863/141‑1 (MQ=255)
                                                                                    tGAGATATAATAAATTGTTAGCGGTTCTTTATCATTTCTATCTCACTCTATTATGATGTGACATTTATTTTACAAACTAATTTGTTTTGAACCTGAAAATAACTTTTTACAACTAAATTGTCGAAAAACAGTGTGTACATg                                       <  2:260197/141‑1 (MQ=255)
                                                                                     gagaTATAATAAATTGTTAGCGGTTCTTTATCATTTCTATCTCTCTCTATTATGATGTGACATTTATTTTACAAACTAATTTGTTTTGAACCTGAAAATAACTTTTTACAACTAAATTGTCGAAAAACAGTGTGTACATGa                                      <  1:365752/141‑1 (MQ=255)
                                                                                        atatAATAAATTGTTAGCGGTTCTTTATCATTTCTATCTCACTCTATTATGATGTGACATTTATTTTACAAACTAATTTGTTTTGAACCTGAAAATAACTTTTTACAACTAAATTGTCGAAAAACAGTGTGTACATGatat                                   >  1:319633/1‑141 (MQ=255)
                                                                                             ataaATTGTTAGCGGTTCTTTATCATTTCTATCTCACTCTATTATGATGTGACATTTATTTTACAAACTAATTTGTTTTGAACCTGAAAATAACTTTTTACAACTAAATTGTCGAAAAACAGTGTGTACATGATAtaataa                              >  2:247669/1‑141 (MQ=255)
                                                                                              taaATTGTTAGCGGGTCTTTATCATTTCTATCTCACTCTATTATGATGTGACATTTATTTTACAAACTAATTTGTTTTGAACCTGACAATAACTTTTTACAACTAAATTGTCGAAAACCAGTGTGTACATGATATAATaaa                             >  2:263164/1‑141 (MQ=255)
                                                                                                   tGTTAGCGGTTCTTTATCATTTCTATCTCACTCTATTATGATGTGACATTTATTTTACAAACTAATTTGTTTTGAACCTGAAAATAACTTTTTACAACTAAATTGTCGAAAAACAGTGTGTACATGATATAATAAATTtat                        <  1:360821/141‑1 (MQ=255)
                                                                                                    gtgaGCGGTTCTTTATCATTTCTATCTCACTCTATTATGATGTGACATTTATTTTACAAACTAATTTGTTTTGAACCTGAAAATAACTTTTTACAACTAAATTGTCGAAAAACAGTGTGTACATGATATAATAAAtatata                       >  2:243008/4‑141 (MQ=255)
                                                                                                      tAGCGGTTCTTTATCATTTCTATCTCACTCTATTATGATGTGACATTTATTTTACAAACTAATTTGTTTTGAACCTGAAAATAACTTTTTACAACTAAATTGTCGAAAAACAGTGTGTACATGATATAATAAATTTATaaa                     >  2:316184/1‑141 (MQ=255)
                                                                                                          ggTTCTTTATCATTTCTATCTCACTCTATTATGATGTGACATTTATTTTACAAACTAATTTGTTTTGAACCTGAAAATAACTTTTTACAACTAAATTGTCGAAAAACAGTGTGTACATGATATAATAAATTTATAAATTGa                 <  2:305746/141‑1 (MQ=255)
                                                                                                          ggTTCTTTATCATTTCTATCTCACTCTATTATGATGTGACATTTATTTTACAAACTAATTTGTTTTGAACCTGAAAATAACTTTTTACAACAAAATAGAAGAAAAACAGTGTGTACATGATATAATAAATTTATAAATTGa                 <  2:131040/141‑1 (MQ=255)
                                                                                                          ggTTCTTTATCATTTCTATCTCACCCTATTATGATGTGACATTTATTTTACAAACTAATTTGTTTTGAACCTGAAAATAACTTTTTACAACTAAATTGTCGAAAAACAGTGTGTACATGATATAATAAATTTATAAATTGa                 <  1:398851/141‑1 (MQ=255)
                                                                                                                  aTCATTTCTATCTCACTCTATTATGATGTGACATTTATTTTACAAACTAATTTGTTTTGAACCTGAAAATAACTTTTTACAACTAAATTGTCGAAAAACAGTGTGTACATGATATAATAAATTTATAAATTGAAAAGAAtt         >  1:279703/1‑141 (MQ=255)
                                                                                                                    taTTTCTATCTCACTCTATTATGATGTGACATTTATTTTACAAACTAATTTGTTTTGAACCTGAAAATAACTTTTTACAACTAAATTGTCGAAAAACAGTGTGTACATGATATAATAAATTTATAAATTGAAAAGAATTCa       <  2:408815/140‑1 (MQ=255)
                                                                                                                         cTATCTCACTCTATTATGATGTGACATTTATTTTACAAACTAATTTGTTTTGAACCTGAAAATAACTTTTTTTTTTTTTTTTTTTGAAAAACAGTGTGAACATAATATAATAAATTTATAAATTGAAAAGAATTCAAAAGa  >  2:394504/1‑141 (MQ=255)
                                                                                                                                  |                                                                                                                                   
AGTAATTTATTGCTTAAGCAACATTTAAGATTTTCAAATTAAAAACGAAGAATTTAAAACATGTGGTGCTAATGTGTAAGAGTCTGAGATATAATAAATTGTTAGCGGTTCTTTATCATTTCTATCTCACTCTATTATGATGTGACATTTATTTTACAAACTAATTTGTTTTGAACCTGAAAATAACTTTTTACAACTAAATTGTCGAAAAACAGTGTGTACATGATATAATAAATTTATAAATTGAAAAGAATTCAAAAGA  >  CP000730/730405‑730666

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