Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 689,983 G→T 8.9% S30Y (TCT→TAT)  USA300HOU_0631 ← hypothetical membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP000730689,9830GT 8.9% 140.6 / 2.1 46S30Y (TCT→TAT) USA300HOU_0631hypothetical membrane protein
Reads supporting (aligned to +/- strand):  ref base G (21/20);  new base T (2/2);  total (24/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.97e-01

TTCGCAAATCTAAATCGTGGTCATACATTTCAAAATGAAAATCTTGATAACCTTTAATCGACTCTAAAAAATAACGCATACGCTTTTTTGATTCTTCTTTTGCATCATAATTTGCTGTTGATATAGACTCAGAACGAACACGATAATACGCCCAAACTTCGCCAGTTTTCGTTAATAACAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTT  >  CP000730/689852‑690118
                                                                                                                                   |                                                                                                                                       
ttCGCAAATCTAAATCGTGGTCATACATTTCAAAATGAAAATCTTGATAACCTTTATTCGACTCTAAAAAATAACGCATACGCTTTTTTGATTCTTCTTTTGCATCATAATTTGCTGTTGATATAGACTCAGAACGAacac                                                                                                                                <  1:218136/141‑1 (MQ=255)
     aaaTCTAAATCGTGGTCATACATTTCAAAATGAAAATCTTGATAACCTTTAATCGACTCTAAAAAATAACGCATACGCTTTTTTGATTCTTCTTTTGCGTCATAATTTGCTGTTGATATAGACTCAGAACGAACACGataa                                                                                                                           >  1:118466/1‑141 (MQ=255)
      aaTCTAAATCGTGGTCATACATTTCAAAATGAAAATCTTGATAACCTTTAATCGACTCTAAAAAATAACGCATACGCTTTTTTGATTCTTCTTTTGCATCATAATTTGCTGTTGATATAGACTCAGAACGAACACGataat                                                                                                                          <  1:403742/141‑1 (MQ=255)
          tAAATCGTGGTCATACATTTCAAAATGAAAATCTTGATAACCTTTAATCGACTCTAAAAAATAACGCATACGCTTTTTTGATTCTTCTTTTGCATCATAATTTGCTGTTGATATAGACTCATAACGAACACGATAATACGc                                                                                                                      >  2:428226/1‑141 (MQ=255)
          tAAATCGTGGTCATACATTTCAAAATGAAAATCTTGATAACCTTTAATCGACTCTAAAAAATAACGCATACGCTTTTTTGATTCTTCTTTTGCATCATAATTTGCTGTTGATATAGACTCATAACGAACACGATAATACGc                                                                                                                      >  2:346595/1‑141 (MQ=255)
            aaTCGTGGTCATACATTTCAAAATGAAAATCTTGATAACCTTTAATCGACTCTAAAAAATAACGCATACGCTTTTTTGATTCTTCTTTTGCATCATAATTTGCTGTTGATATAGACTCAGAACGAACACGATAATACGccc                                                                                                                    <  1:332324/141‑1 (MQ=255)
            aaTCGTGGTCATACATTTCAAAATGAAAATCTTGATAACCTTTAATCGACTCTAAAAAATAACGCATACGCTTTTTTGATTCTTCTTTTGCATCATAATTTGCTGTTGATATAGACTCAGAACGAACACGATAATACGccc                                                                                                                    <  2:245922/141‑1 (MQ=255)
              tCGTGGTCATACATTTCAAAATGAAAATCTTGATAACCTTTAATCGACTCTAAAAAATAACGCATACGCTTTTTTGATTCTTCTTTTGCATCATAATTTGCTGTTGATATAGACTCAGAACGAACACGATAATACGCCCaa                                                                                                                  >  1:46838/1‑141 (MQ=255)
                  ggTCATACATTTCAAAATGAAAATCTTGATAACCTTTAATCGACTCTAAAAAATAACGCATACGCTTTTTTGATTCTTCTTTTGCATCATAATTTGCTGTTGATATAGACTCAGAACGAACACGATAATACGCCCAAACtt                                                                                                              >  1:219210/1‑141 (MQ=255)
                      aTACATTTCAAAATGAAAATCTTGATAACCTTTAATCGACTCTAAAAAATAACGCATACGCTTTTTTGATTCTTCTTTTGCATCATAATTTGCTGTTGATATAGACTCAGAACGAACACGATAATACGCCCAAACTTCGcc                                                                                                          <  1:69701/141‑1 (MQ=255)
                        aCATTTCAAAATGAAAATCTTGATAACCTTTAATCGACTCTAAAAAATAACGCATACGCTTTTTTGATATTTATTTTCCATCATAATTTGCTGTTGATATAGAATTACAACGAACACGATAATACGCCCAAACTTCGCCaa                                                                                                        >  2:354081/1‑140 (MQ=255)
                          aTTTCAAAATGAAAATCTTGATAACCTTTAATCGACTCTAAAAAATAACGCATACGCTTTTTTGATTCTTCTTTTGCATCATAATTTGCTGTTGATATAGACTCAGAACGAACACGATAATACGCCCAAACTTCGCCAGtt                                                                                                      <  2:395049/141‑1 (MQ=255)
                                aaaTGAAAATCTTGATAACCTTTAATCGACTCTAAAAAATAACGCATACGCTTTTTTTTTTTTTTTTTTGCATCATAATTTGCTGTTGATATAGACTCAGAACGAACACGATAATACGCCCAAACTTCGCCAGTTTTCGtt                                                                                                >  2:301811/1‑141 (MQ=255)
                                 aaTGAAAATCTTGATAACCTTTAATCGACTCTAAAAAATAACGCATACGCTTTTTTGATTCTTCTTTTGCATCATAATTTGCTGTTGATATAGACTCAGAACGAACACGATAATACGCCCAAACTTCGCCAGTTTTCGTta                                                                                               >  1:239142/1‑141 (MQ=255)
                                  aTGAAAATCTTGATAACCTTTAATCGACTCTAAAAAATAACGCATACGCTTTTTTGATTCTTCTTTTGCATCATAATTTGCTGTTGATATAGACTCAGAACGAACACGATAATACGCCCAAACTTCGCCAGTTTTCGTtaa                                                                                              <  1:299172/141‑1 (MQ=255)
                                   tGAAAATCTTGATAACCTTTAATCGACTCTAAAAAATAACGCATACGCTTTTTTGATTCTTCTTTTGCATCATAATTTGCTGTTGATATAGACTCAGAACGAACACGATAATACGCCCAAACTTCGCCAGTTTTCGTtaat                                                                                             <  2:84282/141‑1 (MQ=255)
                                    gAAAATCTTGATAACCTTTAATCGACTCTAAAAAATAACGCATACGCTTTTTTGATTCTTCTTTTGCATCATAATTTGCTGTTGATATAGACTCAGAACGAACACGATAATACGCCCAAACTTCGCCAGTTTTCGTtaata                                                                                            <  2:441679/141‑1 (MQ=255)
                                    aaaaaTCTTGATAACCTTTACTCGACTCTAAAAAATAACGCATACGCTTTTTTGATTCTTCTTTTGCATCATAATTTGCTGTTGATATAGACTCAGAACGAACACGATAATACGCCCAAACTTCGCCAGTTTTCGTtaata                                                                                            <  2:382988/140‑1 (MQ=255)
                                     aaaaTCTTGATAACCTTTAATCGACTCTAAAAAATAACGCATACGCTTTTTTGATTCTTCTTTTGCATCATAATTTGCTGTTGATATAGACTCAGAACGAACACGATAATATCCCCAAACTTCGCCAGTTTCAGTtaataa                                                                                           >  2:361417/1‑141 (MQ=255)
                                     aaaaTCTTGATAACCTTTAATCGACTCTAAAAAATAACGCATACGCTTTTTTGATTCTTCTTTTGCATCATAATTTGCTGTTGATATAGACTCAGAACGAACACGATAATACGCCCAAACTTCGCCAGTTTTCGTtaataa                                                                                           >  1:271679/1‑141 (MQ=255)
                                      aaaTCTTGATCACCTTTAACCACCCCTCAAAAATCACGCATACGCTTTTTTGCTTCTTCTCTTGCATCATAATCTGCTGTAGATATCGACTCAGAACGAACACGATAATACGCCCAAACTTCGCCAGTTTTCGTTaacaac                                                                                          <  2:87074/141‑1 (MQ=255)
                                               tAACCTTTAATCGACTCTAAAAAATAACGCATACGCTTTTTTGATTCTTCTTTTGCATCATAATTTGCTGTTGATATAGACTCAGAACGAACACGATAATACGCCCAAACTTCGCCAGTTTTCGTTAATAACAAATTAtct                                                                                 <  2:339775/141‑1 (MQ=255)
                                                   cTTTAATCGACTCTAAAAAATAACGCATACGCTTTTTTGATTCTTCTTTTGCATCATAATTTGCTGTTGATATAGACTCAGAACGAACACGATAATACGCCCAAACTTCGCCAGTTTTCGTTAATAACAAATTATCTCTTa                                                                             >  1:80354/1‑141 (MQ=255)
                                                      tAATCGACTCTAAAAAATAACGCATACGCTTTTTTTTTTTTTTTTTTATATCATAATTTGATTTTGATTTTGATTCAGTACTTTTACGATAATACTCTCACACTTCGCCCTTTTTCGTTTATATTAAATTATCTTTTataa                                                                          >  2:337841/1‑138 (MQ=255)
                                                       aaTCGACTCTAAAAAATAACGCATACGCTTTTTTGATTCTTCTTTTGCATCATAATTTGCTGTTGATATAGACTCAGAACGAACACGATAATACGCCCAAACTTCGCCAGTTTTCGTTAATAACAAATTATCTCTTAAAcc                                                                         >  1:52645/1‑141 (MQ=255)
                                                             ctctAAAAAATAACGCATACGCTTTTTTTTTTTTTTTTTTGCATCATAATTTGCTGTTGATATAGACTCAGAACGAACACGATAATACGCCCAAACTTCGCCAGTTTTCTTTAATAACAAATTATCTCTTAAACCTGCATa                                                                   >  2:412118/1‑141 (MQ=255)
                                                                tAAAAAATAACGCATACGCTTTTTTGATTCTTCTTTTGCATCATAATTTGCTGTTGATATAGACTCAGAACGAACACGATAATACGCCCAAACTTCGCCAGTTTTCGTTAATAACAAATTATCTCTTAAACCTGCATACgg                                                                >  1:368730/1‑141 (MQ=255)
                                                                 aaaaaaTAACGCATACGCTTTTTTGATTCTTCTTTTGCATCATAATTTGCTGTTGATATAGACTCAGAACGAACACGATAATACGCCCAAACTTCGCCAGTTTTCGTTAATAACAAATTATCTCTTAAACCTGCATACGGt                                                               <  2:163293/141‑1 (MQ=255)
                                                                 aaaaaaTAACGCATACGCTTTTTTGATTCTTCTTTTGCATCATAATTTGCTGTTGATATAGACTCAGAACGAACACGATAATACGCCCAAACTTCGCCAGTTTTCGTTAATAACAAATTATCTCTTAAACCTGCATACGGt                                                               <  2:92474/141‑1 (MQ=255)
                                                                 aaaaaaTAACGCATACGCTTTTTTGATTCTTCTTTTGCATCATAATTTGCTGTTGATATAGACTCAGAACGAACACGATAATACGCCCAAACTTCGCCAGTTTTCGTTAATAACAAATTATCTCTTAAACCTGCATACGGt                                                               <  1:126488/141‑1 (MQ=255)
                                                                 aaaaaaTAACGCATACGCTTTTTTGATTCTTCTTTTGCATCATAATTTGCTGTTGATATAGACTCAGAACGAACACGATAATACGCCCAAACTTCGCCAGTTTTCGTTAATAACAAATTATCTCTTAAACCTGCATACGGt                                                               <  2:432636/141‑1 (MQ=255)
                                                                  aaaaaTAACGCATACGCTTTTTTGATTCTTCTTTTGCATCATAATTTGCTGTTGATATAGACTCAGAACGAACACGATAATACGCCCAAACTTCGCCAGTTTTCGTTAATAACAAATTATCTCTTAAACCTGCATACGGTg                                                              >  1:320614/1‑141 (MQ=255)
                                                                  aaaaaTAACGCATACGCTTTTTTGATTCTTCTTTTGCATCATAATTTGCTGTTGATATAGACTCAGAACGAACACGATAATACGCCCAAACTTCGCCAGTTTTCGTTAATAACAAATTATCTCTTAAACCTGCATACGGTg                                                              >  1:410868/1‑141 (MQ=255)
                                                                      aTAACGCATACGCTTTTTTGATTCTTCTTTTGCATCATAATTTGCTGTTGATATAGACTCAGAACGAACACGATAATACGCCCAAACTTCGCCAGTTTTCGTTAATAACAAATTATCTCTTAAACCTGCATACGGTGCttt                                                          <  1:183933/141‑1 (MQ=255)
                                                                         aCGCATACGCTTTTTTTTTTTTTTTTTTGCATCATAATTTGCTGTTGATATAGACTCAGAACGAACACGATAATACGCCCAAACTTCGCCAGTTTTCGTTAATAACAAATTATCTCTTAAACCTTCATACGGTGCTTTTaa                                                       >  2:8606/1‑141 (MQ=255)
                                                                          cGCATACGCTTTTTTGATTCTTCTTTTGCATCATAATTTGCTGTTGATATAGACTCAGAACGAACACGATAATACGCCCAAACTTCGCCAGTTTTCGTTAATAACAAATTATCTCTTAAACCTGCATACGGTGCTTTTaaa                                                      >  1:288963/1‑141 (MQ=255)
                                                                                        tGATTCTTCTTTTGCATCATAATTTGCTGTTGATATAGACTCAGAACGAACACGATAATACGCCCAAACTTCGCCAGTTTTCGTTAATAACAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCt                                        >  2:216310/1‑141 (MQ=255)
                                                                                                     gCATCATAATTTGCTGTTGATATAGACTCAGAACGAACACGATAATACGCCCAAACTTCGCCAGTTTTCGTTAATAACAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATaa                           >  2:122332/1‑141 (MQ=255)
                                                                                                     gCATCATAATTTGCTGTTGATATAGACTCAGAACGAACACGATAATACGCCCAAACTTCGCCAGTTTTCGTTAATAACAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATaa                           <  2:34541/141‑1 (MQ=255)
                                                                                                                  cTGTTGATATAGACTCATAACGAACACGATAATACGCCCAAACTTCGCCAGTTTTCGTTAATAACAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGc              <  1:346595/141‑1 (MQ=255)
                                                                                                                  cTGTTGATATAGACTCATAACGAACACGATAATACGCCCAAACTTCGCCAGTTTTCGTTAATAACAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGc              <  1:428226/141‑1 (MQ=255)
                                                                                                                  cTGTTGATATAGACTCAGAACGAACACGATAATACGCCCAAACTTCGCCAGTTTTCGTTAATAACAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGc              <  1:429472/141‑1 (MQ=255)
                                                                                                                  cTGTTGAGATAGACTCAGAACGAACACGATAATACGCCCAAACTTCGCCAGTTTTCGTTAATAACAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCGCTCCAATAAAAATAACCCCAATAGGGGc              >  1:403123/1‑141 (MQ=255)
                                                                                                                             gACTCAGAACGAACACGATAATACGCCCAAACTTCGCCAGTTTTCGTTAATAACAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGt   >  2:53214/1‑141 (MQ=255)
                                                                                                                             gACTCAGAACGAACACGATAATACGCCCAAACTTCGCCAGTTTTCGTTAATAACAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGt   >  2:263896/1‑141 (MQ=255)
                                                                                                                              aCTCAGAACGAACACGATAATACGCCCAAACTTCGCCAGTTTTCGTTAATAACAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGtt  <  2:54892/141‑1 (MQ=255)
                                                                                                                                   |                                                                                                                                       
TTCGCAAATCTAAATCGTGGTCATACATTTCAAAATGAAAATCTTGATAACCTTTAATCGACTCTAAAAAATAACGCATACGCTTTTTTGATTCTTCTTTTGCATCATAATTTGCTGTTGATATAGACTCAGAACGAACACGATAATACGCCCAAACTTCGCCAGTTTTCGTTAATAACAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTT  >  CP000730/689852‑690118

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