Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 356,720 G→T 7.8% G22V (GGC→GTC)  tdcG‑I → L‑serine dehydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947356,7200GT 7.8% 161.1 / 2.9 51G22V (GGC→GTC) tdcG‑IL‑serine dehydratase
Reads supporting (aligned to +/- strand):  ref base G (23/24);  new base T (2/2);  total (25/26)
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.48e-01

GCTGTTTGAGCCTAGCAGCTACCCCGCCCTGCGCATGGCGCGCAGAGACAAGCCCAGCAAGAGGTGATTCGACAATGGCCATCAGCGTGTTCGACCTGTTCAAGATCGGTGTCGGGCCCTCCAGCTCCCACACCGTCGGCCCCATGCGTGCTGGCGCCCTGTTCGTTCAGGGCCTGCGCGAACGCGGTGAGCTGGAGCGTGTAAAACGGATCGAAGTGCGCCTGTATGGCTCGCTGTCAGCCACCGGTATCGGCCATGGCACCGACAACGCCACCATCATGGGCC  >  NC_002947/356582‑356866
                                                                                                                                          |                                                                                                                                                  
gCTGTTTGAGCCTAGCAGCTACCCCGCCCTGCGCATGGCGCGCAGAGACAAGCCCAGCAAGAGGTGATTCGACAATGGCCATCAGCGTGTTCGACCTGTTCAAGATCGGTGTCGGGCCCTCCAGCTCCCACACCGTCGGCCCCATGCGt                                                                                                                                          >  4:139855/1‑149 (MQ=255)
     ttGAGCCTAGCAGCTACCCCGGCCTGCGCATGGCGCTCTGAGACAAGCCCAGCAAGATGTGATTCGACAATGGCCATCAGCGTGTTCGACCTGTTCAATATCGTTGTCGGGCCCTCCAGCTCCCACAGCGTCGGCCCCATGCGTGCTg                                                                                                                                      <  3:275924/148‑1 (MQ=255)
           cTAGCAGCTACCCCGCCCTGCGCATGGCGCGCAGAGACAAGCCCAGCAAGAGGTGATTCGACAATGGCCATCAGCGTGTTCGACCTGTTCAAGATCGGTGTCGGGCCCTCCAGCTCCCACACCGTCGGCCCCATGCGTGCTGGCGCCCt                                                                                                                               <  4:50406/149‑1 (MQ=255)
            tAGCAGCTACCCCGCCCTGCGCATGGCGCGCAGAGACAAGCCCAGCAAGAGGTGATTCGACAATGGCCATCAGCGTGTTCGACCTGTTCAAGATCGGTGTCGGGCCCTCCAGCTCCCACACCGTCGGCCCCATGCGTGCTGGCGCCCTg                                                                                                                              >  6:239327/1‑149 (MQ=255)
                 gcTACCCCGCCCTGCGCATGGCGCGCAGAGACAAGCCCAGCAAGAGGTGATTCGACAATGGCCATCAGCGTGTTCGACCTGTTCAAGATCGGTGTCGGGCCCTCCAGCTCCCACACCGTCGGCCCCATGCGTGCTGGCGCCCTgttcgt                                                                                                                         >  7:95927/1‑149 (MQ=255)
                     ccccGCCCTGCGCATGGCGCGCAGAGACAAGCCCAGCAAGAGGTGATTCGACAATGGCCATCAGCGTGTTCGACCTGTTCAAGATCGGTGTCGGGCCCTCCAGCTCCCACACCGTCGGCCCCATGCGTGCTGGCGCCCTGTTCGTTCAg                                                                                                                     >  5:189739/1‑149 (MQ=255)
                      cccGCCCTGCGCATGGCGCGCAGAGACAAGCCCAGCAAGAGGTGATTCGACAATGGCCATCAGCGTGTTCGACCTGTTCAAGATCGGTGTCGGGCCCTCCAGCTCCCACACCGTCGGCCCCATGCGTGCTGGCGCCCTGTTCGTTCAgg                                                                                                                    >  1:120882/1‑149 (MQ=255)
                      cccGCCCTGCGCATGGCGCGCAGAGACAAGCCCAGCAAGAGGTGATTCGACAATGGCCATCAGCGTGTTCGACCTGTTCAAGATCGGTGTCGGGCCCTCCAGCTCCCACACCGTCGGCCCCATGCGTGCTGGCGCCCTGTTCGTTCAgg                                                                                                                    <  3:139855/149‑1 (MQ=255)
                        cGCCCTGCGCATGGCGCGCAGAGACAAGCCCAGCAAGAGGTGATTCGACAATGGCCATCAGCGTGTTCGACCTGTTCAAGATCGGTGTCGGGCCCTCCAGCTCCCACACCGTCGGCCCCATGCGTGCTGGCGCCCTGTTCGTTCAGGGc                                                                                                                  >  2:59503/1‑149 (MQ=255)
                             tGCGCATGGCGCGCAGAGACAAGCCCAGCAAGAGGTGATTCGACAATGGCCATCAGCGTGTTCGACCTGTTCAAGATCGGTGTCGGGCCCTCCAGCTCCCACACCGTCGGCCCCATGCGTGCTGGCGCCCTGTTCGTTCAGGGCCTgcg                                                                                                             >  7:112598/1‑149 (MQ=255)
                                 cATGGCGCGCAGAGACAAGCCCAGCAAGAGGTGATTCGACAATGGCCATCAGCGTGTTCGACCTGTTCAAGATCGGTGTCGGGCCCTCCAGCTCCCACACCGTCGGCCCCATGCGTGCTGGCGCCCTGTTCGTTCAGGGCCTGCGCGaa                                                                                                         <  6:189739/149‑1 (MQ=255)
                                     gcgcgcAGAGACAAGCCCAGCAAGAGGTGATTCGACAATGGCCATCAGCGGGTTCGACCTGGTCAAGATCGGTGTCGGGCCCTCCAGCGCCGCCACGGTCGGGCCCCTGCGTGCTGGCGCGCTGTGCGTGAAGGGGCTGGCGGAAcgcg                                                                                                     >  4:262715/1‑149 (MQ=255)
                                        cgcAGAGACAAGCCCAGCAAGAGGTGATTCGACAATGGCCATCAGCGTGTTCGACCTGTTCAAGATCGGTGTCGGGCCCTCCAGCTCCCACACCGTCGGCCCCATGCGTGCTGGCGCCCTGTTCGTTCAGGGCCTGCGCGAACGCGGTg                                                                                                  <  1:4405/149‑1 (MQ=255)
                                                 aaGCCCAGCAAGAGGTGATTCGACAATGGCCATCAGCGTGTTCGACCTGTTCAAGATCGGTGTCGGGCCCTCCAGCTCCCACACCGTCGGCCCCATGCGTGCTGGCGCCCTGTTCGTTCAGGGCCTGCGCGAACGCGGTGAGCTGGAGc                                                                                         >  5:51873/1‑149 (MQ=255)
                                                   gCCCAGCAAGAGGTGATTCGACAATGGCCATCAGCGTGTTCGACCTGTTCAAGATCGGTGTCGGGCCCTCCAGCTCCCACACCGTCGTCCCCATGCGTGCTGGCGCCCTGTTCGTTCAGGGCCTGCGCGAACGCGGTGAGCTGGAGCgt                                                                                       >  8:238352/1‑149 (MQ=255)
                                                   gCCCAGCAAGAGGTGATTCGACAATGGCCATCAGCGTGTTCGACCTGTTCAAGATCGGTGTCGGGCCCTCCAGCTCCCACACCGTCGTCCCCATGCGTGATGTCGCCCTGTTCGTTCAGGGACTGCGCTAACGCGGTGAGCTGGAGCgt                                                                                       >  8:233952/1‑149 (MQ=255)
                                                   gCCCAGCAAGAGCTGATTCCCCCAGGGCCATCAGCGTGTTCGCCCTGTTCAAGATTGGTGTCGGGCCCTCCAGCTCCCACACCGTCGGCCGCCTGCGTGCTGGCGCCCTGCTCGTTCAGGGCCTGCGCGATCGCGGTGAGCTGGAGCgt                                                                                       <  4:235128/149‑1 (MQ=255)
                                                      cAGCAAGAGGTGATTCGACAATGGCCATCAGCGTGTTCGACCTGTTCAAGATCGGTGTCGGGCCCTCCAGCTCCCACACCGTCGGCCCCATGCGTGCTGGCGCCCTGTTCGTTCAGGGCCTGCGCGAACGCGGTGAGCTGGAGCgtgt                                                                                     >  5:145440/1‑148 (MQ=255)
                                                      cAGCAAGAGGTGATTCGACAATGGCCATCAGCGTGTTCGACCTGTTCAAGATCGGTGTCGGGCCCTCCAGCTCCCACACCGTCGGCCCCATGCGTGCTGGCGCCCTGTTCGTTCAGGGCCTGCGCGAACGCGGTGAGCTGGAGCgtgt                                                                                     >  7:215077/1‑148 (MQ=255)
                                                               gTGATTCGACAATGGCCATCAGCGTGTTCGACCTGTTCAAGATCGGTGTCGGGCCCTCCAGCTCCCACACCGTCGGCCCCATGCGTGCTGGCGCCCTGTTCGTTCAGGGCCTGCGCGAACGCGGTGAGCTGGAGCGTGTAAAACGGATc                                                                           >  4:104501/1‑149 (MQ=255)
                                                                  aTTCGACAATGGCCATCAGCGTGTTCGACCTGTTCAAGATCGGTGTCGGGCCCTCCAGCTCCCACACCGTCGGCCCCATGCGTGCTGGCGCCCTGTTCGTTCAGGGCCTGCGCGAACGCGGTGAGCTGGAGCGTGTAAAACGGATCGaa                                                                        <  5:14713/149‑1 (MQ=255)
                                                                     cGACAATGGCCATCAGCGTGTTCGACCTGTTCAAGATCGGTGTCGGGCCCTCCAGCTCCCACACCGTCGGCCCCATGCGTGCTGGCGCCCTGTTCGTTCAGGGCCTGCGCGAACGCGGTGAGCTGGAGCGTGTAAAACGGATCGAAGt                                                                      <  1:59503/148‑1 (MQ=255)
                                                                      gACAATGGCCATCAGCGTGTTCGACCTGTTCAAGATCGGTGTCGGGCCCTCCAGCTCCCACACCGTCGGCCCCATGCGTGCTGGCGCCCTGTTCGTTCAGGGCCTGCGCGAACGCGGTGAGCTGGAGCGTGTAAAACGGATCGAAGTgc                                                                    >  2:266490/1‑149 (MQ=255)
                                                                       aCAATGGCCATCAGCGTGTTCGACCTGTTCAAGATCGGTGTCGGGCCCTCCAGCTCCCACACCGTCGGCCCCATGCGTGCTGGCGCCCTGTTCGTTCAGGGCCTGCGCGAACGCGGTGAGCTGGAGCGTGTAAAACGGATCGAAGTgcg                                                                   >  2:288819/1‑149 (MQ=255)
                                                                        cAATGGCCATCAGCGTGTTCGACCTGTTCAAGATCGGTGTCGGGCCCTCCAGCTCCCACACCGTCGGCCCCATGCGTGCTGGCGCCCTGTTCGTTCAGGGCCTGCGCGAACGCGGTGAGCTGGAGCGTGTAAAACGGATCGAAGTgcgc                                                                  <  2:169162/149‑1 (MQ=255)
                                                                                    gCGTGTTCGACCTGTTCAAGATCGGTGTCGGGCCCTCCAGCTCCCACACCGTCGGCCCCATGCGTGCTGGCGCCCTGTTCGTTCAGGGCCTGCGCGAACGCGGTGAGCTGGAGCGTGTAAAACGGATCGAAGTGCGCCTGTATGGCTCg                                                      <  2:179957/149‑1 (MQ=255)
                                                                                      gtgtTCGACCTGTTCAAGATCGGTGTCGGGCCCTCCAGCTCCCACACCGTCGGCCCCATGCGTGCTGGCGCCCTGTTCGTTCAGGGCCTGCGCGAACGCGGTGAGCTGGAGCGTGTAAAACGGATCGAAGTGCGCCTGTATGGCTCGCt                                                    <  8:95927/149‑1 (MQ=255)
                                                                                       tgtTCGACCTGTTCAAGATCGGTGTCGGGCCCTCCAGCTCCCACACCGTCGGCCCCATGCGTGCTGGCGCCCTGTTCGTTCAGGGCCTGCGCGAACGCGGTGAGCTGGAGCGTGTAAAACGGATCGAAGTGCGCCTGTATGGCTCGCTg                                                   >  8:170930/1‑149 (MQ=255)
                                                                                        gtTCGACCTGTTCAAGATCGGTGTCGGGCCCTCCAGCTCCCACACCGTCGGCCCCATGCGTGCTGGCGCCCTGTTCGTTCAGGGCCTGCGAGAACGCGGTGAGCTGGAGCGTGTAAAACGGATCGAAGTGCGCCTGTATGGCTCGCTGt                                                  >  6:50782/1‑149 (MQ=255)
                                                                                         ttCGACCTGTTCAAGATCGGTGTCGGGCCCTCCAGCTCCCACACCGTCGGCCCCATGCGTGCTGGCGCCCTGTTCGTTCAGGGCCTGCGCGAACGCGGTGAGCTGGAGCGTGTAAAACGGATCGAAGTGCGCCTGTATGGCt                                                        >  5:249611/1‑142 (MQ=255)
                                                                                         ttCGACCTGTTCAAGATCGGTGTCGGGCCCTCCAGCTCCCACACCGTCGGCCCCATGCGTGCTGGCGCCCTGTTCGTTCAGGGCCTGCGCGAACGCGGTGAGCTGGAGCGTGTAAAACGGATCGAAGTGCGCCTGTATGGCt                                                        <  6:249611/142‑1 (MQ=255)
                                                                                           cGACCTGTTCAAGATCGGTGTCGGGCCCTCCAGCTCCCACACCGTCGGCCCCTTGCGTGCTGGCGCCCTGTTCGTTCAGGGCCTGCGCGAACGCGGTGAGCTGGAGCGTGTAAAACGGATCGAAGTGCGCCTGTATGGCTCGCTGTCa                                                <  1:266490/148‑1 (MQ=255)
                                                                                            gACCTGTTCAAGATCGGTGTCGGGCCCTCCAGCTCCCACACCGTCGGCCCCATGCGTGCTGGCGCCCTGTTCGTTCAGGGCCTGCGCGAACGCGGTGAGCTGGAGCGTGTAAAACGGATCGAAGTgcgc                                                                  <  5:63032/129‑1 (MQ=255)
                                                                                            gACCTGTTCAAGATCGGTGTCGGGCCCTCCAGCTCCCACACCGTCGGCCCCATGCGTGCTGGCGCCCTGTTCGTTCAGGGCCTGCGCGAACGCGGTGAGCTGGAGCGTGTAAAACGGATCGAAGTgcgc                                                                  >  6:63032/1‑129 (MQ=255)
                                                                                                  ttCAAGATCGGTGTCGGGCCCTCCAGCTCCCACACCGTCGGCCCCATGCGTGCTGGCGCCCTGTTCGTTCAGGGCCTGCGCGAACGCGGTGAGCTGGAGCGTGTAAAACGGATCGAAGTGCGCCTGTATGGCTCGCTGTCAGCCACCgg                                        <  6:51873/149‑1 (MQ=255)
                                                                                                  ttCAAGATCGGGGTCGGGCCCTCCAGCTCCCACACCGGCGGCCCCATGCGTGCTGGCGCCCTGTTCGTTCAGGGCCTGCGCGAACGCGGTGAGCTGGAGCGTGTAAAACGGATCGAAGTGCGCCTGTATGGCTCGCTGTCAGCCACCg                                         <  8:112598/148‑1 (MQ=255)
                                                                                                   tCAAGATCGGTGTCGGGCCCTCCAGCTCCCACACCGTCGGCCCCATGCGTGCTGGCGCCCTGTTCGTTCAGGGCCTGCGCGAACGCGGTGAGCTGGAGCGTGTAAAACGGATCGAAGTGCGCCTGTATGGCTCGCTGTCAGCCACCgg                                        <  6:145440/148‑1 (MQ=255)
                                                                                                     aaGATCGGTGTCGGGCCCTCCAGCTCCCACACCGTCGGCCCCATGCGTGCTGGCGCCCTGTTCGTTCAGGGCCTGCGCGAACGCGGTGAGCTGGAGCGTGTAAAACGGATCGAAGTGCGCCTGTATGGCTCGCTGTCAGCCACCGGTa                                      <  2:120882/148‑1 (MQ=255)
                                                                                                        aTCGGTGTCGGGCCCTCCAGCTCCCACACCGTCGGCCCCATGCGTGCTGGCGCCCTGTTCGTTCAGGGCCTGCGCGAACGCGGTGAGCTGGAGCGTGTAAAACGGATCGAAGTGCGCCTGTATGGCTCGCTGTCAGCCACCGGTATCgg                                  <  4:20544/149‑1 (MQ=255)
                                                                                                           ggTGTCGGGCCCTCCAGCTCCCACACCGTCGGCCCCATGCGTGCTGGCGCCCTGTTCGTTCAGGGCCTGCGCGAACGCGGTGAGCTGGAGCGTGTAAAACTGATCGAAGTGCGCCTGTATGGCTCGCTGTCAGCCACCGGTATCGGCCa                               <  5:50782/149‑1 (MQ=255)
                                                                                                                 gggCCCTCCAGCTCCCACACCGTCGGCCCCATGCGTGCTGGCGCCCTGTTCGTTCAGGGCCTGCGCGAACGCGGTGAGCTGGAGCGTGTAAAACGGATCGAAGTGCGCCTGTATGGCTCGCTGTCAGCCACCGGTATCGGCCATGGCAc                         >  8:21497/1‑149 (MQ=255)
                                                                                                                     ccTCCAGCTCCCACACCGTCGGCCCCATGCGTGCTGGCGCCCTGTTCGTTCAGGGCCTGCGCGAACGCGGTGAGCTGGAGCGTGTAAAACGGATCGAAGTGCGCCTGTATGGCTCGCTGTCAGCCACCGGTATCGGCCATGGCACCGAc                     <  1:280749/149‑1 (MQ=255)
                                                                                                                      cTCCAGCTCCCACACCGTCGGCCCCATGCGTGCTGGCGCCCTGTTCGTTCAGGGCCTGCGCGAACGCGGTGAGCTGGAGCGTGTAAAACGGATCGAAGTGCGCCTGTATGGCTCGCTGTCAGCCACCGGTATCGGCCATGGCACCGACa                    >  7:261845/1‑149 (MQ=255)
                                                                                                                          aGCTCCCACACCGTCGTCCCCATGCGTGCTGGCGCCCTGTTCGTTCAGGGCCTGCGCGAACGCGGTGAGCTGGAGCGTGTAAAACGGATCGAAGTGCGCCTGTATGGCTCGCTGTCAGCCACCGGTATCGGCCATGGCACCGACAACGc                <  7:233952/149‑1 (MQ=255)
                                                                                                                          aGCTCCCACACCGTCGTCCCCATGCGTGCTGGCGCCCTGTTCGTTCAGGGCCTGCGCGAACGCGGTGAGCTGGAGCGTGTAAAACGGATCGAAGTGCGCCTGTATGGCTCGCTGTCAGCCACCGGTATCGGCCATGGCACCGACAACGc                <  7:238352/149‑1 (MQ=255)
                                                                                                                             tCCCACACCGTCGGCCCCATGCGTGCTGGCGCCCTGTTCGTTCAGGGCCTGCGCGAACGCGGTGAGCTGGAGCGTGTAAAACGGATCGAAGTGCGCCTGTATGGCTCGCTGTCAGCCACCGGTATCGGCCATGGCACCGACAACGccac             >  5:233888/1‑149 (MQ=255)
                                                                                                                                 acacCGTCGGCCCCATGCGTGCTGGCGCCCTGTTCGTTCAGGGCCTGCGCGAACGCGGTGAGCTGGAGCGTGTAAAACGGATCGAAGTGCGCCTGTATGGCTCGCTGTCAGCCACCGGTATCGGCCATGGCACCGACAACGCCACcat          >  1:122482/1‑148 (MQ=255)
                                                                                                                                     cGTCGGCCCCATGCGTGCTGGCGCCCTGTTCGTTCAGGGCCTGCGCGAACGCGGTGAGCTGGAGCGTGTAAAACGGATCGAAGTGCGCCTGTATGGCTCGCTGTCAGCCACCGGTATCGGCCATGGCACCGACAACGCCACCATCATgg     <  8:139876/149‑1 (MQ=255)
                                                                                                                                     cGTCGGCCCCATGCGTGCTGGCGCCCTGTTCGTTCAGGGCCTGCGCGAACGCGGTGAGCTGGAGCGTGTAAAACGGATCGAAGTGCGCCTGTATGGCTCGCTGTCAGCCACCGGTATCGGCCATGGCACCGACAACGCCACCATCATgg     >  7:139876/1‑149 (MQ=255)
                                                                                                                                     cGTCGGCCCCATGCGTGCTGGCGCCCTGTTCGTTCAGGGCCTGCGCGAACGCGGTGAGCTGGAGCGTGTAAAACGGATCGAAGTGCGCCTGTATGGCTCGCTGTCAGCCACCGGTATCGGCCATGGCACCGACAACGCCACCATCATgg     <  3:104501/149‑1 (MQ=255)
                                                                                                                                        cGGCCCCATGCGTGCTGGCGCCCTGTTCGTTCAGGGCCTGCGCGAACGCGGTGAGCTGGAGCGTGTAAAACGGATCGAAGTGCGCCTGTATGGCTCGCTGTCAGCCACCGGTATCGGCCATGGCACCGACAACGCCACCATCATGGGcc  <  1:288819/149‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                                  
GCTGTTTGAGCCTAGCAGCTACCCCGCCCTGCGCATGGCGCGCAGAGACAAGCCCAGCAAGAGGTGATTCGACAATGGCCATCAGCGTGTTCGACCTGTTCAAGATCGGTGTCGGGCCCTCCAGCTCCCACACCGTCGGCCCCATGCGTGCTGGCGCCCTGTTCGTTCAGGGCCTGCGCGAACGCGGTGAGCTGGAGCGTGTAAAACGGATCGAAGTGCGCCTGTATGGCTCGCTGTCAGCCACCGGTATCGGCCATGGCACCGACAACGCCACCATCATGGGCC  >  NC_002947/356582‑356866

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