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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947908,1320CG 8.0% 149.6 / 2.9 49A28P (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

ACCAGGTGCTCAGGGCGATGCTGGTCGAGGTAGACCGGGTGCAGGCCAAGCGCGGCGAACAGGCGCGGCTCGGCGCAGGCCAGTTCCCACACCCGCTGGAAATTCGCCTGGTACACCCCCAGCACCACCATTCGCTCCAGCCCACGCGCCGCCGCGTTGGCCAGCAGGCGCGGGCGGTCGGCGTCGAAGTCGGGGAAGTCCAGGTGGGTGTGGGTGTCGATCAGGCGCATGTTCAGGCCACTGAGATGCGTTGCTTGAAGGTGCGGCCAACGGCGTGCACGCCGG  >  NC_002947/907984‑908268
                                                                                                                                                    |                                                                                                                                        
aCCAGGTGCTCAGGGCGATGCTGGTCGAGGTAGACCGGGTGCAGGCCAAGCGCGGCGAACAGGCGCGGCTCGGCGCAGGCCAGTTCCCACACCCGCTGGAAATTCGCCTGGTACACCCCCAGCACCACCATTCGCTCCAGCCCACGcgcc                                                                                                                                         >  1:18318/1‑150 (MQ=255)
        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                                                                                                                            >  2:451288/1‑131 (MQ=255)
                              tAGACCGGGTGCAGGCCAAGCGCGGCGAACAGGCGCGGCTCGGCGCAGGCCAGTTCCCACACCCGCTGGAAATTCGCCTGGTACACCCCCAGCACCACCATTCGCTCCAGCCAACGCGGCGGCGCGTGGGctg                                                                                                                            <  1:451288/133‑3 (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:612488/150‑1 (MQ=255)
                                                               cgcgGCTCGGCGCAGGCCAGTTCCCACACCCGCTGGAAATTCGCCTGGTACACCCCCAGCACCACCATTCGCTCCAGCCCACGCGCCGCCGCGTTGGCCAGCAGGCGCGGGCGGTCGGCGTCGAAGTCGGGGAAGTCCAGgtgggtgtgg                                                                          <  1:251729/150‑1 (MQ=255)
                                                               cgcgGCTCGGCGCAGGCCAGTTCCCACACCCGCTGGAAATTCGCCTGGTACACCCCCAGCACCACCATTCGCTCCAGCCCACGCGCCGCCGCGTTGGCCAGCAGGCGCGGGCGGTCGGCGTCGAAGTCGGGGAAGTCCAGgtgggtgtgg                                                                          <  2:462340/150‑1 (MQ=255)
                                                                         cgcAGGCCAGTTCCCACACCCGCTGGAAATTCGCCTGGTACACCCCCAGCACCACCATTCGCTCCAGCCCACGCGCCGCCGCGTTGGCCAGCAGGCGCGGGCGGTCGGCGTCGAAGTCGGGGAAGTCCAGGTGGGTGTGGGTGTCGATCa                                                                >  1:146890/1‑150 (MQ=255)
                                                                                cAGTTCCCACACCCGCTGGAAATTCGCCTGGTACACCCCCAGCACCACCATTCGCTCCAGCCCACGCGCCGCCGCGTTGGCCAGCAGGCGCGGGCGGTCGGCGTCGAAGTCGGGGAAGTCCAGGTGGGTGTGGGTGTCGATCAgg                                                              >  2:265040/1‑145 (MQ=255)
                                                                                cAGTTCCCACACCCGCTGGAAATTCGCCTGGTACACCCCCAGCACCACCATTCGCTCCAGCCCACGCGCCGCCGCGTTGGCCAGCAGGCGCGGGCGGTCGGCGTCGAAGTCGGGGAAGTCCAGGTGGGTGTGGGTGTCGATCAgg                                                              <  1:265040/145‑1 (MQ=255)
                                                                                cAGTTCCCACACCCGCTGGAAATTCGCCTGGTACACCCCCAGCACCACCATTCGCTCCAGCCCACGCGCCGCCGCGTTGGCCAGCAGGCGCGGGCGGTCGGCGTCGAAGTCGGGGAAGTCCAGGTGGGTGTGGGTGTCGATCAGGCGCAt                                                         >  1:34987/1‑150 (MQ=255)
                                                                                 aGTTCCCACACCCGCTGGAAATTCGCCTGGTACACCCCCAGCACCACCATTCGCTCCAGCCCACGCGGCGGCGCGtgg                                                                                                                                <  1:458241/78‑3 (MQ=255)
                                                                                 aGTTCCCACACCCGCTGGAAATTCGCCTGGTACACCCCCAGCACCACCATTCGCTCCAGCCCACGCGGCGGCGCGtgg                                                                                                                                >  2:458241/1‑76 (MQ=255)
                                                                                   ttCCCACACCCGCTGGAAATTCGCCTGGTACACCCCCAGCACCACCATTCGCTCCAGCCCACGCGCCGCCGCGTTGGCCAGCAGGCGCGGGCGGTCGGCGTCGAAGTCGGGGAAGTCCAGGTGGGTGTGGGTGTCGATCAg                                                               >  2:595925/1‑141 (MQ=255)
                                                                                   ttCCCACACCCGCTGGAAATTCGCCTGGTACACCCCCAGCACCACCATTCGCTCCAGCCCACGCGCCGCCGCGTTGGCCAGCAGGCGCGGGCGGTCGGCGTCGAAGTCGGGGAAGTCCAGGTGGGTGTGGGTGTCGATCAg                                                               <  1:595925/141‑1 (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                                                                                  <  1:573826/86‑1 (MQ=255)
                                                                                                                       cAGCACCACCATTCGCTCCAGCCCACGCGCCGCCGCGTTGGCCAGCAGGCGCGGGCGGTCGGCGTCGAAGTCGGGGAAGTCCAggt                                                                                  >  2:573826/1‑86 (MQ=255)
                                                                                                                              accaTTCGCTCCAGCCCACGCGCCGCCGCGTTGGCCAGCAGGCGCGGGCGGTCGGCGTCGAAGTCGGGGAAGTCCAGGTGGGTGTGGGTGTCGATCAGGCGCATGTTCAGGCCACTGAGATGCGTTGCTTGAAGGTGCGGCCAACGGCGt           <  2:146173/150‑1 (MQ=255)
                                                                                                                                       tCCAGCCCACGCGCCGCCGCGTTGGCCAGCAGGCGCGGGCGGTCGGCGTCGAAGTCGGGGAAGTCCAGGTGGGTGTGGGTGTCGATCAGGCGCATGTTCAGGCCACTGAGATGCGTTGCTTGAAGGTGAGGCCAACGGCGTGAACGTCgg  >  1:299520/1‑150 (MQ=255)
                                                                                                                                              cACGCGCCGCCGCGTTGGCCAGCAGGCGCGGGCGGTCGGCGTCGAAGTCGGGGAAGTCCAggtgg                                                                                <  2:714350/65‑1 (MQ=255)
                                                                                                                                              cACGCGCCGCCGCGTTGGCCAGCAGGCGCGGGCGGTCGGCGTCGAAGTCGGGGAAGTCCAggtgg                                                                                >  1:714350/1‑65 (MQ=255)
                                                                                                                                                    |                                                                                                                                        
ACCAGGTGCTCAGGGCGATGCTGGTCGAGGTAGACCGGGTGCAGGCCAAGCGCGGCGAACAGGCGCGGCTCGGCGCAGGCCAGTTCCCACACCCGCTGGAAATTCGCCTGGTACACCCCCAGCACCACCATTCGCTCCAGCCCACGCGCCGCCGCGTTGGCCAGCAGGCGCGGGCGGTCGGCGTCGAAGTCGGGGAAGTCCAGGTGGGTGTGGGTGTCGATCAGGCGCATGTTCAGGCCACTGAGATGCGTTGCTTGAAGGTGCGGCCAACGGCGTGCACGCCGG  >  NC_002947/907984‑908268

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: