Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 1,929,804 A→G 45.2% intergenic (‑887/‑281) 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,8040AG45.2% ‑1.4 / 46.0 40intergenic (‑887/‑281)USA300HOU_1794/USA300HOU_1796hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (10/12);  new base G (8/10);  total (18/22)
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.01e-01

AATAAAATCGCCTACTTTCTTAAATTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑GT‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCT  >  CP000730/1929665‑1929940
                                                                                                                                           |                                                                                                                                            
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                                               <  2:119294/141‑1 (MQ=255)
                                                                                          tCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑GT‑GATGCTTATTTTCGTTCACATCATAGc                                               <  1:148237/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                           <  1:4168/131‑1 (MQ=255)
                                                                                                             gaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TT‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgtg                           >  2:75261/11‑141 (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)
                                                                                                                            gtatattcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TT‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtattata            >  1:299753/12‑138 (MQ=255)
                                                                                                                            gTATATTCACAGATGGTATAATAATAGCGTCGACTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGCGGTGCATTAC‑T‑TT‑GCTCATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtattata            <  2:19705/141‑7 (MQ=255)
                                                                                                                            gTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TT‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtattatc            >  2:148393/1‑138 (MQ=255)
                                                                                                                               tattcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑TT‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACtt         >  2:221547/9‑141 (MQ=255)
                                                                                                                                attcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑CTTT‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTa        >  2:374983/8‑141 (MQ=255)
                                                                                                                                 ttcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TT‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTAt       >  2:248475/7‑141 (MQ=255)
                                                                                                                                     cagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TT‑GCTAATCATTTTCGTTCACATCAGAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGc   >  1:49341/3‑141 (MQ=255)
                                                                                                                                      cGATAATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGTCGCAGAGAGGAGTTGTA‑TAA‑G‑GT‑GATGCTTATTTTCTTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGc   <  2:307426/141‑1 (MQ=255)
                                                                                                                                      agatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑TT‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCt  >  2:146035/2‑141 (MQ=255)
                                                                                                                                      agatggtATAATAATAGAGTCGCCGATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TT‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCt  >  1:268162/2‑141 (MQ=255)
                                                                                                                                       gatgatATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGTGGTGTA‑TAA‑G‑GT‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCt  <  2:144472/141‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                            
AATAAAATCGCCTACTTTCTTAAATTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑GT‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCT  >  CP000730/1929665‑1929940

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