Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 1,929,946 G→A 31.1% intergenic (‑1029/‑139) 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,9460GA31.1% 48.6 / 22.8 41intergenic (‑1029/‑139)USA300HOU_1794/USA300HOU_1796hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (14/14);  new base A (6/7);  total (20/21)
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.78e-01

AATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCTAATTTTTCTTTTTGTAAAATAT  >  CP000730/1929808‑1930086
                                                                                                                                            |                                                                                                                                            
aataGTAGAGGCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCAACAGTCCTCAGTGGCTGTGCTGTTGCGTATTATACTTATCGGCTTTGTAAACg                                                                                                                                           >  2:349105/1‑141 (MQ=255)
 ataataGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGc                                                                                                                                          <  1:48139/141‑1 (MQ=255)
  taatagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa                                                                                                                                         <  2:182433/138‑1 (MQ=255)
   aatagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaa                                                                                                                                        <  1:73059/139‑1 (MQ=255)
    gtagAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCaa                                                                                                                                        <  2:372521/141‑1 (MQ=255)
       gagTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAtgc                                                                                                                                     >  2:109714/1‑139 (MQ=255)
         gTCGCCTATCGCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaata                                                                                                                                  >  2:180323/1‑135 (MQ=255)
          tCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGccacaactaa                                                                                                                                 >  2:26099/1‑132 (MQ=255)
           cggcTATCTCTCAGGCGTAAAATTAGACTCAGAGAGGTGGTGTATAAGGTGATGCTTACTTTCGTTCACATCATAGCACCAGTCATCAGTGGCGGTGCCATTGCGTGTTTTTCTTATTGGCTAAGTAGACGCAATACaaaa                                                                                                                                 <  2:86103/138‑1 (MQ=255)
            gCCTATCTCTCAGGCTTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatt                                                                                                                               <  2:53395/141‑10 (MQ=255)
            gCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatt                                                                                                                               <  2:194111/141‑10 (MQ=255)
             ccTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAAtg                                                                                                                               >  1:100377/1‑140 (MQ=255)
              cTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaa                                                                                                                             >  1:50601/1‑130 (MQ=255)
               tATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAgg                                                                                                                             >  1:337437/1‑141 (MQ=255)
               tATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaat                                                                                                                            >  1:37091/1‑129 (MQ=255)
                aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaata                                                                                                                           <  2:336622/141‑14 (MQ=255)
                aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaata                                                                                                                           >  1:102414/1‑128 (MQ=255)
                 tctctcAGGCGTCAATT‑‑GACGAAGAGAGTAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTTGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaatag                                                                                                                          <  2:48917/141‑15 (MQ=255)
                          cGTCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAACTCGGTGACatttcgc                                                                                                                  >  2:239887/1‑136 (MQ=255)
                          cGTCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGc                                                                                                                  <  1:298543/141‑1 (MQ=255)
                                                  gtgtATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACt                                                                                          >  1:238469/1‑141 (MQ=255)
                                                         aaGGTGATGCTTATTTTCGTTCACATCATTGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCa                                                                                    <  2:181922/141‑1 (MQ=255)
                                                          aGGTGATGCTTATTTTCGTTAACATCATAGCACTAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCaa                                                                                   <  1:15621/141‑1 (MQ=255)
                                                           ggTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTatttgcgg                                                                                  <  1:90535/141‑8 (MQ=255)
                                                            gTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAAt                                                                                 >  2:23976/1‑141 (MQ=255)
                                                               gATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAg                                                                               >  1:313288/1‑141 (MQ=255)
                                                                  gCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGa                                                                            <  1:315181/141‑1 (MQ=255)
                                                                  gCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAAAAGTGa                                                                            >  2:7877/1‑141 (MQ=255)
                                                                     tATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAggg                                                                         >  2:160267/1‑141 (MQ=255)
                                                                      aTTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATGGGTTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAgggg                                                                        <  1:251299/141‑1 (MQ=255)
                                                                        tttCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGAt                                                                      <  1:68082/141‑1 (MQ=255)
                                                                         ttCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGAtt                                                                     >  1:204214/1‑141 (MQ=255)
                                                                            gTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTGTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGgt                                                                  <  2:94523/141‑1 (MQ=255)
                                                                               cacaTCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGta                                                               >  2:21745/1‑141 (MQ=255)
                                                                                acaTCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGtat                                                              <  1:170420/141‑1 (MQ=255)
                                                                                 ttttttAGCACCAGTCATCAGTGGCTGTGCAATTGCTTTTTTTTGTTATTGGCTGAGTAGACTCCCTACTAATTTTTTTTTTTTTTTCCGCACCAATAACGATCCGCTCTCTACCGCAAATCGTTAGGGGATTTGTGtatg                                                             >  2:31876/6‑140 (MQ=255)
                                                                                     atAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAACAGTGAGGAGACTGGTGTATAAGTa                                                         >  2:129588/1‑141 (MQ=255)
                                                                                           ccAGTCATCAGTGGCTGGGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGGGAGGGGATTGGTGTATAAGTAAATACt                                                   >  1:198418/1‑141 (MQ=255)
                                                                                                 atcaGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTAtttt                                             <  2:352187/141‑1 (MQ=255)
                                                                                                                         ttttCTTAGTGGCTAAGTAGTCGGAATACAGAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCTAAt                     <  1:67403/141‑1 (MQ=255)
                                                                                                                                         gTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCTAATTTTTCTTTTTGTaaaa     <  1:109714/141‑1 (MQ=255)
                                                                                                                                            gACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAAAAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCTAATTTTTCTTTTTGTAAAatat  >  1:182525/1‑141 (MQ=255)
                                                                                                                                            |                                                                                                                                            
AATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCTAATTTTTCTTTTTGTAAAATAT  >  CP000730/1929808‑1930086

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