Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,236,191 G→T 14.6% G415C (GGT→TGT) ‡ PP_2834 → D‑galactarate dehydratase/altronate dehydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,236,1910GT14.6% 106.6 / 7.9 41G415C (GGT→TGT) ‡PP_2834D‑galactarate dehydratase/altronate dehydratase
Reads supporting (aligned to +/- strand):  ref base G (20/15);  new base T (2/4);  total (22/19)
Fisher's exact test for biased strand distribution p-value = 3.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.28e-01

TGACCACTATCCTCGAGAAATCCCTGGGCGCAGTGGCCAAGGCAGGCTCCAGCAACCTGGTCGATGTGTATCAGTACGCGCAAGCCGTGCGTGCCCAGGGCCTGGTGTTCATGGACACCCCCGGCTACGACCCGGTCTCGGCCACCGGTCAGGTGGCCGGCGGGGCCAACCTGATCGCCTTCACCACCGGCCGCGGCTCGGCCTACGGCTGTGCCCCGGCGCCATCGATCAAGCTGGCAACCAACAACCGGGTGTTCGAGCACCAGCAAGAGGACATGGACGTCAATTGCGG  >  NC_002947/3236045‑3236336
                                                                                                                                                  |                                                                                                                                                 
tGACCACTATCCTCGAGAAATCCCTGGGCGCAGTGGCCAAGGCAGGCTCCAGCAACCTGGTCGATGTGTATCAGTACGCGCAAGCCGTGCGTGCCCAGGGCCTGGTGTTCATGGACACCCCCGGCTACGACCCGGTCTCGGCCACCGGTc                                                                                                                                                <  1:377368/150‑1 (MQ=255)
        aTCCTCGAGAAATCCCTGGGCGCAGTGGCCAAGGCAGGCTCCAGCAACCTGGTCGATGTGTATCAGTACGCGCAAGCCGTGCGTGCCCAGGGCCTGGTGTTCATGGACACCCCCGGCTACGACCCGGTCTCGGCCACCTGACCGGTGGcc                                                                                                                                        <  2:879538/150‑1 (MQ=255)
            tCGAGAAATCCCTGGGCGCAGTGGCCAAGGCAGGCTCCAGCAACCTGGTCGATGTGTATCAGTACGCGCACGCCGTCCGTGCCCAGGGCCTGGTGTTCATGGACACCCCCGCCTACCACCCGGTCTCGGCCACCGGGCAGGTGGCcggcg                                                                                                                                    >  2:618255/1‑150 (MQ=255)
            tCGAGAAATCCCTGGGCGCAGTGGCCAAGGCAGGCTCCAGCAACCTGGTCGATGTGTATCAGTACGCGCAAGCCGTGCGTGCCCAGGGCCTGGTGTTCATGGACACCCCCGGCTACGACCCGGTCTCGGCCACCTGACCGGTGGCcgaga                                                                                                                                    >  1:1070635/1‑147 (MQ=255)
            tCGAGAAATCCCTGGGCGCAGTGGCCAAGGCAGGCTCCAGCAACCTGGTCGATGTGTATCAGTACGCGCAAGCCGTGCGTGCCCAGGGCCTGGTGTTCATGGACACCCCCGGCTACGACCCGGTCTCGGCCACCTGACCGGTGGCcgaga                                                                                                                                    <  2:1070635/150‑4 (MQ=255)
             cGAGAAATCCCTGGGCGCAGTGGCCAAGGCAGGCTCCAGCAACCTGGTCGATGTGTATCAGTACGCGCAAGCCGTGCGTGCCCAGGGCCTGGTGTTCATGGACACCCCCGGCTACGACCCGGTCTCGGCCACCGGTCAGGTGGCcggcgg                                                                                                                                   <  1:327792/150‑1 (MQ=255)
                gaAAGCCCTGGGCGCAGTGACCAAGGCAGGTTCCAGCAACCTGGTCGATGTGTATCAGTACGCGCAAGCCGTGCGTGCCCAGGGCCTGGTGTTCATGGACACCCCCGGCTACGACCCGGTCCCGGCCACCTGGCCTGTGGCcgagaccgg                                                                                                                                <  2:685306/150‑8 (MQ=255)
                                     cAAGGCAGGCTCCAGCAACCTGGTCGATGTGTATAAGTAAGCGCACGCGCTGCGTGCCCAGGGCCTGGTGTGCATGGACACCCCCGGCTACGACCCGGTCTCGGCCACCGGTCAGGTGGCCGGCGGGGCCAACCTGATCGCCTTcaccac                                                                                                           >  1:1315883/1‑150 (MQ=255)
                                      aaGGCAGGCTCCAGCAACCTGGTCGATGTGTATCAGTACGCGCAAGCCGTGCGTGCCCAGGGCCTGGTGTTCAGGGACACCCCCGGCTACGACCCGGTCTCGGCCACCGGTCAGGTGGCCGGCGGGGCCAACCTGATCGCCTTCaccacc                                                                                                          >  1:885601/1‑150 (MQ=255)
                                       aggcaggcTCCAGCAACCTGGTCGATGTGTATCAGTACGCGCAAGCCGTGCGTGCCCAGGGCCTGGTGTTCATGGACACCCCCGGCTACGACCCGGTCTCGGCCACCGGTCAGGTGGCCGGCGGGGCCAACCTGATCGCCTTCACCACCg                                                                                                         <  1:149239/150‑1 (MQ=255)
                                                   gCAACCTGGTCGATGTGTATCAGTACGCGCAAGCCGTGCGTGCCCAGGGCCTGGTGTTCATGGACACCCCCGGCTACGACCCGGTCTCGGCCACCTGACCGGTGGcc                                                                                                                                        <  2:788424/107‑1 (MQ=255)
                                                   gCAACCTGGTCGATGTGTATCAGTACGCGCAAGCCGTGCGTGCCCAGGGCCTGGTGTTCATGGACACCCCCGGCTACGACCCGGTCTCGGCCACCTGACCGGTGGcc                                                                                                                                        >  1:788424/1‑107 (MQ=255)
                                                       ccTGGTCGATGTGTATCAGTACGCGCAAGCCGTGCGTGCCCAGGGCCTGGTGTTCATGGACACCCCCGGCTACGACCCGGTCTCGGCCACCGGTCAGGTGGc                                                                                                                                         >  1:1159655/1‑102 (MQ=255)
                                                       ccTGGTCGATGTGTATCAGTACGCGCAAGCCGTGCGTGCCCAGGGCCTGGTGTTCATGGACACCCCCGGCTACGACCCGGTCTCGGCCACCGGTCAGGTGGc                                                                                                                                         <  2:1159655/102‑1 (MQ=255)
                                                              tatctttATCATTACTCTCACGCCGTGCGTGCCCAGGGCCTGGTCTTCATCGCCACCCCCGTCTACGACCCGGTCTCGGCCACCGGGCAGGTGGCCGGTCGGGCCATCCGGATCGCCTTCACCACCGGCCGCGGCTCGGCCTACGGCtgt                                                                                  <  2:595832/144‑1 (MQ=255)
                                                                      tCAGTACGCGCAAGCCGTGCGTGCCCAGGGCCTGGTGTTCATGGACACCCCCGGCTACGACCCGGTCTCGGCCACCGGTCAGGTGGCCGGCGGGGCCAACCTGATCGCCTTCACCACCGGcc                                                                                                      <  2:1249199/122‑1 (MQ=255)
                                                                      tCAGTACGCGCAAGCCGTGCGTGCCCAGGGCCTGGTGTTCATGGACACCCCCGGCTACGACCCGGTCTCGGCCACCGGTCAGGTGGCCGGCGGGGCCAACCTGATCGCCTTCACCACCGGcc                                                                                                      >  1:1249199/1‑122 (MQ=255)
                                                                               gcAAGCCGTGCGTGCCCAGGGCCTGGTGTTCATGGACACCCCCGGCTACGACCCGGTCTCGGCCACCGGTCAGGTg                                                                                                                                           >  2:780206/1‑76 (MQ=255)
                                                                               gcAAGCCGTGCGTGCCCAGGGCCTGGTGTTCATGGACACCCCCGGCTACGACCCGGTCTCGGCCACCGGTCAGGTg                                                                                                                                           <  1:780206/76‑1 (MQ=255)
                                                                                 aaGCCGTGCGTGCCCAGGGCCTGGTGTTCATGGACACCCCCGGCTACGACCCGGTCTCGGCCACCGGTCAGGTgg                                                                                                                                          <  1:1183735/75‑1 (MQ=255)
                                                                                 aaGCCGTGCGTGCCCAGGGCCTGGTGTTCATGGACACCCCCGGCTACGACCCGGTCTCGGCCACCGGTCAGGTgg                                                                                                                                          >  2:1183735/1‑75 (MQ=255)
                                                                                 aaGCCGTGCGTGCCCAGGGCCTGGTGTTCATGGACACCCCCGGCTACGACCCGGTCTCGGCCACCGGTCAGGTGGACGGCGGGGCCAACCTGATCGCCTTCACCACCGGCCGCGGCTCGGCCTACGGCTGTGCCCCGGCGCCATCGATCa                                                               >  1:229421/1‑150 (MQ=255)
                                                                                  aGCCGTGCGTGCCCAGGGCCTGGTGTTCATGGGCAACCCCGGCTACGACTCGTTCTCGGCCACCGGTCAGTTGGCCGGCGGGAGCAACCTGATCGCCTTCACCACCGGCCGCGGCTCGGCCTACGGCTGTGCCCCGGCGCCATCGATCaa                                                              >  1:1033245/1‑150 (MQ=255)
                                                                                               cAGGGCCTGGTGTTCATGGACACCCCCGGCTACGACCCGGTCTCGGCCACCGGTCAGGTGGCCGGCGGGGCCAACCTGATCGCCTTCACCACCGGCCGCGGCTCGGCCTACGGCTGTGCCCCGGCGCCATCGATCAAGCTGGcaaccaac                                                 >  1:1419884/1‑150 (MQ=255)
                                                                                                     cTGGTGTTCATGGACACCCCCGGCTACGACCCGGTCTCGGCCACCGGTCAGGTGGCCGGCGGGGCCAACCTGATCGCCTTCACCACCGGCCGCGGCTCGGCCTACGGCTGTGCCCCGGCGCCATCGATCAAGCTGGCAACCAACAACCgg                                           <  1:167317/150‑1 (MQ=255)
                                                                                                     cTGGTGTTCATGGACACCCCCGGCTACGACCCGGTCTCGGCCACCGGTCAGGTGGCCGGCGGGGCCAACCTGATCGCCTTCACCACCGGCCGCGGCTCGGCCTACGGCTGTGCCCCGGCGCCATCGATCAAGCTGACAACCAACAACCgg                                           <  2:6312/150‑1 (MQ=255)
                                                                                                      tGGTGTTCATGGACACCCCCGGCTACGACCCGGTCTCGGCCACCGGTCAGGTGGCCGGCGGGGCCAACCTGATCGCCTTCACCACCGGCCGCGGCTCGGCCTACGGCTGTGCCCCGGCGCCATCGATCAAGCTGGCAACCAACAACCggg                                          >  1:354035/1‑150 (MQ=255)
                                                                                                       ggTGTTCATGGACACCCCCGGCTACGACCCGGTCTCGGCCACCGGTCAGGTGGCCGGCGGGGCCAACCTGATCGCCTTCACCACCGGCCGCGGCTCGGCCTACGGCTGTGCCCCGGCGCCATCGATCAAGCTGGCAACCAACAACCGGgt                                         <  2:1033245/150‑1 (MQ=255)
                                                                                                              aTGGACACCCCCGGCTACGACCCGGTCTCGGCCACCGGTCAGGTGGCCGGCGGGGCCAACCTGATCGCCTTCACCACCGGCCGCGGCTCGGCCTACGGCTGTGCCCCGGCGCCATCGATCAAGCTGGCAACCAACAACCGGGTGTTCGAg                                  >  1:1189211/1‑150 (MQ=255)
                                                                                                               tGGACACCCCCGGCTACGACCCGGTCTCGGCCACCGGTCAGGTGGCCGGCGGGGCCAACCTGATCGCCTTCACCACCGGCCGCGGCTCGGCCTACGGCTGTGCCCCGGCGCCATCGATCAAGCTGGCAACCAACAACCGGGTGTTCGAGc                                 >  1:340035/1‑150 (MQ=255)
                                                                                                                       cccGGCTACGACCCGGTCTCGGCCACCGGTCAGGTGGCCGGCGGGGCCAACCTGATCGCCTTCACCACCGGCCGCGGCTCGGCCTACGGCTGTGCCCCGGCGCCATCGATCAAGCTGGCAAACAGCAACCGGGTGTTCGAGCACCAGCaa                         >  2:1039682/1‑150 (MQ=255)
                                                                                                                                gACCCGGTCTCGGCCACCGGTCAGGTGGCCGGCGGGGCCAACCTGATCGCCTTCACCACCGGCCGCGGCTCGGCCTACGGCTGTGCCCCGGCGCCATCGATCAAGCTGGCAACCAACAACCggg                                          >  1:785380/1‑124 (MQ=255)
                                                                                                                                gACCCGGTCTCGGCCACCGGTCAGGTGGCCGGCGGGGCCAACCTGATCGCCTTCACCACCGGCCGCGGCTCGGCCTACGGCTGTGCCCCGGCGCCATCGATCAAGCTGGCAACCAACAACCggg                                          <  2:785380/124‑1 (MQ=255)
                                                                                                                                gACCCGGTCTCGGCCACCGGTCAGGTGGCCGGCGGGGCCAACCTGATCGCCTTCACCACCGGCCGCGGCTCGGCCTACGGCTGTGCCCCGGCGCCATCGATCAAGCTGGCAACCAACAACCGGGTGTTCGAGCACCAGCAAGAGGACATg                >  2:697972/1‑150 (MQ=255)
                                                                                                                                gACCCGGTCTCGGCCACCGGTCAGGTGGCCGGCGGGGCCAACCTGATCGCCTTCACCACCGGCCGCGGCTCGGCCTACGGCTGTGCCCCGGCGCCATCGATCAAGCTGGCAACCAACAACCGGGTGTTCGAGCACCAGCAAGAGGACATg                <  2:885601/150‑1 (MQ=255)
                                                                                                                                  cccGGTCTCGGCCACCGGTCAGGTGGCCGGCGGGGCCAACCTGATCGCCTTCACCACCGGCCGCGGCTCGGCCTACGGCTGTGCCCCGGCGCCATCGATCAAGCTGGCAACCAACAACCGGGTGTTCGAGCACCAGCAAGAGGACATGGa              >  1:545885/1‑150 (MQ=255)
                                                                                                                                    cGGTCTCGGCCACCGGTCAGGTGGCCGGCGGGGCCAACCTGATCGCCTTCACCACCGGCCGCGGCTCGGCCTACGGCTGTGCCCCGGCGCCATCGATCAAGCTGGCAACCAACAACCGGGTGTTCGAGCACCAGCAAGAGGACATGGACg            >  2:996609/1‑150 (MQ=255)
                                                                                                                                      gTCTCGGCCACCGGTCAGGTGGCCGGCGGGGCCAACCTGATCGCCTTCACCACCGGCCGCGGCTCGGCCTACGGCTGTGCCCCGGCGCCATCGATCAAGCTGGCAACCAACAACCGGGTGTTCGAGCACCAGCAAGAGGACATGGACGTc          <  2:1189211/150‑1 (MQ=255)
                                                                                                                                       tctcGGCCACCGGTCAGGTGGCCGGCGGGGCCAACCTGATCGCCTTCACCACCGGCCGCGGCTCGGCCTACGGCTGTGCCCCGGCGCCATCGATCAAGCTGGCAACCAACAACCGGGTGTTCGAGCACCAGCAAGAGGACATGGACGTCa         >  1:1250443/1‑150 (MQ=255)
                                                                                                                                           ggCCACCGGTCAGGTGGCCGGCGGGGCCAACCTGATCGCCTTCACCACCGGCCGCGGCTCGGCCTACGGCTGTGCCCCGGCGCCATCGATCAAGCTGGCAACCAACAACCGGGTGTTCGAGCACCAGCAAGAGGACATGGACGTCAATTg     >  2:1078929/1‑150 (MQ=255)
                                                                                                                                              cACCGGTCAGGTGGCCGGCGGGGCCAACCTGATCGCCTTCACCACCGGCCGCGGCTCGGCCTACGGCTGTGCCCCGGCGCCATCGATCAAGCTGGCAACCAACAACCGGGTGTTCGAGCACCAGCAAGAGGACATGGACGTCAATTgcgg  <  2:1250443/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
TGACCACTATCCTCGAGAAATCCCTGGGCGCAGTGGCCAAGGCAGGCTCCAGCAACCTGGTCGATGTGTATCAGTACGCGCAAGCCGTGCGTGCCCAGGGCCTGGTGTTCATGGACACCCCCGGCTACGACCCGGTCTCGGCCACCGGTCAGGTGGCCGGCGGGGCCAACCTGATCGCCTTCACCACCGGCCGCGGCTCGGCCTACGGCTGTGCCCCGGCGCCATCGATCAAGCTGGCAACCAACAACCGGGTGTTCGAGCACCAGCAAGAGGACATGGACGTCAATTGCGG  >  NC_002947/3236045‑3236336

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