Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,560,699 C→T 9.3% S137S (AGC→AGT)  tdcG‑III → L‑serine dehydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,560,6990CT 9.3% 116.7 / 3.3 43S137S (AGC→AGT) tdcG‑IIIL‑serine dehydratase
Reads supporting (aligned to +/- strand):  ref base C (23/16);  new base T (2/2);  total (25/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.79e-01

CAGGTCCTGCAAGAGCACTGCCTGATGCTCGACGGCAGCCATCCGCTGGCCTTCGACTACGGCCGTGACATGCTGCTGCTCGACGAGAGCCTGCCTTATCACCCCAATGCCATGAGCCTGGAAGCCATGGACGGGCAGCGCAGCCTGTTGAGCCAGACCTATTTTTCGGTGGGCGGCGGCTTCATCGTCGAGCAAGCCGAGATCGATGCGCCGGCAAGCGAAATGAATGCCGTGATGCTGCCTTATGAGTTCGACAGTGGTGCGCAACTGCTGGTGTTGTGCAAGGCCC  >  NC_002947/3560556‑3560844
                                                                                                                                               |                                                                                                                                                 
cAGGTCCTGCAAGAGCACTGCCTGATGCTCGACGGCATCCATCCGCTTGCCTTCGACTACGGCCGTGCCATGCTGCTGCTCGACGAGAGCCTGCCTTATCACCCCAATGACATGAGCCTCGAAGCCATGGACGGGCAGCGCAGCCTGTTg                                                                                                                                             >  1:498932/1‑150 (MQ=255)
      cTGCAAGAGCACTGCGTGATGCTCGACGGCAGCCATCCGCTGGCCTTCGACTACGGCCGTGACATGCTGCTGCTCAACGAGAGCTTGCCGTAGCGCCCCAATGCCATGAACCTGGAAGCCATGGACGGGCAGCGCAGCCTGTTGAGCCAg                                                                                                                                       <  1:26700/150‑1 (MQ=255)
       tGCAAGAGCACTGCCTGATGCTCGACGGCAGCCATCCGCTGGCCTTCGACTACGGCCGTGACATGCTGCTGCTCGACGAGAGCCTGCCTTATCACCCCAATGCCATGAGCCTGGAAGCCATGGACGGGCAGCGCAGCCTGTTGAGCCAGa                                                                                                                                      >  1:179228/1‑150 (MQ=255)
            gagCACTGCCTGATGCTCGACGGCAGCCATCCGCTGGCCTTCGACTACGGCCGTGACATGCTGCTGCTCGACGAGAGCCTGCCTTATCACCCCAATGCCATGAGCCTGGAAGCCATGGACGGGCAGCGCAGCCTGTTGAGCCAGACCTAt                                                                                                                                 >  1:247134/1‑150 (MQ=255)
                    ccTGATGCTCGACGGCAGCCATCCGCTGGCCTTCGACTACGGCCGTGACATGCTGCTGCTCGCCGAGGGCCGGCCTTCTCACCCCAATGCCATGAGCCTGGAAGCCATGGACGGGCAGCGCAGCCTGTTGAGCCAGACCTATTTTTCGGt                                                                                                                         <  2:943803/150‑1 (MQ=255)
                            tCGACGGCAGCCATCCGCTGGCCTTCGACTACGGCCGTGACATGCTGCTGCTCGACGAGAGCCTGCCTTATCACCCCAATGCCATGAGCCTGGAAGCCATGGACGGGCAGCGCAGCCTGTTGAGCCAGACCTATTTTTCGGTGggcggcg                                                                                                                 <  1:398963/150‑1 (MQ=255)
                              gACGGCAGCCATCCGCTGGCCTTCGACTACGGCCGTGACATGCTGCTGCTCGACGAGAGCCTGCCTTATCACCCCAATGCCATGAGCCTGGAAGCCATGGACGGGCAGCGCAGCCTGTTGAGCCAGACCTATTTTTCGGTGggcggcggc                                                                                                               >  1:954168/1‑150 (MQ=255)
                                       cATCCGCTGGCCTTCGACTACGGCCGTGACATGCTGCTGCTCGACGAGAGCCTGCCTTATCACCCCAATGCCATGAGCCTGGAAGCCATGGACGGGCAGCGCAGCCTGTTGAGCCAGACCTATTTTTCGGTGGGCGGCGGCTTCAtcgtc                                                                                                      >  1:1074493/1‑150 (MQ=255)
                                       cATCCGCTGGCCTTCGACTACGGCCGTGACATGCTGCTGCTCGACGAGAGCCTGCCTTATCACCCCAATGCCATGAGCCTGGAAGCCATGGACGGGCAGCGCAGCCTGTTGAGCCAGACCTATTTTTCGGTGGGCGGCGGCTTCAtcgtc                                                                                                      >  1:1066468/1‑150 (MQ=255)
                                            gCTGGCCTTCGACTACGGCCGTGACATGCTGCTGCTCGACGAGAGCCTGCCTTATCACCCCAATGCCATGAGCCTGGAAGCCATGGACGGGCAGCGCAGCCTGTTGAGCCAGACCTATTTTTCGGTGGGCGGCGGCTTCATCGTCGAGCa                                                                                                 >  1:436236/1‑150 (MQ=255)
                                                  cTTCGACTACGGCCGTGACATGCTGCTGCTCGACGAGAGCCTGCCTTATCACCCCAATGCCATGAGCCTGGAAGCCATGGACGGGCAGCGCAGCCTGTTGAGCCAGACCTATTTTTCGGTGGGCGGCGGCTTCATCGTCGAGCAAGCCga                                                                                           >  2:35364/1‑150 (MQ=255)
                                                     cGACTACGGCCGTGACATGCTGCTGCTCGACGAGAGCCTGCCTTATCACCCCAATGCCATGAGCCTGGAAGCCATGGACGGGCAGCGCAGCCTGTTGAGCCAGACCTATTTTTCGGTGGGCGGCGGCTTCATCGTCGAGCAAGCCGAGAt                                                                                        <  2:1074493/150‑1 (MQ=255)
                                                     cGACTACGGCCGTGACATGCTGCTGCTCGACGAGAGCCTGCCTTATCACCCCAATGCCATGAGCCTGGAAGCCATGGACGGGCAGCGCAGCCTGTTGAGCCAGACCTATTTTTCGGTGGGCGGCGGCTTCATCGTCGAGCAAGCCGAGAt                                                                                        <  1:144530/150‑1 (MQ=255)
                                                      gACTACGGCCGTGACATGCTGCTGCTCGACGAGAGCCTGCCTTATCACCCCAATGCCATGAGCCTGGAAGCCATGGACGGGCAGCGCAGCCTGTTGAGCCAGACCTATTTTTTTTTTTTTTTCGGCTTCATCGTCGAGCAAGCCGAGATc                                                                                       >  2:805661/1‑150 (MQ=255)
                                                            ggCCGTGACATGCTGCTGCTCGACGAGAGCCTGCCTTATCACCCCAATGCCATGAGCCTGGAAGCCATGGACGGGCAGCGCAGCCTGTTGAGCCAGACCTATTTTTCGGTGGGCGGCGGCTTCATCGTCGAGCAAGCCGAGATCGATgcg                                                                                 >  1:320792/1‑150 (MQ=255)
                                                               cGTGACATGTTGCTGCTCGACGAGAGCCTGCCTTATCACCCCAATGCCATGAGCCTGGAAGCCATGGACGGGCAGCGCAGCCTGTTGAGCCAGACCTATTTTTCGGTGGGCGGCGGCTTCATCGTCGAGCAAGCCGAGATCGATGCGCCg                                                                              >  2:261495/1‑150 (MQ=255)
                                                                gTGACATGCTGCTGCTCGACGAGAGCCTGCCTTATCACCCCAATGCCATGAGCCTGGAAGCCATGGACGGGCAGCGCAGCCTGTTGAGc                                                                                                                                          <  1:233829/89‑1 (MQ=255)
                                                                gTGACATGCTGCTGCTCGACGAGAGCCTGCCTTATCACCCCAATGCCATGAGCCTGGAAGCCATGGACGGGCAGCGCAGCCTGTTGAGc                                                                                                                                          >  2:233829/1‑89 (MQ=255)
                                                                  gACATGCTGCTGCTCGACGAGAGCCTGCCTTATCACCCCAATGCCATGAGCCTGGAAGCCATGGACGGGCAGCGCAGCCTGTTGAGCCAGACCTATTTTTCGGTGGGCGGCGGCTTCATCGTCGAGCAAGCCGAGATCGATGCGCCGGCa                                                                           >  1:755007/1‑150 (MQ=255)
                                                                     aTGCTGCTGCTCGACGAGAACCTGCATTATCACCCGAATGCGTTGAGCCTGGGAGCCATGGACGGGCAGCGCAGCCTGTTGAGCCAGACCTATTTTTCGGTGGGCGGCGGCTTCATCGTCGAGCAAGCCGAGATCGATGCGCCGGCAAGc                                                                        >  1:933042/1‑150 (MQ=255)
                                                                                 gacgaGAGCCTGCCTTATCACCCCAATGCCATGAGCCTGGAAGCCATGGACGGGCAGCGCAGCCTGTTGAGCCAGACCTATTTTTCGGTGGGCGGCGGCTTCATCGTCGAGCAAGCCGAGATCGATGCGCCGGCAAGCGAAATGAATGcc                                                            <  1:261495/150‑1 (MQ=255)
                                                                                      gagCCTGCCTTATCACCCCAATGCCATGAGCCTGGAAGCCATGGACGGGCAGCGCAGCCTGTTGAGCCAGACCTATTTTTCGGTGGGCGGCGGCTTCATCGTCGAGCAAGCCGAGATCGATGCGCCGGCAAGCGAAATGAATGCCGTGAt                                                       <  1:1207694/150‑1 (MQ=255)
                                                                                          ctgcctTATCACCCCAATGCCATGAGCCTGGAAGCCATGGACGGGCAGCGCAGCCTGTTGAGCCAGACCTATTTTTCGGTGGGCGGCGGCTTCATCGTCGAGCAAGCCGAGATCGATGCGCCGGCAAGCGAAATGAATGCCGTGAtgctg                                                   >  1:9786/1‑150 (MQ=255)
                                                                                           tgcctTATCACCCCAATGCCATGAGCCTGGAAGCCATGGACGGGCAGCGCAGCCTGTTGAGCCAGACCTATTTTTCGGTGGGCGGCGGCTTCATCGTCGAGCAAGCCGAGATCGATGCGCCGGCAAGCGAAATGAATGCCGTGAtgctgc                                                  >  1:935690/1‑150 (MQ=255)
                                                                                            gcctTATCACCCCAATGCCATGAGCCTGGAAGCCATGGACGGGCAGCGCAGCCTGTTGAGCCAGACCTATTTTTCGGTGGGCGGCGGCTTCATCGTCGAGCAAGCCGAGATCGATGCGCCGGCAAGCGAAATGAATGCCGTGATGCTGcc                                                 >  1:577183/1‑150 (MQ=255)
                                                                                             cctTATCCACCCAATGCCATGGGCCTGGAAGCCCTGGACCTGCAGCGCAGCCTTGTGAGCCAGACCTATTTTTCGGTGGGCGGCGGCTTCATCGTCGAGCAAGCCGAGATCGATGCGCCGGCAAGCGAAATGAATGCCGTGATGCTGCCt                                                <  2:436236/150‑1 (MQ=255)
                                                                                                 aTCACCCCAATGCCATGAGCCTGGAAGCCATGGACGGGCAGCGCAGCCTGTTGAGCCAGACCTATTTTTCGGTGGGCGGCGGCTTCATCGTCGAGCAAGCCGAGATCGATGCGCCGGCAAGCGAAATGAATGCCGTGATGATGCCTTATg                                            >  2:557558/1‑150 (MQ=255)
                                                                                                   cACCCCAATGCCATGAGCCTGGAAGCCATGGACGGGCAGCGCAGCCTGTTGAGCCAGACCTATTTTTCGGTGGGCGGCGGCTTCATCGTCGAGCAAGCCGAGATCGATGCGCCGGCAAGCGAAATGAATGCCGTGATGCTGCCTTATGAg                                          <  2:577183/150‑1 (MQ=255)
                                                                                                       ccAATGCCATGAGCCTGGAAGCCATGGACGGGCAGCGCAGCCTGTTGAGCCAGACCTATTTTTCGGTGGGCGGCGGCTTCATCGTCGAGCAAGCCGAGATCGATGCGCCGGCAAGCGAAATGAATGCCGTGATGCTGCCTTATGAGTTCg                                      >  2:955218/1‑150 (MQ=255)
                                                                                                         aaTGCCATGAGCCTGGAAGCCATGGACGGGCAGCGCAGCCTGTTGAGCCAGACCTATTTTTCGGTGGGCGGCGGCTTCATCGTCGAGCAAGCCGAGATCGATGCGCCGGCAAGCGAAATGAATGCCGTGATGCTGCCTTATGAGTTCGAc                                    >  2:1113394/1‑150 (MQ=255)
                                                                                                             ccATGAGCCTGGAAGCCATGGACGGGCAGCGCAGCCTGTTGAGCCAGACCTATTTTTCGGTGGGCGGCGGCTTCATCGTCGAGCAAGCCGAGATCGATGCGCCGGCAAGCGAAATGAATGCCGTGATGCTGCCTTATGAGTTCGACAgtg                                <  1:1113394/150‑1 (MQ=255)
                                                                                                             ccATGAGCCTGGAAGCCATGGACGGGCAGCGCAGCCTGTTGAGCCAGACCTATTTTTCGGTGGGCGGCGGCTTCATCGTCGAGCAAGCCGAGATCGATGCGCCGGCAAGCGAAATGAATGCCGTGATGCTGCCTTATGAGTTCGACAgtg                                <  2:935690/150‑1 (MQ=255)
                                                                                                                 gAGCCTGGAAGCCATGGACGGGCAGCGCAGCCTGTTGAGCCAGACCTATTTTTCGGTGGGCGGCGGCTTCATCGTCGAGCAAGCCGAGATCGATGCGCCGGCAAGCGAAATGAATGCCGTGATGCTGCCTTATGAGTTCGACAGTGGTgc                            <  2:179228/150‑1 (MQ=255)
                                                                                                                 gAGCCTGGAAGCCATGGACGGGCAGCGCAGCCTGTTGAGCCAGACCTATTTTTCGGTGGGCGGCGGCTTCATCGTCGAGCAAGCCGAGATCGATGCGCCGGCAAGCGAAATGAATGCCGTGATGCTGCCTTATGAGTTCGACAGTGGTgc                            >  1:685751/1‑150 (MQ=255)
                                                                                                                    ccTGGAAGCCATGGACGGGCAGCGCAGCCTGTTGAGCCAGACCTATTTTTCGGTGGGCGGCGGCTTCATCGTCGAGCAAGCCGAGATCGATGCGCCGGCAAGCGAAATGAATGCCGTGATGCTGCCTTATGAGTTCGACAGTGGTGCGCa                         >  1:433451/1‑150 (MQ=255)
                                                                                                                       ggAAGCCATGGACGGGCAGCGCAGCCTGTTGAGCCAGACCTATTTTTCGGTGGGCGGCGGCTTCATCGTCGAGCAAGCCGAGATCGATGCGCCGGCAAGCGAAATGAATGCCGTGATGCTGCCTTATGAGTTCGACAGTGGTGCGCAAct                      <  2:685751/150‑1 (MQ=255)
                                                                                                                        gAAGCCATGGACGGGCAGCGCAGTCTGTTGAGCCAGACCTATTTTTCGGTGGGCGGCGGCTTCATCGTCGAGCAAGCCGAGATCGATGCGCCGGCAAGCGAAATGAATGCCGTGATGCTGCCt                                                >  2:905903/1‑123 (MQ=255)
                                                                                                                        gAAGCCATGGACGGGCAGCGCAGTCTGTTGAGCCAGACCTATTTTTCGGTGGGCGGCGGCTTCATCGTCGAGCAAGCCGAGATCGATGCGCCGGCAAGCGAAATGAATGCCGTGATGCTGCCt                                                >  2:905571/1‑123 (MQ=255)
                                                                                                                        gAAGCCATGGACGGGCAGCGCAGTCTGTTGAGCCAGACCTATTTTTCGGTGGGCGGCGGCTTCATCGTCGAGCAAGCCGAGATCGATGCGCCGGCAAGCGAAATGAATGCCGTGATGCTGCCt                                                <  1:905903/123‑1 (MQ=255)
                                                                                                                        gAAGCCATGGACGGGCAGCGCAGTCTGTTGAGCCAGACCTATTTTTCGGTGGGCGGCGGCTTCATCGTCGAGCAAGCCGAGATCGATGCGCCGGCAAGCGAAATGAATGCCGTGATGCTGCCt                                                <  1:905571/123‑1 (MQ=255)
                                                                                                                               tGGACGGGCAGCGCAGCCTGTTGAGCAAGACCTATTTTTCGGTGGGCGGCGGCTTCATCGTCGAGCAAGCCGAGATCGATGCGCCGGCAAGCGAAATGAATGCCGTGATGCTGCCTTATGAGTTCGACAGTGGTGCGCAACTGCTGgtgt              <  2:755007/150‑1 (MQ=255)
                                                                                                                                ggACGGGCAGCGCAGCCTGTTGAGCCAGACCTATTTTTCGGTGGGCGGCGGCTTCATCGTCGAGCAAGCCGAGATCGATGCGCCGGCAAGCGAAATGAATGCCGTGATGCTGCCTTATGAGTTCGACAGTGGTGCGCAACTGCTGGtgtt             >  1:210489/1‑150 (MQ=255)
                                                                                                                                ggACGGGCAGCGCAGCCTGTTGAGCCAGACCTATTTTTCGGTGGGCGGCGGCTTCATCGTCGAGCAAGCCGAGATCGATGCGCCGGCAAGCGAAATGAATGCCGTGATGCTGCCTTATGAGTTCGACAGTGGTGCGCAACTGCTGGtgtt             <  2:1066468/150‑1 (MQ=255)
                                                                                                                                      gcagcgcagcCTGTTGAGCCAGACCTATTTTTCGGTGGGCGGCGGCTTCATCGTCGAGCAAGCCGAGATCGATGCGCCGGCAAGCGAAATGAATGCCGTGATGCTGCCTTATGAGTTCGACAGTGGTGCGCAACTGCTGGTGTTGTGCaa       >  1:817258/1‑150 (MQ=255)
                                                                                                                                        agggcagcTTGTTGAGCCAGACCTATTTTTCGGTGGGCGGCGGCTTCATCGTCGAGCAAGCCGAGATCGATGCGCCGGCAAGCGAAATGAATGCCGTGATGCTGCCTTATGAGTTCGACAGTGGTGCGCAACTGCTGGTGTTGTGCAAgg     <  1:805661/147‑1 (MQ=255)
                                                                                                                                        agcgcagcCTGTTGAGCCAGACCTATTTTTCGGTGGGCGGCGGCTTCATCGTCGAGCAAGCCGAGATCGATGCGCCGGCAAGCGAAATGAATGCCGTGATGCTGCCTTATGAGTTCGACAGTGGTGCGCAACTGCTGGTGTTGTGCAAgg     <  2:433451/150‑1 (MQ=255)
                                                                                                                                           gcatcCGGTTGAGCCAGACCTATTTTTCGGTGGGCGGCGGCTTCATCGTCGAGCAAGCCGAGATCGATGCGCCGGCAAGCGAAATGAATGCCGTGATGCTGCCTTATGAGTTCGACAGTGGTGCGCAACTGCTGGTGTTGTGCAAGGccc  <  2:555163/150‑1 (MQ=255)
                                                                                                                                           gcagcCTGTTGAGCCAGACCTATTTTTCGGTGGGCGGCGGCTTCATCGTCGAGCAAGCCGAGATCGATGCGCCGGCAAGCGAAATGAATGCCGTGATGCTGCCTTATGAGTTCGACAGTGGTGCGCAACTGCTGGTGTTGTGCAAGGccc  >  1:326690/1‑150 (MQ=255)
                                                                                                                                               |                                                                                                                                                 
CAGGTCCTGCAAGAGCACTGCCTGATGCTCGACGGCAGCCATCCGCTGGCCTTCGACTACGGCCGTGACATGCTGCTGCTCGACGAGAGCCTGCCTTATCACCCCAATGCCATGAGCCTGGAAGCCATGGACGGGCAGCGCAGCCTGTTGAGCCAGACCTATTTTTCGGTGGGCGGCGGCTTCATCGTCGAGCAAGCCGAGATCGATGCGCCGGCAAGCGAAATGAATGCCGTGATGCTGCCTTATGAGTTCGACAGTGGTGCGCAACTGCTGGTGTTGTGCAAGGCCC  >  NC_002947/3560556‑3560844

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