Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,574,303 C→A 9.8% R73S (CGC→AGC)  pcaC → 4‑carboxymuconolactone decarboxylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,574,3030CA 9.8% 111.0 / 3.0 41R73S (CGC→AGC) pcaC4‑carboxymuconolactone decarboxylase
Reads supporting (aligned to +/- strand):  ref base C (23/14);  new base A (2/2);  total (25/16)
Fisher's exact test for biased strand distribution p-value = 6.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.68e-01

GACCGCAGCCTGGAGAAGCTCAACGACTTCAATGGCGAGTTCCAGGAAATGATCACCCGCCACGCCTGGGGTGACATCTGGACCCGCCCGGGGCTGCCGCGCCATACCCGCAGCCTGATCACCATCGCCATGCTGATTGG‑CATGAACCGCAACGACGAGCTGAAGCTGCACCTGCGTGCGGCGGCCAACAATGGCGTGACCCGCGACGAGATCAAGGAAGTGCTGATGCAGAGCGCGATCTACTGCGGCATTCCGGCGGCCAATGCCACGTTCCACCTGGCTGAGTCGGTGTGGGAT  >  NC_002947/1574156‑1574452
                                                                                                                                                    |                                                                                                                                                     
gACCGCAGCCTGGAGAAGCTCAACGACTTCAAGGGCGAGTTCCAGGAAATGATCACCCGCCACGCCTGGGGTGACATCTGGACCCGCCCGGGGCTGCCGCGCCATACCCGCAGCCTGATCACCATCGCCATGCTGATTGG‑CATGAACCGc                                                                                                                                                     >  2:625628/1‑150 (MQ=255)
           ggAGAAGCTCAACGACTTCAATGGCGAGTTCCAGGAAATGATCACCCGCCACGCCTGGGGTGACATCTGGACCCGCCCGGGGCTGCCGCGCCATACCCGCAGCCTGATCACCATCGCCATGCTGATTGG‑CATGAACCGCAACGACGAGCt                                                                                                                                          <  1:65023/150‑1 (MQ=255)
                   tCAACGACTTCAATGGCGAGTTCCAGGAAATGATAACCCGCCACGCCTGGGGGGACATCTGGACCCGCCCGGGGCTGCCGCGCCATACCCGCAGCCTGATCACCATCGCCATGCTGATTGG‑CATGAACCGCCACGACGAGCGGAAGCTGc                                                                                                                                  >  2:73427/1‑150 (MQ=255)
                              aaTGGCGAGTTCCAGGAAATGATCACCCGCCACGCCTGGGGTGACATCTGGACCCGCCCGGGGCTGCCGCGCCATACCCGCAGCCTGATCACCATCGCCATGCTGATTGG‑CAGGAATCGCTACGACGAGCTGAAGCTGCACCtgcgtgcg                                                                                                                       >  2:819124/1‑150 (MQ=255)
                                           aGGAAATGATCACCCGCCACGCCTGGGGTGACATCTGGACCCGCCCGGGGCTGCCGCGCCATACCCGCAGCCTGATCACCATCGCCATGCTGATTGG‑CATGAACCGCAACGACGAGCTGAAGCTGCACCTGCGTGCGGCGGCCAACAATg                                                                                                          <  1:149046/150‑1 (MQ=255)
                                              aaaTGATCACCCGCCACGCCTGGGGTGACATCTGGACCCGCCCGGGGCTGCCGCGCCATACCCGCAGCCTGATCACCATCGCCATGCTGATTGG‑CATGAACCGCAACGACGAGCTGAAGCTGCACCTGCGTGCGGCGGCCAACAATGGCg                                                                                                       >  1:974886/1‑150 (MQ=255)
                                                 tGATCACCCGCCACGCCTGGGGTGACATCTGGACCCGCCCGGGGCTGCCGCGCCATACCCGCAGCCTGATCACCATCGCCATGCTGATTGG‑CATGAACCGCAACGACGAGCTGAAGCTGCACCTGCGTGCGGCGGCCAACAATGGCGTGa                                                                                                    >  1:985297/1‑150 (MQ=255)
                                                    tCACCCGCCACGCCTGGGGTGACATCTGGACCCGCCCGGGGCTGCCGCGCCATACCCGCAGCCTGATCACCATCGCCATGCTGATTGG‑CATGAACCGCAACGACGAGCTGAAGCTGCACCTGCGTGCGGCGGCCAACAATGGCGTGAccc                                                                                                 >  2:969787/1‑150 (MQ=255)
                                                    tCACCCGCCACGCCTGGGGTGACATCTGGACCCGCCCGGGGCTGCCGCGCCATACCCGCAGCCTGATCACCATCGCCATGCTGATTGG‑CATGAACCGCAACGACGAGCTGAAGCTGCACCTGCGTGCGGCGGCCAACAATGGCGTGAccc                                                                                                 >  2:329991/1‑150 (MQ=255)
                                                       cccGCCACGCCTGGGGTGACATCTGGACCCGCCCGGGGCTGCCGCGCCATACCCGCAGCCTGATCACCATCGCCATGCTGATTGG‑CATGAACCGCAACGACGAGCTGAAGCTGCACCTGCGTGCGGCGGCCAACAATGGCGTGACCcgcg                                                                                              >  2:423750/1‑150 (MQ=255)
                                                        ccGCCACGCCTGGGGTGACATCTGGACCCGCCCGGGGCTGCCGCGCCATACCCGCAGCCTGATCACCATCGCCATGCTGATTGG‑CATGAACCGCAACGACGAGCTGAAGCTGCACCTGCGTGCGGCGGCCAACAATGGCGTGACCCGcga                                                                                             >  1:994328/1‑150 (MQ=255)
                                                         cGCCACGCCTGGGGTGACATCTGGACCCGCCCGGGGCTGCAGCGCCATACCCGCAGCCTGATCACCATCGCCTTGCTGATTGG‑CATGAACCGCAACGACGAGCTGACGCTGCACCTGCGTGCGGCGGCCAACAATGGCGTGACCCGcgac                                                                                            <  1:74034/150‑1 (MQ=255)
                                                                  tGGGGTGACATCTGGACCCGCCCGGGGCTGCCGCGCCATACCCGCAGCCTGATCACCATCGCCATGCTGATTGG‑CATGAACCGCAACGACGAGCTGAAGCTGCACCTGCGTGCGGCGGCCAACAATGGCGTGACCCGCGACGAGATCAAg                                                                                   <  1:361272/150‑1 (MQ=255)
                                                                     ggTGACATCTGGACCCGCCCGGGGCTGCCGCGCCATACCCGCAGCCTGATCACCATCGCCATGCTGATTGG‑CATGAACCGCAACGACGAGCTGAAGCTGCACCTGCGTGCGGCGGCCAACAATGGCGTGACCCGCGACGAGATCAAGGaa                                                                                >  1:102952/1‑150 (MQ=255)
                                                                           aTCTGGACCCGCCCGGGGCTGCCGCGCCATACCCGCAGCCTGATCACCATCGCCATGCTGATTGG‑CATGAACCGCAACGACGAGCTGAAGCTGCACCTGCGTGCGGCGGCCAACAATGGCGTGACCCGCGTCGAGATCAAGGAAGTGCTg                                                                          >  1:335518/1‑150 (MQ=255)
                                                                            tCTGGACCCGCCCGGGGCTGCCGCGCCATACCCGCAGCCTGATCACCATCGCCATGCTGATTGG‑CATGAACCGCAACGACGAGCTGAAGCTGCACCTGCGTGCGGCGGCCAACAATGGCGTGACCCGCGACGAGATCAAGGAAGTGCTGa                                                                         >  2:335233/1‑150 (MQ=255)
                                                                             cTGGACCCGCCCGGGGCTGCCGCGCCATACCCGCAGCCTGATCACCATCGCCATGCTGATTGG‑CATGAACCGCAACGACGAGCTGAAGCTGCACCTGCGTGCGGCGGCCAACAATGGCGTGACCCGCGACGAGATCAAGGAAGTGCTGAt                                                                        >  1:716036/1‑150 (MQ=255)
                                                                              tGGACCCGCCCGGGGCTGCCGCGCCATACCCGCAGCCTGATCACCATCGCCATGCTGATTGG‑CATGAACCGCAACGACGAGCTGAAGCTGCACCTGCGTGCGGCGGCCAACAATGGCGTGACCCGCGACGAGATCAAGGAAGTGCTGATg                                                                       <  2:335518/150‑1 (MQ=255)
                                                                              tGGACCCGCCCGGGGCGCCCGCGCCCTACCCGCAGCCTGATCACCATCGCCATGCTGATTGG‑CATGAACCCCAACGACGAGCTGAAGCTGCACCTGCGTGCGGCGGCCAACAATGGCGTGACCCGCGACGAGATCAAGGAAGTGCTGATg                                                                       <  1:885742/150‑1 (MQ=255)
                                                                                gACCCGCCCGGGGCTGCCGCGCCATACCCGCAGCCTGATCACCATCGCCATGCTGATTGG‑CATGAACCGCAACGACGAGCTGAAGCTGCACCTGCGTGCGGCGGCCAACAATGGCGTGACCCGCGACGAGATCAAGGAAGTGCTGATGCa                                                                     <  2:994328/150‑1 (MQ=255)
                                                                                 aCCCGCCCGGGGCTGCCGCGCCATACCCGCAGCCTGATCACCATCGCCATGCTGATTGG‑CATGAACCGCAACGACGAGCTGAAGCTGCACCTGCGTGCCGCGGCCAACAATGGCGTGACCCGCGACGAGATCAAGGAAGTGCTGATGCag                                                                    >  2:1050815/1‑150 (MQ=255)
                                                                                  cccgcccgGGGCTGCCGCGCCATACCCGCAGCCTGATCACCATCGCCATGCTGATTGG‑CATGAACCGCAACGACGAGCTGAAGCTGCACCTGCGTGCGGCGGCCAACAATGGCGTGACCCGCGACGAGATCAAGGAAGTGCTGATGCaga                                                                   >  1:565156/1‑150 (MQ=255)
                                                                                   ccgcccgGGGCTGCCGCGCCATACCCCCAGCCTGATCACCATCGCCATGAT‑ATTGGTCATGAACCGCAACGACGAGCTGAAGCTGCACCTGCGTGCGGCTGCCAACAAGGGCGTGACCCGCGACCAGATCCAGGAAGTGCCGAATCagag                                                                  >  1:357258/1‑150 (MQ=255)
                                                                                               gCCGCGCCATACCCGCAGCCTGATCACCATCGCCATGCTGATTGG‑CATGAACCGCAACGACGAGCTGAAGCTGCACCTGCGTGCGGCGGCCAACAATGGCGTGACCCGCGACGAGATCAAGGAAGTGCTGATGCAGAGCGCGATCTACTg                                                      <  2:974886/150‑1 (MQ=255)
                                                                                                 cgcgCCATACCCGCAGCCTGATCACCATCGCCATGCTGATTGG‑CATGAACCGCAACGACGAGCTGAAGCTGCACCTGCGTGCGGCGGCCAACAATGGCGTGACCCGCGACGAGATCAAGGAAGTGCTGATGCAGAGCGCGATCTACTGCg                                                    >  2:52811/1‑150 (MQ=255)
                                                                                                               aGCCTGATCACCATCGCCATGCTGATTGG‑CATGAACCGCAACGACGAGCTGAAGCTGCACCTGCGTGCGGCGGCCAACAATGGCGTGACCCGCGACGAGATCAAGGAAGTGCTGATGCAGAGCGCGATCTACTGCGGCATTCCGGCGGcc                                      <  2:102952/150‑1 (MQ=255)
                                                                                                                    gATCACCATCGCCATGCTGATTGG‑CATGAACCGCAACGACGAGCTGAAGCTGCACCTGCGTGCGGCGGCCAACAATGGCGTGACCCGCGACGAGATCAAGGAAGTGCTGATGCAGAGCGCGATCTACTGCGGCATTCCGGCGGCCAATGc                                 >  1:819272/1‑150 (MQ=255)
                                                                                                                        aCCATCGCCATGCTGATTGG‑CATGAACCGCAACGACGAGCTGAAGCTGCACCTGCGTGCGGCGGCCAACAATGGCGTGACCCGCGACGAGATCAAGGAAGTGCTGATGCAGAGCGCGATCTACTGCGGCa                                                 >  1:1004305/1‑130 (MQ=255)
                                                                                                                        aCCATCGCCATGCTGATTGG‑CATGAACCGCAACGACGAGCTGAAGCTGCACCTGCGTGCGGCGGCCAACAATGGCGTGACCCGCGACGAGATCAAGGAAGTGCTGATGCAGAGCGCGATCTACTGCGGCa                                                 <  2:1004305/130‑1 (MQ=255)
                                                                                                                         ccATCGCCATGCTGATTGG‑CATGAACCGCAACGACGAGCTGAAGCTGCACCTGCGTGCGGCGGCCAACAATGGAGTTTCCCGCGACGCGATCAACGGAGTCATAGTGCGCAGCGCTATCTACTGCGGCATTCCGGCGGCCAATGCCACGt                            >  1:805741/1‑150 (MQ=255)
                                                                                                                          cATCGCCATGCTGATTGG‑CATGAACCGCAACGACGAGCTGAAGCTGCACCTGCGTGCGGCGGCCAACAATGGCGTGACCCGCGACGAGATCAAGGAAGTGCTGATGCAGAGCGCGATCTACTGCGGCATTCCGGCGGCCAATGCCACGtt                           >  2:650378/1‑150 (MQ=255)
                                                                                                                            tCGCCATGCTGATTGG‑CATGAACCGCAACGACGAGCTGAAGCTGCACCTGCGTGCGGCGGCCAACAATGGCGTGACCCGCGACGAGATCAAGGAAGTGCTGATGCAGAGCGCGATCTACTGCGGCATTCCGGCGGCCAATGCCACGTTcc                         <  1:969787/150‑1 (MQ=255)
                                                                                                                               ccATGCTGATTGG‑CATGAACCGCAACGACGAGCTGAAGCTGCACCTGCGTGCGGCGGCCAACAATGGCGTGACCCGCGACGAGATCAAGGAAGTGCTGATGCAGAGCGCGATCTACTGCGGCATTCCGGCGGCCAATGCCACGTTCCAcc                      <  1:219648/150‑1 (MQ=255)
                                                                                                                               ccATGCTGATTGG‑CATGAACAGCAACGACGAGCTGAAGCTGCACCTGCGTGCGGCGGCCAACAATGGCGTGACCCGCGACGAGATCAAGGAAGTGCt                                                                           >  1:441327/1‑97 (MQ=255)
                                                                                                                               ccATGCTGATTGG‑CATGAACAGCAACGACGAGCTGAAGCTGCACCTGCGTGCGGCGGCCAACAATGGCGTGACCCGCGACGAGATCAAGGAAGTGCt                                                                           <  2:441327/97‑1 (MQ=255)
                                                                                                                               ccATGCTGATTGG‑CATGAACAGCAACGACGAGCTGAAGCTGCACCTGCGTGCGGCGGCCAACAATGGCGTGACCCGCGACGAGATCAAGGAAGTGCt                                                                           <  2:441595/97‑1 (MQ=255)
                                                                                                                               ccATGCTGATTGG‑CATGAACAGCAACGACGAGCTGAAGCTGCACCTGCGTGCGGCGGCCAACAATGGCGTGACCCGCGACGAGATCAAGGAAGTGCt                                                                           >  1:441595/1‑97 (MQ=255)
                                                                                                                                cATGCTGATTGG‑CATGAACCGCAACGACGAGCTGAAGCTGCACCTGCGTGCGGCGGCCAACAATGGCATGACCCGCGACGAGATCAAGGAAGTGCTGATGCAGAGCGCGATCTACTGCGGCATTCCGCCGGCCAATGCCACGTTCCCCCt                     >  2:313748/1‑150 (MQ=255)
                                                                                                                                              aTGAACCGCAACGACGAGCTGAAGCTGCACCTGCGTGCGGCGGCCAACAATGGCGTGACCCGCGACGAGATCAAGGAAGTGCTGATGCAGAGCGCGATCTACTGCGGCATTCCGGCGGCCAATGCCACGTTCCACCTGGCTGAGTCGgtg        >  1:1125166/1‑150 (MQ=255)
                                                                                                                                               tGAACCGCAACGACGAGCTGAAGCTGCACCTGCGTGCGGCGGCAAAAAATGGCGTGACCCGCGACGAGATCAAAGACGTGCTGATGCAGAGCGCGATCTACTGCGGCATTCCGGCGGCCAATGCCACGTTCCACCTGGCTGAGTCGgtgt       <  1:625628/150‑1 (MQ=255)
                                                                                                                                                gAACCGCAACGACGAGCTGAAGCTGCACCTGCGTGCGGCGGCCAACAATGGCGTGACCCGCGACGAGATCAAGGAAGTGCTGATGCAGAGCGCGATCTACTGCGGCATTCCGGCGGCCAATGCCACGTTCCACCTGGCTGAGTCGGtgtg      <  2:985297/150‑1 (MQ=255)
                                                                                                                                                    cGCAACGACGAGCTGAAGCTGCACCTGCGTGCGGCGGCCAACAATGGCGTGACCCGCGACGAGATCAAGGAAGTGCTGATGCAGAGCGCGATCTACTGCGGCATTCCGGCGGCCAATGCCACGTTCCACCTGGCTGAGTCGGTGTGGGAt  <  1:52811/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GACCGCAGCCTGGAGAAGCTCAACGACTTCAATGGCGAGTTCCAGGAAATGATCACCCGCCACGCCTGGGGTGACATCTGGACCCGCCCGGGGCTGCCGCGCCATACCCGCAGCCTGATCACCATCGCCATGCTGATTGG‑CATGAACCGCAACGACGAGCTGAAGCTGCACCTGCGTGCGGCGGCCAACAATGGCGTGACCCGCGACGAGATCAAGGAAGTGCTGATGCAGAGCGCGATCTACTGCGGCATTCCGGCGGCCAATGCCACGTTCCACCTGGCTGAGTCGGTGTGGGAT  >  NC_002947/1574156‑1574452

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