Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 1,682,458 A→G 11.5% intergenic (‑203/+101) phoH ← / ← USA300HOU_1573 phosphate starvation‑inducible protein PhoH/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007301,682,4580AG11.5% 78.1 / 4.9 35intergenic (‑203/+101)phoH/USA300HOU_1573phosphate starvation‑inducible protein PhoH/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (14/17);  new base G (2/2);  total (16/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.88e-01

ATCATGCAATTGATTGATTTAAACATCTTTGCATTGACTTACTCATTTTTCTTTAATTTAAAATAATCATTTATACTTTTCACCTTGTCATAAAAGCAGATTCAATTTTAGACTAATCAATTTTTTGCTGTAATTTAAATAACTAATATTTCTCGAAAAATTCTCTTTGTTTTTTAATTTTTCTTAACGTAAAAAAACAAGCAACTACCATTACGATAGTGCTTGTTTTCAAAAAATTTTATAATTGTTTTGGTTTAGCTAAAATT  >  CP000730/1682321‑1682586
                                                                                                                                         |                                                                                                                                
atCATGCAATTGATTGCTTTAAACATCTTTGCATTGACTTACTCATTTTTCTTTAATTTAAAATAATCATTTATACTTTTCACCTTGTCATAAAAGCAGATTCAATTTTAGACTAATCAATTTTTTTTTGTTATTTAATTa                                                                                                                               >  2:274343/1‑141 (MQ=255)
         ttgattgaTTTAAACATCTTTGCATTGACTTACTCATTTTTCTTTAATTTAAAATAATCATTTATACTTTTCACCTTGTCATAAAAGCAGATTCAATTTTAGACTAATCAATTTTTTGCTGTAATTTAGATAACTAATAtt                                                                                                                      >  1:42199/1‑141 (MQ=255)
         ttgattgaTTTAAACATCTTTGCATTGACTTACTCATTTTTCTTTAATTTAAAATAATCATTTATACTTTTCACCTTGTCATAAAAGCAGATTCAATTTTAGACTAATCAATTTTTTGCTGTAATTTAGATAACTAATAtt                                                                                                                      >  1:261406/1‑141 (MQ=255)
           gattgattTAAACATCTTTGCATTGACTTACTCATTTTTCTTTAATTTAAAATAATCATTTATACTTTTCACCTTGTCATAAAAGCAGATTCAATTTTAGACTAATCAATTTTTTGCTGTAATTTAAATAACTAATATTtc                                                                                                                    <  1:247924/141‑1 (MQ=255)
            attgattTAAACATCTTTGCATTGACTTACTCATTTTTCTTTAATTTAAAATAATCAATTATACTTTTCACCTTGTCATAAAAGCAGATTCAATTTTAGACTAATCAATTTTTTGCTGTAATTTAAATAACTaaaaaaaaa                                                                                                                   <  2:255496/141‑8 (MQ=255)
              tgattTAAACATCTTTGCATTGACTTACTCATTTTTCTTTAATTTAAAATAATCATTTATACTTTTCACCTTGTCATAAAAGCAGATTCAATTTTAGACTAATCAATTTTTTGCTGTAATTTAGATAACTAATATTTCTCg                                                                                                                 <  2:42199/141‑1 (MQ=255)
              tgattTAAACATCTTTGCATTGACTTACTCATTTTTCTTTAATTTAAAATAATCATTTATACTTTTCACCTTGTCATAAAAGCAGATTCAATTTTAGACTAATCAATTTTTTGCTGTAATTTAGATAACTAATATTTCTCg                                                                                                                 <  2:261406/141‑1 (MQ=255)
                    aaaCATCTTTGCATTGACTTACTCATTTTTCTTTAATTTAAAATAATCATTTATACTTTTCACCTTGTCATAAAAGCAGATTCAATTTTAGACTAATCAATTTTTTGCTGTAATTTAAATAACTAATATTTCTCGAAAAAt                                                                                                           >  1:222698/1‑141 (MQ=255)
                                                 tCTTTAATTTAAAATAATCATTTATACTTTTCACCTTGTCATAAAAGCAGATTCAATTTTAGACTAATCAATTTTTTGCTGTAATTTAAATAACTAATATTTCTCGAAAAATTCTCTTTGTTTTTTAATTTTTCTTAACGt                                                                              >  2:17536/1‑141 (MQ=255)
                                                   tttaatttaaAAAAATCATTTACACTTTTCACCTTGTCATAAAATCAGATTCAATTTTAGACTAATCAATTTTTTGCTGTAATTTAAATAACTAATATTTCTCGAAAAATTCTCTTTGTTTTTTAATTTTTCTTAACGTaa                                                                            <  2:9928/141‑1 (MQ=255)
                                                       atttaaAATAATCATTTATGCTTTGCACCTTGTCTTAAAAGCAGATTCAATTTTAGACTAATCAATTTTTTGCTGTAATTTAAATAACTAATATTTCTCGAAAAATTCTCTTTGTTGTTTAATTTTTCTTAACGTaaaaaa                                                                        <  2:133746/141‑1 (MQ=255)
                                                        tttagAATAATCATTTATACTTTTCCCCTTGTCATAAAAGCAGAGTCAATTTTAGGCTAATCAATTTTTTGCTGTAATTTAAATAACCAATATGTCTCGAAAAATTCTCGTTGTTTTTTAATTTTACTTTAAGTaaaaaaa                                                                       <  2:249259/141‑1 (MQ=255)
                                                                 atCATTTATACTTTTCACCTTGTCATAAAAGCAGATTCAATTTTAGACTAATCAATTTTTTGCTGTAATTTAAATAACTAATATTTCTCGAAAAATTCTCTTTGTTTTTTAATTTTTCTTAACGTAAAAAAACAAGCAACt                                                              <  2:222698/141‑1 (MQ=255)
                                                                            ttttCACCTTGTCATAAAAGCAGATTCAATTTTAGACTAATCAATTTTTTGCTGTAATTTAAATAACTAATATTTCTCGAAAAATTCTCTTTGTTTTTTAATTTTTCTTAACGTAAAAAAACAAGCAACTACCATTACGAt                                                   <  2:222842/141‑1 (MQ=255)
                                                                             tttCACCTTGTCATAAAAGCAGATTCAATTTTAGACTAATCAATTTTTTTCTTTTATTTAAATAACTAATATTTCTCGAACAATTCTCTTTGTTTTTTAATTTTTCTTAAGGTAAAAAACCAAGCAACTCCCAtttctttt                                                  >  2:207103/1‑135 (MQ=255)
                                                                             tttCACCTTGTCATAAAAGCAGATTCAATTTTAGACTAATCAATTTTTTGCTGTAATTTAAATAACTAATATTTCTCGAAAAATTCTCTTTGTTTTTTAATTTTTCTTAACGTAAAAAAACAAGCAACTACCATTACGATa                                                  >  1:36854/1‑141 (MQ=255)
                                                                                cACCTTGTCATAAAAGCAGATTCAATTTTAGACTAATCAATTTTTTGCTGTAATTTAAATAACTAATATTTCTCGCAAAATTCTCTTTGTTTTTTAATTTTTCTTAACGTAAAAAAACAAGCAACTACCATTACGATAGTg                                               <  2:215198/141‑1 (MQ=255)
                                                                                 aCCTTGTCATAAAAGCAGATTCAATTTTAGACTAATCAATTTTTTTCTTTTTTTTAAATAACTAATATTTCTCGAAAAATTCTCTTTGTTTTTTAATTTTTCTTATCGTAAAAAAACAAGCAACTACCATTACGATAGTgt                                              >  2:283017/1‑140 (MQ=255)
                                                                                  ccTTGTCATAAAAGCAGATTCAATTTTAGACTAATCAATTTTTTGCTGTAATTTAAATAACTAATATTTCTCGAAAAATTCTCTTTGTTTTTTATTTTTTCTTAACGTAAAAAAACAAGCAACTACCATTACGATAGTGCt                                             >  1:264513/1‑141 (MQ=255)
                                                                                     tGTCATAAAAGCAGATTCAATTTTAGACTAATCAATTTTTTGCTGTAATTTAAATAACTAATATTTCTCGAAAAATTCTCTTTGTTTTTTAATTTTTCTTAACGTAACAAAACAAGCAACTCCCATTCCGAAAGTGCttct                                          >  2:126883/1‑139 (MQ=255)
                                                                                      tttaaaaatGCAGATTCAATTTTAGACTAATCAATTTTTTTCTGTAATTTAAATAACTAATATTTCTCGAAAAATTCTCTTTGTTTTTTAATTTTTCTTAACGTAAAAAAACAAGCAACTACCATTACGATAGTGCTTGtt                                         <  2:241575/136‑1 (MQ=255)
                                                                                        cATAAAAGCAGATTCAATTTTAGACTAATCAATTTTTTGCTGTAATTTAAATAACTAATATTTCTCGAAAAATTCTCTTTGTTTTTTAATTTTTCTTAACGTAAAAAAACAAGCAACTACCATTACGATAGTGCTTGtttt                                       <  1:274343/141‑1 (MQ=255)
                                                                                         aTAAAAGCAGATTCAATTTTAGACTAATCAATTTTTTGCTGTAATTTAAATAACTAATATTTCTCGAAAAATTCTCTTTGTTTTTTAATTTTTCTTAACGTAAAAAAACAAGCAACTACCATTACGATAGTGCTTGTTTTc                                      <  1:221015/141‑1 (MQ=255)
                                                                                           aaaaGCAGATTCAATTTTAGACTAATCAATTTTTTGCTGTAATTTAAATAACTAATATTTCTCGAAAAATTCTCTTTGTTTTTTAATTTTTCTTAACGTAAAAAAACAAGCAACTACCATTACGATAGTGCTTGTTTTCaa                                    <  1:30056/141‑1 (MQ=255)
                                                                                                  gATTCAATTTTAGACTAATCAATTTTTTGCTGTAATTTAAATAACTAATATTTCTCGAAAAATTCTCTTTGTTTTTTAATTTTTCTTAACGTAAAAAAACAAGCAACTACCATTACGATAGTACTTGTTTTCAAAAAAttt                             >  1:107216/1‑141 (MQ=255)
                                                                                                    ttCAATTTTAGACTAATCAATTTTTTGCTGTAATTTAAATAACTAATATTTCTCGAAAAATTCTCTTTGTTTTTTAATTTTTATTAACATCAACAAACGGGCAGCTGCCATTACGATAGTGCTTGTTTTCAAAAAATTTta                           >  2:212269/1‑141 (MQ=255)
                                                                                                       aaTTTTAGACTAATCAATTTTTTGCTGTAATTTAAATAACTAATATTTCTCGAAAAATTCTCTTTGTTTTTTAATTTTTCTTAACGTAAAAAAACAAGCAACTACCATTACGATAGTGCTTGTTTTCAAAAAATTTTATaa                        >  1:150354/1‑141 (MQ=255)
                                                                                                         ttttAGACTAATCAATTTTTTGCTGTAATTTAAATAACTAATATTTCTCGAAAAATTCTCTTTGTTTTTTAATTTTTCTTAACGTAAAAAAACAAGCAACTACCATTACGATAGTGCTTGTTTTCAAAAAATTTTATAAtt                      <  1:35462/141‑1 (MQ=255)
                                                                                                          tttAGAATAATCTATTTTTTGCTGTAATTTAAATAACTAATATTTCTCGAAATATTCTCTTTGTTTTTTAATTTTTCTTAACGTAAAAAAACAAGCAACTACCATTACGATAGTGCTTGTTTTCAAAAAATTTTATAATTg                     <  2:154521/141‑1 (MQ=255)
                                                                                                                   aTCAATTTTTTGCTGTAATTTAAATAACTAATATTTCTCGAAAAATTCTCTTTGTTTTTTAATTTTTCTTAACGTAAAAAAACAAGCAACTACCATTACGATAGTGCTTGTTTTCAAAAAATTTTATAATTGTTTTGGttt            <  2:40818/141‑1 (MQ=255)
                                                                                                                     cAATTTTTTGCTGTAATTTAAATAACTAATATTTCTCGAAAAATTCTCTTTGTTTTTTAATTTTTCTTAACGTAAAAAAACAAGCAACTACCATTACGATAGTGCTTGTTTTCAAAAAATTTTATAATTGTTTTGGTTTAg          <  1:17536/141‑1 (MQ=255)
                                                                                                                           tttGCTGTAATTTAAATAACTAATATTTCTCGAAAAATTCTCTTTGTTTTTTAATTTTTCTTAACGTAAAAAAACAAGCAACTACCATTACGATAGTGCTTGTTTTCAAAAAATTTTATAATTGTTTTGGTTTAGCTaaaa    <  1:33690/141‑1 (MQ=255)
                                                                                                                            ttGCTGTAATTTAAATAACTAATATTTCTCGAAAAATTCTCTTTGTTTTTTAATTTTTCTTAACGTAAAAAAACAAGCAACTACCATTACGATAGTGCTTGTTTTCAAAAAATTTTATAATTGTTTTGGTTTAGCTAAAAt   >  1:160246/1‑141 (MQ=255)
                                                                                                                             tGCTGTAATTTAAATAACTAATATTTCTCGAAAAATTCTCTTTGTTTTTTAATTTTTCTTAACGTAAAAAAACAAGCAACTACCATTACGATAGTGCTTGTTTTCAAAAAATTTTATAATTGTTTTGGTTTAGCTAAAAtt  >  2:206361/1‑141 (MQ=255)
                                                                                                                             tGCTGTAATTTAAATAACTAATATTTCTCGAAAAATTCTCTTTGTTTTTTAATTTTTCTTAACGTAAAAAAACAAGCAACTACCATTACGATAGTGCTTGTTTTCAAAAAATTTTATAATTGTTTTGGTTTAGCTAAAAtt  >  1:243607/1‑141 (MQ=255)
                                                                                                                                         |                                                                                                                                
ATCATGCAATTGATTGATTTAAACATCTTTGCATTGACTTACTCATTTTTCTTTAATTTAAAATAATCATTTATACTTTTCACCTTGTCATAAAAGCAGATTCAATTTTAGACTAATCAATTTTTTGCTGTAATTTAAATAACTAATATTTCTCGAAAAATTCTCTTTGTTTTTTAATTTTTCTTAACGTAAAAAAACAAGCAACTACCATTACGATAGTGCTTGTTTTCAAAAAATTTTATAATTGTTTTGGTTTAGCTAAAATT  >  CP000730/1682321‑1682586

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