Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 892,892 C→T 9.5% A276T (GCA→ACA)  pta ← phosphate acetyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947892,8920CT 9.5% 135.8 / 4.1 52A276T (GCA→ACA) ptaphosphate acetyltransferase
Reads supporting (aligned to +/- strand):  ref base C (25/22);  new base T (2/3);  total (27/25)
Fisher's exact test for biased strand distribution p-value = 6.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

GTTGGCGGTGTCGTACGAGCCGGTGCTGACCGAAAGGATCGGCAGGCCGGCCTGCAGGGCACCGCGGCACAGGCCAAGGATACGCACGTCAGGCTTGCTGTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCGGGACCGTGTTGGCCACCGTGCGTGCACAGATGATGATCTTGCTCATGCGCCGCTGCTCGTAATCGCCGGCATT  >  NC_002947/892746‑893036
                                                                                                                                                  |                                                                                                                                                
gTTGGCGGTGTCGTACGAGCCGGTGCTGACCGAAAGGATCGGCAGGCCGGCCTGCAGGGCACCGCGGCACAGGCCAAGGATACGCACGTCAGGCTTGCTGTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAg                                                                                                                                               >  1:1119460/1‑150 (MQ=255)
 ttGGCGGTGTCGTACGAGCCGGTGCTGACCGAAAGGATCGGCAGGCCGGCCTGCAGGGCACCGCGGCACAGGCCAAGGATACGCACGTCAGGCTTGCTGTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAgg                                                                                                                                              <  2:1098170/150‑1 (MQ=255)
   ggCGGTGTCGTACGAGCCGGTGCTGACCGAAAGGATCGGCAGGCCGGCCTGCAGGGCACCGCGGCACAGGCCAAGGATACGCACGTCAGGCTTGCTGTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCt                                                                                                                                            <  1:587110/150‑1 (MQ=255)
   ggCGGTGTCGTACGAGCCGGTGCTGACCGAAAGGATCGGCAGGCCGGCCTGCAGGGCACCGCGGCACAGGCCAAGGATACGCACGTCAGGCTTGCTGTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCt                                                                                                                                            <  1:488811/150‑1 (MQ=255)
     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)
                                                  ccTGCAGGGCACCGCGGCACAGGCCAAGGATACGCACGTCAGGCTTGCTGTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCa                                                                                                                                                <  2:391325/99‑1 (MQ=255)
                                                  ccTGCAGGGCACCGCGGCACAGGCCAAGGATACGCACGTCAGGCTTGCTGTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCa                                                                                                                                                >  1:391325/1‑99 (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                                                                                               >  2:1148351/1‑87 (MQ=255)
                                                                                                               cagcaAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAg                                                                                               <  1:1148351/87‑1 (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                                                                 >  1:794453/1‑88 (MQ=255)
                                                                                                                                            tCGCTGTCAGCCTGGCAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCAGGACCGTGtt                                                                 <  2:794453/88‑1 (MQ=255)
                                                                                                                                            tCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCGGGACCGTGTTGGCCACCGTGCGTGCACAGATGATGATCTTGCTCATGCGCCGCTGCTCGTAATCGCCGGCAt   >  2:843643/1‑150 (MQ=255)
                                                                                                                                             cGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCGGGACCGTGTTGGCCACCGTGCGTGCACAGATGATGATCTTGCTCATGCGCCGCTGCTCGTAATCGCCGGCAtt  <  1:832764/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
GTTGGCGGTGTCGTACGAGCCGGTGCTGACCGAAAGGATCGGCAGGCCGGCCTGCAGGGCACCGCGGCACAGGCCAAGGATACGCACGTCAGGCTTGCTGTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCGGGACCGTGTTGGCCACCGTGCGTGCACAGATGATGATCTTGCTCATGCGCCGCTGCTCGTAATCGCCGGCATT  >  NC_002947/892746‑893036

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: ‑