Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,561,078 G→A 9.3% A264T (GCC→ACC)  tdcG‑III → L‑serine dehydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,561,0780GA 9.3% 112.4 / 3.3 43A264T (GCC→ACC) tdcG‑IIIL‑serine dehydratase
Reads supporting (aligned to +/- strand):  ref base G (16/23);  new base A (2/2);  total (18/25)
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.42e-01

TCTGGGCGGCCATGCGTGAGTGCGTGGCGAACGGCCTGCGCAACGAGGGCATCCTGCCGGGCGGCTTGCACGTCAAGCGCCGCGCCGCGCGGTTGCATCGCAGCCTGCAGGAGCTGGGCAAGCCCAACGTGATCGGCTCGACCCTTGGTGCCATGGAGTGGGTCAACCTGTTCGCCTTGGCGGTCAATGAAGAAAATGCTGCCGGTGGGCGTATGGTTACCGCGCCCACCAACGGGGCAGCCGGCATTATCCCGGCGGTGTTGCACTACTACATGAAGTTCAACCC  >  NC_002947/3560929‑3561214
                                                                                                                                                     |                                                                                                                                        
tCTGGGCGGCCATGCGTGAGTGCGTGGCGAACGGCCTGCGCAACGAGGGCATCCTGCCGGGCGGCTTGCACGTCAAGCGCCGCGCCGCGCGGTTGCATCGCAGCCTGCAGGAGCTGGGCAAGCCCAACGTGATCGGCTCGACCCTTGGTg                                                                                                                                          <  2:926099/150‑1 (MQ=255)
    ggcggcCATGCGTGAGTGCGTGGCGAACGGCCTGCGCAACGAGGGCATCCTGCCGGGCGGCTTGCACGTCAAGCGCCGCGCCGCGCGGTTGCATCGCAGCCTGCAGGAGCTGGGCAAGCCCAACGTGATCGGCTCGACCCTTGGTGCCAt                                                                                                                                      <  1:768352/150‑1 (MQ=255)
    ggcggcCATGCGTGAGTGCGTGGCGAACGGCCTGCGCAACGAGGGCATCCTGCCGGGCGGCTTGCACGTCAAGCGCCGCGCCGCGCGGTTGCATCGCAGCCTGCAGGAGCTGGGCAAGCCCAACGTGATCGGCTCGACCCTTGGTGCCAt                                                                                                                                      <  2:722774/150‑1 (MQ=255)
          cATGCGTGAGTGCGTGGCGAACGGCCTGCGCAACGAGGGCATCCTGCCGGGCGGCTTGCACGTCAAGCGCCGCGCCGCGCGGTTGCATCGCAGCCTGCAGGAGCTGGGCAAGCCCAACGTGATCGGCTCGACCCTTGGTGCCATGGAGTg                                                                                                                                >  2:978826/1‑150 (MQ=255)
           aTGCGTGAGTGCGTGGCGAACGGCCTGCGCAACGAGGGCATCCTGCCGGGCGGCTTGCACGTCAAGCGCCGCGCCGCGCGGTTGCATCGCAGCCTGCAGGAGCTGGGCAAGCCCAACGTGATCGGCTCGACCCTTGGTGCCATGGAGTgg                                                                                                                               >  1:718692/1‑150 (MQ=255)
           aTGCGTGAGTGCGTGGCGAACGGCCTGCGCAACGAGGGCATCCTGCCGGGCGGCTTGCACGTCAAGCGCCGCGCCGCGCGGTTGCATCGCAGCCTGCAGGAGCTGGGCAAGCCCAACGTGATCGGCTCGACCCTTGGTGCCATGGAGTgg                                                                                                                               >  1:1078604/1‑150 (MQ=255)
                tGAGTGCGTGGCGAACGGCCTGCGCAACGAGGGCATCCTGCCGGGAGGCTTCCTCATCAATCGCCTCTCCGCGCGGTTGCATCGCAGCCTGCAGGAGCTGGGCAAGCCCAACGTGATCGGCTCGACCCTTGGTGCCATGGAGTGGGTCaa                                                                                                                          <  1:407270/150‑1 (MQ=38)
                     gCGTGGCGAACGGCCTGCGCAACGAGGGCATCCTGCCGGGCGGCTTGCACGTCAAGCGCCGCGCCGCGCGGTTGCATCGCAGCCTGCAGGAGCTGGGCAAGCCCAACGTGATCGGCTCGACCCTTGGTGCCATGGAGTGGGTCAACCTGt                                                                                                                     <  2:96809/150‑1 (MQ=255)
                                     gcgcAACGAGGGCATCCTGCCGGGCGGCTTGCACGTCAAGCGCCGCGCCGCGCGGTTGCATCGCAGCCTGCAGGAGCTGGGCAAGCCCAACGTGATCGGCTCGACCCTTGGTGCCATGGAGTGGGTCAACCTGTTCGCCTTGGCGGTCaa                                                                                                     <  1:343126/150‑1 (MQ=255)
                                        cAACGAGGCTACCGTGCTGGGCGGCTTGCACGAAAAGCGCCTCGCCGCGCGGTTGCATCGCAGCGTGCAGGAGCGGGGCAAGCCCAACGTGATCGGCTCGACCCTTCGTGCTATGTACTGGGTCGACCTGTTCGCCTTGGCGGTCAATga                                                                                                  <  2:405028/150‑1 (MQ=21)
                                                gcagtcTGCCGGGCGGGGTGCACGGCAAGCGCGGGGCCGCGCGGTTGCATCGCAGCCTGCAGGAGCTGGGCAAGCCCAACGTGATCGGCTCGACCCTTGGTGCCATGGAGTGGGTCAACCTGTTCGCCTTGGAGGTTAATGAAGAAAAtg                                                                                          <  2:401322/145‑1 (MQ=21)
                                                  aTCCTGCCGGGCGGCTTGCACGTCAAGCGCCGCGCCGCGCGGTTGCATCGCAGCCTGCAGGAGCTGGGCAAGCCCAACGTGATCGGCTCGACCCTTGGTGCCATGGAGTGGGTCAACCTGTTCGCCTTGGCGGTCAATGAAGAAAAtgct                                                                                        >  1:960036/1‑150 (MQ=255)
                                                      tGCCGGGCGGCTTGCACGTCAAGCGCAGCGCCGCGCGGTTGCATCGCAGCCTGCAGGAGCTGGGCAAGCCCAACGTGATCGGCTCGACCCTTGGTGCCATGGCGTGGGTCAACCTGTCCGCCTTGGCGGTCAAGGAAGAAAATGCTGCCg                                                                                    >  2:559678/1‑150 (MQ=255)
                                                       gCCGGGCGGCTTGCACGTCAAGCGCCGCGCCGCGCGGTTGCATCGCAGCCTGCAGGGGCTGGGCAAGCCCAACGTGATCGGCTCGACCCTTGGTGCCATGGAGTGGGTCAACCTGTTCGCCTTGGCGGTCAATGAAGAAAATGCTGCCgg                                                                                   >  1:1037155/1‑150 (MQ=255)
                                                          gggCGGCTTGCACGTCAAGCGCCGCGCCGCGCGGTTGCATCGCAGCCTGCAGGAGCTGGGCAAGCCCAACGTGATCGGCTCGACCCTTGGTGCCATGGAGTGGGTCAACCTGTTCGCCTTGGCGGTCAATGAAGAAAATGCTGCCGGTgg                                                                                <  2:939761/150‑1 (MQ=255)
                                                          gggCGGCTTGCACGTCAAGCGCCGCGCCGCGCGGTTGCATCGCAGCCTGCAGGAGCTGGGCAAGCCCAACGTGATCGGCTCGACCCTTGGTGCCATGGAGTGGGTCAACCTGTTCGCCTTGGCGGTCAATGAAGAAAATGCTGCCGGTgg                                                                                <  1:17738/150‑1 (MQ=255)
                                                          gggCGGCTTGCACGTCAAGCGCCGCGCCGCGCGGTTGCATCGCAGCCTGCAGGAGCTGGGCAAGCCCAACGTGATCGGCTCGACCCTTGGTGCCATGGAGTGGGTCAACCTGTTCGCCTTGGCGGTCAATGAAGAAAATGCTGCCGGTgg                                                                                <  2:313691/150‑1 (MQ=255)
                                                               gcTTGCACGTCAAGCGCCGCGCCGCGCGGTTGCATCGCAGCCTGCAGGAGCTGGGCAAGCCCAACGTGATCGGCTCGACCCTTGGTACCATGGAGTGGGTCAACCTGTTCGCCTTGGCGGTCAATGAAGAAAATGCTGCCGGTGGGCGTa                                                                           >  1:571363/1‑150 (MQ=255)
                                                               gcTTGCACGTCAAGCGCCGCGCCGCGCGGTTGCATCGCAGCCTGCAGGAGCTGGGCAAGCCCAACGTGATCGGCTCGACCCTTGGTACCATGGAGTGGGTCAACCTGTTCGCCTTGGCGGTCAATGAAGAAAATGCTGCCGGTGGGCGTa                                                                           <  2:571363/150‑1 (MQ=255)
                                                               gcTTGCACGTCAAGCGCCGCGCCGCGCGGTTGCATCGCAGCCTGCAGGAACTGGGCAAGCCCAACGTGATCGGCTCGACCCTTGGTACCATGGAGTGGGTCAACCTGTTCGCCTTGGCGGTCAATGAAGAAAATGCTGCCGGTGGGCGTa                                                                           >  1:571227/1‑150 (MQ=255)
                                                               gcTTGCACGTCAAGCGCCGCGCCGCGCGGTTGCATCGCAGCCTGCAGGAACTGGGCAAGCCCAACGTGATCGGCTCGACCCTTGGTACCATGGAGTGGGTCAACCTGTTCGCCTTGGCGGTCAATGAAGAAAATGCTGCCGGTGGGCGTa                                                                           <  2:571227/150‑1 (MQ=255)
                                                                cTTGCACGTCAAGCGCCGCGCCGCGCGGTTGCATCGCAGCCTGCAGGAGCTGGGCAAGCCCAACGTGATCGGCTCGACCCTTGGTGCCATGGAGTGGGTCAACCTGTTCGCCTTGGCGGTCAATGAAGAAAATGCTGCCGGTGGGCGTAt                                                                          <  1:141904/150‑1 (MQ=255)
                                                                                     cgcgcgGTTGCATCGCAGCCTGCAGGAGCTGGGCAAGCCCAACGTGATCGGCTCGACCCTTGGTGCCATGGAGTGGGTCAACCTGTTCGCCTTGGCGGTCAATGAAGAAAATGCTGCCGGTGGGCGTATGGTTACCGCGCCCACCAACgg                                                     <  2:1078604/150‑1 (MQ=255)
                                                                                      gcgcgGTTGCATCGCAGCCTGCAGGAGCTGGGCAAGCCCAACGTGATCGGCTCGACCCTTGGTGCCATGGAGTGGGTCAACCTGTTCGCCTTGGCGGTCAATGAAGAAAATGCTGCCGGTGGGCGTATGGTTACCGCGCCCACCAACggg                                                    <  2:960036/150‑1 (MQ=255)
                                                                                               cATCGCAGACTGCAGGAGCTGGGCAAGCCCAACGTGATCGGCTCGACCCTTGGTGCCATGGAGTGGGTCAACCTGTTCGCCTTGGCGGTCAATga                                                                                                  >  1:194050/1‑95 (MQ=255)
                                                                                               cATCGCAGACTGCAGGAGCTGGGCAAGCCCAACGTGATCGGCTCGACCCTTGGTGCCATGGAGTGGGTCAACCTGTTCGCCTTGGCGGTCAATga                                                                                                  <  2:194050/95‑1 (MQ=255)
                                                                                                aTCGCAGCCTGCAGGAGCTGGGCAAGCCCAACGTGATCGGCTCGACCCTTGGTGCCATGGAGTGGGTCAACCTGTTCGCCTTGGCGGTCAATGAAGAAAATGCTGCCGGTGGGCGTATGGTTACCGCGCCCACCAACGGGGCAGCCGGCa                                          >  2:447946/1‑150 (MQ=255)
                                                                                                aTCGCAGCCTGCAGGAGCTGGGCAAGCCCAACGTGATCGGCTCGACCCTTGGTGCCATGGAGTGGGTCAACCTGTTCGCCTTGGCGGTCAATGAAGAAAATGCTGCCGGTGGGCGTATGGTTACCGCGCCCACCAACGGGGCAGCCGGCa                                          >  1:1020494/1‑150 (MQ=255)
                                                                                                         tGCAGGAGCTGGGCAAGCCCAACGTGATCGGCTCGACCCTTGGTGCCATGGAGTGGGTCAACCTGTTCGCCTTGGCGGTCAATGAAGAAAATGCTGCCGGTGGGCGTATGGTTACCGCGCCCACCAACGGGGCAGCCGGCATTATCCcgg                                 <  1:103648/150‑1 (MQ=255)
                                                                                                          gCAGGAGCTGGGCAAGCCCAACGTGATCGGCTCGACCCTTGGTGCCATGGAGTGGGTCAACCTGTTCGCCTTGGCGGTCAATGAAGAAAATGCTGCCGGTGGGCGTATGGTTACCGCGCCCACCAACGGGGCAGCCGGCATTATCCcggc                                >  1:1089179/1‑150 (MQ=255)
                                                                                                           cAGGAGCTGGGCAAGCCCAACGTGATCGGCTCGACCCTTGGTGCCATGGAGTGGGTCAACCTGTTCGCCTTGGCGGTCAATGAAGAAAATGCTGCCGGTGGGCGTATGGTTACCGCGCCCACCAACGGGGCAGCCGGCATTATCCcggcg                               <  1:163502/150‑1 (MQ=255)
                                                                                                           cAGGAGCTGGGCAAGCCCAACGTGATCGGCTCGACCCTTGGTGCCATGGAGTGGGTCAACCTGTTCGCCTTGGCGGTCAATGAAGAAAATGCTGCCGGTGGGCGTATGGTTACCGCGCCCACCAACGGGGCAGCCGGCATTATCCcggcg                               <  1:163408/150‑1 (MQ=255)
                                                                                                             ggAGCTGGGCAAGCCCAACGTGATCGGCTCGACCCTTGGTGCCATGGAGTGGGTCAACCTGTTCGCCTTGGCGGTCAATGAAGAAAATGCTGCCGGTGGGCGTATGGTTACCGCGCCCACCAACGGGGCAGCCGGCATTATCCCGGCGgt                             <  2:718692/150‑1 (MQ=255)
                                                                                                                   gggCAAGCCCAACGTGATCGGCTCGACCCTTGGTGCCATGGAGTGGGTCAACCTGTTCGCCTTGGCGGTCAATGAAGAAAATGCTGCCGGTGGGCGTATGGTTACCGCGCCCACCAACGGGGCAGCCGGCATTATCCCGGCGGTGTTGCa                       <  2:1037155/150‑1 (MQ=255)
                                                                                                                     gCAAGCCCAACGTGATCGGCTCGACCCTTGGTGCCATGGAGTGGGTCAACCTGTTCGCCTTGGCGGTCAATGAAGAAAATGCTGCCGGTGGGCATATGGTTACCGTGCCCACCAACGGGCCAGCCGGCAGAATCCCGGCGGGGTTGCact                     >  2:67520/1‑150 (MQ=21)
                                                                                                                      cAAGCCCAACGTGATCGGCTCGACCCTTGGTGCCATGGAGTGGGTCAACCTGTTCGCCTTGGCGGTCAATGAAGAAAATGCTGCCGGTGGGCGTATGGTTACCGCGCCCACCAACGGGGCAGCCGGCATTATCCCGGCGGTGTTGCacta                    >  1:510239/1‑150 (MQ=255)
                                                                                                                       aaGCCCAACGTGATCGGCTCGACCCTTGGTGCCATGGAGTGGGTCAACCTGTTCGCCTTGGCGGTCAATGAAGAAAATGCTGCCGGTGGGCGTATGGTTACCGCGCCCACCAACgg                                                     >  1:95281/1‑116 (MQ=255)
                                                                                                                       aaGCCCAACGTGATCGGCTCGACCCTTGGTGCCATGGAGTGGGTCAACCTGTTCGCCTTGGCGGTCAATGAAGAAAATGCTGCCGGTGGGCGTATGGTTACCGCGCCCACCAACgg                                                     <  2:95281/116‑1 (MQ=255)
                                                                                                                        aGCCCAACGTGATCGGCTCGACCCTTGGTGCCATGGAGTGGGTCAACCTGTTCGCCTTGGCGGTCAATGAAGAAAATGCTGCCGGTGGGCGTATGGTTACCGCGCCCACCAACGGGGCAGCCGGCATTATCCCGGCGGTGTTGCa                       <  1:913307/145‑1 (MQ=255)
                                                                                                                        aGCCCAACGTGATCGGCTCGACCCTTGGTGCCATGGAGTGGGTCAACCTGTTCGCCTTGGCGGTCAATGAAGAAAATGCTGCCGGTGGGCGTATGGTTACCGCGCCCACCAACGGGGCAGCCGGCATTATCCCGGCGGTGTTGCa                       >  2:913307/1‑145 (MQ=255)
                                                                                                                                     cGGCTCGACCCTTGGTGCCATGGAGTGGGTCAACCTGTTCGCCTTGGCGGTCAATGAAGAAAATGCTGCCGGTGGGCGTATGGTTACCGCGCCCACCAACGGGGCAGCCGGCATTATCCCGGCGGTGTTGCACTACTACATGAAGTTCaa     <  2:1020494/150‑1 (MQ=255)
                                                                                                                                        cTCGACCCTTGGTGCCATGGAGTGGGTCAACCTGTTCGCCTTGGCGGTCAATGAAGAAAATGCTGCCGGTGGGCGTATGGTTACCGCGCCCACCAACGGGGCAGCCGGCATTATCCCGGCGGTGTTGCACTACTACATGAAGTTCAAccc  >  1:376067/1‑150 (MQ=255)
                                                                                                                                         tCGACCCTTGGTGCCATGGAGTGGGTCAACCTGTTCGCCTTGGCGGTCAATGAAGAAAATGCTGCCGGTGGGCGTATGGTTACCGCGCCCACCAACGGGGCAGc                                               <  2:6874/104‑1 (MQ=255)
                                                                                                                                         tCGACCCTTGGTGCCATGGAGTGGGTCAACCTGTTCGCCTTGGCGGTCAATGAAGAAAATGCTGCCGGTGGGCGTATGGTTACCGCGCCCACCAACGGGGCAGc                                               >  1:6874/1‑104 (MQ=255)
                                                                                                                                                     |                                                                                                                                        
TCTGGGCGGCCATGCGTGAGTGCGTGGCGAACGGCCTGCGCAACGAGGGCATCCTGCCGGGCGGCTTGCACGTCAAGCGCCGCGCCGCGCGGTTGCATCGCAGCCTGCAGGAGCTGGGCAAGCCCAACGTGATCGGCTCGACCCTTGGTGCCATGGAGTGGGTCAACCTGTTCGCCTTGGCGGTCAATGAAGAAAATGCTGCCGGTGGGCGTATGGTTACCGCGCCCACCAACGGGGCAGCCGGCATTATCCCGGCGGTGTTGCACTACTACATGAAGTTCAACCC  >  NC_002947/3560929‑3561214

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