Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 1,034,902 G→T 12.6% intergenic (‑1319/‑6387) USA300HOU_0983 ← / → USA300HOU_0984 competence protein/glycosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007301,034,9020GT12.6% 126.7 / 5.0 43intergenic (‑1319/‑6387)USA300HOU_0983/USA300HOU_0984competence protein/glycosyltransferase
Reads supporting (aligned to +/- strand):  ref base G (14/23);  new base T (4/2);  total (18/25)
Fisher's exact test for biased strand distribution p-value = 2.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTCAGTGAATTTTTGGAAATTGTGGAAACTGTTACCATCATCATCTGTTATAAAACTAAAGTTTAAATTGCCAGTATCATGATAAACAGCGCCACCACCAGAAATTCTTCTTACTACATCGATGTTGTGAGCATCGATATATGTCTGATTTACTTCCTCTATCGTATTTTGATTCTTTCCAACAATGATAGATGGTCTATTTATGTAAAATAAAAAGTAACTTTCTTCTGCTGGTAAATTTTTTAAAACATATTCTTCCATTGCTAAG  >  CP000730/1034767‑1035034
                                                                                                                                       |                                                                                                                                    
ttCAGTGAATTTTTGGAAATTGTGGAAACTGTTACCATCATCATCTGTTATAAAACTAAAGTTTAAATTGCCAGTATCATGATAAACAGCGCCACCACCAGAAATTCTTCTTACTACATCGATGTTGTGAGCATCGatata                                                                                                                                 <  2:76114/141‑1 (MQ=255)
         tttttGGAAATTGTGGAAACTGTTACCATCATCATCTGTTATAAAACTAAAGTTTAAATTGCCAGTATCATGATAAACAGCGCCACCACCAGAAATTCTTCTTACTACATCGATGTTGTGAGCATCGATATATGTCTGAtt                                                                                                                        >  2:125254/1‑141 (MQ=255)
         tttttGGAAATTGTGGAAACTGTTACCATCATCATCTGTTATAAAACTAAAGTTTAAATTGCCAGTATCATGATAAACAGCGCCACAACCAGAAATTCTTCTTACTACATCGATGTTGTGAGCATCGATATATGTCTGAtt                                                                                                                        <  1:436399/141‑1 (MQ=255)
           tttGGAAATTGTGGAAACTGTTACCATCATCATCTGTTATAAAACTAAAGTTTAAATTGCCAGTATCATGATAAACAGCGCCACCACCAGAAATTCTTCTTACTACATCGATGTTGTGAGCATCGATATATGTCTGATTTa                                                                                                                      <  1:325503/141‑1 (MQ=255)
             tGGAAATTGTGGAAACTGTTACCATCATCATCTGTTATAAAACTAAAGTTTAAATTGCCAGTATCATGATAAACAGCGCCACCACCAGAAATTCTTCTTACTACATCGATGTTGTGAGCATCGATATATGTCTGATTTACt                                                                                                                    <  2:320195/141‑1 (MQ=255)
                  aTTGTGGAAACTGTTACCATCATCATCTGTTATAAAACTAAAGTTTAAATTGCCAGTATCATGATAAACAGCGCCACCACCAGAAATTCTTCTTACTACATCGATGTTGTGAGCATCGATATATGTCTGATTTACTTCctc                                                                                                               <  2:322342/141‑1 (MQ=255)
                        gAAACTGTTACCATCATCATCTGTTATACAACTAAAGTTTAAATTGCCAGTATCATGATAAACAGCGCCACCACCAGAAATTCTTCTTACTACATCGATGTTGTGAGCATCGATATATGTCTGATTTACTTCCTCTATCGt                                                                                                         <  1:20225/141‑1 (MQ=255)
                        gAAACTGTTACCATCATCATCTGTTATAATACTAAAGTTTAAATTGCCAGTATCATGATAAACAGCGCCACCACCAGAAATTCTTCTTACTACATCGATGTTGTGAGCATCGATATATGTCTGATTTACTTCCTCTATCGt                                                                                                         <  1:20223/141‑1 (MQ=255)
                        gAAACTGTTACCATCATCATCTGTTATAAAACTAAAGTTTAAATTGCCAGTATCATGATAAACAGCGCCACCACCAGAAATTCTTCTTACTACATCGATGTTGTGAGCATCGATATATGTCTGATTTACTTCCTCTATCGt                                                                                                         <  2:45042/141‑1 (MQ=255)
                           aCTGTTACCATCATCATCTGTTATAAAACTAAAGTTTAAATTGCCAGTATCATGATAAACAGCGCCACCACCAGAAGTTCTTCTTACTACATCGATGTTGTGAGCATCGATATATGTCTGATTTACTTCCTCTATCGTAtt                                                                                                      >  1:99375/1‑141 (MQ=255)
                             tGTTACCATCATCATCTGTTATAAAACTAAAGTTTAAATTGCCAGTATCATGATAAACAGCGCCACCACCAGAAATTCTTCTAACTACATCGATGTTATGAGCATCTATATATGTCTGATTTACTTCCTCTATCGTAtttt                                                                                                    <  1:128717/141‑1 (MQ=255)
                                              gTTATAAAACTAAAGTTTAAATTGCCAGTATCATGATAAACAGCGCCACCACCAGAAATTCTTCTTACTACATCGATGTTGTGAGCATCGATATATGTCTGATTTACTTCCTCTATCGTATTTTGATTCTTTCCAACAATg                                                                                   >  1:240211/1‑141 (MQ=255)
                                                    aaaCTAAAGTTTAAATTGCCAGTATCATGATAAACAGCGCCACCACCAGAAATTCTTCTTACTACATCGATGTTGTGAGCATCGATATATGTCTGATTTACTTCCTCTATCGTATTTTGATTCTTTCCAACAATGATAGAt                                                                             <  1:429/141‑1 (MQ=255)
                                                          aaGTTTAAATTGCCAGTATCATGATAAACAGCGCCACCACCAGAAATTCTTCTTACTACATCGATGTTGTGAGCATCGATATATGTCTGATTTACTTCCTCTATCGTATTTTGATTCTTTCCAACAATGATAGATGGTCTa                                                                       >  1:121171/1‑141 (MQ=255)
                                                                aaaTTGCCAGTATCATGATAAACAGCGCCACCACCAGAAATTCTTCTTACTACATCGATGTTGTGAGCATCGATATATGTCTGATTTACTTCCTCTATCGTATTTTGATTCTTTCCAACAATGATAGATGGTCTATTTATg                                                                 <  1:190337/141‑1 (MQ=255)
                                                                 aaTTGCCAGTATCATGATAAACAGCGCCACTACCAGAAATTCTTCTTACTACTTCTATGTTGTGAGCATCGATATATGTCTGATTTACTTTCTCTATCGTATTTTGATTCTTTCCAACTATGATAGATAGTCTATTTATGt                                                                <  2:274331/141‑1 (MQ=255)
                                                                 aaTTGCCAGTATCATGATAAACAGCGCCACCACCAGAAATTCTTCTTACTACATCGATGTTGTGAGCATCGATATATGTCTGATTTACTTCCTCTATCGTATTTTGATTCTTTCCAACAATGATAGATGGTCTATTTATGt                                                                >  1:323754/1‑141 (MQ=255)
                                                                     gCCAGTATCATGATAAACAGCGCCACCACCAGAAATTCTTCTTACTACATCGATGTTGTGAGCATCGATATATGTCTGATTTACTTCCTCTATCGTATGTTGATTCTTTCCAACAATGATAGATGGTCTATTTATGtaaaa                                                            <  2:31810/141‑1 (MQ=255)
                                                                             cATGATAAACAGCGCCACCACCAGAAATTCTTCTTACTACATCGATGTTGTGAGCATCGATATATGTCTGATTTACTTCCTCTATCGTATTTTGATTCTTTCCAACAATGATAGATGGTCTATTTATGTAAAATAAAAAGt                                                    <  2:26581/141‑1 (MQ=255)
                                                                             cATGATAAACAGCGCCACCACCAGAAATTCTTCTTACTACATCGATGTTGTGAGCATCGATATATGTCTGATTTACTTCCTCTATCGTATTTTGATTCTTTCCAACAATGATAGATGGTCTATTTATGTAAAATAAAAAGt                                                    <  2:408751/141‑1 (MQ=255)
                                                                                 aTAAACAGCGCCACCACCAGAAATTCTTCTTACTACATCGATGATGTGAGCATCGATATATGTCTGATTTACTTCCTCTATCGTATTTTGATTCTTTCCAACAATGATAGATGGTCTATTTATGTAAAATAAAAAGTAACt                                                <  1:64665/141‑1 (MQ=255)
                                                                                   aaaCAGCGCCACCACCAGAAATTCTTCTTACTACATCGATGTTGTGAGCATCGATATATGTCTGATTTACTTCCTCTATCGTATTTTGATTCTTTCCAACAATGATAGATGGTCTATTTATGTAAAATAAAAAGTAACttt                                              <  2:438448/141‑1 (MQ=255)
                                                                                   aaaCAGCGCCACCACCAGAAATTCTTCTTACTACATCGATGTTGTGAGCATCGATATATGTCTGATTTACTTCCTCTATCGTATTTTGATTCTTTCCAACAATGATAGATGGTCTATTTATGTAAAATAAAAAGTAACttt                                              >  2:80803/1‑141 (MQ=255)
                                                                                   aaaCAGCGCCACCACCAGAAATTCTTCTTACTACATCGATGTTGTGAGCATCGATATATGTCTGATTTACTTCCTCTATCGTATTTTGATTCTTTCCAACAATGATAGATGGTCTATTTATGTAAAATAAAAAGTAACttt                                              >  2:80812/1‑141 (MQ=255)
                                                                                   aaaCAGCGCCACCAACAGAAATTCTTCTTACTACATCGATGTTGTGAGCATCGATATATGTCTGATTTACTTCCTCTATCGTATTTTGATTCTTTCCAACAATGATAGATGGTCTATTTATGTAAAATAAAAAGTAACttt                                              <  2:438446/141‑1 (MQ=255)
                                                                                      cAGCGCCACCACCAGAAATTCTTCTTACTACATCGATGTTGTGAGCATCGATATATGTCTGATTTACTTCCTCTATCGTATTTTGATTCTTTCCAACAATGATAGATGGTCTATTTATGTAAAATAAAAAGTAACTttctt                                           <  1:290547/141‑1 (MQ=255)
                                                                                        gcgcCACCACCAGAAATTCTTCTTACTACATCGATGTTGTGAGCATCGATATATGTCTGATTTACTTCCTCTATCGTATTTTGATTCTTTCCAACAATGATAGATGGTCTATTTATGTAAAATAAAAAGTAACTTtcttct                                         >  1:61616/1‑141 (MQ=255)
                                                                                            caccaccaGAAATGCTTCTAACTACATCGATGTTATGAGCATCTATATATGTCTGATTTACTTCCTCTATCGTATTTTGATTCTTTCCAACAATGATAGATGGTCTATTTATGTAAAATACAAAGTAACTTTCTTctgctg                                     >  2:455876/1‑141 (MQ=255)
                                                                                            caccaccaGAAATGCTTCTAACTACATCGATGTTATGAACATCTATATATGTCTGATTTACTTCCTCTATCGTATTTTGATTCTTTCCAACAATGATAGATGGTCTATTTATGTAAAATAGAAAGTAACTTTCTTctgctt                                     >  2:455982/1‑140 (MQ=255)
                                                                                             accaccaGAAATTCTTCTTACTACATCGATGTTGTGAGCATCGATATATGTCTGATTTACTTCCTCTATCGTATTTTGATTCTTTCCAACAATGATAGATGGTCTATTTATGTAAAATAAAAAGTAACTTTCTTCTGCTgg                                    <  1:231155/141‑1 (MQ=255)
                                                                                              ccaccaGAAATTCTTCTTACTACATCGATGTTGTGAGCATCTATATATGTCTGATTTACTTCCTCTATCGTATTTTGATTCTTTCCAACAATGATAGATGGTCTATTTATGTAAAATAAAAAGTAACTTTCTTCTGCTGGt                                   >  1:447527/1‑141 (MQ=255)
                                                                                                  caGAAATTCTTCTTACTACATCGATGTTGTGAGCATCGATATATGTCTGATTTACTTCCTCTATCGTATTTTGATTCTTTCCAACAATGATAGATGGTCTATTTATGTAAAATAAAAAGTAACTTTCTTCTGCTGGTAAAt                               >  1:87225/1‑141 (MQ=255)
                                                                                                     aaaTTCTTCTTACTACATCGATGTTGTGAGCATCGATATATGTCTGATTTACTTCCTCTATCGTATTTTGATTCTTTCCAACAATGATAGATGGTCTATTTATGTAAAATAAAAAGTAACTTTCTTCTGCTGGTAAAtttt                            >  2:133966/1‑141 (MQ=255)
                                                                                                      aaTTCTTCTTACTACATCGATGTTGTGAGCATCGATATATGTCTGATTTACTTCCTCTATCGTATTTTGATTCTTTCCAACAATGATAGATGGTCTATTTATGTAAAATAAAAAGTAACTTTCTTCTGCTAGTAAAttttt                           <  1:202757/141‑1 (MQ=255)
                                                                                                         tcttctTACTACATCGATGTTGTGAGCATCGATATATGTCTGATTTACTTCCTCTATCGTATTTTGATTCTTTCCAACAATGATAGATGGTCTATTTATGTAAAATAAAAAGTAACTTTCTTCTGCTGGTAAATTTTTTaa                        >  1:191896/1‑141 (MQ=255)
                                                                                                          cttcttACTACATCGATGTTGTGAGCATCGATATATGTCTGATTTACTTCCTCTATCGTATTTTGATTCTTTCCAACAATGATAGATGGTCTATTTATGTAAAATAAAAAGTAACTTTCTTCTGCTGGTAAATTTTTTaaa                       >  2:221275/1‑141 (MQ=255)
                                                                                                          cttctaACTACATCGATGTTATGAGCATCTATATATGTCTGATTTACTTCCTCTATCGTATTTTGATTCTTTCca                                                                                         >  1:104032/1‑75 (MQ=255)
                                                                                                          cttctaACTACATCGATGTTATGAGCATCTATATATGTCTGATTTACTTCCTCTATCGTATTTTGATTCTTTCca                                                                                         <  2:104032/75‑1 (MQ=255)
                                                                                                                        gATGTTGTGAGCATCGATATATGTCTGATTTACTTCCTCTATCGTATTTTGATTCTTTCCAACAATGATAGATGGTCTATTTATGTAAAATAAAAAGTAACTTTCTTCTGCTGGTAAATTTTTTAAAACATATTCTTCCAt         >  1:194190/1‑141 (MQ=255)
                                                                                                                         aTGTTGTGAGCATCGATATATGTCTGATTTACTTCCTCTATCGTATTTTGATTCTTTCCAACAATGATAGATGGTCTATTTATGTAAAATAAAAAGTAACTTTCTTCTGCTGGTAAATTTTTTAAAACATATTCTTCCAtt        <  1:64679/141‑1 (MQ=255)
                                                                                                                            ttgtGAGCATCGATATATGTCTGATTTACTTCCTCTATCGTATTTTGATTCTTTCCAACAATGATAGATGGTCTATTTATGTAAAATAAAAAGTAACTTTCTTCTGCTGGTAAATTTTTTAAAACATATTCTTCCATTGCt     <  1:417770/141‑1 (MQ=255)
                                                                                                                              gtgAGCATCGATATATGTCTGATTTACTTCCTCTATCGTATTTTGATTCTTTCCAACAATGATAGATGGTCTATTTATGTAAAATAAAAAGTAACTTTCTTCTGCTGGTAAATTTTTTAAAACATATTCTTCCATTGCTaa   >  2:258541/1‑141 (MQ=255)
                                                                                                                               tgAGTATCGATATATGTCTGATTTACTTCCTCTATCGTATTTTGATTCTTTCCACCAATGATAGATGGACTATTTATGTAAAATAAAAAGTAACTTTCTTCTGCTGGTAAATTTTTTAAAACATATTCTTCCATTGCTAAg  <  2:18547/141‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                    
TTCAGTGAATTTTTGGAAATTGTGGAAACTGTTACCATCATCATCTGTTATAAAACTAAAGTTTAAATTGCCAGTATCATGATAAACAGCGCCACCACCAGAAATTCTTCTTACTACATCGATGTTGTGAGCATCGATATATGTCTGATTTACTTCCTCTATCGTATTTTGATTCTTTCCAACAATGATAGATGGTCTATTTATGTAAAATAAAAAGTAACTTTCTTCTGCTGGTAAATTTTTTAAAACATATTCTTCCATTGCTAAG  >  CP000730/1034767‑1035034

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