Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 908,135 C→G 8.1% A27P (GCG→CCG)  yjjV ← putative hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947908,1350CG 8.1% 152.2 / 3.6 49A27P (GCG→CCG) yjjVputative hydrolase
Reads supporting (aligned to +/- strand):  ref base C (22/23);  new base G (2/2);  total (24/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

CTCAGGGCGATGCTGGTCGAGGTAGACCGGGTGCAGGCCAAGCGCGGCGAACAGGCGCGGCTCGGCGCAGGCCAGTTCCCACACCCGCTGGAAATTCGCCTGGTACACCCCCAGCACCACCATTCGCTCCAGCCCACGCGCCGCCGCGTTGGCCAGCAGGCGCGGGCGGTCGGCGTCGAAGTCGGGGAAGTCCAGGTGGGTGTGGGTGTCGATCAGGCGCATGTTCAGGCCACTGAGATGCGTTGCTTGAAGGTGCGGCCAACGGCGTGCACGCCGGGTTCATAGCGTTTTTC  >  NC_002947/907992‑908284
                                                                                                                                               |                                                                                                                                                     
cTCAGGGCGATGCTGGTCGAGGTAGACCGGGTGCAGGCCAAGCGCGGCGAACAGGCGCGGCTCGGCGCAGGCCAGTTCCCACACCCGCTGGCAATTCGCCTGGTACACCCCCAGCACCACCATTCGCTCCAGCCCACGCGCCGCCGCGtt                                                                                                                                                 <  2:630197/150‑1 (MQ=255)
 tCAGGGCGATGCTGGTCGAGGTAGACCGGGTGCAGGCCAAGCGCGGCGAACAGGCGCGGCTCGGCGCAGGCCAGTTCCCACACCCGCTGGAAATTCGCCTGGTACACCCCCAGCACCACCATTCGCTCCAGCCCATGCGCCGCCGCGTTg                                                                                                                                                >  2:169786/1‑150 (MQ=255)
 tCAGGGCGATGCTGGTCGAGGTAGACCGGGTGCAGGCCAAGCGCGGCGAACAGGCGCGGCTCGGCGCAGGCCAGTTCCCACACCCGCTGGAAATTCGCCTGGTACACCCCCAGCACCACCATTCGCTCCAGCCCACGCGCCGCCGCGTTg                                                                                                                                                >  1:198367/1‑150 (MQ=255)
   aGGGCGATGCTGGTCGAGGTAGACCGGGTGCAGGCCAAGCGCGGCGAACAGGCGCGGCTCGGCGCAGGCCAGTTCCCACACCCGCTGGAAATTCGCCTGGTACACCCCCAGCACCACCATTCGCTCCAGCCCACGCGCCGCCGCGTTGGc                                                                                                                                              >  1:462340/1‑150 (MQ=255)
            cTGGTCGAGGTAGACCGGGTGCGGGCCAAGCGCGGCGAACAGGCGCGGCTCGGCGCAGGCCAGTTCCCACACCCGCTGGAAATTCGCCTGGTACCCCCCCAGCACCACCATTCGCTCCAGCCCACGCGCCGCCGCGTTGGCCAGCAGgcg                                                                                                                                     <  1:410023/150‑1 (MQ=255)
                    ggTAGACCGGGTGCAGGCCAAGCGCGGCGAACAGGCGCGGCTCGGCGCAGGCCAGTTCCCACACCCGCTGGAAATTCGCCTGGTACACCCCCAGCACCACCATTCGCTCCAGCCCACGCGCCGCCGCGTTGGCCAGCAGGCGCGGGCGGt                                                                                                                             >  1:470709/1‑150 (MQ=255)
                      tAGACCGGGTGCAGGCCAAGCGCGGCGAACAGGCGCGGCTCGGCGCAGGCCAGTTCCCACACCCGCTGGAAATTCGCCTGGTACACCCCCAGCACCACCATTCGCTCCAGCCCACGCGCCGCCGCGTTGGCCAGCAGGCGCGGGCGGTCg                                                                                                                           <  2:466982/150‑1 (MQ=255)
                      tAGACCGGGTGCAGGCCAAGCGCGGCGAACAGGCGCGGCTCGGCGCAGGCCAGTTCCCACACCCGCTGGAAATTCGCCTGGTACACCCCCAGCACCACCATTCGCTCCAGCCAACGCGGCGGCGCGTGGGctg                                                                                                                                            <  1:451288/133‑3 (MQ=255)
                      tAGACCGGGTGCAGGCCAAGCGCGGCGAACAGGCGCGGCTCGGCGCAGGCCAGTTCCCACACCCGCTGGAAATTCGCCTGGTACACCCCCAGCACCACCATTCGCTCCAGCCAACGCGGCGGCGCGTGGGctg                                                                                                                                            >  2:451288/1‑131 (MQ=255)
                           cGGGTGCAGGCCAAGCGCGGCGAACAGGCGCGGCTCGGCGCAGGCCAGTTCCCACACCCGCTGGAAATTCGCCTGGTACACCCCCAGCACCACCATTCGCTCCAGCCCACGCGCCGCCGCGTTGGCCAGCAGGCGCGGGCGGTCGGCGTc                                                                                                                      >  1:484135/1‑150 (MQ=255)
                                gCAGGCCAAGCGCGGCGAACAGGCGCGGCTCGGCGCAGGCCAGTTCCCACACCCGCTGGAAATTCGCCTGGTACACCCCCAGCACCACCATTCGCTCCAGCCCACGCGCCGCCGCGTTGGCCAGCAGGCGCGGGCGGTCGGCGTCGAAGt                                                                                                                 >  1:696/1‑150 (MQ=255)
                                  aGGCCAAGCGCGGCGAACAGGCGCGGCTCGGCGCAGGCCAGTTCCCACACCCGCTGGAAATTCGCCTGGTACACCCCCAGCACCACCATTCGCTCCAGCCCACGCGCCGCCGCGTTGGCCAGCAGGCGCGGGCGGTCGGCGTCGAAGTCg                                                                                                               >  1:96744/1‑150 (MQ=255)
                                   ggCCAAGCGCGGCGAACAGGCGCGGCTCGGCGCAGGCCAGTTCCCACACCCGCTGGAAATTCGCCTGGTACACCCCCAGCACCACCATTCGCTCCAGCCCATGCGCCGCCGCGTTGGCCAGCAGGCGCGGGCGGTCGGCGTCGAAGTCgg                                                                                                              <  1:169786/150‑1 (MQ=255)
                                    gCCAAGCGCGGCGAACAGGCGCGGCTCGGCGCAGGCCAGTTCCCACACCCGCTGGAAATTCGCCTGGTACACCCCCAGCACCACCATTCGCTCCAGCCCACGCGCCGCCGCGTTGGCCAGCAGGCGCGGGCGGTCGGCGTCGAAGTCggg                                                                                                             <  1:599556/150‑1 (MQ=255)
                                    gCCAAGCGCGGCGAACAGGCGCGGCTCGGCGCAGGCCAGTTCCCACACCCGCTGGAAATTCGCCTGGTACACCCCCAGCACCACCATTCGCTCCAGCCCACGCGCCGCCGCGTTGGCCAGCAGGCGCGGGCGGTCGGCGTCGAAGTCggg                                                                                                             <  2:96744/150‑1 (MQ=255)
                                     ccAAGCGCGGCGAACAGGTCCGGCGCGGCGCAGGCCAGTTCCCACACCCGCTGGCAATTCGCCTGCTACATCCCCAGCACCACCATTCGCTCCATCCCACGCGCCGCCGCGTTGGCCAGCAGGCGCGGGCGGTCGGCGTCGAAGTCgggg                                                                                                            <  1:228529/150‑1 (MQ=255)
                                                 aaCAGGCGCGGCTCGGCGCAGGCCAGTTCCCACACCCGCTGGAAATTCGCCTGGTACACCCCCAGCACCACCATTCGCTCCAGCCCACGCGCCGCCGCGTTGGCCAGCAGGCGCGGGCGGTCGGCGTCGAAGTCGGGGAAGTCCAggtgg                                                                                                <  1:250880/150‑1 (MQ=255)
                                                   cAGGCGCGGCTCGGCGCAGGCCAGTTCCCACACCCGCTGGAAATTCGCCTGGTACACCCCCAGCACCACCATTCGCTCCAGCCCACGCGCCGCCGCGTTGGCCAGCAGGCGCGGGCGGTCGGCGTCGAAGTCGGGGAAGTCCAggtgggt                                                                                              <  1:622673/150‑1 (MQ=255)
                                                    aGGCGCGGCTCGGCGCAGGCCAGTTCCCACACCCGCTGGAAATTCAACTGCTAGAGCCCCAGCACCACCATTCGCTCCAGCCCACGCGCCGCCGCGTTGGCCAGCAGGCGCGGGCGGTCGGTGTCGAAGTCGGGGAAGTCCAggtgggtg                                                                                             <  2:395571/150‑1 (MQ=255)
                                                      gcgcGGCTCGGCGCAGGCCAGTTCCCACACCCGCTGGAAATTCGCCTGGTACACCCCCAGCACCACCATTCGCTCCAGCCCACGCGCCGCCGCGTTGGCCAGCAGGCGCGGGCGGTCGGCGTCGAAGTCGGGGAAGTCCAGgtgggtgtg                                                                                           >  1:146173/1‑150 (MQ=255)
                                                       cgcgGCTCGGCGCAGGCCAGTTCCCACACCCGCTGGAAATTCGCCTGGTACACCCCCAGCACCACCATTCGCTCCAGCCCACGCGCCGCCGCGTTGGCCAGCAGGCGCGGGCGGTCGGCGTCGAAGTCGGGGAAGTCCAGgtgggtgtgg                                                                                          <  2:462340/150‑1 (MQ=255)
                                                       cgcgGCTCGGCGCAGGCCAGTTCCCACACCCGCTGGAAATTCGCCTGGTACACCCCCAGCACCACCATTCGCTCCAGCCCACGCGCCGCCGCGTTGGCCAGCAGGCGCGGGCGGTCGGCGTCGAAGTCGGGGAAGTCCAGgtgggtgtgg                                                                                          <  2:612488/150‑1 (MQ=255)
                                                       cgcgGCTCGGCGCAGGCCAGTTCCCACACCCGCTGGAAATTCGCCTGGTACACCCCCAGCACCACCATTCGCTCCAGCCCACGCGCCGCCGCGTTGGCCAGCAGGCGCGGGCGGTCGGCGTCGAAGTCGGGGAAGTCCAGgtgggtgtgg                                                                                          <  1:251729/150‑1 (MQ=255)
                                                                 cgcAGGCCAGTTCCCACACCCGCTGGAAATTCGCCTGGTACACCCCCAGCACCACCATTCGCTCCAGCCCACGCGCCGCCGCGTTGGCCAGCAGGCGCGGGCGGTCGGCGTCGAAGTCGGGGAAGTCCAGGTGGGTGTGGGTGTCGATCa                                                                                >  1:146890/1‑150 (MQ=255)
                                                                        cAGTTCCCACACCCGCTGGAAATTCGCCTGGTACACCCCCAGCACCACCATTCGCTCCAGCCCACGCGCCGCCGCGTTGGCCAGCAGGCGCGGGCGGTCGGCGTCGAAGTCGGGGAAGTCCAGGTGGGTGTGGGTGTCGATCAgg                                                                              <  1:265040/145‑1 (MQ=255)
                                                                        cAGTTCCCACACCCGCTGGAAATTCGCCTGGTACACCCCCAGCACCACCATTCGCTCCAGCCCACGCGCCGCCGCGTTGGCCAGCAGGCGCGGGCGGTCGGCGTCGAAGTCGGGGAAGTCCAGGTGGGTGTGGGTGTCGATCAgg                                                                              >  2:265040/1‑145 (MQ=255)
                                                                        cAGTTCCCACACCCGCTGGAAATTCGCCTGGTACACCCCCAGCACCACCATTCGCTCCAGCCCACGCGCCGCCGCGTTGGCCAGCAGGCGCGGGCGGTCGGCGTCGAAGTCGGGGAAGTCCAGGTGGGTGTGGGTGTCGATCAGGCGCAt                                                                         >  1:34987/1‑150 (MQ=255)
                                                                         aGTTCCCACACCCGCTGGAAATTCGCCTGGTACACCCCCAGCACCACCATTCGCTCCAGCCCACGCGGCGGCGCGtgg                                                                                                                                                >  2:458241/1‑76 (MQ=255)
                                                                         aGTTCCCACACCCGCTGGAAATTCGCCTGGTACACCCCCAGCACCACCATTCGCTCCAGCCCACGCGGCGGCGCGtgg                                                                                                                                                <  1:458241/78‑3 (MQ=255)
                                                                           ttCCCACACCCGCTGGAAATTCGCCTGGTACACCCCCAGCACCACCATTCGCTCCAGCCCACGCGCCGCCGCGTTGGCCAGCAGGCGCGGGCGGTCGGCGTCGAAGTCGGGGAAGTCCAGGTGGGTGTGGGTGTCGATCAg                                                                               <  1:595925/141‑1 (MQ=255)
                                                                           ttCCCACACCCGCTGGAAATTCGCCTGGTACACCCCCAGCACCACCATTCGCTCCAGCCCACGCGCCGCCGCGTTGGCCAGCAGGCGCGGGCGGTCGGCGTCGAAGTCGGGGAAGTCCAGGTGGGTGTGGGTGTCGATCAg                                                                               >  2:595925/1‑141 (MQ=255)
                                                                               cacaCCCGCTGGAAATTCGCCTGGTACACCCCCAGCACCACCATTCGCTCCAGCCCACGCGCCGCCGCGTTGGCCAGCAGGCGCGGGCGGTCGGCGTCGAAGTCGGGGAAGTCCAGGTGGGTGTGGGTGTCGATCAGGCGCATGTTCAgg                                                                  <  1:6506/150‑1 (MQ=255)
                                                                                    ccGCTGGAAATTCGCCTGGTACACCCCCAGCACCACCATTCGCTCCAGCCCACGCGCCGCCGCGTTGGCCAGCAGGCGCGGGCGGTCGGCGTCGAAGTCGGGGAAGTCCAGGTGGGTGTGGGTGTCGATCAGGCGCATGTTCAGGCCACt                                                             >  1:518985/1‑150 (MQ=255)
                                                                                               tCGCCTGGTACACCCCCAGCACCACCATTCGCTCCAGCCCACGCGCCGCCGCGTTGGCCAGCAGGCGCGGGCGGTCGGCGTCGAAGTCGGGGAAGTCCAGGTGGGTGTGGGTGTCGATCAGGCGAATGTTCAGGCCACTGAGATGCGTTg                                                  <  2:198367/150‑1 (MQ=255)
                                                                                                cGCCTGGTACACCCCCAGCACCACCATTCGCTCCAGCCCACGCGCCGCCGCGTTGGCCAGCAGGCGCGGGCGGTCGGCGTCGAAGTCGGGGAAGTCCAGGTGGGTGTGGGTGTCGATCAGGCGCATGTTCAGGCCACTGAGATGCGTTGc                                                 <  1:420760/150‑1 (MQ=255)
                                                                                                  ccTGGTACACCCCCAGCACCACCATTCGCTCCAGCCCACGCGCCGCCGCGTTGGCCAGCAGGCGCGGGCGGTCGGCGTCGAAGTCGGGGACGTCCAGGTGGGTGTGGGTGTCGATCAGGCGCATGTTCAGGCCACTGAGATGCGTTGCtt                                               >  1:384815/1‑150 (MQ=255)
                                                                                                   cTGGTACACCCCCAGCACCACCATTCGCTCCAGCCCACGCGCCGCCGCGTTGGCCAGCAGGCGCGGGCGGTCGGCGTCGAAGTCGGGGAAGTCCAGGTGGGTGTGGGTGTCGATCAGGCGCATGTTCAGGCCACTGAGATGCGTTGCTTg                                              >  2:564647/1‑150 (MQ=255)
                                                                                                      gTACACCCCCAGCACCACCATTCGCTCCAGCCCACGCGCCGCCGCGTTGGCCAGCAGGCGCGGGCGGTCGGCGTCGAAGTCGGGGAAGTCCAGGTGGGTGTGGGTGTCGATCAGGCGCATGTTCAGGCCACTGAGATGCGTTGCTTGAAg                                           <  2:436646/150‑1 (MQ=255)
                                                                                                           cccccAGCACCACCATTCGCTCCAGCCCACGCGCCGCCGCGTTGGCCAGCAGGCGCGGGCGGTCGGCGTCGAAGTCg                                                                                                               >  2:224059/1‑77 (MQ=255)
                                                                                                           cccccAGCACCACCATTCGCTCCAGCCCACGCGCCGCCGCGTTGGCCAGCAGGCGCGGGCGGTCGGCGTCGAAGTCg                                                                                                               <  1:224059/77‑1 (MQ=255)
                                                                                                             cccAGCACCACCATTCGCTCCAGCCCACGCGCCGCCGCGTTGGCCAGCAGGCGCGGGCGGTCGGCGTCGAAGTCGGGGAAGTCCAGGTGGGTGTGTGTGTCGATCAGGCGCATGTTCAGGCCACTGAGATGCGTTGCTTGAAGGTGCGGc                                    >  2:356873/1‑150 (MQ=255)
                                                                                                             cccAGCACCACCATTCGCTCCAGCCCACGCGCCGCCGCGTTGGCCAGCAGGCGCGGGCGGTCGGCGTCGAAGTCGGGGAAGTCCAGGTGGGTGTGGGTGTCGATCAGGCGCATGTTCAGGCCACTGAGATGCGTTACTTGAAGGTGCGGc                                    >  1:569218/1‑150 (MQ=255)
                                                                                                              ccAGCACCACCATTCGCTCCAGCCCACGCGCCGCCGCGTTGGCCAGCAGGCGCGGGCGGTCGGCGTCGAAGTCGGGGAAGTCCAGGTGGGTGTGGGTGTCGATCAGGCGCATGTTCAGGCCACTGAGATGCGTTGCTTGAAGGTGCGGcc                                   >  2:416456/1‑150 (MQ=255)
                                                                                                               cAGCACCACCATTCGCTCCAGCCCACGCGCCGCCGCGTTGGCCAGCAGGCGCGGGCGGTCGGCGTCGAAGTCGGGGAAGTCCAggt                                                                                                  >  2:573826/1‑86 (MQ=255)
                                                                                                               cAGCACCACCATTCGCTCCAGCCCACGCGCCGCCGCGTTGGCCAGCAGGCGCGGGCGGTCGGCGTCGAAGTCGGGGAAGTCCAggt                                                                                                  <  1:573826/86‑1 (MQ=255)
                                                                                                                      accaTTCGCTCCAGCCCACGCGCCGCCGCGTTGGCCAGCAGGCGCGGGCGGTCGGCGTCGAAGTCGGGGAAGTCCAGGTGGGTGTGGGTGTCGATCAGGCGCATGTTCAGGCCACTGAGATGCGTTGCTTGAAGGTGCGGCCAACGGCGt                           <  2:146173/150‑1 (MQ=255)
                                                                                                                               tCCAGCCCACGCGCCGCCGCGTTGGCCAGCAGGCGCGGGCGGTCGGCGTCGAAGTCGGGGAAGTCCAGGTGGGTGTGGGTGTCGATCAGGCGCATGTTCAGGCCACTGAGATGCGTTGCTTGAAGGTGAGGCCAACGGCGTGAACGTCgg                  >  1:299520/1‑150 (MQ=255)
                                                                                                                                      cACGCGCCGCCGCGTTGGCCAGCAGGCGCGGGCGGTCGGCGTCGAAGTCGGGGAAGTCCAggtgg                                                                                                >  1:714350/1‑65 (MQ=255)
                                                                                                                                      cACGCGCCGCCGCGTTGGCCAGCAGGCGCGGGCGGTCGGCGTCGAAGTCGGGGAAGTCCAggtgg                                                                                                <  2:714350/65‑1 (MQ=255)
                                                                                                                                       gggggcggccgcGTTGGCCAGCAGGCGCGGGCGGTCGGCGTCGAAGTCGGGGAAGTCCAGGTGGGTGTGGGTGTCGATCAGGCGCATGTTCAGGCCACTGAGATGCGTTGCTTGAAGGTGCGGCCAACGGCGTGCACGCCGGGTTCATAg          >  2:200366/8‑150 (MQ=255)
                                                                                                                                               ccgcGTTGGCCAGCAGGCGCGGGCGGTCGGCGTCGAAGTCGGGGAAGTCCAGGTGGGTGTGGGTGTCGATCAGGCGCATGTTCAGGCCACTGAGATGCGTTGCTTGAAGGTGCGGCCAACGGCGTGCACGCCGGGTTCATAGCGTTTTTc  >  2:48040/1‑150 (MQ=255)
                                                                                                                                               ccgcGTTGGCCAGCAGGCGCGGGCGGTCGGCGTCGAAGTCGGGGAAGTCCAGGTGGGTGTGGGTGTCGATCAGGCGCATGTTCAGGCCACTGAGATGCGTTGCTTGAAGGTGCGGCCAACGGCGTGCACGCCGGGTTCATAGCGTTTTTc  >  1:22386/1‑150 (MQ=255)
                                                                                                                                               |                                                                                                                                                     
CTCAGGGCGATGCTGGTCGAGGTAGACCGGGTGCAGGCCAAGCGCGGCGAACAGGCGCGGCTCGGCGCAGGCCAGTTCCCACACCCGCTGGAAATTCGCCTGGTACACCCCCAGCACCACCATTCGCTCCAGCCCACGCGCCGCCGCGTTGGCCAGCAGGCGCGGGCGGTCGGCGTCGAAGTCGGGGAAGTCCAGGTGGGTGTGGGTGTCGATCAGGCGCATGTTCAGGCCACTGAGATGCGTTGCTTGAAGGTGCGGCCAACGGCGTGCACGCCGGGTTCATAGCGTTTTTC  >  NC_002947/907992‑908284

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