Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 643,642 T→A 21.5% intergenic (‑290/+259) yjjP_1 ← / ← USA300HOU_RS03070 Inner membrane protein YjjP/Putative heme‑dependent peroxidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE643,6420TA21.5% 45.1 / 11.4 28intergenic (‑290/+259)yjjP_1/USA300HOU_RS03070Inner membrane protein YjjP/Putative heme‑dependent peroxidase
Reads supporting (aligned to +/- strand):  ref base T (13/9);  new base A (3/3);  total (16/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.75e-01

TCTCTTGAAGTATAAAAAGGAAAAATTATTGAAAAAAATTCGACAAAAATATGTTTCAAATATTTCAAATAGGTTAAAATAATATTAGTTAACCATTACAAAAATTGTATAGAGTAGCGACTGTATAATTTCTATTGAGGTTAACGTTTATATGTAGTGATAGTAGTTAAAGTTCTCCCAAGGAAGACTACTCGGGTACACTTTGCTATGAGCAAAGTGTACTTTGTTATTGATAATACATTAGCACATATATATAAGTTTAAACATACAGATTTC  >  USA300TCH1516_ALE/643507‑643782
                                                                                                                                       |                                                                                                                                            
tctcTTGAAGTATAAAAAGGAAAAATTATTGAAAAAAATTCGACAAAAATATGTTTCAAATATTTCAAATAGGTTAAAATAATATTAGTTAACCATTACAAAAATTGTATAGAGTAGCGACTGTATAATTTCTATTGAGGt                                                                                                                                         <  1:61685/141‑1 (MQ=255)
      gAAGTATAAAAAGGAAAAATTATTGAAAAAAATTCGACAAAAATATGTTTCAAATATTTCGAATAGGTTAAAATAATATTAGTTAACCATTACAAAAATTGTATAGAGTAGCGACTGTATAATTTCTATTGAGGTTAACGt                                                                                                                                   >  1:2216/1‑141 (MQ=255)
           ataAAAAGGAAAAATTATTGAAAAAAATTCGACAAAAATATGTTTCAAATATTTCAAATAGGTTAAAATAATATTAGTTAACCATTACAAAAATTGTATAGAGTAGCGACTGTATAATTTCTATTGAGGTTAACGTTtata                                                                                                                              >  1:169905/1‑141 (MQ=255)
           ataAAAAGGAAAAATTATTGAAAAAAATTCGACAAAAATATGTTTCAAATATTTCAAATAGGTTAAAATAATATAAGTAAACGATTACAAAAATTGTATAGAGTAGCGACTGTATAATTTCTATTGAGGTTAACGTTtata                                                                                                                              >  2:213802/1‑141 (MQ=255)
                   gAAAAATTATTGAAAAAAATTCGACAAAAATATGTTTCAAATATTTCAAATAGGTTAAAATAATATTAGTTAACCATTACAAAAAATGTATAGAGTAGCGACTGTATAATTTCTATTGAGGTTAACGTTTATATGTAGTGa                                                                                                                      >  2:239532/1‑141 (MQ=255)
                     aaaaTTATTGAAAAAAATTCGACAAAAATATGTTTCAAATATTTCAAATAGGTTAAAATAATATTAGTTAACCATTACAAAAATTGTATAGAGTAGCGACTGTATAATTTCTATTGAGGTTAACGATTATATGTAGTGAta                                                                                                                    >  2:236085/1‑141 (MQ=255)
                             tGAAAAAAATTCGACAAAAATATGTTTCAAATATTTCAAATAGGTTAAAATAATATTAGTTAACCATTACAAAAATTGTATAGAGTAGCGACTGTATAATTTCTATTGAGGTTAACGTTTATATGTAGTGATAGTAGTTaa                                                                                                            <  1:212729/141‑1 (MQ=255)
                               aaaaaaaTTCGACAAAAATATGTTTCAAATATTTCAAATAGGTTAAAATAATATTAGTTAACCATTACAAAAATTGTATAGAGTAGCGACTGTATAATTTCTATTGAGGTTAACGTTTATATGTAGTGATAGTAGTTAAAg                                                                                                          <  1:67050/141‑1 (MQ=255)
                                                     tttcaaatatttcaaataGGTTAAAATAATATTAGTTAACCATTACGAAGCTTGTCACCCTTCCCGACTGTATAATTTTTGTAATGGTTAACGTTTATATGTAGTGATAGTAGTTAAAGTTCTCCCAAGGAAGACTACTCg                                                                                    <  1:116104/141‑1 (MQ=255)
                                                       tcaaatatttcaaataGGTTAAAATAATATTAGTTAACCATTACAAAAATTGTATAGAGTAGCGACTGTATAATTTTTGTAATGGTTAACGTTTATATGTAGTGATAGTAGTTAAAGTTCTCCCAAGGAAGACTACTCggg                                                                                  <  1:107723/141‑1 (MQ=255)
                                                           atatttcaaataGGTTAAAATAATATTAGTTAACCATTACAAAAATTGTATAGAGTAGCGACTGTATAATTTCTATTGAGGTTAACGTTTATATGTAGTGATAGTAGTTAAAGTTCTCCCAAGGAAGACTACTCGGGTaca                                                                              <  2:169905/141‑1 (MQ=255)
                                                             atttcaaataGGTTAAAATAATATTAGTTAACCATTACAAAAATTGTATAGAGTAGCGACTGTATAATTTCTATTGAGGTTAACGTTTATATGTAGTGATAGTAGTTAAAGTTCTCCCAAGGAAGACTACTAGAGTaaaat                                                                            >  2:84627/1‑137 (MQ=255)
                                                              tttcaaataGGTTAAAATAATATTAGTTAACCATTACAAAAATTGTATAGAGTAGCGACTGTATAATTTCTATTGAGGTTAACGTTTATATGTAGTGATAGTAGTTAAAGTTCTCCCAAGGAAGACTACTCGGGTACACtt                                                                           >  2:105870/1‑141 (MQ=255)
                                                                            aaaTAATATTAGTTAACCATTACAAAAATTGTATAGAGTAGCGACTGTATAATTTTTGTAATGGTTAACGTTTATATGTAGTGATAGTAGTTAAAGTTCTCCCAAGGAAGACTACTCGGGTACACTTTGCTATGAGCAAAg                                                             <  2:48191/141‑1 (MQ=255)
                                                                                     tAGTTAACCATTACAAAAATTGTATAGAGTAGCGACTGTATAATTTCTATTGAAGTTAACGTTTATATGTAGTGATAGTAGTTAAAGTTCTCCCAAGGAAGACTACTCGCGTACACTTTGCTATGAGCAAAGTGTACtttc                                                    >  2:71748/1‑140 (MQ=255)
                                                                                      aGTTAACCATTACAAAAATTGTATAGAGTAGCGACTGTATAATTTTTGTAATGGTTAACGTTTATATGTAGTGATAGTAGTTAAAGTTCTCCCAAGGAAGAATACTCGGGTACACTTTGATATGAGCAAAGTGTACTTTGt                                                   >  2:223308/1‑141 (MQ=255)
                                                                                             cATTACAAAAATTTTATAGAGTAGTGATTGTATAATTTCTATTGAGTTTAATGTTTATATGTAGTGATAGTAGTTAAAGTTCTCCCAAGGAAGACTACTCGGGTACACTTTGCTATGAGCAAAGTGTACTTTGTTATTGat                                            <  2:242674/141‑1 (MQ=255)
                                                                                             cATTACAAAAATTGTATAGAGTAGCGACTGTATAATTTTTGTAATGGTTAACGTTTATATGTAGTGATAGTAGTTAAAGTTCTCCCAAGGAAGACTACTCGGGTACACTTTGCTATGAGCAAAGTGTACTTTTTTATTGat                                            >  2:126882/1‑141 (MQ=255)
                                                                                             cATTACAAAAATTGTATAGAGTAGCGACTGTATAATTTCTATTGAGGTTAACGTTTATATGTAGTGATAGTAGTTAAAGTTCTCCCAAGGAAGACTACTCGGGTACACTTTGCTATGAGCAAAGTGTACTTTGTTATTGat                                            <  2:122812/141‑1 (MQ=255)
                                                                                                 aCAAAAATTGTATAGAGTATCTGCTGTATAATTTCGATTGAGGTTAACGTTTATATGTAGTGATAGTAGTTAAAGTTCTCCCAAGGAAGACTACTCGGGTACACTTTGCTATGAGCAAAGTGTACTTTGTTATTGataata                                        >  1:183379/1‑141 (MQ=255)
                                                                                                 aCAAAAATGGTATAGAGTAGCGACAGTATAATTTTTGTAATGGTTAACGTCTATATGTAGTGATAGTAGTTAAAGTTCTCCCAAGAAAGACTACTCGCGTACACTTTGCTAAGAGCAAAGTCTCCTTTGTTATTGataata                                        >  2:48237/1‑141 (MQ=255)
                                                                                                  aaaaatTTTTATGGAGTAGCGACTGTATAATTTCTATTGAGGTTAACGTTTATATGTAGTGATAGTAGTTAAAGTTCTCCCAAGGAAGACTACTCGGGTACACTTTGCTATGAGCAAAGTGTACTTTGTTATTGATAATAc                                       <  2:34112/140‑1 (MQ=255)
                                                                                                           tataGAGTAGCGACTGTATAATTTCTATTGAGGTTAACGTTTATATGTAGTGATAGTAGTTAAAGTTCTCCCAAGGAAGACTACTCGGGTACACTTTGCTATGAGCAAAGTGTACTTTGTTATTGATAATACATTAGcaca                              >  1:172518/1‑141 (MQ=255)
                                                                                                            ataGAGTAGCGACTGTATAATTTCTATTGAGGTTAACGTTTATATGTAGTGATAGTAGTTAAAGTTCTCCCAAGGAAGACTACTCGGGTACACTTTGCTATGAGCAAAGTGTACTTTGTTATTGATAATACATTAGCACat                             >  1:139177/1‑141 (MQ=255)
                                                                                                                            ataATTTCTATTGAGGTTAACGTTTATATGTAGTGATAGTAGTTAAAGTTCTCCCAAGGAAGACTACTCGGGTACACTTTGCTATGAGCAAAGTGTACTTTGTTATTGATAATACATTAGCACAAATATATAAGTTTaacc             >  2:127323/1‑139 (MQ=255)
                                                                                                                               aTTTCTATTGAGGTTAACGTTTATATGTAGTGATAGTAGTTAAAGTTCTCCCAAGGAAGACTACTCGGGTACACTTTGCTATGAGCAAAGTGTACTTTGTTATTGATAATACATTAGCACATATATATAAGTTTAAACATa          >  1:233119/1‑141 (MQ=255)
                                                                                                                                tttCTATTGAGGTTAACGTTTATATGTAGTGATAGTAGTTAAAGTTCTCCCAAGGAAGACTACTCGGGTACACTTTGCTATGAGCAAAGTGTACTTTGTTATTGATAATACATTAGCACATATATATAAGTTTAAACATAc         <  1:13579/141‑1 (MQ=255)
                                                                                                                                   cTATTGAGGTTAACGTTTATATGTAGTGATAGTAGTTAAAGTTCTCCCAAGGAAGACTACTCGGGTACACTTTGCTATGAGCAAAGTGTACTTTGTTATTGATAATACATTAGCACATATATATAAGTTTAAACATACAGa      <  1:207919/141‑1 (MQ=255)
                                                                                                                                       tGAGGTTAACGTTTATATGTAGTGATAGTAATTAAAGTTCTCCCAAGGAAGACTACTCGGGTACACTTTGCTATGAGCAAAGTGTACTTTGTTATTGATAATACATTAGCACATATATATAAGTTTAAACATACAGATTTc  <  1:93349/141‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                            
TCTCTTGAAGTATAAAAAGGAAAAATTATTGAAAAAAATTCGACAAAAATATGTTTCAAATATTTCAAATAGGTTAAAATAATATTAGTTAACCATTACAAAAATTGTATAGAGTAGCGACTGTATAATTTCTATTGAGGTTAACGTTTATATGTAGTGATAGTAGTTAAAGTTCTCCCAAGGAAGACTACTCGGGTACACTTTGCTATGAGCAAAGTGTACTTTGTTATTGATAATACATTAGCACATATATATAAGTTTAAACATACAGATTTC  >  USA300TCH1516_ALE/643507‑643782

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: