Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 1,034,893 G→A 12.2% intergenic (‑1310/‑6396) 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,8930GA12.2% 104.8 / 5.2 41intergenic (‑1310/‑6396)USA300HOU_0983/USA300HOU_0984competence protein/glycosyltransferase
Reads supporting (aligned to +/- strand):  ref base G (15/21);  new base A (3/2);  total (18/23)
Fisher's exact test for biased strand distribution p-value = 6.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.97e-01

ACAATTGGTTCAGTGAATTTTTGGAAATTGTGGAAACTGTTACCATCATCATCTGTTATAAAACTAAAGTTTAAATTGCCAGTATCATGATAAACAGCGCCACCACCAGAAATTCTTCTTACTACATCGATGTTGTGAGCATCGATATATGTCTGATTTACTTCCTCTATCGTATTTTGATTCTTTCCAACAATGATAGATGGTCTATTTATGTAAAATAAAAAGTAACTTTCTTCTGCTGGTAAATTTTTTAAAACATATTCTTCCATTGCTAA  >  CP000730/1034759‑1035033
                                                                                                                                      |                                                                                                                                            
aCAATTGGTTCAGTGAATTTTTGGAAATTGTGGAAACTGTTACCATCATCATCTGTTATAAAACTAAAGTTTAAATTGCCAGTATCATGATAAACAGCGCCACCACCAGAAATTCTTCTTACTACATCGATGTTGTGAGCa                                                                                                                                        >  2:301007/1‑141 (MQ=255)
        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:408751/141‑1 (MQ=255)
                                                                                     cATGATAAACAGCGCCACCACCAGAAATTCTTCTTACTACATCGATGTTGTGAGCATCGATATATGTCTGATTTACTTCCTCTATCGTATTTTGATTCTTTCCAACAATGATAGATGGTCTATTTATGTAAAATAAAAAGt                                                   <  2:26581/141‑1 (MQ=255)
                                                                                         aTAAACAGCGCCACCACCAGAAATTCTTCTTACTACATCGATGATGTGAGCATCGATATATGTCTGATTTACTTCCTCTATCGTATTTTGATTCTTTCCAACAATGATAGATGGTCTATTTATGTAAAATAAAAAGTAACt                                               <  1:64665/141‑1 (MQ=255)
                                                                                           aaaCAGCGCCACCACCAGAAATTCTTCTTACTACATCGATGTTGTGAGCATCGATATATGTCTGATTTACTTCCTCTATCGTATTTTGATTCTTTCCAACAATGATAGATGGTCTATTTATGTAAAATAAAAAGTAACttt                                             >  2:80812/1‑141 (MQ=255)
                                                                                           aaaCAGCGCCACCACCAGAAATTCTTCTTACTACATCGATGTTGTGAGCATCGATATATGTCTGATTTACTTCCTCTATCGTATTTTGATTCTTTCCAACAATGATAGATGGTCTATTTATGTAAAATAAAAAGTAACttt                                             >  2:80803/1‑141 (MQ=255)
                                                                                           aaaCAGCGCCACCACCAGAAATTCTTCTTACTACATCGATGTTGTGAGCATCGATATATGTCTGATTTACTTCCTCTATCGTATTTTGATTCTTTCCAACAATGATAGATGGTCTATTTATGTAAAATAAAAAGTAACttt                                             <  2:438448/141‑1 (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)
                                                                                                                                      |                                                                                                                                            
ACAATTGGTTCAGTGAATTTTTGGAAATTGTGGAAACTGTTACCATCATCATCTGTTATAAAACTAAAGTTTAAATTGCCAGTATCATGATAAACAGCGCCACCACCAGAAATTCTTCTTACTACATCGATGTTGTGAGCATCGATATATGTCTGATTTACTTCCTCTATCGTATTTTGATTCTTTCCAACAATGATAGATGGTCTATTTATGTAAAATAAAAAGTAACTTTCTTCTGCTGGTAAATTTTTTAAAACATATTCTTCCATTGCTAA  >  CP000730/1034759‑1035033

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: