Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,236,195 A→C 16.7% 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,1950AC16.7% 76.0 / 11.4 36Q416P (CAG→CCG) PP_2834D‑galactarate dehydratase/altronate dehydratase
Reads supporting (aligned to +/- strand):  ref base A (16/14);  new base C (3/3);  total (19/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 = 9.57e-01

TCGAGAAATCCCTGGGCGCAGTGGCCAAGGCAGGCTCCAGCAACCTGGTCGATGTGTATCAGTACGCGCAAGCCGTGCGTGCCCAGGGCCTGGTGTTCATGGACACCCCCGGCTACGACCCGGTCTCGGCCACCGGTCAGGTGGCCGGCGGGGCCAACCTGATCGCCTTCACCACCGGCCGCGGCTCGGCCTACGGCTGTGCCCCGGCGCCATCGATCAAGCTGGCAACCAACAACCGGGTGTTCGAGCACCAGCAAGAGGACATGGACGTCAATTGCGGC  >  NC_002947/3236057‑3236337
                                                                                                                                          |                                                                                                                                              
tCGAGAAATCCCTGGGCGCAGTGGCCAAGGCAGGCTCCAGCAACCTGGTCGATGTGTATCAGTACGCGCAAGCCGTGCGTGCCCAGGGCCTGGTGTTCATGGACACCCCCGGCTACGACCCGGTCTCGGCCACCGGTCAGGTGGCcggcg                                                                                                                                     >  1:944936/1‑150 (MQ=255)
tCGAGAAATCCCTGGGCGCAGTGGCCAAGGCAGGCTCCAGCAACCTGGTCGATGTGTATCAGTACGCGCAAGCCGTGCGTGCCCAGGGCCTGGTGTTCATGGACACCCCCGGCTACGACCCGGTCTCGGCCACCGGCCAGGTGGCcggcg                                                                                                                                     >  2:204470/1‑150 (MQ=255)
 cGAGAAATCCCTGGGCGCAGTGGCCAAGGCAGGCTCCAGCAACCTGGTCGATGTGTATCAGTACGCGCAAGCCGTGCGTGCCCAGGGCCTGGTGTTCATGGACACCCCCGGCTACGACCCGGTCTCGGCCACCTGACCGGTGGCcgagac                                                                                                                                    <  2:889876/150‑5 (MQ=255)
 cGAGAAATCCCTGGGCGCAGTGGCCAAGGCAGGCTCCAGCAACCTGGTCGATGTGTATCAGTACGCGCAAGCCGTGCGTGCCCAGGGCCTGGTGTTCATGGACACCCCCGGCTACGACCCGGTCTCGGCCACCGGTCAGGTGGCcggcgg                                                                                                                                    <  2:361753/150‑1 (MQ=255)
         cccTGGGCGCAGTGGCTAAGGCCGGCTCAAGCAACCTGGTAGATGTGTATCAGTACGCGCAACCCGTGCGTGCACAGGGCCTGGTGTTCATAGAAACCCACGATTACTCCCCGGTCTCGGCCACCGGTCAGGTCGCCGGCGGGGCCAAcc                                                                                                                            <  2:472031/150‑1 (MQ=255)
                 gcAGTGGCCAAGGCAGGCTCCAGCAACCTGGTCGATGTGTATCAGTACGCGCAAGCCGTGCGTGCCCAGGGCCTGGTGTTCATGGACACCCCCGGCTACGACCCGGTCTCGGCCACCGGTCAGGTGGCCGGCGGGGCCAACCTGATCGcc                                                                                                                    <  1:1165019/150‑1 (MQ=255)
                       gCCACGGCAGGCTCCAGCAACCTGGTCGATGTGTATCAGTACGCGCAAGCCGTGCGTGCCCAGGGCCTGGTGTTCATGGACACCCCCGGCTACGACCCGGTCTCGGCCACCGGTCAGGTGGCcggcgg                                                                                                                                    <  2:801107/128‑1 (MQ=255)
                       gCCAAGGCAGGCTCCAGCAACCTGGTCGATGTGTATCAGTACGCGCAAGCCGTGCGTGCCCAGGGCCTGGTGTTCATGGACACCCCCGGCTACGACCCGGTCTCGGCCACCGGTCAGGTGGCcggcgg                                                                                                                                    >  1:801107/1‑128 (MQ=255)
                                  cTCCAGCAACCTGGTCGATGTGTATCAGTACGCGCAAGCCGTGCGTGCCCAGGGCCTGGTGTTCATGGACACCCCCGGCTACGACCCGGTCTCGGCCACCGGTCAGGTGGCCGGCGGGGCCAACCTGATCGCCTTCACCACCGGCCGCgg                                                                                                   <  1:1030957/150‑1 (MQ=255)
                                         aaCCTGGTCGATGTGTATCAGTACGCGCAAGCCGTGCGTGCCCAGGGCCTGGTGTTCATGGACACCCCCGGCTACGACCCGGTCCCGGCCACCTGACCGGTGGCcgaga                                                                                                                                     >  2:760311/1‑106 (MQ=255)
                                         aaCCTGGTCGATGTGTATCAGTACGCGCAAGCCGTGCGTGCCCAGGGCCTGGTGTTCATGGACACCCCCGGCTACGACCCGGTCCCGGCCACCTGACCGGTGGCcgaga                                                                                                                                     <  1:760311/109‑4 (MQ=255)
                                           ccTGGTCGATGTGTATCAGTACGCGCAAGCCGTGCGTGCCCAGGGCCTGGTGTTCATGGACACCCCCGGCTACGACCCGGTCTCGGCCACCGGTCAGGTGGCCGGCGGGGCCAACCTGATCGCCTTCACCACCGGCCGCGGCTTGGCCTa                                                                                          >  1:963268/1‑150 (MQ=255)
                                                        tATCAGTACGCGCAAGCCGTGCGTGCCCAGGGCCTGGTGTTCATGGACACCCCCGGCTACGACCCGGTCTCGGCCACCGGTCAGGTGGCCGGCGGGGCCAACCTGATCGCCTTcac                                                                                                               >  1:843248/1‑116 (MQ=255)
                                                        tATCAGTACGCGCAAGCCGTGCGTGCCCAGGGCCTGGTGTTCATGGACACCCCCGGCTACGACCCGGTCTCGGCCACCGGTCAGGTGGCCGGCGGGGCCAACCTGATCGCCTTcac                                                                                                               <  2:843248/116‑1 (MQ=255)
                                                              tACGCGCAAGCCGTGCGTGCCCAGGGCCTGGTGTTCATGGACACCCCCGGCTACGACCCGGTCTCGGCCACCGGTCAGGTGGCCGGCGGGGCCAACCTGATCGCCTTCACCACCGGCCGCGGCTCGGCCTACGGCTGTGCTCCGGCGCCa                                                                       >  1:227754/1‑150 (MQ=255)
                                                                   gcAAGCCGTGCGTGCCCAGGGCCTGGTGTTCATGGACACCCCCGGCTACGACCCGGTCTCGGCCACCGGTCAGGTGGCCGGCGGGGCCAACCTGATc                                                                                                                       >  1:480220/1‑97 (MQ=255)
                                                                   gcAAGCCGTGCGTGCCCAGGGCCTGGTGTTCATGGACACCCCCGGCTACGACCCGGTCTCGGCCACCGGTCAGGTGGCCGGCGGGGCCAACCTGATc                                                                                                                       <  2:480220/97‑1 (MQ=255)
                                                                         cgtgcgtgCCCAGGGCCTGGTGTTCATGGACACCCCCGGCTACGACCCGGTCTCGGCCACCGGTCAGGTGGCCGGCGGGGCCAACCTGATCGCCTTCACCACCGGCCGCGGCTTGGCCTACGGCTGTGCCCCGGCGCCATCGATCAAGCt                                                            <  2:963268/150‑1 (MQ=255)
                                                                                  ccAGCGCCTGGGGTTCATGGACACCCCCGGCTACGACCCGGACTCGGCCACCGGTCAGGTGGCCGGCGGCGCCAACCTGATCGCCTTCACCGCCGGCCGCGGCTCGGCCTACGGCTGTGCCCCGGCGCCATCGATCAAGCTGGcaaccaa                                                   <  2:1032131/150‑1 (MQ=255)
                                                                                      ggCCTGGTGTTCATGGACACCCCCGGCTACGACCCGGTCTCGGCCACCGGTCAGGTGGCCGGCGGGGCCAACCTGATCGCCTTCACCACCGGCCGCGGCTCGGCCTACGGCTGTGCCCCGGCGCCATCGATCAAGCTGGCAACCaacaac                                               <  2:944936/150‑1 (MQ=255)
                                                                                          tGGTGTTCATGGACACCCCCGGCTACGACCCGGTCTCGGCCACCGGTCAGGTGGCCGGCGGGGCCAACCTGATCGCCTTCACCACCGGCCGCGGCTCGGCCTACGGCTGTGCCCCGGCGCCATCGATCAAGCTGGCAACCAACAACCggg                                           >  2:256032/1‑150 (MQ=255)
                                                                                                 cATGGACACCCCCGGCTACGACCCGGTCTCGGCCACCGGTCAGGTGGCCGGCGGGGCCAACCTGATCGCCTTCACCACCGGCCGCGGCTCGGCCTACGGCTGTGCCCCGGCGCCATCGATCAAGCTGGCAACCAACAACCGGGTGTTCGa                                    <  1:256032/150‑1 (MQ=255)
                                                                                                      acacCCCCGGCTACGACCCGGTCTCGGCCACCGGTCAGGTGGCCGGCGGGGCCAACCTGATCGCCTTCACCACCGGCCGCGGCTCGGCCTACGGCTGTGCCCCGGCGCCATCGATCAAGCTGGCAACCAACAACCgg                                            <  2:1107133/137‑1 (MQ=255)
                                                                                                      acacCCCCGGCTACGACCCGGTCTCGGCCACCGGTCAGGTGGCCGGCGGGGCCAACCTGATCGCCTTCACCACCGGCCGCGGCTCGGCCTACGGCTGTGCCCCGGCGCCATCGATCAAGCTGGCAACCAACAACCgg                                            >  1:1107133/1‑137 (MQ=255)
                                                                                                          ccccGGCTACGACCCGGTCTCGGCCACCGGTCAGGTGGCCGGCGGGGCCAACCTGATCGCCTTCACCACCGGCCGCGGCTCGGCCTACGGCTGTGCCCCGGCGCCATCGATCAAGCTGGCAACCAACAACCGGGTGTTCGAGCACCAGCa                           >  2:836798/1‑150 (MQ=255)
                                                                                                            ccGGCTACGACCCGGTCTCGGCCACCGGTCAGGTGGCCGGCGGGGCCAACCTGATCGCCTTCACCACCGGCCGCGGCTCGGCCTACGGCTGTGCCCCGGCGCCATCGATCAAGCTGGCAACCAACAACCGGGTGTTCGAGCACCAGCAag                         >  2:1088928/1‑150 (MQ=255)
                                                                                                                 tACGACCCGGTCTCGGCCACCGGTCAGGTGGCCGGCGGGGCCAACCTGATCGCCTTCACCACCGGCCGCGGCTCGGCCTACGGCTGTGCCCCGGCGCCATCGATCAAGCTGGCAACCAACAACCGGGTGTTCGAGCACCAGCAAGAGGAc                    <  1:836798/150‑1 (MQ=255)
                                                                                                                       ccGGTCTCGGCCACCGGTCAGGTGGCCGGCGGGGCCAACCTGATCGCCTTCACCACCGGCCGCGGCTCGGCCTACGGCTGTGCCCCGGCGCCATCGATCAAGCTGGCAACCAACAACCGGGTGTTCGAGCACCAGCAAGAGGACATGGAc              >  2:274489/1‑150 (MQ=255)
                                                                                                                          gTCTCGGCCACCGGTCAGGTGGCCGGAAGGGCCAACCTGATCGCCTTCACCACCGGCCGCGGCTCGGTCTACGGCTGTGCCCCGGCGCCATCGATCAAGCTGGCAACCAACAACCGGGTGTTCGAGCACCAGCAAGAGGACATGGACGTc           >  1:939294/1‑150 (MQ=255)
                                                                                                                          ccgccGGCCACCTGACCGGTGGCCGGCGGGGCCAACCTGATCGCCTTCACCACCGGCCGCGGCTCGGCCTACGGCTGTGCCCCGGCGCCATCGa                                                                   <  2:282875/90‑1 (MQ=38)
                                                                                                                          ccgccGGCCACCTGACCGGTGGCCGGCGGGGCCAACCTGATCGCCTTCACCACCGGCCGCGGCTCGGCCTACGGCTGTGCCCCGGCGCCATCGa                                                                   >  1:282875/5‑94 (MQ=38)
                                                                                                                           tctcGGCCACCGGTCAGGTGGCCGGCGGGGCCAACCTGATCGCCTTCACCACCGGCCGCGGCTCGGCCTACGGCTGTGCCCCGGCGCCATCGATCAAGCTGGCAACCAACAACCGAGCGTGCGAGTGCCAGCAAGAGGCCATGGACGTCa          >  2:1056815/1‑150 (MQ=255)
                                                                                                                              cGGCCACCTGACCGGTGGCCGGCGGGGCCAACCTGATCGCCTTCACCACCGGCCGCGGCTCGGCCTACGGCTGTGCCCCGGCGCCATCGATCAAGCTGGCAACCAACAACCGGGTGTTCGAGCACCAGCAAGAGGACATGGACGTCAAtt       >  2:297236/1‑150 (MQ=255)
                                                                                                                               ggCCACCGGTCAGGTGGCCGGCGGGGCCAACCTGATCGCCTTCACCACCGGCCGCGGCTCGGCCTACGGCTGTGCCCCGGCGCCATCGATCAAGCTGGCAACCAACAACCGGGTGTTCGAGCACCAGCAAGAGGACATGGACGTCAATTg      >  2:364622/1‑150 (MQ=255)
                                                                                                                                   aCCGGTCAGGTGGCCGGCGGGGCCAACCTGATCGCCTTCACCACCGGCCGCGGCTCGGCCTACGGCTGTGCCCCGGCGCCATCGATCAAGCTGGCAACCAACAACCGGGTGTTCGAGCACCAGCAAGAGGACATGGACGTCAATTgcggc  <  1:1088928/150‑1 (MQ=255)
                                                                                                                                   aCCGGTCAGGTGGCCGGCGGGGCCAACCTGATCGCCTTCACCACCGGCCGCGGCTCGGCCTACGGCTGTGCCCCGGCGCCATCGATCAAGCTGACAACCAACAACCGGGTGTTCGAGCACCAGCAAGAGGACATGGACGTCAATTgcggc  >  2:119501/1‑150 (MQ=255)
                                                                                                                                          |                                                                                                                                              
TCGAGAAATCCCTGGGCGCAGTGGCCAAGGCAGGCTCCAGCAACCTGGTCGATGTGTATCAGTACGCGCAAGCCGTGCGTGCCCAGGGCCTGGTGTTCATGGACACCCCCGGCTACGACCCGGTCTCGGCCACCGGTCAGGTGGCCGGCGGGGCCAACCTGATCGCCTTCACCACCGGCCGCGGCTCGGCCTACGGCTGTGCCCCGGCGCCATCGATCAAGCTGGCAACCAACAACCGGGTGTTCGAGCACCAGCAAGAGGACATGGACGTCAATTGCGGC  >  NC_002947/3236057‑3236337

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: ‑