Predicted mutation
evidence seq id position mutation annotation gene description
RA CP000730 2,026,955 G→A T331I (ACA→ATA)  vraS ← sensor histidine kinase VraS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007302,026,9550GA96.3% 81.7 / NA 27T331I (ACA→ATA) vraSsensor histidine kinase VraS
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (18/8);  minor base C (1/0);  total (19/8)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACTACTTCAATATCACTTTGCGTTGATAGATAACTTGAAATTCCTATACGTACCATTTCATGATCATCCACAAACAATACTTTAATCGTCATACGAATCCTCCTTATTTAAAGGTGCTTTCACCTCGATACGTGTACCTGAATCTGGCAATGATACAATATGGAACGTTGCACCAATTTCCAAAGCTCTTTCACGCATATTTTTAAGTCCATAACTTTGTTCTAATTTTTCATCAACATTAAAACCTTTACCATTATCTTGAATTCT  >  CP000730/2026822‑2027088
                                                                                                                                     |                                                                                                                                      
actactTCAATATCACTTTGCGTTGATAGATAACTTGAAATTCCTATACGTACCATTTCATGATCATCCACAAACAATACTTTAATCGTCATACGAATCCTCCTTATTTAAAGGTGCTTTCACCTCGATACGTATACCTGa                                                                                                                                 <  1:20455/141‑1 (MQ=255)
      tCAATATCACTTTGCGTTGATAGATAACTTGAAATTCCTATACGTACCATTTCATGATCATCCACAAACAATACTTTAATCGTCATACGAATCCTCCTTATTTAAAGGTGCTTTCACCTCGATACGTATACCTGAATCTgg                                                                                                                           <  2:170541/141‑1 (MQ=255)
      tCAATATCACTTTGCGTTGATAGATAACTTGAAATTCCTATACGTACCATTTCATGATCATCCACAAACAATACTTTAATCGTCATACGAATCCTCCTTATTTAAAGGTGCTTTCACCTCGATACGTATACCTGAATCTgg                                                                                                                           >  1:205686/1‑141 (MQ=255)
                     gTTGATAGATAACTTGAAATTCCTATACGTACCATTTCATGATCATCCACAAACAATACTTTAATCGTCATACGAATCCTCCTTATTTAAAGGTGCTTTCACCTCGATACGTATACCTGAATCTGGCAATGATACAATATg                                                                                                            >  1:158931/1‑141 (MQ=255)
                     gTTGATAGATAACTTGAAATTCCTATACGTACCATTTCATGATCATCCACAAACAATACTTTAATCGTCATACGAATCCTCCTTATTTAAAGGTGCTTTCACCTCGATACGTATACCTGAATCTGGCAATGATACAATATg                                                                                                            >  1:99014/1‑141 (MQ=255)
                      ttGATAGATAACTTGAAATTCCTATACGTACCATTTCATGATCATCCACAAACAATACTTTAATCGTCATACGAATCCTCCTTATTTAAAGGTGCTTTCACCTCGATACGTATACCTGAATCTGGCAATGATACAATATgg                                                                                                           <  1:273736/141‑1 (MQ=255)
                               aaCTTGAAATTCCTATACGTACCATTTCATGATCATCCACAAACAATACTTTAATCGTCATACGAATCCTCCTTATTTAAAGGTGCTTTCACCTCGATACGTATACCTGAATCTGGCAATGATACAATATGGAACGTTGCa                                                                                                 <  2:198803/141‑1 (MQ=255)
                                  ttGAAATTCCTATACGTACCATTTCATGATCATCCACAAACAATACTTTAATCGTCATACGAATCCTCCTTATTTAAAGGTGCTTTCACCTCGATACGTATACCTGAATCTGGCAATGATACAATATGGAACGTTGCACCa                                                                                              >  1:157810/1‑141 (MQ=255)
                                     aaaTTCCTATACGTACCATTTCATGATCATCCACAAACAATACTTTAATCGTCATACGAATCCTCCTTATTTAAAGGTGCTTTCACCTCGATACGTATACCTGAATCTGGCAATGAAACATCTGGTAACGTTGCACCAAtt                                                                                           >  1:177937/1‑141 (MQ=255)
                                      aaTTCCTATACGTACCATTTCATGATCATCCACAAACAATACTTTAATCGTCATACGAATCCTCCTTATTTAAAGGTGCTTTCACCTCGATACGTATACCTGAATCTGGCAATGATACAATATGGAACGTTGCACCAAttt                                                                                          >  1:153670/1‑141 (MQ=255)
                                                cGTACCATTTCATGATCATCCACAAACAATACTTTAATCGTCATACGAATCCTCCTTATTTAAAGGTGCTTTCACCTCGATACGTATACCt                                                                                                                                   <  2:94999/91‑1 (MQ=255)
                                                cGTACCATTTCATGATCATCCACAAACAATACTTTAATCGTCATACGAATCCTCCTTATTTAAAGGTGCTTTCACCTCGATACGTATACCt                                                                                                                                   >  1:94999/1‑91 (MQ=255)
                                                    ccATTTCATGATCATCCACAAACAATACTTTAATCGTCATACGAATCCTCCTTATTTTATGTTGTTCTGGCCCGGATTCGGATGCCTGAATCTGGCAATGATACAATATGGAACGTTGCACCAATTTCCAAAGCTCTTTCa                                                                            >  1:155659/1‑141 (MQ=255)
                                                    ccATTTCATGATCATCCACAAACAATACTTTAATCGTCATACGAATCCTCCTTATTTAAAGGTGCTTTCACCTCGATACGTATACCTGAATCTGGCAATGATACAATATGGAACGTTGCACCAATTTCCAAAGCTCTTTCa                                                                            >  1:321891/1‑141 (MQ=255)
                                                      aTTTCATGATCATCCACAAACAATACTTTAATCGTCATACGAATCCTCCTTATTTAAAGGTGCTTTCACCTCGATACGTATACCTGAATCTGGCAATGATACAATATGGAACGTTGCACCAATTTCCAAAGCTCTTTCACg                                                                          <  1:66523/141‑1 (MQ=255)
                                                          cATGATCATCCACAAACAATACTTTAATCGTCATACGAATCCTCCTTATTTAAAGGTGCTTTCACATCGATACGTATACCTGAATCTGGCATTGATACAATATGGATCGTTGCACCAATTTCTAAAGCTCTTTCACGCcta                                                                      >  1:286785/1‑141 (MQ=255)
                                                           aTGATCATCCACAAACAATACTTTAATCGTCATACGAATCCTCCTTATTTAAAGGTGCTTTCACCTCGATACGTATACCTGAATCTGGCAATGATACAATATGGAACGTTGCACCAATTTCCAAAGCTCTTTCACGCatat                                                                     >  1:278020/1‑141 (MQ=255)
                                                                              aCTTTAATAGACATACGAATCCTCCTTATTTAAAGGTGATTTCACATCGCTCCGTCTACCTGAATCTGGCAATGATACAATATGGAACGTTGCACCAATTTCCAAAGCTCTTTCACGCATATTTTTAAGTCAATCACttta                                                  >  2:293379/1‑140 (MQ=255)
                                                                                    aTCGTCATACGAATCCTCCTTATTTAAAGGTGCTTTCACCTCGATACGTATACCTGAATCTGGCAATGATACAATATGGAACGTTGCACCAATTTCCAAAGCTCTTTCACGCATATTTTTAAGTCCATAACTTTGTTCTaa                                            >  1:254755/1‑141 (MQ=255)
                                                                                         cATACGAATCCTCCTTATTTAAAGGTGCTTTCACCTCGATACGTATACCTGAATCTGGCAATGATACAATATGGAACGTTGCACCAATTTCCAAAGCTCTTTCACGCATATTTTTAAGTCCATAACTTTGTTCTAAttttt                                       >  1:10628/1‑141 (MQ=255)
                                                                                         cATACGAATCCTCCTTATTTAAAGGTGCTTTCACCTCGATACGTATACCTGAATCTGGCAATGATACAATATGGAACGTTGCACCAATTTCCAAAGCTCTTTCACGCATATTTTTAAGTCCATAACTTTGTTCTAAttttt                                       >  1:246420/1‑141 (MQ=255)
                                                                                            aCGAATCCTCCTTATTTAAAGGTGCTTTCACGTCGATACGTATACCTGAATCTGGCAATGATACAATATGGAACGTTGCACCAATTTCCAAAGCTCTTTCACGCATCTTTTTAAGTCCATAACTTTGTTCTAATTTTtcaa                                    >  2:174479/1‑140 (MQ=255)
                                                                                               agtcctccTTGTTTAAATGTGCGTTCTCCTCGATACGTATACCTGAATCGCGAAATGATACAATATGGAACCTTGCACCAATTTCCCAAGCCCTTTCACGCATATTTTTAAGTCCATAACTATTTTCTAATTTTtcatcca                                 <  2:137103/139‑3 (MQ=255)
                                                                                                     cctTATTTAAAGGTGCTTTCACCTCGATACGTATACCTGAATCTGGCAATGATACAATATGGAACGTTGCACCAATTTCCAAAGCTCTTTCACGCATATTTTTAAGTCCATAACTTTGTTCTAATTTTTCATCAACATTaa                           <  1:342059/141‑1 (MQ=255)
                                                                                                             aaaGGTGCTTTCACCTCGATACGTATACCTGAATCTGGCAATGATACAATATGGAACGTTGCACCAATTTCCAAAGCTCTTTCACGCATATTTTTAAGTCCATAACTTTGTTCTAATTTTTCATCAACATTAAAACCTTTa                   >  1:58008/1‑141 (MQ=255)
                                                                                                                   gCTTTCACCTCGATACGTATACCTGAATCTGGCAATGATACAATATGGAACGTTGCACCAATTTCCAAAGCTCTTTCACGCATATTTTTAAGTCCATAACTTTGTTCTAATTTTTCATCAACATTAAAACCTTTACCATTa             >  1:14856/1‑141 (MQ=255)
                                                                                                                              gATACGTATACCTGAATCTGGCAATGATACAATATGGAACGTTGCACCAATTTCCAAAGCTCTTTCACGCATATTTTTAAGTCCATAACTTTGTTCTAATTTTTCATCAACATTAAAACCTTTACCATTATCTTGAATtct  >  2:266706/1‑141 (MQ=255)
                                                                                                                                     |                                                                                                                                      
ACTACTTCAATATCACTTTGCGTTGATAGATAACTTGAAATTCCTATACGTACCATTTCATGATCATCCACAAACAATACTTTAATCGTCATACGAATCCTCCTTATTTAAAGGTGCTTTCACCTCGATACGTGTACCTGAATCTGGCAATGATACAATATGGAACGTTGCACCAATTTCCAAAGCTCTTTCACGCATATTTTTAAGTCCATAACTTTGTTCTAATTTTTCATCAACATTAAAACCTTTACCATTATCTTGAATTCT  >  CP000730/2026822‑2027088

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: