Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 1,929,797 G→A 43.6% intergenic (‑880/‑288) USA300HOU_1794 ← / → USA300HOU_1796 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007301,929,7970GA43.6% 19.9 / 41.2 42intergenic (‑880/‑288)USA300HOU_1794/USA300HOU_1796hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (9/14);  new base A (8/10);  total (18/24)
Fisher's exact test for biased strand distribution p-value = 7.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTGTCATAATAAAATCGCCTACTTTCTTAAATTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑GT‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTT  >  CP000730/1929658‑1929930
                                                                                                                                           |                                                                                                                                         
ttGTCATAATAAAATCGCCTACTTTCTTAAATTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGc                                                                                                                                          <  1:200802/141‑1 (MQ=255)
    cATAATAAAATCGCCTACTTTCTTAAATTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCgat                                                                                                                                      <  1:167988/141‑1 (MQ=255)
       aataaAATCGCCTACTTTCTTAAATTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCgatgat                                                                                                                                   <  1:134826/141‑1 (MQ=255)
                         tttAAATTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCt                                                                                                                 <  2:79189/140‑1 (MQ=255)
                               ttAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTAtctct                                                                                                           >  2:331103/1‑141 (MQ=255)
                                tAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATGGTAGAGTCGCCTGactctc                                                                                                          >  1:181922/1‑141 (MQ=255)
                                  aCAATCTCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAg                                                                                                        <  2:304180/141‑1 (MQ=255)
                                   cAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAgg                                                                                                       >  1:374347/1‑141 (MQ=255)
                                    aaTATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGc                                                                                                      <  1:323341/141‑1 (MQ=255)
                                                cctataaTTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTctaaattgac                                                                                          <  2:333373/136‑10 (MQ=255)
                                                    gactttGTACAATTAAGTGTAGAGGACTCAATCTATAAGCTAGAATCATATAAGTTATCAGAAGCAAGTGACTGTGCTTGTATATTTTACGATGAGCAAAAATNAGAGTCGACGATCTCTTAGTCTTCAAGTTAAACGCag                                                                                      >  2:626/4‑141 (MQ=255)
                                                    gAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGAGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCag                                                                                      >  1:372521/1‑141 (MQ=255)
                                                       tttGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCagaga                                                                                   <  1:82290/141‑1 (MQ=255)
                                                                                 aaTCTATCAGCTAGAATCATATAACTTATCCGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑GT‑GATGCTTAttt                                                     <  1:45392/141‑1 (MQ=255)
                                                                                                 tCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑GT‑GATGCTTATTTTCGTTCACATCATAGc                                     <  1:148237/141‑1 (MQ=255)
                                                                                                 tCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑GT‑GATGCTTATTTTCGTTCACATCATAGc                                     <  2:119294/141‑1 (MQ=255)
                                                                                                    aataACCTATCAGAAGCAACTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑GT‑GATGCTTATTTTCGTTCACATCATAGCAcc                                  <  2:371816/140‑1 (MQ=255)
                                                                                                             tCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑GT‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAg                         <  1:331103/141‑1 (MQ=255)
                                                                                                             tCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGGCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑GT‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAg                         >  1:352187/1‑141 (MQ=255)
                                                                                                                gAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑GT‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCCGTgg                      >  2:228990/1‑141 (MQ=255)
                                                                                                                atcggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TT‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGc                     <  1:258370/127‑1 (MQ=255)
                                                                                                                atcggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGGGCATTAC‑T‑TT‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGc                     >  1:309528/15‑141 (MQ=255)
                                                                                                                 tcggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑TT‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCt                    <  2:159485/128‑1 (MQ=255)
                                                                                                                 tcggaagaagatatgCTTGTATATTCACAGATGGTATAATAATATAGTCGCCTATCTCTCAGGCGTGAATT‑‑GATGACGAGAGGAGGGGCAGTAC‑G‑GT‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTTGcg                    >  2:262988/14‑140 (MQ=255)
                                                                                                                 aaGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑GT‑GATGCTTATTTTCGTTCCCATCATAGCACCAGTCATCAGTGGc                     >  1:286233/1‑141 (MQ=255)
                                                                                                                  agccagcgtcggtgCTTGTCTATTTGCTGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGGCTATTTAGCCGCAGAGAGGAGGTGTA‑TAA‑T‑GT‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCt                    <  1:240800/130‑1 (MQ=255)
                                                                                                                    gaagtagttatgCTTGGATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGCCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑CTTT‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgtg                 <  1:4203/131‑1 (MQ=255)
                                                                                                                    gaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TT‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgtg                 >  2:75261/11‑141 (MQ=255)
                                                                                                                    gaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TT‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgtg                 <  1:4168/131‑1 (MQ=255)
                                                                                                                    gaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCAATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TT‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgtg                 <  1:161606/131‑1 (MQ=255)
                                                                                                                      agtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTTTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TT‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGct               <  2:320757/133‑2 (MQ=255)
                                                                                                                      agtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TT‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGct               <  1:36955/133‑2 (MQ=255)
                                                                                                                       gtagatatggttGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑TT‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctg              <  2:157222/131‑3 (MQ=255)
                                                                                                                       gtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGGA‑TTA‑C‑TTGGCTAATCATTTTCGTTCACATTATAGTACCAGTCATCCGTGGCTGTGctg              >  2:325281/8‑139 (MQ=255)
                                                                                                                       gtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑TT‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctg              >  1:356263/8‑139 (MQ=255)
                                                                                                                       gtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTATCAGGCGTTAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TT‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctg              >  1:254871/8‑139 (MQ=255)
                                                                                                                        tagagatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TT‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgt             <  1:327937/135‑4 (MQ=255)
                                                                                                                        tGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑GT‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGccc              >  2:113043/1‑140 (MQ=255)
                                                                                                                          gatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGCCCATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑CTTT‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCGGTGctgttg           >  2:133121/5‑136 (MQ=255)
                                                                                                                               gCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑GT‑GATGCTTATTTTCGTTCACATCATAGCACCAGTAATAAGTGGCTGTGCCATTGCGtt       >  2:170420/1‑141 (MQ=255)
                                                                                                                                   gTATATTCACAGATGGTATAATAATAGCGTCGACTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGCGGTGCATTAC‑T‑TT‑GCTCATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtattata  <  2:19705/141‑7 (MQ=255)
                                                                                                                                   gTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TT‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtattatc  >  2:148393/1‑138 (MQ=255)
                                                                                                                                           |                                                                                                                                         
TTGTCATAATAAAATCGCCTACTTTCTTAAATTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑GT‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTT  >  CP000730/1929658‑1929930

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: ‑