Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 1,929,804 A→G 40.8% 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,8040AG40.8% 8.1 / 51.0 47intergenic (‑887/‑281)USA300HOU_1794/USA300HOU_1796hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (15/13);  new base G (7/12);  total (22/25)
Fisher's exact test for biased strand distribution p-value = 3.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.23e-01

TAAAATCGCCTACTTTCTTAAATTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTA  >  CP000730/1929667‑1929941
                                                                                                                                         |                                                                                                                                                
taaAATCGCCTACTTTCTTAAATTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGatat                                                                                                                                               <  1:239638/141‑1 (MQ=255)
  aaaTCGCCTACTTTCTTAAATTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATCTAACTCATCTGAAGCAAGTGACTGTGCTTGTATATTTGCCgatgatcgaa                                                                                                                                             >  1:29218/1‑137 (MQ=255)
   aaTCGCCTACTTTCTTAAATTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCCGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATataat                                                                                                                                            >  1:241423/1‑141 (MQ=255)
   aaTCGCCTACTTTCTTAAATTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGATTGTATATTTGCCGataatataat                                                                                                                                            >  1:254164/1‑141 (MQ=255)
     tCGCCTACTTTCTTAAATTAACAATATCTATTCTCCTAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAAtag                                                                                                                                          >  1:155225/1‑141 (MQ=255)
              ttCTTAAATTAACAATATCTATTCTCGTAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGCCTGTGATGGTATATTTGCCGTAGAGATAATGTTagagcggc                                                                                                                                 >  1:122341/1‑137 (MQ=25)
                         aCAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCAAGAATCATATAACTTATCAGAAGCAAGTGACGGTGCTTGTATATTTGCCGATGATATAATAGTAGAGGCGCCTATCTCTCAg                                                                                                                      >  1:351220/1‑141 (MQ=255)
                          cAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAgg                                                                                                                     >  1:257350/1‑141 (MQ=255)
                                          aGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCAGATGAGATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCa                                                                                                     <  1:337362/141‑1 (MQ=255)
                                              ttggTCCAAGTAAGTGTAGACGATTCAATCGATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTAAATTTGCCGATAATATAATAGTAGAGTCGCCTATCTCTCAGGCATCAATTTAGACGCagaga                                                                                                 <  2:224567/138‑1 (MQ=255)
                                               ttGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCagagag                                                                                                <  1:141453/141‑1 (MQ=255)
                                                 gTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGAGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAgcggg                                                                                              >  2:66977/1‑140 (MQ=255)
                                                   ccAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGTCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGaggagg                                                                                            >  2:249357/1‑141 (MQ=255)
                                                                    aTTCAATCTATCAGCTAGAATCGTATAACTGATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGCtt                                                                       <  1:272708/141‑1 (MQ=255)
                                                                    aTTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTCGAGTCGCCTATCTCTCAGTCGTCAATTTAGACGCAGAATGGAGATGTATAAGGTGCTGCtt                                                                       >  2:205052/1‑141 (MQ=255)
                                                                      tCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTAt                                                                     <  1:151651/141‑1 (MQ=255)
                                                                        aaTCCATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAAGAGGTGTATAAGGTGATGCTTAt‑‑tt                                                                 <  1:206777/141‑1 (MQ=255)
                                                                          tCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTAT‑‑TTTc                                                               <  2:257350/141‑1 (MQ=255)
                                                                                 gCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTAAATTTAGACGAAGAGAGGAGGTGTAGNANGNGATGCTTAT‑‑TTTAGTTcaca                                                       >  1:265429/1‑141 (MQ=255)
                                                                                          atatAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTAT‑‑TTTCGTTCACATCATAGCAc                                              <  1:289258/141‑1 (MQ=255)
                                                                                          atatAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGAATAAGGTGATGCTTCT‑‑TTTCGTTCACATCTTAGCAc                                              >  1:65610/1‑141 (MQ=255)
                                                                                             taACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTAT‑‑TTTCGTTCACATCATAGCACCAg                                           <  2:287930/141‑1 (MQ=255)
                                                                                                   aTCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTAT‑‑TTTCGTTCACATCATAGCACCAGtcatca                                     >  2:148616/1‑141 (MQ=255)
                                                                                                    tCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAg                                    >  1:106868/1‑141 (MQ=255)
                                                                                                    tCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAg                                    >  1:134483/1‑141 (MQ=255)
                                                                                                    tCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAg                                    >  1:150063/1‑141 (MQ=255)
                                                                                                       atcggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGc                                <  1:110454/127‑1 (MQ=255)
                                                                                                        tcggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCAT‑‑TTTCGTTCACATCATATCACCAGTCATCAGTGGCt                               >  2:106091/14‑141 (MQ=255)
                                                                                                          ggaagtagataggCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgt                             <  2:14923/129‑1 (MQ=255)
                                                                                                           gaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCAT‑‑TTTCGTTCACATCATAGCACTAGTCATCAGTGGCtgtg                            >  1:346576/11‑141 (MQ=255)
                                                                                                           gaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgtg                            <  1:127911/131‑1 (MQ=255)
                                                                                                            aagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGc                           <  1:5700/132‑1 (MQ=255)
                                                                                                            aagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGc                           >  2:251848/10‑141 (MQ=255)
                                                                                                             aGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGc                           >  2:119275/1‑141 (MQ=255)
                                                                                                             aGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGc                           <  1:57279/141‑1 (MQ=255)
                                                                                                               tagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgt                        <  2:205543/135‑4 (MQ=255)
                                                                                                                 gatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTTTCAGGTGTCAA‑‑TCGACGAAGAGAGGAGGTGCATTACTTTGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgttg                      <  1:67338/137‑6 (MQ=255)
                                                                                                                  atatgCTTGTATATTCACAGATGGTATAATAATAGCGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgttgc                     <  2:184533/138‑7 (MQ=255)
                                                                                                                    ttgCTTGTATTCTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCg                    <  2:245553/140‑1 (MQ=255)
                                                                                                                    atgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GCCGAAGAGAGGAGGTGCATTACTTTGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGt                   >  1:291399/2‑141 (MQ=255)
                                                                                                                    atgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGt                   >  1:189774/2‑141 (MQ=255)
                                                                                                                    atgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGt                   >  1:163342/2‑141 (MQ=255)
                                                                                                                    atgCTTGTATATTCACAGAGGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTTAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGt                   >  1:232664/2‑141 (MQ=255)
                                                                                                                       cTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtatt                <  2:216324/141‑4 (MQ=255)
                                                                                                                       cTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtatt                <  2:275754/141‑4 (MQ=255)
                                                                                                                        ttGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtatta               <  2:173931/141‑5 (MQ=255)
                                                                                                                        ttGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtatta               <  1:251848/141‑5 (MQ=255)
                                                                                                                           tatattcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtattatac            >  1:261721/11‑137 (MQ=255)
                                                                                                                           tataTTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtattatac            <  1:96084/141‑8 (MQ=255)
                                                                                                                             tattcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTT‑‑TGCTAATGGTTACTTATCACATCATAGCACCAGTCATCAGGGGCTGGGCTGTTGCGTATTATACtt          >  1:99374/9‑141 (MQ=255)
                                                                                                                                 cacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTg      >  2:197643/5‑141 (MQ=255)
                                                                                                                                  acagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTTTTGCGTATTATACTTATTgg     <  1:217272/138‑1 (MQ=255)
                                                                                                                                   ccGATGATATAATAGTAGAGGCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTAT‑‑TTTCGTTCACATCATAGCACCCGTTATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTgg     <  2:58816/141‑1 (MQ=255)
                                                                                                                                   cagatggtATAATAATAGAGGCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGc    <  1:175042/139‑1 (MQ=255)
                                                                                                                                    cGAGGATATAATAGGAGAGTCGCCGATGTCTCAGGCGGCAATGTAGGCGCAGAGAGGAGGTGTATAAGGTGATGCTTAT‑‑TTTCGTTCACATCAGAGCACCAGTCATCAGTGGCTGTGCCACTGCGTTTTTTTCTTATTGGc    <  1:270609/141‑1 (MQ=255)
                                                                                                                                    agatggtATAATAAAAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCt   >  2:243199/2‑141 (MQ=255)
                                                                                                                                     gatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCtt  >  1:87064/1‑140 (MQ=255)
                                                                                                                                         |                                                                                                                                                
TAAAATCGCCTACTTTCTTAAATTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTA  >  CP000730/1929667‑1929941

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: