Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 892,900 A→G 10.4% V273A (GTC→GCC)  pta ← phosphate acetyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947892,9000AG10.4% 125.1 / 4.7 48V273A (GTC→GCC) ptaphosphate acetyltransferase
Reads supporting (aligned to +/- strand):  ref base A (25/18);  new base G (2/3);  total (27/21)
Fisher's exact test for biased strand distribution p-value = 6.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.57e-01

CGGTGTCGTACGAGCCGGTGCTGACCGAAAGGATCGGCAGGCCGGCCTGCAGGGCACCGCGGCACAGGCCAAGGATACGCACGTCAGGCTTGCTGTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCGGGACCGTGTTGGCCACCGTGCGTGCACAGATGATGATCTTGCTCATGCGCCGCTGCTCGTAATCGCCGGCATTGAGCACCTG  >  NC_002947/892751‑893045
                                                                                                                                                     |                                                                                                                                                 
cGGTGTCGTACGAGCCGGTGCTGACCGAAAGGATCGGCAGGCCGGCCTGCAGGGCACCGCGGCACAGGCCAAGGAGACGCACGTCAGGCTTGCTGTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGa                                                                                                                                                   <  1:1385141/150‑1 (MQ=255)
 ggTGTCGTACGAGCCGGTGCTGACCGAAAGGATCGGCAGGCCGGCCTGCAGGGCACCGCGGCACAGGCCAAGGATACGCACGTCAGGCTTGCTGTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGAc                                                                                                                                                  >  2:1123898/1‑150 (MQ=255)
      cGTACGAGCCGGTGCTGACCGAAAGGATCGGCAGGCCGGCCTGCAGGGCACCGCGGCACAGGCCAAGGATACGCACGTCAGGCTTGCTGTCGCTGGTCCGCATCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGa                                                                                                                                             <  1:1105480/150‑1 (MQ=255)
      cGTACGAGCCGGTGCTGACCGAAAGGATCGGCAGGCCGGCCTGCAGGGCACCGCGGCACAGGCCAAGGATACGCACGTCAGGCTTGCTGTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGTCAGCCTGGCAGCGa                                                                                                                                             <  1:773832/150‑1 (MQ=255)
      cGTACGAGCCGGTGCTGACCGAAAGGATCGGCAGGCCGGCCTGCAGGGCACCGCGGCACAGGCCAAGGATACGCACGTCAGGCTTGCTGTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGa                                                                                                                                             <  2:669123/150‑1 (MQ=255)
          cGAGCCGGTGCTGACCGAAAGGATCGGCAGGCCGGCCTGCAGGGCACCGCGGCACAGGCCAAGGATACGCACGTCAGGCTTGCTGTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGgat                                                                                                                                         <  1:1123898/150‑1 (MQ=255)
          agcgCCGGTGCTGACCGAAAGGATCGGCAGGGCGGCCTGCAGGGCACCGCGGCAGAGGCCAAGGATACGCCCGTCAGGCTTGCTGTTGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCGAGGCTGACAGCGAGgat                                                                                                                                         <  2:469309/147‑1 (MQ=255)
           gAGCCGGTGCTGACCGAAAGGATCGGCAGGCCGGCCTGCAGGGCACCGCGGCACAGGCCAAGGATACGCACGTCAGGCTTGCTGTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGgatg                                                                                                                                        <  1:443962/150‑1 (MQ=255)
             gCCGGTGCTGACCGAAAGGATCGGCAGGCCGGCCTGCAGGGCACCGCGGCACAGGCCAAGGATACGCACGTCAGGCTTGCTGTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGgatgat                                                                                                                                      <  2:549953/150‑1 (MQ=255)
             gCCGGTGCTGACCGAAAGGATCGGCAGGCCGGCCTGCAGGGCACCGCGGCACAGGCCAAGGATACGCACGTCAGGCTTGCTGTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCAGTGATCGCTGCCAGGCTGACAGCGAGgatgat                                                                                                                                      <  1:1402999/150‑1 (MQ=255)
                         cGAAAGGATCGGCAGGCCGGCCTGCAGGGCACCGCGGCACAGGCCAAGGATACGCACGTCAGGCTTGCTGTCGCTGGTCTTCGATTGGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTc                                                                                                                          >  1:1156510/1‑150 (MQ=255)
                           aaaGGATCGGCAGGCCGGCCTGCAGGGCACCGCGGCACAGGCCAAGGATACGCACGTCAGGCTTGCTGTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGc                                                                                                                        >  2:491172/1‑150 (MQ=255)
                                aTCGGCAGGCCGGCCTGCAGGGCACCGCGGCACAGGCCAAGGATACGCACGTCAGGCTTGCTGTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCggg                                                                                                                   <  1:1145334/150‑1 (MQ=255)
                                 cctgCAGGCCGGCCTGCAGGGCACCGCGGCACAGGCCAAGGATACGCACGTCAGGCTTGCTGTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCgggg                                                                                                                  >  2:635280/4‑150 (MQ=255)
                                        gccggccTGCAGGGCACCGCGGCACAGGCCAAGGATACGCACGTCAGGCTTGCTGTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTaccac                                                                                                           >  1:452556/1‑150 (MQ=255)
                                          cggccTGCAGGGCACCGCGGCACAGGCCAAGCATACGCACGTCAGGCTTGCTGTCGCTGGTCAGCAGCAAGCCGGACAGGGGGACGCCATTGATCGCTGCCAGGCTCACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTAccacca                                                                                                         >  2:646355/1‑150 (MQ=255)
                                               tGCAGGGCACCGCGGCACAGGCCAAGGATACGCACGTCAGGCTTGCTGTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTg                                                                                                    >  2:1104316/1‑150 (MQ=255)
                                                    ggCACCGCGGCACAGGCCAAGGATACGCACGTCAGGCTTGCTGTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGa                                                                                               <  2:720556/150‑1 (MQ=255)
                                                         cgcgGCACAGGCCAAGGATACGCACGTCAGGCTTGCTGTCGCTGGAAAGCGTGGTACGGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCa                                                                                          <  2:1107251/150‑1 (MQ=255)
                                                                      aaGGATACGCACGTCAGGCTTGCTGTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCGGGACCgt                                                                             >  2:201919/1‑150 (MQ=255)
                                                                              gCACGTCAGGCTTGCTGTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCGGGACCGTGTTGGCCa                                                                     >  1:1149306/1‑150 (MQ=255)
                                                                                  gTCAGGCTTGCTGTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCGGGACCGTGTTGGCCACcgt                                                                 >  2:1165095/1‑150 (MQ=255)
                                                                                   tCAGGCTTGCTGTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCGGGACCGTGTTGGCCACcgtg                                                                >  1:214829/1‑150 (MQ=255)
                                                                                    cAGGCTTGCTGTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCGGGACCGTGTTGGCCACcgtgc                                                               <  2:1119460/150‑1 (MQ=255)
                                                                                      ctctggCTGTCGCGGGTCAGCAG‑‑‑GCGGGGCAGTGGGACGCCATTGATCGCTGCGAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCGGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCGGGACCGTGTTGGCCACCGTGCGTGca                                                          <  2:82747/145‑1 (MQ=255)
                                                                                         ttGCTGTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCGGGACCGTGTTGGCCACCGTGCGTGca                                                          >  1:682831/1‑150 (MQ=255)
                                                                                         ttGCTGTCGCTGGT‑‑‑CAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCGGGACCGTGTTGGCCACCGTGCGTGCACAg                                                       >  1:82747/1‑150 (MQ=255)
                                                                                              gTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCGGGACCGTGTTGGCCACCGTGCGTGCACAgat                                                     >  1:264031/1‑150 (MQ=255)
                                                                                              gTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCGGGACCGTGTTGGCCACCGTGCGTGCACAgat                                                     >  2:40225/1‑150 (MQ=255)
                                                                                                cGCCGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCGGGACCGTGTTGGCCACCGTGCGTGCACAgatga                                                   >  1:49607/1‑150 (MQ=255)
                                                                                                       cagcagCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCGGGACCGTGTTGGCCACCGTGCGTGCACAGATGATGATCtt                                            <  1:478849/150‑1 (MQ=255)
                                                                                                        agcagcAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCGGGACCGTGTTGGCCACCGTGCGTGCACAGATGATGATCTTg                                           >  2:832764/1‑150 (MQ=255)
                                                                                                         gcagcaAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCGGGACCGTGTTGGCCACCGTGCGTGCACAGATGATGATCTTGc                                          >  1:527975/1‑150 (MQ=255)
                                                                                                          cagcaAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAg                                                                                                        <  1:1148351/87‑1 (MQ=255)
                                                                                                          cagcaAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAg                                                                                                        >  2:1148351/1‑87 (MQ=255)
                                                                                                                      aaGGGGGGCGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGGCGTCGCGGTCGCCCGAGGTTACCACCAGGGAGCCCGAGGCCAGTAGCGGGCCCGTATTGGCCACCGTGCGTGCACAGATGATGATCTTGCTCATGCGCCgctg                             >  2:89487/2‑150 (MQ=255)
                                                                                                                               gCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCGGGACCGTGTTGGCCACCGTGCGTGCACAGATGATGATCTTGCTCATGCGCCGCTGCTCGTAATc                    >  2:476873/1‑150 (MQ=255)
                                                                                                                                ccATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCGGGACCGTGTTGGCCACCGTGCGTGCACAGATGATGATCTTGCTCATGCGCCGCTGCTCGTAATCg                   >  2:681061/1‑150 (MQ=255)
                                                                                                                                  aTTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTAACACCAGGGTGCCCGAGGTCAGTAGCGGGACCGTGTTGGCCACCGTGCGTGCACAGATGATGATCTTGCTCATGCGCCGCTGCTCGTAATCGcc                 <  1:491172/150‑1 (MQ=255)
                                                                                                                                      aTCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCGGGACCGTGTTGGCCACCGTGCGTGCACAGATGATGATCTTGCTCATGCGCCGCTGCTCGTAATCGCCGGCa             <  2:527975/150‑1 (MQ=255)
                                                                                                                                       tCGCTGTCAGCCTGGCAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCGGGACCGTGTTGGCCACCGTGCGTGCACAgatgatgat                                               >  1:33396/1‑115 (MQ=255)
                                                                                                                                       tCGCTGTCAGCCTGGCAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCGGGACCGTGTTGGCCACCGTGCGTGCACAgatgatgat                                               <  2:33396/115‑1 (MQ=255)
                                                                                                                                       tCGCTGTCAGCCTGGCAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCAGGACCGTGtt                                                                          <  2:794453/88‑1 (MQ=255)
                                                                                                                                       tCGCTGTCAGCCTGGCAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCAGGACCGTGtt                                                                          >  1:794453/1‑88 (MQ=255)
                                                                                                                                       tCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCGGGACCGTGTTGGCCACCGTGCGTGCACAGATGATGATCTTGCTCATGCGCCGCTGCTCGTAATCGCCGGCAt            >  2:843643/1‑150 (MQ=255)
                                                                                                                                        cGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCGGGACCGTGTTGGCCACCGTGCGTGCACAGATGATGATCTTGCTCATGCGCCGCTGCTCGTAATCGCCGGCAtt           <  1:832764/150‑1 (MQ=255)
                                                                                                                                                ggCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCGGGACCGTGTTGGCCACCGTGCGTGCACAGATGATGATCTTGCTCATGCGCCGCTGCTCGTAATCGCCGGCAtt           <  1:1213666/142‑1 (MQ=255)
                                                                                                                                                ggCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCGGGACCGTGTTGGCCACCGTGCGTGCACAGATGATGATCTTGCTCATGCGCCGCTGCTCGTAATCGCCGGCAtt           >  2:1213666/1‑142 (MQ=255)
                                                                                                                                                 gCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCGGTAGCGGGACCGTGTTGGCCACCGTGCGTGCACAGATGATGATCTTGCTCATGCGCCGCTGCTCGTAATCGCCGGCATTGAGCACCTg  >  2:1262634/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
CGGTGTCGTACGAGCCGGTGCTGACCGAAAGGATCGGCAGGCCGGCCTGCAGGGCACCGCGGCACAGGCCAAGGATACGCACGTCAGGCTTGCTGTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCGGGACCGTGTTGGCCACCGTGCGTGCACAGATGATGATCTTGCTCATGCGCCGCTGCTCGTAATCGCCGGCATTGAGCACCTG  >  NC_002947/892751‑893045

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