Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,572,891 T→G 9.1% L338R (CTG→CGG) ‡ pcaB → 3‑carboxy‑cis,cis‑muconate cycloisomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,572,8910TG 9.1% 160.1 / 7.9 55L338R (CTG→CGG) ‡pcaB3‑carboxy‑cis,cis‑muconate cycloisomerase
Reads supporting (aligned to +/- strand):  ref base T (25/25);  new base G (3/2);  total (28/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.37e-01

GTGTTGATCGGTGCCGCGACCCGCGTGCCGGGCCTGCTGTCGACGCTGTTCGCAGCCATGCCTCAGGAGCACGAACGCAGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCCGCGCGCATGCGCCGTAACCTCGACCTGACCCAAGGCCTGGTGCTGGCCGAAGCGGTGAGCATCGTCCTCGCCCAGCGTCTGGGTCGCGACCGTGCCCACC  >  NC_002947/1572746‑1573034
                                                                                                                                                  |                                                                                                                                                
gtgtTGATCGGTGCCGCGACCCGCGTGCCGGGCCTGCTGTCGACGCTGTTCGCAGCCATGCCTCAGGAGCACGAACGCAGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTgcg                                                                                                                                               >  2:152071/1‑149 (MQ=255)
      aTCGGTGCCGCGACCCGCGTGCCGGGCCTGCTGTCGACGCTGTTCGCAGCCATGCCTCAGGAGCACGAACGCAGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCAGGGCGCCAga                                                                                                                                          <  4:48085/148‑2 (MQ=255)
          gTGCCGCGACCCGCGTGCCGGGCCTGCTGTCGACGCTGTTCGCAGCCATGCCTCAGGAGCACGAACGCAGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCaa                                                                                                                                     <  1:152071/149‑1 (MQ=255)
                cgACCCGCGTGCCGGGCCTGCTGTCGACGCTGTTCGCAGCCATGCCTCAGGAGCACGAACGCAGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGAtt                                                                                                                               >  4:167350/1‑149 (MQ=255)
                 gACCCGCGTGCCGGGCCTGCTGTCGACGCTGTTCGCAGCCATGCCTCAGGAGCACGAACGCAGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTg                                                                                                                              >  2:58205/1‑149 (MQ=255)
                    ccGCGTGCCGGGCCTGCTGTCGACGCTGTTCGCAGCCATGCCTCAGGAGCACGAACGCAGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCAGGGCGCCa                                                                                                                                            <  5:215825/132‑1 (MQ=255)
                    ccGCGTGCCGGGCCTGCTGTCGACGCTGTTCGCAGCCATGCCTCAGGAGCACGAACGCAGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCAGGGCGCCa                                                                                                                                            >  6:215825/1‑132 (MQ=255)
                      gcgTGCCGGGCCTGCTGTCGACGCTGTTCGCAGCCATGCCTCAGGAGCACGAACGCAGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAg                                                                                                                         >  6:194854/1‑149 (MQ=255)
                       cgTGCCGGGCCTGCTGTCGACGCTGTTCGCAGCCATGCCTCAGGAGCACGAACGCAGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAgg                                                                                                                        <  4:297175/149‑1 (MQ=255)
                        gTGCCGGGCCTGCTGTCGACGCTGTTCGCAGCCATGCCTCAGGAGCACGAACGCAGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAggg                                                                                                                       <  4:285764/149‑1 (MQ=255)
                             gggCCTGCTGTCGACGCTGTTCGCAGCCATGCCTCAGGAGCACGAACGCAGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATgg                                                                                                                  <  3:61716/149‑1 (MQ=255)
                                ccTGCTGTCGACGCTGTTCGCAGCCATGCCTCAGGAGCACGAACGCAGCCTGGGCCTATGGCATGCCGAGTGGGAACCCCTGCCGGATATCTGATGCCTGGACTCTGCCGCCCTGCGCCAGGCTCAATTGATTGCCGAGGGCATGGAgg                                                                                                               >  3:67949/1‑149 (MQ=255)
                                 ctgctgTCGACGCTGTTCGCAGCCATGCCTCAGGAGCACGAACGCAGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGt                                                                                                              <  1:58205/149‑1 (MQ=255)
                                     tgTCGACGCTGTTCGCAGCCATGCCTCAGGAGCACGAACGCAGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAg                                                                                                                         >  6:258116/1‑134 (MQ=255)
                                     tgTCGACGCTGTTCGCAGCCATGCCTCAGGAGCACGAACGCAGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAg                                                                                                                         <  5:258116/134‑1 (MQ=255)
                                       tCGACGCTGTTCGCAGCCATGCCTCAGGAGCACGAACGCAGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGAt                                                                                                          <  5:194854/147‑1 (MQ=255)
                                       tCGACGCTGTTCGCAGCCATGCCTCAGGAGCACGAACGCAGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGc                                                                                                        >  7:102593/1‑149 (MQ=255)
                                       tCGACGCTGTTCGCAGCCATGCCTCACGAGCACGAACGCAGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATg                                                                                                         <  2:176630/148‑1 (MQ=255)
                                         gACGCTGTTCGCAGCCATGCCTCAGGAGCACGAACGCAGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAg                                                                                                                         <  1:261979/130‑1 (MQ=255)
                                         gACGCTGTTCGCAGCCATGCCTCAGGAGCACGAACGCAGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAg                                                                                                                         >  2:261979/1‑130 (MQ=255)
                                         gACGCTGTTCGCAGCCATGCCTCAGGAGCACGAACGCAGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCcg                                                                                                      >  1:220274/1‑149 (MQ=255)
                                          aCGCTGTTCGCAGCCATGCCTCAGGAGCACGAACGCAGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCcgc                                                                                                     <  2:220274/149‑1 (MQ=255)
                                                ttCGCAGCCATGCCTCAGGAGCACGAACGCAGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCCGCGCGCAt                                                                                               >  4:102578/1‑149 (MQ=255)
                                                   gCAGCCATGCCTCAGGAGCACGAACGCAGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCCGCGCGCATgc                                                                                             <  8:102593/148‑1 (MQ=255)
                                                            ccTCAGGAGCACGAACGCAGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGa                                                                                                           <  6:9572/125‑1 (MQ=255)
                                                            ccTCAGGAGCACGAACGCAGCCTGGGACTATGGCATGCCGAGTGCGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGa                                                                                                           >  5:9572/1‑125 (MQ=255)
                                                              tCAGGAGCACGAACGCAGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCCGCGCGCATGCGCCGTAACCTc                                                                                  >  2:217984/1‑148 (MQ=255)
                                                                    gcacgcacgcacCCTGGGCCGAAGGCATGGCGCGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCCGCGCGCATGCGCCGTAACCTCGACCTGa                                                                           <  4:67949/149‑1 (MQ=255)
                                                                     cACGAACGCAGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCCGCGCGCATGCGCCGTAACCTCGACCTGAc                                                                          >  3:89720/1‑149 (MQ=255)
                                                                        gAACGCAGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCCGCGCGCATgc                                                                                             <  2:98560/127‑1 (MQ=255)
                                                                        gAACGCAGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCAAGGCTCAAGTGATTGCCGAGGGCATGGAGGGGGAGGCCGCGCGCAGGCGCCGTCTCCTCGAcaggtccc                                                                        >  8:255990/1‑141 (MQ=255)
                                                                              aGCCTGGGCCTATGGCATGCCGAGTGAGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCCGCGCGCATGCGcc                                                                                          >  1:3328/1‑125 (MQ=255)
                                                                              aGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCCGCGCGCATGCGcc                                                                                          <  2:3328/124‑1 (MQ=255)
                                                                              aGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCCGCGCGCATGCGCCGTAACCTCGACCTGACCCAAGGCCt                                                                 >  2:342725/1‑149 (MQ=255)
                                                                                       cTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCCGCGCGCATGCGCCGTAACCTCGACCTGACCCAAGGCCTGGTGCTgg                                                         <  5:145145/148‑1 (MQ=255)
                                                                                       cTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCCGCGCGCATGCGCCGTAACCTCGACCTGACCCAAGGCCTGGTGCTgg                                                         >  6:145145/1‑148 (MQ=255)
                                                                                                 ccGCGTGGGAAACCCTGCCGGAAATCTGCTGCCTGGGCCCGGGCGCCCTGCGGCATGCGCAAGTGAATGCCGATGGCATGGATGTGGATGCCGCGCGCAAGCGCCGTAACCTCGACCTGACCCAAGGCCTGGTGCTGGCCGAAGCGGTg                                              <  5:379197/149‑1 (MQ=255)
                                                                                                 ccGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCCGCGCGCATGCGCCGTAACCTCGACCTGACCCAAGGCCTGGTGCTGGCCGAAGCGGt                                               <  4:89720/148‑1 (MQ=255)
                                                                                                   gAGTGGGAACCCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCCGCGCGCATGCGCCGTAACCTCGACCTGACCCAAGGCCTGGTGCTGGCCGAAGCGGTGAg                                            >  2:264000/1‑149 (MQ=255)
                                                                                                    aGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCCGCGCGCATGCGCCGTAACCTCGACCTGACCCAAGGCCTGGTGCTGGCCGAAGCGGTGAGc                                           <  3:167350/149‑1 (MQ=255)
                                                                                                         ggAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCCGCGCGCATGCGCCGTAACCTCGACCTGACCCAAGGCCTGGTGCTGGCCGAAGCGGTGAGCATCg                                       >  6:222769/1‑149 (MQ=255)
                                                                                                             aCCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCCGCGCGCATGCGCCGTAACCTCGACCTGACCCAAGGCCTGGTGCTGGCCGAAGCGGTGAGCATCGTCCt                                   >  5:148788/1‑149 (MQ=255)
                                                                                                               ccTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCCGCGCGCATGCGCCGTAACCTCGACCTGACCCAAGGCCTGGTGCTGGCCGAAGCGGTGAGCATCGTCCTCg                                 >  8:81987/1‑149 (MQ=255)
                                                                                                                     ggATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCCGCGCGCATGCGCCGTAACCTCGACCTGACCCAAGGCCTGGTGCTGGCCGAAGCGGTGAGCATCGTCCTCGCCCAGc                           >  7:207317/1‑149 (MQ=255)
                                                                                                                        tatCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCCGCGCGCATGCGCCGTAACCTCGACCTGACCCAAGGCCTGGTGCTGGCCGAAGCGGTGAGCATCGTCCTCGCCCAGCGTc                        >  4:285711/1‑149 (MQ=255)
                                                                                                                          tCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCCGCGCGCATGCGCCGTAACCTCGACCTGACCCAAGGCCTGGTGCTGGCCGAAGCGGTGAGCATCGTCCTCGCCCAGCGTCTg                     >  1:149696/1‑149 (MQ=255)
                                                                                                                              ctgcCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCCGCGCGCATGCGCCGTAACCTCGACCTGACCCAAGGCCTGGTGCTGGCCGAAGCGGTGAGCATCGTCCTCGCCCAGCGTCTGGGt                  <  3:102578/148‑1 (MQ=255)
                                                                                                                              ctgcCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCCGCGCGCATGCGCCGTAACCTCGACCTGACCCAAGGCCTGGTGCTGGCCGAAGCGGTGAGCATCGTCCTCGCCCAGCGTCTGGGt                  <  1:217984/148‑1 (MQ=255)
                                                                                                                              ctgcCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCCGCGCGCATGCGCCGTAACCTCGACCTGACCCAAGGCCTGGTGCTGGCCGAAGCGGTGAGCATCGTCCTCGCCCAGCGTCTGGGTc                 >  5:85537/1‑149 (MQ=255)
                                                                                                                              ctgcCTGGTCGCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCCGCGCGCATGCGCCGTAACCTCGACCTGACCCAAGGCCTGGTGCTGGCCGAAGCGGTGAGCATCGTCCTCGCCCAGCGTCTgg                    <  7:255990/146‑1 (MQ=255)
                                                                                                                                gcCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGGCGAGGGCATGGAGGTGGATGCCGCGCGCATGCGCCGTAACCTCGACCTGACCCAAGGCCTGGTGCTGGCCGAAGCGGTGAGCATCGTCCTCGCCCAGCGTCTGGGTcgc               <  4:54164/149‑1 (MQ=255)
                                                                                                                                         ctGGCGCAGGGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCCGCGCGCATGCGCCGTAACCTCGACCTGACCCAAGGCCTGGTGCTGGCCGAAGCGGTGAGCATCGTCCTCGCCCAGCGTCATAGATCGCGACCGTGc       >  1:312928/1‑149 (MQ=255)
                                                                                                                                         ctGGCGCAGGGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCCGCGCGCATGCGCCGTAACCTCGACCTGACCCAAGGCCTGGTGCTGGCCGAAGCGGTGAGCATCGTCCTCGCCCAGCGTCTGGGTCGCGACCGTGcc      >  3:23867/1‑149 (MQ=255)
                                                                                                                                            gcgcCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCCGCGCGCATGCGCCGTAACCTCGACCTGACCCAAGGCCTGGTGCTGGCCGAAGCGGTGAGCATCGTCCTCGCCCAGCGTCTGGGTCGCGACCGTGCccac   <  5:222769/149‑1 (MQ=255)
                                                                                                                                             cgcCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCCGCGCGCATGCGCCGTAACCTCGACCTGACCCAAGGCCTGGTGCTGGCCGAAGCGCAGAGCATCGTCCTCGCCCAGCGTCTGGGTCGCGACCGTGCccacc  <  1:264000/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
GTGTTGATCGGTGCCGCGACCCGCGTGCCGGGCCTGCTGTCGACGCTGTTCGCAGCCATGCCTCAGGAGCACGAACGCAGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCCGCGCGCATGCGCCGTAACCTCGACCTGACCCAAGGCCTGGTGCTGGCCGAAGCGGTGAGCATCGTCCTCGCCCAGCGTCTGGGTCGCGACCGTGCCCACC  >  NC_002947/1572746‑1573034

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