Predicted mutation
evidence seq id position mutation annotation gene description
RA CP000730 2,027,751 C→G V66L (GTA→CTA)  vraS ← sensor histidine kinase VraS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007302,027,7510CG96.0% 82.2 / ‑2.9 25V66L (GTA→CTA) vraSsensor histidine kinase VraS
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base G (14/10);  total (15/10)
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.87e-01

AGTGTATGGTATAAATCTAACGTTTCACTATATATTTCTATATTTTGATCATTAATGCCAACTGTTTCGCCTTCCATTGAACGCTCAATTTGCGTCTTAATCCAATCATTTTGCTGATTGATTTTGTAAGCGAGTACCGAACCAACAATAATACACAATAATATGATGATGAGATTTAAAAATAAAAAAACTGGTATTCCGAATATTTGTGTATAAAACATACCTTGAAAATAGATGATATTTACAAAAACTTTATCGATGAACAGAAATG  >  CP000730/2027615‑2027885
                                                                                                                                        |                                                                                                                                      
aGTGTATGGTATAAATCTAACGTTTCACTATATATTTCTATATTTTGATCATTAATGCCAACTGTTTCGCCTTCCATTGAACGCTCAATTTGCGTCTTAATCCAATCATTTTGCTGATTGATTTTGTAATCGAGTAGCGaa                                                                                                                                    >  2:282605/1‑141 (MQ=255)
   gtATGGTATAAATCTAACGTTTCACTATATATTTCTATATTTTGATCATTAATGCCAACTGTTTCGCCTTCCATTGAACGCTCAATTTGCGTCTTAATCCAATCATTTTGCTGATTGATTTTGTAAGCGAGTAGCGAACca                                                                                                                                 <  2:316341/141‑1 (MQ=255)
                cTAACGTTTCACTATATATTTCTATATTTTGATCATTAATGCCAACTGTTTCGCCTTCCATTGAACGCTCAATTTGCGTCTTAATCCAATCATTTTGCTGATTGATTTTGTAAGCGAGTAGCGAACCAACAATAATAcaca                                                                                                                    <  1:259771/141‑1 (MQ=255)
                  aaCGTTTCACTATATATTTCTATATTTTGATCATTAATGCCAACTGTTTCGCCTTCCATTGAACGCTCAATTTGCGTCTTAATCCAATCATTTTGCTGATTGATTTTGTAAGCGAGTAGCGAACCAACAATAATACACaat                                                                                                                  >  2:310207/1‑141 (MQ=255)
                    cGTTTCACTATATATTTCTATATTTTGATCATTAATGCCAACTGTTTCGCCTTCCATTGAACGCTCAATTTGCGTCTTAATCCAATCATTTTGCTGATTGATTTTGTAAGCGAGTAGCGAACCAACAATAATACACaataa                                                                                                                <  1:244856/141‑1 (MQ=255)
                       ttCACTATATATTTCTATATTTTGATCATTAATGCCAACTGTTTCGCCTTCCATTGAACGCTCAATTTGCGTCTTAATCCAATCATTTTGCTGATTGATTTTGTAAGCGAGTAGCGAACCAACAATAATACACAATAatat                                                                                                             >  2:225276/1‑141 (MQ=255)
                       ttCACTATATATTTCTATATTTTGATCATTAATGCCAACTGTTTCGCCTTCCATTGAACGCTCAATTTGCGTCTTAATCCAATCATTTTGCTGATTGATTTTGTAAGCGAGTAGCGAACCAACAATAATACACAATAatat                                                                                                             >  1:105118/1‑141 (MQ=255)
                             atatatTTCTATATTTTGATCATTAATGCCAACTGTTTCGCCTTCCATTGAACGCTCAATTTGCGTCTTAATCCAATCATTTTGCTGATTGATTTTGTAAGCGAGTAGCGAACCAACAATAATACACAATAATatgatgat                                                                                                       <  1:65108/141‑1 (MQ=255)
                                 atTTCTATATTTTGATCATTAATGCCAACTGTTTCGCCTTCCATTGAACGCTCAATTTGCGTCTTAATCCAATCATTTTGCTGATTGATTTTGTAAGCGAGTAGCGAACCAACAATAATACACAATAATATGATGATgaga                                                                                                   >  1:176705/1‑141 (MQ=255)
                                   ttCTATATTTTGATCATTAATGCCAACTGTTTCGCCTTCCATTGAACGCTCAATTTGCGTCTTAATCCAATCATTTTGCTGATTGATTTTGTAAGCGAGTAGCGAACCAACAATAATACACAATAATATGATGATGAGAtt                                                                                                 >  2:34134/1‑141 (MQ=255)
                                   ttCTATATTTTGATCATTAATGCCAACTGTTTCGCCTTCCATTGAACGCTCAATTTGCGTCTTAATCCAATCATTTTGCTGATTGATTTTGTAAGCGAGTAGCGAACCAACAATAATACACAATAATATGATGATGAGAtt                                                                                                 <  2:30210/141‑1 (MQ=255)
                                                   ttAATGCCAACTGTTTCGCCTTCCATTGAACGCTCAATTTGCGTCTTAATCCAATCATTTTGCTGATTGATTTTGTAAGCGAGTAGCGAACCAACAATAATACACAATAATATGATGATGAGATTTAAAAATAAAAAAACt                                                                                 >  1:332651/1‑141 (MQ=255)
                                                             cTGTTTCGCCTTCCATTGAACGCTCAATTTGCGTCTTAATCCAATCATTTTGCTGATTGATTTTGTAAGCGAGTAGCGAACCAACAATAATACACAATAATATGATGATGAGATTTAAAAATAAAAAAACTGGTATTCCGa                                                                       >  1:87133/1‑141 (MQ=255)
                                                                                aCGCTCAATTTGCGTCTTAATCCAATCATTTTGCTGATTGATTTTGTATGCGAGTACCGAACCAACAATAATACACAATAATATGATGATGAGATTTAAAAATAAACAGACTGGTATTCCGAATATTTGTGTATAAAACAt                                                    >  2:283875/1‑141 (MQ=255)
                                                                                aCGCTCAATTTGCGTCTTAATCCAATCATTTTGCTGATTGATTTTGTAAGCGAGTAGCGAACCAACAATAATACACAATAATATGATGATGAGATTTAAAAATAAAAAAACTGGTATTCCGAATATTTGTGTATAAAACAt                                                    <  1:135636/141‑1 (MQ=255)
                                                                                        tttGCGTCTTAATCCAATCATTTTGCTGATTGATTTTGTAAGCGAGTAGCGAACCAACAATAATACACAATAATATGATGATGAGATTTAAAAATAAAAAAACTGGTATTCCGAATATTTGTGTATAAAACATACCTTGaa                                            >  1:342829/1‑141 (MQ=255)
                                                                                        tttGCGTCTTAATCCAATCATTTTGCTGATTGATTTTGTAAGCGAGTAGCGAACCAACAATAATACACAATAATATGATGATGAGATTTAAAAATAAAAAAACTGGTATTCCGAATATTTGTGTATAAAACATACCTTGaa                                            >  2:211422/1‑141 (MQ=255)
                                                                                                  aaTCCAATCATTTTGCTGATTGATTTTGTAAGCGAGTAGCGAACCAACAATAATACACAATAATATGATGATGAGATTTAAAAATAAAAAAACTGGTATTCCGAATATTTGTGTATAAAACATACCTTGAAAATAgatgat                                  >  2:97158/1‑141 (MQ=255)
                                                                                                  aaTCCAATCATTTTGCTGATTGATTTTGTAAGCGAGTAGCGAACCAACAATAATACACAATAATATGATGATGAGATTTAAAAATAAAAAAACTGGTATTCCGAATATTTGTGTATAAAACATACCTTGAAAATAgatgat                                  >  2:118931/1‑141 (MQ=255)
                                                                                                      cAATCATTTTGCTGATTGATTTTGTAAGCGAGTAGCGAACCAACAATAATACACAATAATATGATGATGAGATTTAAAAATAAAAAAACTGGTATTCCGAATATTTGTGTATAAAACATACCTTGAAAATAGATGATAttt                              >  1:138887/1‑141 (MQ=255)
                                                                                                          cATTTTGCTGATTGATTTTGTAAGCGAGTAGCGAACCAACAATAATACACAATAATATGATGATGAGATTTAAAAATAAAAAAACTGGTATTCCGAATATTTGTGTATAAAACATACCTTGAAAATAGATGATATTTACaa                          <  1:143033/141‑1 (MQ=255)
                                                                                                                gCTGATTGATTTTGTAAGCGAGTAGCGAACCAACAATAATACACAATAATATGATGATGAGATTTAAAAATAAAAAAACTGGTATTCCGAATATTTGTGTATAAAACATACCTTGAAAATAGATGATATTTACAAAAACtt                    <  2:69173/141‑1 (MQ=255)
                                                                                                                          tttGTAAGCGAGTAGCGAACCAACAATAATACACAATAATATGATGATGAGATTTAAAAATAAAAAAACTGGTATTCCGAATATTTGTGTATAAAACATACCTTGAAAATAGATGATATTTACAAAAACTTTATCGATGaa          <  1:310207/141‑1 (MQ=255)
                                                                                                                           ttGTAAGCGAGTAGCGAACCAACAATAATACACAATAATATGATGATGAGATTTAAAAATAAAAAAACTGGTATTCCGAATATTTGTGTATAAAACATACCTTGAAAATAGATGATATTTACAAAAACTTTATCGATGAAc         <  1:211422/141‑1 (MQ=255)
                                                                                                                                  cGAGTAGCGAACCAACAATAATACACAATAATATGATGATGAGATTTAAAAATAAAAAAACTGGTATTCCGAATATTTGTGTATAAAACATACCTTGAAAATAGATGATATTTACAAAAACTTTATCGATGAACAGAAATg  >  1:243591/1‑141 (MQ=255)
                                                                                                                                        |                                                                                                                                      
AGTGTATGGTATAAATCTAACGTTTCACTATATATTTCTATATTTTGATCATTAATGCCAACTGTTTCGCCTTCCATTGAACGCTCAATTTGCGTCTTAATCCAATCATTTTGCTGATTGATTTTGTAAGCGAGTACCGAACCAACAATAATACACAATAATATGATGATGAGATTTAAAAATAAAAAAACTGGTATTCCGAATATTTGTGTATAAAACATACCTTGAAAATAGATGATATTTACAAAAACTTTATCGATGAACAGAAATG  >  CP000730/2027615‑2027885

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