Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 76,889 G→A 22.8% intergenic (+947/‑187) USA300HOU_0072 → / → USA300HOU_0075 universal stress protein/lysophospholipase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP00073076,8890GA22.8% 73.2 / 11.5 35intergenic (+947/‑187)USA300HOU_0072/USA300HOU_0075universal stress protein/lysophospholipase
Reads supporting (aligned to +/- strand):  ref base G (12/15);  new base A (2/6);  total (14/21)
Fisher's exact test for biased strand distribution p-value = 4.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

TATCAAAATATCAATACAGCAAAGAATACTTTAAAAGGTATTGAATGTATTTACGCTCTATATAAAAAGAACCGCAGGTCTCTTCAGATCTACGGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTA  >  CP000730/76751‑77027
                                                                                                                                          |                                                                                                                                          
tATCAAAATATCAATACGGCAAAGAATACTTTAAAAGGTATTGATTGTATTTACGCTCTATATAAAAAGAACCGTAGGTCTCTTCAGATCTACGGATTTTCCCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACg                                                                                                                                          <  2:262835/141‑1 (MQ=11)
       atatCAATACAGCAAAGAATACTTTAAAAGGTATTGAATGTATTTACGCTCTATATAAAAAGAACCGCAGGTCTCTTCAGATCTACGGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTg                                                                                                                                   >  1:230292/1‑141 (MQ=11)
              aacggCAAAGAATACTTTAAAAGGTATTGAATGTATTTACGCTTTAAATGGTATGAACCGTAGGTCTCTTCAGATCTACGGATTTTCCCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACTCTGACATGAt                                                                                                                            <  1:288281/137‑1 (MQ=1)
                               tacaaGGGATTGAATGTACTTACGCTATTTATAAAAAGAAACGTAGGTCTCTTCAGATCTACGGATTTTCCCCATGACACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCta                                                                                                           <  1:119266/138‑1 (MQ=1)
                                    ggtcttgggtggaTTTACGCTCTATATAAAAAGAACCGCAGGTCTCTTCAGATCTACGGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATAttt                                                                                                      >  2:344883/13‑141 (MQ=1)
                                           agtgcaTTTACGCTCTATATAAAAAGAACCGCACGTCTCTTCAGATCTACGGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATAttttttactt                                                                                               <  1:20127/136‑5 (MQ=1)
                                                 ttgaCGCTCTATATAAAAAGAACCGCAGGTCTCTTCAGATCTACGGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATATTTTTTACTTTGCAAc                                                                                        >  2:111397/4‑141 (MQ=1)
                                                        tcgatataAAAATAACCGCAGGTCTCTTCAGATCTACGGATTTTCGCCATGCCACGAAATGAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATATTTTTTACTTTGCAACAGAAcca                                                                                 <  2:225532/138‑2 (MQ=1)
                                                        tccatctaAAAAGAACCGCAGGTCTCTTCACATCTTCGGATTTTGGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAAcc                                                                                  <  2:217049/135‑1 (MQ=1)
                                                        cctacagaaaaaGAACCGTAGGTCTCTTCAGATCTACGGATTTTCCCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATATTTTTTACTTTGCAACAGAAcca                                                                                 <  1:218019/134‑2 (MQ=1)
                                                          tatataAAAAGAACCGCAGGTCTCTTCAGATCTACGGATTTTCGCCATGCCACGAAATTAGCATCAGCCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCta                                                                                >  2:162075/1‑141 (MQ=11)
                                                           atataAAAAGAACCGCAGGTCTCTTCAGATCTACGGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCtat                                                                               >  1:145919/1‑141 (MQ=31)
                                                             ataAAAAGAACCGCAGGTCTCTTCAGATCTACGGATTTTCGCCATGCCACGAAATGAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATaa                                                                             <  2:21387/141‑1 (MQ=17)
                                                               aaaaaGAACCGCAGGTCTCTTCAGATCTACGGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATaaaa                                                                           <  1:237099/141‑1 (MQ=34)
                                                               aaaaaGAACCGCAGGTCTCATCAGATCTACGGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATaaaa                                                                           >  1:238324/1‑141 (MQ=21)
                                                                aaaaGAAACGCAGCTCTTTTCAGATCTACGGATGTTCGCTATGCCACGAAATTAGCACCATGATAGCAAGTTAAGCGAACACTGACATTATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATaaaaa                                                                          <  2:86919/141‑1 (MQ=11)
                                                                      aCCGCAGGCCTCTTCAGATCTACGGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATgg                                                                    <  2:347105/141‑1 (MQ=255)
                                                                        cGCAGGTCTCTTCAGATCTACGGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGAGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAACCATGGGc                                                                  >  1:303927/1‑141 (MQ=255)
                                                                         gCAGGTCTCTTCAGATCTACGGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAAAACCTACAAAAAACATGGGCt                                                                 <  1:213476/141‑1 (MQ=255)
                                                                         gCAGGTCTCTTCAGATCTACGGACTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCt                                                                 <  2:107346/141‑1 (MQ=255)
                                                                             gTCTCTTCAGATCTACGGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTaaa                                                             >  2:115110/1‑141 (MQ=255)
                                                                                  ttCAGATCTACGGATTTCTGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGtt                                                        <  2:18655/141‑1 (MQ=255)
                                                                                  ctCAGAGCTACGGATTTTCGCCATGCCACGCAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATCTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGtt                                                        <  2:256252/140‑1 (MQ=255)
                                                                                           gcgtgttttCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTtatat                                               <  2:195641/136‑1 (MQ=255)
                                                                                           aCGGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTtatat                                               >  1:122281/1‑141 (MQ=255)
                                                                                            cGGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTtatata                                              <  1:163463/141‑1 (MQ=255)
                                                                                               aTTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAGTa                                           >  1:171136/1‑141 (MQ=255)
                                                                                               aTTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAGTa                                           >  2:321990/1‑141 (MQ=255)
                                                                                                ttttCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTAAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAGTAt                                          <  2:50317/141‑1 (MQ=255)
                                                                                                    cGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAGTATTACt                                      >  2:398822/1‑141 (MQ=255)
                                                                                                             aCGAAAGTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAGTATTACTGCAAACAtt                             <  1:96047/141‑1 (MQ=255)
                                                                                                                        atcatGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCa                  >  1:334940/1‑141 (MQ=255)
                                                                                                                         tcatGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAGTATTACTGCAAACAGTGATTTGGGTCaa                 <  2:279530/141‑1 (MQ=255)
                                                                                                                                       tAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATACAAAACATGGGCTTAAAGGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGAGGTAtt   >  2:1362/1‑141 (MQ=255)
                                                                                                                                        aaGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTa  <  2:122281/141‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                          
TATCAAAATATCAATACAGCAAAGAATACTTTAAAAGGTATTGAATGTATTTACGCTCTATATAAAAAGAACCGCAGGTCTCTTCAGATCTACGGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTA  >  CP000730/76751‑77027

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: