Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 643,640 A→G 21.3% intergenic (‑288/+261) 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,6400AG21.3% 49.2 / 9.5 28intergenic (‑288/+261)yjjP_1/USA300HOU_RS03070Inner membrane protein YjjP/Putative heme‑dependent peroxidase
Reads supporting (aligned to +/- strand):  ref base A (13/9);  new base G (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 = 9.81e-01

TCTCTTGAAGTATAAAAAGGAAAAATTATTGAAAAAAATTCGACAAAAATATGTTTCAAATATTTCAAATAGGTTAAAATAATATTAGTTAACCATTACAAAAATTGTATAGAGTAGCGACTGTATAATTTCTATTGAGGTTAACGTTTATATGTAGTGATAGTAGTTAAAGTTCTCCCAAGGAAGACTACTCGGGTACACTTTGCTATGAGCAAAGTGTACTTTGTTATTGATAATACATTAGCACATATATATAAGTTTAAACATACAGA  >  USA300TCH1516_ALE/643507‑643778
                                                                                                                                     |                                                                                                                                          
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)
                                                                                                                                     |                                                                                                                                          
TCTCTTGAAGTATAAAAAGGAAAAATTATTGAAAAAAATTCGACAAAAATATGTTTCAAATATTTCAAATAGGTTAAAATAATATTAGTTAACCATTACAAAAATTGTATAGAGTAGCGACTGTATAATTTCTATTGAGGTTAACGTTTATATGTAGTGATAGTAGTTAAAGTTCTCCCAAGGAAGACTACTCGGGTACACTTTGCTATGAGCAAAGTGTACTTTGTTATTGATAATACATTAGCACATATATATAAGTTTAAACATACAGA  >  USA300TCH1516_ALE/643507‑643778

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: