Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 76,915 G→A 15.2% intergenic (+973/‑161) 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,9150GA15.2% 93.1 / 7.1 38intergenic (+973/‑161)USA300HOU_0072/USA300HOU_0075universal stress protein/lysophospholipase
Reads supporting (aligned to +/- strand):  ref base G (12/20);  new base A (2/4);  total (14/24)
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.90e-01

TAAAAGGTATTGAATGTATTTACGCTCTATATAAAAAGAACCGCAGGTCTCTTCAGATCTACGGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTA  >  CP000730/76782‑77052
                                                                                                                                     |                                                                                                                                         
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                                                                    >  2:321990/1‑141 (MQ=255)
                                                                aTTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAGTa                                                                    >  1:171136/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)
                                                                                                                  acgACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATCTATATAGTATTACTGCAAACATTGATTTGGGTCCAACATTAGGGGTATTAAAAAcatga                  <  2:395580/139‑1 (MQ=255)
                                                                                                                          acaaaaTAGTGGTTAGAAATAATAAGATACGTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAGTATTACTGTAAACATAAAAAAAAAAAAAAAATTAGGGGTATTAAAAAAATGAATGAATCa          <  2:197230/139‑1 (MQ=255)
                                                                                                                            aaaTTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATa        <  1:198534/141‑1 (MQ=255)
                                                                                                                            aaaTTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATa        <  1:318955/141‑1 (MQ=255)
                                                                                                                                 aGTGGTTAGCTATATTTTGTTACTTTGCAACAGAACCTATAAAAAACATGGGCGTAAAGGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAAtttt   <  1:45340/141‑1 (MQ=255)
                                                                                                                                  gTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTa  >  1:93082/1‑141 (MQ=255)
                                                                                                                                     |                                                                                                                                         
TAAAAGGTATTGAATGTATTTACGCTCTATATAAAAAGAACCGCAGGTCTCTTCAGATCTACGGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTA  >  CP000730/76782‑77052

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: