Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,236,195 A→C 11.9% Q416P (CAG→CCG)  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,1950AC11.9% 82.4 / 5.6 35Q416P (CAG→CCG) PP_2834D‑galactarate dehydratase/altronate dehydratase
Reads supporting (aligned to +/- strand):  ref base A (16/14);  new base C (2/2);  total (18/17)
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.22e-01

GACCACTATCCTCGAGAAATCCCTGGGCGCAGTGGCCAAGGCAGGCTCCAGCAACCTGGTCGATGTGTATCAGTACGCGCAAGCCGTGCGTGCCCAGGGCCTGGTGTTCATGGACACCCCCGGCTACGACCCGGTCTCGGCCACCGGTCAGGTGGCCGGCGGGGCCAACCTGATCGCCTTCACCACCGGCCGCGGCTCGGCCTACGGCTGTGCCCCGGCGCCATCGATCAAGCTGGCAACCAACAACCGGGTGTTCGAGCACCAGCAAGAGGACATGGACGTCAATTGCGGCGGCATCG  >  NC_002947/3236046‑3236344
                                                                                                                                                     |                                                                                                                                                     
gACCACTATCCTCGAGAAATCCCTGGGCGCAGTGGCCAAGGCAGGCTCCAGCAACCTGGTCGATGTGTATCAGTACGCGCAAGCCGTGCGTGCCCAGGGCCTGGTGTTCATGGACACCCCCGGCTACGACCCGGTCTCGGCCACCGGTCa                                                                                                                                                       >  1:242978/1‑150 (MQ=255)
 aCCACTATCCTCGAGAAATCCCTGGGCGCAGTGGCCAAGGCAGGCTCCAGCAACCTGGTCGATGTGTATCAGTACGCGCAAGCCGTGCGTGCCCAGGGCCTGGTGTTCATGGACACCCCCGGCTACGACCCGGTCTCGGCCACCGGTCAg                                                                                                                                                      <  1:15022/150‑1 (MQ=255)
   cACTATCCTCGAGAAATCCCTGGGCGCAGTGGCCAAGGCAGGCTCCAGCAACCTGGTCGATGTGTATCAGTACGCGCAAGCCGTGCGTGCCCAGGGCCTGGTGTTCATGGACACCCCCGGCTACGACCCGGTCTCGGCCACCGGTCAGGt                                                                                                                                                    >  2:634108/1‑150 (MQ=255)
          cTCGAGAAATCCCTGGGCGCAGTGGCCAAGGCAGGCTCCAGCAACCTGGTCGATGTGTATCAGTACGCGCAAGCCGTGCGTGCCCAGGGCCTGGTGTTCATGGACACCCCCGGCTACGACCCGGTCTCGGCCACCGGTCAGGTGGCcggc                                                                                                                                             >  1:579342/1‑150 (MQ=255)
           tCGAGAAATCCCTGGGCGCAGTGGCCAAGGCAGGCTCCAGCAACCTGGTCGATGTGTATCAGTACGCGCAAGCCGTGCGTGCCCAGGGCCTGGTGTTCATGGACACCCCCGGCTACGACCCGGTCTCGGCCACCGGTCAGGTGGCcggcg                                                                                                                                            <  2:579342/150‑1 (MQ=255)
            cGAGACATCCCTGGGCGCAGTGGCCAAGGCAGGCTCCAGCAACCTGGTCGATGTGTATCAGTACGCGCAAGCCGTGCGTGCCCAGGGCCTGGTGTTCATGGACACCCCCGGCTACGACCCGGTCTACGCCACCGGTCTGGTAGCcggcgg                                                                                                                                           <  2:402592/150‑1 (MQ=255)
               gaAATCCCTGGGCGCAGTGGCCAAGGCAGGCTCCAGCAACCTGGTCGATGTGTATCAGTACGCGCAAGCCGTGCGTGCCCAGGGCCTGGTGTTCATGGACACCCCCGGCTACGACCCGGTCTCGGCCACCGGTCAGGTGGCCGGCGGGGc                                                                                                                                        <  2:209795/150‑1 (MQ=255)
                      cTGGGCGCAGTGGCCAAGGCAGGCTCCAGCAACCTGGTCGATGTGTATCAGTACGCGCAAGCCGTGCGTGCCCAGGGCCTGGTGTTCATGGACACCCCCGGCTACGACCCGGTCTCGGCCACCGGTCAGGTGGCCGGCGGGGCCAACCTg                                                                                                                                 >  1:213484/1‑150 (MQ=255)
                      cTGGGCGCAGTGGCCAAGGCAGGCTCCAGCAACCTGGTCGATGTGTATCAGTACGCGCAAGCCGTGCGTGCCCAGGGCCTGGTGTTCATGGACACCCCCGGCTACGACCCGGTCCCGGCCACCTGACCGGTGGCcga                                                                                                                                              >  2:12572/1‑136 (MQ=255)
                      cTGGGCGCAGTGGCCAAGGCAGGCTCCAGCAACCTGGTCGATGTGTATCAGTACGCGCAAGCCGTGCGTGCCCAGGGCCTGGTGTTCATGGACACCCCCGGCTACGACCCGGTCCCGGCCACCTGACCGGTGGCcga                                                                                                                                              <  1:12572/137‑2 (MQ=255)
                                                cAGCAACCTGGTCGATGTGTATCAGTACGCGCAAGCCGTGCGTGCCCAGGGCCTGGTGTTCATGGACACCCCCGGCTACGACCCGGTCTCGGCCACCGGTCAGGTGGCCGGCGGGGCCAACCTGATCGCCTTCACCACCg                                                                                                                 <  2:674667/140‑1 (MQ=255)
                                                cAGCAACCTGGTCGATGTGTATCAGTACGCGCAAGCCGTGCGTGCCCAGGGCCTGGTGTTCATGGACACCCCCGGCTACGACCCGGTCTCGGCCACCGGTCAGGTGGCCGGCGGGGCCAACCTGATCGCCTTCACCACCg                                                                                                                 >  1:674667/1‑140 (MQ=255)
                                                cAGCAACCTGGTCGATGTGTATCAGTACGCGCAAGCCGTGCGTGCCCAGGGCCTGGTGTTCATGGACACCCCCGGCTACGACCCGGTCTCGGCCACCGGTCAGGTGGCCGGCGGGGCCAACCTGATCGCCTTCACCACCGGCCGCGGCTc                                                                                                       >  1:442445/1‑150 (MQ=255)
                                                cAGCAACCTGGTCGATGTGTATCAGTACGCGCAAGCCGTGCGTGCCCAGGGCCTGGTGTTCATGGACACCCCCGGCTACGACCCGGTCTCGGCCACCGGTCAGGTGGCCGGCGGGGCCAACCTGATCGCCTTCACCACCGGCCGCGGCTc                                                                                                       <  2:213484/150‑1 (MQ=255)
                                                      ccTGGTCGATGTGTATCAGTACGCGCAAGCCGTGCGTGCCCAGGGCCTGGTGTTCATGGACACCCCCGGCTACGACCCGGTCTCGGCCACCGGTCAGGTGGCCGGCGGGGCCAACCTGATCGCCTTCACCACCGGCCGCGGCTCGGCCTa                                                                                                 >  1:655334/1‑150 (MQ=255)
                                                                   tATCAGTACGCGCAAGCCGTGCGTGCCCAGGGCCTGGTGTTCATGGACACCCCCGGCTACGACCCGGTCTCGGCCACCGGTCAGGTGGCCGGCGGGGCCAACCTGATCGCCTTCACCACCGGCCGCGGCTCGGCCTACGGCTGTGCCCCg                                                                                    >  2:229836/1‑150 (MQ=255)
                                                                     tCAGTACGCGCAAGCCGTGCGTGCCCAGGGCCTGGTGTTCATGGACACCCCCGGCTACGACCCGGTCTCGGCCACCGGTCAGGTGGCCGGCGGGGCAAACCCGATCGCCTTCACCACCGGCCGCGGCTCGGCCTACGGCTGTGCCCCGgc                                                                                  >  1:188977/1‑150 (MQ=255)
                                                                        gTACGCGCAAGCCGTGCGTGCCCAGGGCCTGGTGTTCATGGACACCCCCGGCTACGACCCGGTCTCGGCCACCAGTCAGGTGGCCGGCGGGGCCAACCTGATCGCCTTCACCACCGGCCGCGGCTCGGCCTACGGCTGTGCCCCGGCGcc                                                                               >  2:231553/1‑150 (MQ=255)
                                                                             cgcAAGCCGTGCGTGCCCAGGGCCTGGTGTTCATGGACACCCCCGGCTACGACCCGGTCTCGGCCACCGGTCAGGTGGCCGGCGGGGCCAACCTGATCGCCTTCACCACCGGCCGCGGCTCGGCCTACGGCTGTGCCCCGGCGCCATCGa                                                                          <  2:242978/150‑1 (MQ=255)
                                                                                    cgtgcgtgCCCAGGGCCTGGTGTTCATGGACACCCCCGGCTACGACCCGGTCTCGGCCACCGGTCAGGTGGCCGGCGGGGCCAACCTGATCGCCTTCACCACCGGCCGCGGCTCGGCCTACGGCTGTGCCCCGGCGCCATCGATCAAGCt                                                                   <  1:634108/150‑1 (MQ=255)
                                                                                      tgcgtgcCCAGGGCCTGGTGTTCATGGACACCCCCGGCTACGACCCGGTCTCGGCCACCGGTCAGGTGGCCGGCGGGGCCAACCTGATCGCCTTCACCACCGGCCGCGGCTCGGCCTACGGCTGTGCCCCGGCGCCATCGATCAAGCTgg                                                                 <  2:655334/150‑1 (MQ=255)
                                                                                                gggCCTGGTGTTCATGGACACCCCCGGCTACGACCCGGTCTCGGCCACCGGTCAGGTGGCCGGCGGGGCCAACCTGATCGCCTTCACCACCGGCCGCGGCTCGGCCTACGGCTGTGCCCCGGCGCCATCGATCAAGCTGGCAACcaacaa                                                       >  2:2797/1‑150 (MQ=255)
                                                                                                                   acCCCCGGCTACGACCCGGTCTCGGCCACCGGTCAGGTGGCCGGCGGGGCCAACCTGATCGCCTTCACCACCGGCCGCGGCTCGGCCTACGGCTGTGCCCCGGCGCCATCGATCAAGCTGGCAACCAACAACCgg                                                   >  1:631398/1‑135 (MQ=255)
                                                                                                                   acCCCCGGCTACGACCCGGTCTCGGCCACCGGTCAGGTGGCCGGCGGGGCCAACCTGATCGCCTTCACCACCGGCCGCGGCTCGGCCTACGGCTGTGCCCCGGCGCCATCGATCAAGCTGGCAACCAACAACCgg                                                   <  2:631398/135‑1 (MQ=255)
                                                                                                                     ccccGGCTACGACCCGGTCTCGGCCACCGGTCAGGTGGCCGGCGGGGCCAACCTGATCGCCTTCACCACCGGCCGCGGCTCGGCCTACGGCTGTGCCCCGGCGCCATCGATCaa                                                                      <  1:553956/114‑1 (MQ=255)
                                                                                                                     ccccGGCTACGACCCGGTCTCGGCCACCGGTCAGGTGGCCGGCGGGGCCAACCTGATCGCCTTCACCACCGGCCGCGGCTCGGCCTACGGCTGTGCCCCGGCGCCATCGATCaa                                                                      >  2:553956/1‑114 (MQ=255)
                                                                                                                       ccGGCTACGACCCGGTCTCGGCCACCGGTCAGGTGGCCGGCGGGGCCAACCTGATCGCCTTCACCACCGGCCGCGGCTCGGCCTACGGCTGTGCCCCGGCGCCATCGATCAAGCTGGCAACCAACAACCGGGTGTTCGAGCACCAGCAag                                >  2:376220/1‑150 (MQ=255)
                                                                                                                          gCTACGACCCGGTCTCGGCCACCGGTCAGGTGGCCGGCGGGGCCAACCTGATCGCCTTCACCACCGGCCGCGGCTCGGCCTACGGCTGTGCCCCGGCGCCATCGATCAAGCTGGCAACCAACAACCGGGTGTTCGAGCACCAGCAAGAgg                             <  1:2797/150‑1 (MQ=255)
                                                                                                                           cTACGACCCGGTCTCGGCCACCAGTCAGGTGGCCGGCGGGGCCAACCTGATCGCCTTCACCACCGGCCGCGGCTCGGCCTACGGCTGTGCCCCGGCGCCATCGATCAAGCTGGCAACCAACAACCGGGTGTTCGAGCACCAGCAAGAGGa                            <  1:231553/150‑1 (MQ=255)
                                                                                                                               gACCCGGTCTCGGCCACCGGTCAGGTGGCCGGCGGGGCCAACCTGATCGCCTTCACCACCGGCCGCGGCTCGGCCTACGGCTGTGCCCCGGCGCCATCGATCAAGCTGGCAACCAACAACCGGGTGTTCGAGCACCAGCaa                                 <  1:737314/141‑1 (MQ=255)
                                                                                                                               gACCCGGTCTCGGCCACCGGTCAGGTGGCCGGCGGGGCCAACCTGATCGCCTTCACCACCGGCCGCGGCTCGGCCTACGGCTGTGCCCCGGCGCCATCGATCAAGCTGGCAACCAACAACCGGGTGTTCGAGCACCAGCaa                                 >  2:737314/1‑141 (MQ=255)
                                                                                                                                     ccgccGGCCACCTGACCGGTGGCCGGCGGGGCCAACCTGATCGCCTTCACCACCGGCCGCGGCTCGGCCTACGGCTGTGCCCCGGCGCCATCGATCAAGCTGGCAACCAACAACCGGGTGTTCGAGCACCAGCAaga                               <  2:273158/133‑1 (MQ=255)
                                                                                                                                     ccgccGGCCACCTGACCGGTGGCCGGCGGGGCCAACCTGATCGCCTTCACCACCGGCCGCGGCTCGGCCTACGGCTGTGCCCCGGCGCCATCGATCAAGCTGGCAACCAACAACCGGGTGTTCGAGCACCAGCAaga                               >  1:273158/5‑137 (MQ=255)
                                                                                                                                          ggCCACCGGTCAGGTGGCCGGCGGGGCCAACCTGATCGCCTTCACCACCGGCCGCGGCTCGGCCTACGGCTGTGCCCCGGCGCCATCGATCAAGCTGGCAACCAACAACCGGGTGTTCGAGCACCAGCAAGAGGACATGGACGTCAATTg             >  2:556218/1‑150 (MQ=255)
                                                                                                                                             cACCGGTCAGGTGGCCGGCGGGGCCAACCTGATCGCCTTCACCACCGGCCGCGGCTCGGCCTACGGCTGTGCCCCGGCGCCATCGATCAAGCTGGCAACCAACAACCGGGTGTTCGAGCACCAGCAAGAGGACATGGACGTCAATTgcgg          >  2:588519/1‑150 (MQ=255)
                                                                                                                                                cGGTCAGGTGGCCGGCGGGGCCAACCTGATCGCCTTCACCACCGGCCGCGGCTCGGCCTACGGCTGTGCCCCGGCGCCATCGATCAAGCTGGCAACCAACAACCGGGTGTTCGAGCACCAGCAAGAGGACATGGACGTCAATTGcggcgg       <  1:229836/150‑1 (MQ=255)
                                                                                                                                                     aGGTGGCCGGCGGGGCCAACCTGATCGCCTTCACCACCGGCCGCGGCTCGGCCTACGGCTGTGCCCCGGCGCCATCGATCAAGCTGGCAACCAACAACCGGGTGTTCGAGCACCAGCAAGAGGACATGGACGTCAATTGCGGCGGCATCg  >  1:653533/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GACCACTATCCTCGAGAAATCCCTGGGCGCAGTGGCCAAGGCAGGCTCCAGCAACCTGGTCGATGTGTATCAGTACGCGCAAGCCGTGCGTGCCCAGGGCCTGGTGTTCATGGACACCCCCGGCTACGACCCGGTCTCGGCCACCGGTCAGGTGGCCGGCGGGGCCAACCTGATCGCCTTCACCACCGGCCGCGGCTCGGCCTACGGCTGTGCCCCGGCGCCATCGATCAAGCTGGCAACCAACAACCGGGTGTTCGAGCACCAGCAAGAGGACATGGACGTCAATTGCGGCGGCATCG  >  NC_002947/3236046‑3236344

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