Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 1,929,812 G→A 46.2% intergenic (‑895/‑273) 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,8120GA46.2% 11.3 / 47.4 39intergenic (‑895/‑273)USA300HOU_1794/USA300HOU_1796hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (10/11);  new base A (10/8);  total (20/19)
Fisher's exact test for biased strand distribution p-value = 7.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.63e-01

GCCTACTTTCTTAAATTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑GTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACG  >  CP000730/1929674‑1929949
                                                                                                                                          |                                                                                                                                            
gCCTACTTTCTTAAATTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAAtagta                                                                                                                                            <  1:130671/141‑1 (MQ=255)
gCCTACTTTCTTAAATTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAAtagta                                                                                                                                            <  1:70890/141‑1 (MQ=255)
              aTTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTAtctc                                                                                                                              >  2:6300/1‑141 (MQ=255)
                    aaTATCTATTCTCATAGATTTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGc                                                                                                                        <  1:224594/141‑1 (MQ=255)
                    aaTATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCTATCTAGCAGCTAGAATCATATAACTTATCAGGAGCAACTGACTGTGCTTGTATAGTTGCCGATGATAGAATAGTAGAGGCGCCTATCTCTCAGGc                                                                                                                        <  1:23639/141‑1 (MQ=38)
                    aaTATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGc                                                                                                                        <  1:52326/141‑1 (MQ=255)
                      tatCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGt                                                                                                                      <  2:88820/141‑1 (MQ=255)
                      tatCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGt                                                                                                                      <  2:254269/141‑1 (MQ=255)
                         cTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCaa                                                                                                                   >  1:17149/1‑141 (MQ=255)
                         cTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTAGTATATTTGCCGAAGATATCATAGTAGAGTCGCCTATCTCTCAGGCGTCaa                                                                                                                   >  1:206243/1‑141 (MQ=255)
                           aTTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAAtt                                                                                                                 <  1:47824/141‑1 (MQ=255)
                           aTTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAAtt                                                                                                                 <  1:182440/141‑1 (MQ=255)
                                   aGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAGGCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGGCGCCTATCTCTCAGGCGTCAATTTAGACGCa                                                                                                         >  1:92591/1‑141 (MQ=255)
                                                              ttCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑GTGATGCTTa                                                                           >  2:168695/1‑141 (MQ=255)
                                                              ttCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑GTGATGCTTa                                                                           >  2:168798/1‑141 (MQ=255)
                                                               tCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑GTGATGCTTAt                                                                          >  1:168045/1‑141 (MQ=255)
                                                                             aGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑GTGATGCTTATTTTCGTTCAcatca                                                            <  1:250419/141‑1 (MQ=255)
                                                                             aGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑GTGATGCTTATTTTCGTTCAcatca                                                            <  1:228917/141‑1 (MQ=255)
                                                                              gAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑GTGATGCTTATTTTCGTTCAcatcat                                                           >  1:149889/1‑141 (MQ=255)
                                                                                            aTCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑GTGATGCTTATTTTCGTTCACATCATAGCACCAGtcatca                                             >  2:30236/1‑141 (MQ=255)
                                                                                            aTCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑GTGATGCTTATTTTCGTTCACATCATAGCACCAGtcatca                                             <  1:40376/141‑1 (MQ=255)
                                                                                                atcggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGc                                        <  2:123401/127‑1 (MQ=255)
                                                                                                atcggaagtaagtatgCTTGTATATTCACAGAGGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGc                                        <  2:174121/127‑1 (MQ=255)
                                                                                                 tcggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCt                                       <  1:90029/128‑1 (MQ=255)
                                                                                                 aaGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑GTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGc                                        >  1:219071/1‑141 (MQ=255)
                                                                                                      agtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGct                                  >  1:136894/9‑140 (MQ=255)
                                                                                                      agtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑CTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGct                                  >  1:110942/9‑140 (MQ=255)
                                                                                                       gtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctg                                 <  2:39487/134‑3 (MQ=255)
                                                                                                            tatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgttgcg                            <  2:113243/139‑8 (MQ=255)
                                                                                                                 ttGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtatta                       >  1:45114/1‑140 (MQ=255)
                                                                                                                 ttGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑C‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtatta                       >  2:203814/1‑137 (MQ=255)
                                                                                                                 ttGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtatta                       >  1:214784/1‑137 (MQ=255)
                                                                                                                   gTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtattata                     >  1:131453/1‑138 (MQ=255)
                                                                                                                         ttGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑GTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTAt                <  1:107900/141‑1 (MQ=255)
                                                                                                                          cacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTg              >  1:174630/5‑141 (MQ=255)
                                                                                                                          cacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTg              >  2:3331/5‑141 (MQ=255)
                                                                                                                              gatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCtt          <  2:134521/141‑2 (MQ=255)
                                                                                                                                tggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCttag        <  1:108433/138‑4 (MQ=255)
                                                                                                                                tggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGCTGCATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCttag        <  2:207811/138‑4 (MQ=255)
                                                                                                                                  gtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTa      >  1:88667/2‑141 (MQ=255)
                                                                                                                                   tataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaa     >  1:119775/1‑140 (MQ=255)
                                                                                                                                      aataatAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACg  >  1:4037/1‑141 (MQ=255)
                                                                                                                                          |                                                                                                                                            
GCCTACTTTCTTAAATTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑GTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACG  >  CP000730/1929674‑1929949

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