Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,425,301 G→C 13.8% R290T (AGG→ACG)  yegS → lipid kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,425,3010GC13.8% 97.1 / 6.7 35R290T (AGG→ACG) yegSlipid kinase
Reads supporting (aligned to +/- strand):  ref base G (13/17);  new base C (3/2);  total (16/19)
Fisher's exact test for biased strand distribution p-value = 6.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGGGGCTCTGCGGGACCTGTTGGGCGGTGACGGCCTGTTCGTTCGGGCCAGGTTGCCATGGGTCGAGATCAAGAGTTCGCAAGGGTTGGACATCAACCTCGATGGCGAGCCGTTGCAGGCTGCCAACCTGCGCTTTCAGGCAAGGCCGGCGGCGTTGCACCTGCACCTGCCGGCCGGGTCACCGTTGCTCAGTCATCCAGGCTGATGATTTTTTCGCGTACGGCGAACAGCACCAGGCCGGCTACGTCGTGAATCTGCAGGCGCTTCATGATTTGCGAGCGGTGAGTTTCCA  >  NC_002947/2425158‑2425449
                                                                                                                                               |                                                                                                                                                    
gggggCTCTGCGGGACCTGTTGGGCGGTGACGGCCTGTTCGTTCGGGCCAGGTTGCCATGGGTCGAGATCAAGAGTTCGCAAGGGTTGGACATCAACCTCGATGGCGAGCCGTTGCAGGCTGCCAACCTGCGCTTTCAGGCAAGGCcggc                                                                                                                                                <  1:1100155/150‑1 (MQ=255)
      tctGCGGGACGTGTTGGACGGTGACGGCCTGTTCGTTCGGGCCGGGTTGCTATGAGAGGAGGTCAAGAGTTCGCAAGGGTTGGACATCAACGTCGATGGCGCGCCGTTGCAGGCTGCCGGCCTGCGCTTTCAGGCAACGCGGgcggcggt                                                                                                                                          >  2:485393/1‑148 (MQ=255)
          cGGGACCTGTTGGGCGGTGACGGCCTGTTCGTTCGGGCCAGGTTGCCATGGGTCGAGATCAAGAGTTCGCAAGGGTTGGACATCAACCTCGATGGCGAGCCGTTGCAGGCTGCCAACCTGCGCTTTCAGGCAACGCCGCCGGCCTtgcct                                                                                                                                      >  2:905573/1‑148 (MQ=255)
          cGGGACCTGTTGGGCGGTGACGGCCTGTTCGTTCGGGCCAGGTTGCCATGGGTCGAGATCAAGAGTTCGCAAGGGTTGGACATCAACCGCGATGGCGCGCCGTTGCAGGCTGCCAACCTGCGCTTTCAGGCAAGGCCGGCGGCGTtgcac                                                                                                                                      <  1:485393/150‑1 (MQ=255)
              aCCTGTTGGGCGGTGACGGCCTGTTCGTTCGGGCCAGGTTGCCATGGGTCGAGATCAAGAGTTCGCAAGGGTTGGACATCAACCTCGATGGCGAGCCGTTGCAGGCTGCCAACCTGCGCTTTCAGGCAACGCCGCCGGCCTtgcctgaaa                                                                                                                                  <  1:905573/150‑7 (MQ=255)
                                   ggttcgttcGGGCCAGGTTGCCATGGGTCGAGATCAAGAGTTCGCAAGGGTTGGACATCAACCTCGATGGCGAGCCGTTGCAGGCTGCCAACCTGCGCTTTCAGGAAAGGCCGGCGGCGTTGCACCTGCACCTGCCGGCCGGGTCACCGt                                                                                                             <  2:188262/149‑1 (MQ=255)
                                    gttcgttcGGGCCAGGTTGCCATGGGTCGAGATCAAGAGTTCGCAAGGGTTGGACATCAACCTCGATGGCGAGCCGTTGCAGGCTGCCAACCTGCGCTTTCAGGCAAGGCCGGCGGCGTTGCACCTGCACCTGCCGGCCGGGTCACCGtt                                                                                                            <  2:1205488/150‑1 (MQ=255)
                                         ttcgGGCCAGGTTGCCATGGGTCGAGATCAAGAGTTCGAAAGGGTTGGACATCAACCTCGATGGCGAGCCGTTGCAGGCTGCCAACCTGCGCTTTCAGGCAAGGCCGGCGGCGTTGCACCTGCACCTGCCGGCCGGGTCACCGTTGCTCa                                                                                                       >  2:697314/1‑150 (MQ=255)
                                               ccAGGTTGCCATGGGTCGAGATCAAGAGTTCGCAAGGGTTGGACATCAACCTCGATGGCGAGCCGTTGCAGGCTGCCAACCTGCGCTTTCAGGCAAGGCCGGCGGCGTTGCACCTGCACCTGccggccgg                                                                                                                     <  1:136259/130‑1 (MQ=255)
                                               ccAGGTTGCCATGGGTCGAGATCAAGAGTTCGCAAGGGTTGGACATCAACCTCGATGGCGAGCCGTTGCAGGCTGCCAACCTGCGCTTTCAGGCAAGGCCGGCGGCGTTGCACCTGCACCTGccggccgg                                                                                                                     >  2:136259/1‑130 (MQ=255)
                                                            ggTCGAGATCAAGAGTTCGCAAGGGTTGGACATCAACCTCGATGGCGAGCCGTTGCAGGCTGCCAACCTGCGCTTTCAGGCAAGGCCGGCGGCGTTGCACCTGCACCTGCCGGCCGGGTCACCGTTGCTCAGTCATCCAGGCTGATGAtt                                                                                    >  1:557688/1‑150 (MQ=255)
                                                                 agaTCAAGAGTTCGCAAGGGTTGGACATCAACCTCGATGGCGAGCCGTTGCAGGCTGCCAACCTGCGCTTTCAGGCAAGGCCGGCGGCGTTGCACCTGCACCTGCCGGCCGGGTCACCGTTGCTCAGTCATCCAGGCTGATGATTTTTTc                                                                               >  1:48560/1‑150 (MQ=255)
                                                                 agaTCAAGAGTTCGCAAGGGTTGGACATCAACCTCGATGGCGAGCCGTTGCAGGCTGCCAACCTGCGCTTTCAGGCAAGGCCGGCGGCGTTGCACCTGCACCTGCCGGCCGGGTCACCGTTGCTCAGTCATCCAGGCTGATGATTTTTTc                                                                               <  2:690325/150‑1 (MQ=255)
                                                                   aTCAAGAGTTCGCAAGGGTTGGACATCAACCTCGATGGCGAGCCGTTGCAGGCTGCCAACCTGCGCTTTCAGGCAAGGCCGGCGGTGTTGCACCTGCACCTGCCGGCCGGGTCACCGTTGCTCAGTCATCCAGGCTGATGATTTTTTcgc                                                                             <  1:562837/150‑1 (MQ=255)
                                                                      aaGAGTTCACAAGGGTTGGACATCAACCTCGATGGCGAGCCGTTGCAGGCTGCCAACCTGCGCTTTCAGGCAAGGCCGGCGGCGTTGCACCTGCACCTGCCGGCCGGGTCACCGTTGCTCAGTCATCCAGGCTGATGATTTTTTCGCGTa                                                                          >  2:208419/1‑150 (MQ=255)
                                                                                   ggTTGGACATCAACCTCGATGGCGAGCCGTTGCAGGCTGCCAACCTGCGCTTTCAGGCAAGGCCGGCGGCGTTGCACCTGCACCTGCCGGCCGGGTCACCGTTGCTCAGTCATCCAGGCTGATGATTTTTTCGCGTACGGCGAACAGCAc                                                             <  1:208419/150‑1 (MQ=255)
                                                                                          cATCAACCTCGATGGCGAGCCGTTGCAGGCTGCCAACCTGCGCTTTCAGGCAAGTCCGGCGGCGTTGCACCTGCACATGCCGGACGGGTCACCGTTGCTCAGATACCCATGCTGATCATTGTTTCGCGTTCTGCGAACAGCACCACGCCg                                                      >  2:264074/1‑150 (MQ=255)
                                                                                           aTCAACCTCGATGGCGAGCCGTTGCAGGCTGCCAACCTGCGCTTTCAGGCAAGGCCGGCGGCGTTGCACCTGCACCTACCGGCCGGGTCACCGTTGCTCAGTCATCCAGGCTGATGATTTTTTCGCGTACGGCGAACAGCACCAGGCCgg                                                     >  1:997169/1‑150 (MQ=255)
                                                                                                        gcgcgcCGTTGCAGGCTGCCAACCTGCGCTTTCAGGCAAGGCCGGCGGCGTTGCACCTGCACCTGCCGGCCGGGTCACCGTTGCTCAGTCATCCAGGCTGATGATTTTTTCGCGTACGGCGAACAGCACCAGGCCGGCTACGTCGTGAAt                                        <  1:49937/150‑1 (MQ=255)
                                                                                                             ccGTTGCAGGCTGCCAACCTGCGCTTTCAGGCAAGGCCGGCGGCGTTGCACCTGCACCTGCCGGCCGGGTCACCGTTGCTCAGTCATCCAGGCTGATGATTTTTTCGCGTACGGCGAACAGCACCAGGCCGGCTACGTCGTGAATCTGCa                                   <  2:48560/150‑1 (MQ=255)
                                                                                                                         gCCAACCTGCGCTTTCAGGCAAGGCCGGCGGCGTTGCACCTGCACCTGCCGGCCGGGTCACCGTTGCTCAGTCATCCAGGCTGATGATTTTTTCGCGTACGGCGAACAGCACCAGGCCGGCTACGTCGTGAATCTGCAGGCGCTTCATGa                       >  1:114661/1‑150 (MQ=255)
                                                                                                                            aaCCTGCGCTTTCAGGCAAGGCCGGCGGCGTTGCACCTGCACCTGCCGGCCGGGTCACCGTTGCTCAGTCATCCAGGCTGATGATTTTTTCGCGTACGGCGAACAGCACCAGGCCGGCTACGTCGTGAATCTGCAGGCGCTTCAt                         >  2:656226/1‑145 (MQ=255)
                                                                                                                            aaCCTGCGCTTTCAGGCAAGGCCGGCGGCGTTGCACCTGCACCTGCCGGCCGGGTCACCGTTGCTCAGTCATCCAGGCTGATGATTTTTTCGCGTACGGCGAACAGCACCAGGCCGGCTACGTCGTGAATCTGCAGGCGCTTCAt                         <  1:656226/145‑1 (MQ=255)
                                                                                                                            aaCCTGCGCTTTCAGGCAAGGCCGGCGGCGTTGCACCTGCACCTGCCGGCCGGGTCACCGTTGCTCAGTCATCCAGGCTGATGATTTTTTCGCGTACGGCGAACAGCACCAGGCCGGCTACGTCGTGAATCTGCAGGCGCTTCATGAttt                    >  2:930152/1‑150 (MQ=255)
                                                                                                                                tGCGCTTTCAGGCAAGGCCGGCGGCGTTGCACCTGCACCTGCCGGCCGGGTCACCGTTGCTCAGTCATCCAGGCtga                                                                                         >  2:928345/1‑77 (MQ=255)
                                                                                                                                ggcgcTTTCAGGCAAGGCCGGCGGCGTTGCACCTGCACCTGCCGGCCGGGTCACCGTTGCTCAGTCATCCAGGCtga                                                                                         <  1:928345/76‑1 (MQ=255)
                                                                                                                                    cTTTCAGGCAAGGCCGGCGGCGTTGCACCTGCACCTGCCGGCCGGGTCACCGTTGCTCAGTCATCCAGGCTGATGATTTTTTCGCGTACGGCGAACAGCACCAGGCCGGCTACGTCGTGAATCTGCAGGCGCTTCATGATTTGCGAGCgg            <  1:930152/150‑1 (MQ=255)
                                                                                                                                     tttCAGGCAAGGCCGGCGGCGTTGCACCTGCACCTGCCGGCCGGGTCACCGTTGCTCAGTCATCCAGGCTGATGATTTTTTCGCGTACGGCGAACAGCACCAGGCCGGCTACGTCGTGAATCTGCAGGCGCTTCATGATTTGCGAGCGGt           >  1:291537/1‑150 (MQ=255)
                                                                                                                                       tCAGGCAAGGCCGGCGGCGTTGCACCTGCACCTGCCGGCCGGGTCACCGTTGCTCAGTCATCCAGGCTGATGATTTTTTCGCGTACGGCGAACAGCACCAGGCCGGCTACGTCGTGAATCTGCAGGCGCTTCATGATTTGCGAGCGGTGa         <  1:697314/150‑1 (MQ=255)
                                                                                                                                       tCAGGCAAGGCCGGCGGCGTTGCACCTGCACCTGCCGGCCGGGTCACCGTTGCTCAGTCATCCAGGCTGATGATTTTTTCGCGTACGGCGAACAGCACCAGGCCGGCTACGTCGTGAATCTGCAGGCGCTTCATGATTTGCGAGCGGTGa         <  2:557688/150‑1 (MQ=255)
                                                                                                                                        cAGGCAAGGCCGGCGGCGTTGCACCTGCACCTGCCGGCCGGGTCACCGTTGCTCAGTCATCCAGGCTGATGATTTTTTCGCGTACGGCGAACAGCACCAGGCCGGCTACGTCGTGAATCTGCAGGCGCTTCATGATTTGCGAGCGGTGAg        <  2:114661/150‑1 (MQ=255)
                                                                                                                                         aGGCAAGGCCGGCGGCGTTGCACCTGCACCTGCCGGCCGGGTCACCGTTGCTCAGTCATCCAGGCTGATGATTTTTTCGCGTACGGCGAACAGCACCAGGCCGGCTACGTCGTGAATCTGCAGGCGCTTCATGATTTGCGAGCGGTGAGt       >  2:1054695/1‑150 (MQ=255)
                                                                                                                                           gCAACGCCGCCGGCCTTGCACCTGCACCTGCCGGCCGGGTCACCGTTGCTCAGTCATCCAGGCTGATGATTTTTTCGCGTACGGCGAACAGCACCAGGCCGGCTACGTCGTGa                                          <  1:845346/113‑1 (MQ=255)
                                                                                                                                           gCAACGCCGCCGGCCTTGCACCTGCACCTGCCGGCCGGGTCACCGTTGCTCAGTCATCCAGGCTGATGATTTTTTCGCGTACGGCGAACAGCACCAGGCCGGCTACGTCGTGa                                          >  2:845346/1‑113 (MQ=255)
                                                                                                                                              aGGCCGGCGGCGTTGCACCTGCACCTGCCGGCCGGGTCACCGTTGCTCAGTCATCCAGGCTGATGATTTTTTCGCGTACGGCGAACAGCACCAGGCCGGCTACGTCGTGAATCTGCAGGCGCTTCATGATTTGCGAGCGGTGAGTTTCCa  <  1:1054695/150‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                    
GGGGGCTCTGCGGGACCTGTTGGGCGGTGACGGCCTGTTCGTTCGGGCCAGGTTGCCATGGGTCGAGATCAAGAGTTCGCAAGGGTTGGACATCAACCTCGATGGCGAGCCGTTGCAGGCTGCCAACCTGCGCTTTCAGGCAAGGCCGGCGGCGTTGCACCTGCACCTGCCGGCCGGGTCACCGTTGCTCAGTCATCCAGGCTGATGATTTTTTCGCGTACGGCGAACAGCACCAGGCCGGCTACGTCGTGAATCTGCAGGCGCTTCATGATTTGCGAGCGGTGAGTTTCCA  >  NC_002947/2425158‑2425449

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