Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 1,233,504 T→A 100% V470D (GTT→GAT)  USA300HOU_1157 → non‑specific serine/threonine protein kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007301,233,5040TA100.0% 143.9 / ‑6.5 46V470D (GTT→GAT) USA300HOU_1157non‑specific serine/threonine protein kinase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base A (27/18);  minor base C (1/0);  total (28/18)
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GAACGTGTTGAACGTGGTGACAGTGTTGATGTTGTTATATCAAAGGGCCCTGAAAAGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGTTACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACTGAAATCGCTATTCATGATTCTAATATTAAACTATATGAATCTTTAGGCATTAAGCAAGTT  >  CP000730/1233371‑1233637
                                                                                                                                     |                                                                                                                                     
gAACGTGTTGAACGTGGTGACAGTGTTGATGTTGTTATATCAAAGGGCCCTGAAAAGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGAtt                                                                                                                                >  1:1176/1‑141 (MQ=255)
  aCGTGTTGAACGTGGTGACAGTGTTGATGTTGTTATATCAAAGGGCCCTGAAAAGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGa                                                                                                                              >  1:541483/1‑141 (MQ=255)
   cGTGTTGAACGTGGTGAAAGTGTTGATGTTGTTATATCAAAGGGCGCAGAAAAGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGGAGCCTTGCAGAAATTAAAAGCGAGAGGACATAGAGATGCTACggaggaa                                                                                                                             >  2:340419/1‑135 (MQ=255)
              tggtgACAGTGTTGATGTTGTTATATCAAAGGGCCCTGAAAAGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTCTctaa                                                                                                                  >  1:367023/1‑141 (MQ=255)
                   aCAGTGTTGATGTTGTTAAATCAAAGGGCCCTGAAAAGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATc                                                                                                             <  2:54445/141‑1 (MQ=255)
                                   tataTCAAAGGGCCCTGAAAAGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATa                                                                                             <  2:18395/141‑1 (MQ=255)
                                    atatCAAAGGGCCCTGAATAGGTTAAAATGCCCAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATAc                                                                                            <  2:52649/141‑1 (MQ=255)
                                    atatCAAAGGGCCCTGAAAAGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATAc                                                                                            >  1:459344/1‑141 (MQ=255)
                                     tcaccaaggaaccTAAAAAGGTAAAAATGACAAAGGACAATGGTTGACCTAACGAGGAAGCCCTGCAAAAATTAAAATCGTTATTTCTTAAAGAAAAAAAAAAAAAAAAAGTAGATAATAATCAAGCGCCAAAAGGATACa                                                                                           <  2:263137/130‑1 (MQ=255)
                                       tCAAGGGGCCCTGAAAAGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACAtt                                                                                         >  2:21450/1‑141 (MQ=255)
                                       tCAAAGGGCCCTGAAAAGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCATTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACAtt                                                                                         >  2:153205/1‑141 (MQ=255)
                                        cAAAGGGCCCTGAAAAGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGATGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATTAAGCGCCACAAGGATACATTg                                                                                        <  1:297490/141‑1 (MQ=255)
                                               cccTGAAAAGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCa                                                                                 >  1:226429/1‑141 (MQ=255)
                                                ccTGAAAAGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGCTACATTGCAAATCaa                                                                                >  2:199743/1‑141 (MQ=255)
                                                ccTGAAAAGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCaa                                                                                >  1:466022/1‑141 (MQ=255)
                                                  tGAAAAGGTTAAAATGGCAAATGTCATTGGTTTAACTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCATCCAAATTACATTGCAAATCAAAg                                                                              <  2:550314/141‑1 (MQ=255)
                                                             aaaTGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAACTTGGTACGGCTGAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCaa                                                                   <  1:131989/141‑1 (MQ=255)
                                                             aaaTGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCaa                                                                   <  1:241096/141‑1 (MQ=255)
                                                             aaaTGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCaa                                                                   >  1:125884/1‑141 (MQ=255)
                                                                  ccAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACt                                                              >  1:344630/1‑141 (MQ=255)
                                                                  ccAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACt                                                              >  1:28600/1‑141 (MQ=255)
                                                                      aggttaTTGGTTTACCTAAGGGGGAAGGCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAACCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACTGaaa                                                          <  2:147016/136‑1 (MQ=255)
                                                                        gTCATTGGTTTACCGAAGGCGGAAGACTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACTGAAATc                                                        >  1:313590/1‑141 (MQ=255)
                                                                           aTTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTCTAACCGAAAATACTGAAATCGCt                                                     >  2:166168/1‑141 (MQ=255)
                                                                             tGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACTGAAATCGCTAt                                                   >  1:237685/1‑141 (MQ=255)
                                                                             tGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACTGAAATCGCTAt                                                   <  1:83873/141‑1 (MQ=255)
                                                                                tttACTTAATGGGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGTTACATTGCAAATCAAAGTGTAACCGCAAATACTGAAATCGCTATTCa                                                <  1:196059/141‑1 (MQ=255)
                                                                                       aaGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACTGAAATCGCTATTCATGATTCt                                         <  1:160991/141‑1 (MQ=255)
                                                                                         ggaggaAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACTGAAATCGCTGTTCATGATTCTaa                                       >  2:56457/1‑141 (MQ=255)
                                                                                          gaggaAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACTGAAATCGCTATTCATGATTCTAat                                      <  2:9002/141‑1 (MQ=255)
                                                                                          gaggaAGCCTTGCAGAAATTAAAAGCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAATGATACATTGCAAACCAAAGTGTAACCGCAAATACTGAAATCGCTATTCATGATTCTAat                                      <  1:250504/141‑1 (MQ=255)
                                                                                              aaGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACTGAAATCGCTATTCATGATTCTAATATTa                                  >  2:506177/1‑141 (MQ=255)
                                                                                                 ccTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACTGAAATCGCTATTCATGATTCTAATATTAAAc                               >  1:186189/1‑141 (MQ=255)
                                                                                                      cAGAAATTAAAACTGTTAGGTCTTAAAGATGATACGAATGAAACAGTATATATAAATCAAGCGCCAAAAGGATACATTTCAAATCAAAGTGTAACCGCAAATACTGAAATCGCTATTCATGATTCTAATATTAAACTatat                          <  2:15808/141‑1 (MQ=255)
                                                                                                       aGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACTGAAATCGCTATTCATGATTCTAATATTAAACTATATg                         <  1:368316/141‑1 (MQ=255)
                                                                                                           aTTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACTGAAATCGCTATTCATGATTCTAATATTAAACTATATGAATc                     >  2:155136/1‑141 (MQ=255)
                                                                                                           aTTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACTGAAATCGCTATTCATGATTCTAATATTAAACAATATGAATc                     >  2:491534/1‑141 (MQ=255)
                                                                                                              aaaaTCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACTGAAATCGCTATTCATGATTCTAATATTCAACTATATGAATCttt                  <  1:339369/141‑1 (MQ=255)
                                                                                                                aaTCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACTGAAATCGCTATTCATGATTCTGATATTAAACTATATGAATATTTAg                >  2:357011/1‑141 (MQ=255)
                                                                                                                  tCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATACTCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACTGAAATAGCTATTCATGATTCTAATATTAAACTGTATGAATCTTTAGGc              >  2:202931/1‑141 (MQ=255)
                                                                                                                  tCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACTGAAATCGCTATTCATGATTCTGATATTAAACTATATGAATCTTTAgca              >  2:7562/1‑139 (MQ=255)
                                                                                                                    gTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACTGAAATCGCTATTCATGATTCTAATATTAAACGATATGAATCTTTAGGCAt            <  1:420897/141‑1 (MQ=255)
                                                                                                                    gTTAGGTCTTAAAGAGGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACTGAAATAGCTATTCATGATTCTAATATTAAACTATATGAATCTTTAGGCAt            >  1:423126/1‑141 (MQ=255)
                                                                                                                        ggTCTTAAAGATGATACGATTGAAAGAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACTGAAATCGCTATTCATGATTCTAATATTAAACTATATGAATCTTTAGGCATTAAg        >  1:21679/1‑141 (MQ=255)
                                                                                                                        ggTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACTGAAATCGCTATTCATGATTCTAATATTAAACTATATGAATCTTTAGGCATTAAg        >  1:12209/1‑141 (MQ=255)
                                                                                                                              aaaGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAGTACTGAAATCGCTATTCATGATTCTAATATTAAACTATATGAATCTTTAGGCATTAAGCAAGtt  >  1:219079/1‑141 (MQ=255)
                                                                                                                                     |                                                                                                                                     
GAACGTGTTGAACGTGGTGACAGTGTTGATGTTGTTATATCAAAGGGCCCTGAAAAGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGTTACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACTGAAATCGCTATTCATGATTCTAATATTAAACTATATGAATCTTTAGGCATTAAGCAAGTT  >  CP000730/1233371‑1233637

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: