Predicted mutation
evidence seq id position mutation annotation gene description
RA CP000730 423,247 A→G V112V (GTA→GTG)  USA300HOU_0394 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP000730423,2470AG92.1% 109.5 / 1.8 38V112V (GTA→GTG) USA300HOU_0394hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (3/0);  new base G (19/16);  total (22/16)
Fisher's exact test for biased strand distribution p-value = 2.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.37e-01

AGTGTTTTTAATAAAGAAGGGAAAGAGACGCAATTAATGTATACTGCTACATTTGATGTTAAACCGCATAGATACTTGAAAATTACACATAAAGGTCATCATGTAGAAACTTTTGAAGAAGTTGAAAAGGAAGAAGTACCAAAGAAAGCATTAGACAAACTGAGTCGATAATAGCATGCTTATATTACATGGTTCATTTATAAAAGGAGTACGAACGAAAGTAACGCATGACGATTAATTTAAAAATATTTGTAATAATTATGAATAA  >  CP000730/423110‑423377
                                                                                                                                         |                                                                                                                                  
aGTGTTTTTAATAAAGAAGGGAAAGAGACGCAATTAATGTATACTGCTACATTTGATGTTAAACCGCATAGATACTTGAAAATTACACATAAAGGTCATCATGTAGAAACTTTTGAAGAAGTTGAAAAGGAAGAAGTGCCa                                                                                                                                 >  2:152208/1‑141 (MQ=255)
   gtTTTTAATAAAGAAGGGAAAGAGACGCAATTAATGTATACTGCTACATTTGATGTTAAACCGCATAGATACTTGAAAATTACACATAAAGGTCATCATGTAGAAACTTTTGAAGAAGTTGAAAAGGAAGAAGTGCCaaag                                                                                                                              <  1:11657/141‑1 (MQ=255)
   gtTTTTAATAAAGAAGGGAAAGAGACGCAATTAATGTATACTGCTACATTTGATGTTAAACCGCATAGATACTTGAAAATTACACATAAAGGTCATCATGTAGAAACTTTTGAAGAAGTTGAAAAGGAAGAAGTGCCaaag                                                                                                                              <  2:131827/141‑1 (MQ=255)
              agaagGGAAAGAGACGCAATTAATGTATACTGCTACATTTGATGTTAAACCGCATAGATACTTGAAAATTACACATAAAGGTCATCATGTAGAAACTTTTGAAGAAGTTGAAAAGGAAGAAGTGCCAAAGAAAGCATTAGa                                                                                                                   >  2:353857/1‑141 (MQ=255)
               gaagGGAAAGAGACGCAATTAATGTATACTGCTACATTTGATGTTAAACCGCATAGATACTTGAAAATTACACATAAAGGTCATCATGTAGAAACTTTTGAAGAAGTTAAAAAGGAAGAAGTGCCAAAGAAAGCATTAGAc                                                                                                                  >  2:65675/1‑141 (MQ=255)
               gaagAGAAAGCGACGCAATTAATGTATACTGATAACTTAGAAGTTAAACCGCATAAAAACTAAAAAATTAAAAATAAAGGTAAAAATGTAAAAAAATTTGAAGAAGTTGAAAAGGAAGAAGAGAGAAAGAAAGCATTAGAc                                                                                                                  <  2:55157/141‑1 (MQ=255)
                     aaaGAGACGCAATTAATGTATACTGCTACATTTGATGTTAAACCGCATAGATACTTGAAAATTACACATAAAGGTCATCATGTAGAAACTTTTGAAGAAGTTGAAAAGGAAGAAGTGCCAAAGAAAGCATTAGACAAACTg                                                                                                            >  2:107347/1‑141 (MQ=255)
                                   aaTGTATACTGCTACATTTGATGTTAAACCGCATAGATACTTGAAAATTACACATAAAGGTCATCATGTAGAAACTTTTGAAGAAGTTGAAAAGGAAGAAGTGCCAAAGAAAGCATTAGACAAACTGAGTCGATAATAGCa                                                                                              <  1:233713/141‑1 (MQ=255)
                                       tataCTGCTACATTTGATGTTAAACCGCATAGATACTTGAAAATTACACATAAAGGTCATCATGTAGAAACTTTTGAAGAAGTTGAAAAGGAAGAAGTGCCAAAGAAAGCATTAGACAAACTGAGTCGATAATAGCATGCt                                                                                          >  1:230684/1‑141 (MQ=255)
                                        ataCTGCTACATTTGATGTTAAACCGCATAGATACTTGAAAATTACACATAAAGGTCATCATGTAGAAACTTTTGAAGAAGTTGAAAAGGAAGAAGTGCCAAAGAAAGCATTAGACAAACTGAGTCGATAATAGCATGCtt                                                                                         >  1:101454/1‑141 (MQ=255)
                                        ataCTGCTACATTTGATGTTAAACCGCATAGATACTTGAAAATTACACATAAAGGTCATCATGTAGAAACTTTTGAAGAAGTTGAAAAGGAAGAAGTGCCAAAGAAAGCATTAGACAAACTGAGTCGATAATAGCATGCtt                                                                                         >  1:365752/1‑141 (MQ=255)
                                         taCTGCTACATTTGATGTTAAACCGCATAGATACTTGAAAATTACACATAAAGGTCATCATGTAGAAACTTTTGAAGAAGTTGAAAAGGAAGAAGTGCCAAAGAAAGCATTAGACAAACCGAGTCGATAATAGCATGCTta                                                                                        <  1:315053/141‑1 (MQ=255)
                                         taCTGCTACATTTGATGTTAAACCGCATAGATACTTGAAAATTACACATAAAGGTCATCATGTAGAAACTTTTGAAGAAGTTGAAAAGGAAGAAGTGCCAAAGAAAGCATTAGACAAACCGAGTCGATAATAGCATGCTta                                                                                        <  1:315045/141‑1 (MQ=255)
                                             gCTACATTTGATGTTAAACCGCATAGATACTTGAAAATTACACATAAAGGTCATCATGTAGAAACTTTTGAAGAAGTTGAAAAGGAAGAAGTACCAAAGAAAGCATTAGACAAACTGAGTCGATAATAGCATGCTTATAtt                                                                                    >  2:302008/1‑141 (MQ=255)
                                              cTACATTTGATGTTAAACCGCATAGATACTTGAAAATTACACATAAAGGTCATCATGTAGAAACTTTTGAAGAAGTTGAAAAGGAAGAAGTACCAAAGAAAGCATTAGACAAACTGAGTCGATAATAGCATGCTTATATTa                                                                                   >  1:268631/1‑141 (MQ=255)
                                              cTACATTTGATGTTAAACCGCATAGATACTTGAAAATTACACATAAAGGTCATAATGTAGAAACTTTTGAAGAAGTTGAAAAGGAAGAAGTACCAAAGAAAGCATTAGACAAACTGAGTCGATAATAGCATGCTTATATTa                                                                                   >  1:268645/1‑141 (MQ=255)
                                                 cATTTGATGTTAAACCGCATAGATACTTGAAAATTACACATAAAGGTCATCATGTAGAAACTTTTGAAGAAGTTGAAAAGGAAGAAGTGCCAAAGAAAGCATTAGACAAACTGAGTCGATAATAGCATGCTTATATTACAt                                                                                >  1:302625/1‑141 (MQ=255)
                                                    ttGATGTTAAACCGCATAGATACTTGAAAATTACACATAAAGGTCATCATGTAGAAACTTTTGAAGAAGTTGAAAAGGAAGAAGTGCCAAAGAAAGCATTAGACAAACTGAGTCGATAATAGCATGCTTATATTACATGGt                                                                             <  1:353857/141‑1 (MQ=255)
                                                        tGTTAAACCGCATAGATACTTGAAAATTACACATAAAGGTCATCATGTAGAAACTTTTGAAGAAGTTGAAAAGGAAGAAGTGCCAAAGAAAGCATTAGACAAACTGAGTCGATAATAGCATGCTTATATTACATGGTTCAt                                                                         <  1:304258/141‑1 (MQ=255)
                                                         gTTAAACCGCATAGATACTTGAAAATTACACATAAAGGTCATCATGTAGAAACTTTTGAAGAAGTTGAAAAGGAAGAAGTGCCAAAGAAAGCATTAGACAAACTGAGTCGATAATAGCATGCTTATATTACATGGTTCAtt                                                                        >  2:61282/1‑141 (MQ=255)
                                                              aCCGCATAGATACTTGAAAATTACACATAAAGGTCATCATGTAGAAACTTTTGAAGAAGTTGAAAAGGAAGAAGTGCCAAAGAAAGCATTAGACAAACTGAGTCGATAATAGCATGCTTATATTACATGGTTCATTTATaa                                                                   >  2:154710/1‑141 (MQ=255)
                                                                   aTAGATACTTGAAAATTACACATAAAGGTCATCATGTAGAAACTTTTGAAGAAGTTGAAAAGGAAGAAGTGCCAAAGAAAGCATTAGACAAACTGAGTCGATAATAGCATGCTTATATTACATGGTTCATTTATAAAAGGa                                                              >  2:314291/1‑141 (MQ=255)
                                                                               aaaTTACACATAAAGGTCATCATGTAGAAACTTTTGAAGAAGTTGAAAAGGAAGAAGTGCCAAAGAAAGCATTAGACAAACTGAGTCGATAATAGCATGCTTATATTACATGGTTCATTTATAAAAGGAGTACGAACGaaa                                                  <  2:218802/141‑1 (MQ=255)
                                                                                  ttACACATAAAGGTCATCATGTAGAAACTTTTGAAGAAGATGAAAAGGAAGAAGTGCCAAAGAAAGCATTAGACAAACTGAGTCGATAATAGCATGCTTATATTACATGGTTCATTTATAAAAGGAGTACGAACGAAAGTa                                               >  2:189421/1‑141 (MQ=255)
                                                                                         aaaaGGTCATCATGAAAAACCTTTTAAAGAAGTTGAAAAGGAAGAAGTGCC‑AAGAAAGCATTAGACAAACTGAGTCGATAATAGCATGCTTATATTACATGGTTCATTTATAAAAGGAGTACGAACGAAAGTAACGCATGa                                       <  1:215606/140‑1 (MQ=255)
                                                                                                 atcatGTAGAAACTTTTGAAGAAGTTGAAAAGGAAGAAGTGCCAAAGAAAGCATTAGACAAACTGAGTCGATAATAGCATGCTTATATTACATGGTTCATTTATAAAAGGAGTACGAACGAAAGTAACGCATGACGATTaa                                >  1:24842/1‑141 (MQ=255)
                                                                                                 atcatGTAGAAACTTTTGAAGAAGTTGAAAAGGAAGAAGTGCCAAAGAAAGCATTAGACAAACTGAGTCGATAATAGCATGCTTATATTACATGGTTCATTTATAAAAGGAGTACGAACGAAAGTAACGCATGACGATTaa                                >  2:250542/1‑141 (MQ=255)
                                                                                                   catGTAGAAACTTTTGAAGAAGTTGAAGAGGAAGAAGTGCCAAAGAAAGCATTAGAGAAACTGATGCGAGAAGAGCATGCATAGATGACATGGTTCATTTATAAAAGTAGTACGAACGAAGGTAACTCATGACGAGTAAtt                              >  1:312845/1‑141 (MQ=255)
                                                                                                     tGTAGAAACTTTTGAAGAAGTTGAAAAGGAAGAAGTGACAAAAAAATAAATAAACAAACTGAGTCGATAATAGCATGCTTATATTACATGGTTCATTTATAAAAAGAGTACGAACGCAAGTAACGCATTAcaaaaaaaaaa                            <  2:201965/141‑11 (MQ=255)
                                                                                                               tttGAAGAAGTTGAAAAGGAAGAAGTGCCAAAGAAAGCATTAGACAAACTGAGTCGATAATAGCATGCTTATATTACATGGTTCATTTATAAAAGGAGTACGAACGAAAGTAACGCATGACGATTAATTTAAAAATATTTg                  >  1:360024/1‑141 (MQ=255)
                                                                                                                ttGAAGAAGTTGAAAAGGAAGAAGTGCCAAAGAAAGCATTAGACAAACTGAGTCGATAATAGCATGCTTATATTACATGGTTCATTTATAAAAGGAGTACGAACGAAAGTAACGCATGACGATTAATTTAAAAATATTTGt                 <  2:254126/141‑1 (MQ=255)
                                                                                                                   aagaagTTGAAAAGGAAGAAGTGCCAAAGAAAGCATTAGACAAACTGAGTCGATAATAGCATGCTTATATTACATGGTTCATTTATAAAAGGAGTACGAACGAAAGTAACGCATGACGATTAATTTAAAAATATTTGtaat              >  2:246288/1‑141 (MQ=255)
                                                                                                                        gTTGAAAAGGAAGAAGTGCCAAAGAAAGCATTAGACAAACTGAGTCGATAATAGCATGCTTATATTACATGGTTCATTTATAAAAGGAGTACGAACGAAAGTAACGCATGACGATTAATTTAAAAATATTTGTAATAATTa         <  1:177884/141‑1 (MQ=255)
                                                                                                                        gTTGAAAAGGAAGAAGTGCCAAAGAAAGCATTAGACAAACTGAGTCGATAATAGCATGCTTATATTACATGGTTCATTTATAAAAGGAGTACGAACGAAAGTAACGCATGACGATTAATTTAAAAATATTTGTAATAATTa         <  2:124504/141‑1 (MQ=255)
                                                                                                                         ttGAAAAGGAAGAAGTGCCAAAGAAAGCATTAGACAAACTGAGTCGATAATAGCATGCTTATATTACATGGTTCATTTATAAAAGGAGTACGAACGAAAGTAACGCATGACGATTAATTTAAAAATATTTGTAATAATTAt        >  1:242084/1‑141 (MQ=255)
                                                                                                                            aaaaGGAAGAAGTGCCAAAGAAAGCATTAGACAAACTGAGTCGATAATAGCATGCTTATATTACATGGTTCATTTATAAAAGGAGTACGAACGAAAGTAACGCATGACGATTAATTTAAAAATATTTGTAATAATTATGaa     >  2:114711/1‑141 (MQ=255)
                                                                                                                             aaaGGAAGAAGTGCCAAAGAAAGCATTAGACAAACTGAGTCGATAATAGCATGCTTATATTACATGGTTCATTTATAAAAGGAGTACGAACGAAAGTAACGCATGACGATTAATTTAAAAATATTTGTAATAATTATGAAt    <  1:196018/141‑1 (MQ=255)
                                                                                                                               aGGAAGAAGTGCCAAAGAAAGCATTAGACAAACTGAGTCGATAATAGCATGCTTATATTACATGGTTCATTTATAAAAGGAGTACGAACGAAAGTAACGCATGACGATTAATTTAAAAATATTTGTAATAATTATGAATaa  <  2:101454/141‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                  
AGTGTTTTTAATAAAGAAGGGAAAGAGACGCAATTAATGTATACTGCTACATTTGATGTTAAACCGCATAGATACTTGAAAATTACACATAAAGGTCATCATGTAGAAACTTTTGAAGAAGTTGAAAAGGAAGAAGTACCAAAGAAAGCATTAGACAAACTGAGTCGATAATAGCATGCTTATATTACATGGTTCATTTATAAAAGGAGTACGAACGAAAGTAACGCATGACGATTAATTTAAAAATATTTGTAATAATTATGAATAA  >  CP000730/423110‑423377

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: ‑