Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,212,492 C→T 8.0% intergenic (+56/+38) PP_4593 → / ← PP_4594 hypothetical protein/cystathionine gamma‑synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,212,4920CT 8.0% 142.0 / 2.2 50intergenic (+56/+38)PP_4593/PP_4594hypothetical protein/cystathionine gamma‑synthase
Reads supporting (aligned to +/- strand):  ref base C (21/25);  new base T (2/2);  total (23/27)
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

CACCAACGACGAAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGA‑GCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTACCGACAACCGCAGCAGGTTCTCGCTGATCCCCATCACC  >  NC_002947/5212346‑5212638
                                                                                                                                                  |                                                                                                                                                   
cACCAACGACGAAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACcgtcg                                                                                                                                                    >  4:147542/1‑148 (MQ=255)
  ccAACGACGAAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAg                                                                                                                                                 <  3:147542/149‑1 (MQ=255)
    aaCGACGAAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAg                                                                                                                                                 >  1:286434/1‑147 (MQ=255)
    aaCGACGAAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAg                                                                                                                                                 >  4:299851/1‑147 (MQ=255)
      cgacgaAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGA‑Gc                                                                                                                                            <  2:286434/149‑1 (MQ=255)
       gacgaAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCCTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGA‑GCg                                                                                                                                           >  4:47321/1‑149 (MQ=255)
        acgaAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACGACGGTGTAGGA‑GCgg                                                                                                                                          <  2:212576/149‑1 (MQ=255)
        acgaAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGA‑GCgg                                                                                                                                          <  5:268045/149‑1 (MQ=255)
        acgaAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGA‑GCgg                                                                                                                                          <  3:299851/149‑1 (MQ=255)
              ccGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGA‑GCGGGTTTAc                                                                                                                                    <  2:185938/149‑1 (MQ=255)
                   aaCCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGA‑GCGGGTTTACCcgc                                                                                                                                >  5:21979/1‑148 (MQ=255)
                        ccTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCCTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGA‑GCGGGTTTACCCGCGAAAgg                                                                                                                          <  3:47321/149‑1 (MQ=255)
                             gcAAGGAAAGATTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAACCCGCTCCTACGACGGTGTAGGA‑GCGGGcttg                                                                                                                                     <  7:37352/133‑2 (MQ=37)
                             gcAAGGAAAGATTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAACCCGCTCCTACGACGGTGTAGGA‑GCGGGcttg                                                                                                                                     >  8:37352/1‑132 (MQ=37)
                               aaGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGA‑GCGGGTTTACCCGCGAAAGGGCCTGCg                                                                                                                   >  8:123957/1‑149 (MQ=255)
                                    aaGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGA‑GCGGGTTTACCCGCGAAAGGGCCTGCGGCCCt                                                                                                              <  8:12289/149‑1 (MQ=255)
                                       gTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGA‑GCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTcac                                                                                                           <  2:283241/149‑1 (MQ=255)
                                        ttGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGA‑GCGGGTTTAc                                                                                                                                    <  5:50136/123‑1 (MQ=255)
                                        ttGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGA‑GCGGGTTTAc                                                                                                                                    >  6:50136/1‑123 (MQ=255)
                                            tGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGA‑GCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGc                                                                                                      >  6:253936/1‑149 (MQ=255)
                                               tCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGA‑GCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCgcc                                                                                                   >  4:222819/1‑147 (MQ=255)
                                                 aaGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGA‑GCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCa                                                                                                 <  7:123957/149‑1 (MQ=255)
                                                      cAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGA‑GCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCt                                                                                            >  4:233573/1‑149 (MQ=255)
                                                      cAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGA‑GCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCt                                                                                            <  7:75610/149‑1 (MQ=255)
                                                      cAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGA‑GCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCt                                                                                            >  6:210027/1‑149 (MQ=255)
                                                        aGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGA‑GCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCttt                                                                                          >  3:62575/1‑149 (MQ=255)
                                                           tCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGACGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTc                                                                                         <  1:86590/148‑1 (MQ=255)
                                                             gAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGA‑GCGGGTTTACCCGCGAAAGGGCCTGCTGCCCTCACACAGCGACCAGCGCTTTCAACg                                                                                     >  1:266362/1‑149 (MQ=255)
                                                              aaCGGATCATTCATCCGCCCAAATCCTGACAGCCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGA‑GCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGc                                                                                    >  1:207150/1‑149 (MQ=255)
                                                                      attcatCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGA‑GCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATcc                                                                             <  6:21979/148‑1 (MQ=255)
                                                                       ttcatCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGA‑GCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAg                                                                           <  3:233573/149‑1 (MQ=255)
                                                                               cccAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGA‑GCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGa                                                                   <  5:136923/149‑1 (MQ=255)
                                                                                ccAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGA‑GCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGAc                                                                  >  7:136589/1‑149 (MQ=255)
                                                                                 cAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGA‑GCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGAc                                                                  >  8:48455/1‑148 (MQ=255)
                                                                                                        aaTCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGA‑GCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATccc                                          >  2:168914/1‑149 (MQ=255)
                                                                                                         aTCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGA‑GCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCt                                         <  5:251682/149‑1 (MQ=255)
                                                                                                         aTCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGA‑GCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCt                                         <  1:168914/149‑1 (MQ=255)
                                                                                                          tCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGA‑GCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAg                                                                           >  5:312737/1‑114 (MQ=255)
                                                                                                          tCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGA‑GCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAg                                                                           <  6:312737/114‑1 (MQ=255)
                                                                                                             aCCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGA‑GCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGGTGTCAACGCCCCAGCCAGGTCCTCGACCAGGTCCGGGGACTCGTCGATCGCTACCg                                     >  8:118115/1‑149 (MQ=37)
                                                                                                              ccGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGA‑GCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTACCGa                                    <  2:207150/149‑1 (MQ=255)
                                                                                                                      ttCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGA‑GCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTACCGa                                    <  4:219635/141‑1 (MQ=255)
                                                                                                                      ttCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGA‑GCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTACCGa                                    >  3:219635/1‑141 (MQ=255)
                                                                                                                      ttCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGA‑GCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTACCGACAACCg                              >  7:125355/1‑147 (MQ=255)
                                                                                                                        cgcgGGCAAGCCCGCTCCTACACCGTCGTAGGA‑GCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCTTCGATCCCTACCGACAACCgca                            <  7:48455/147‑1 (MQ=255)
                                                                                                                        cgcgGGCAAGCCCGCTCCTACACCGTCGTAGGA‑GCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTACCGACAACCgca                            <  3:280322/147‑1 (MQ=255)
                                                                                                                                aGCCCGCTCCTACACCGTCGTAGGA‑GCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTc                                                <  5:141054/119‑1 (MQ=255)
                                                                                                                                aGCCCGCTCCTACACCGTCGTAGGA‑GCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTc                                                >  6:141054/1‑119 (MQ=255)
                                                                                                                                  cccGCTCCTACGACGGTGTAGGA‑GCGGGCTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTACCGACAACCGCAGCAGGTTCTCGc                >  4:274130/1‑149 (MQ=39)
                                                                                                                                    agCTCCTACACCGTCGTAGGA‑GCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCAAGGTAATCCTCGATCCCTACCGACAACCGCAGCAGGTTCTCGCTg              <  5:210027/148‑1 (MQ=255)
                                                                                                                                     gCTCCTACACCGTCGTAGGA‑GCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTACCGACAACCGCAGCAGGTTCTCGCTGa             <  5:253936/149‑1 (MQ=255)
                                                                                                                                                gtcgtAGGA‑GCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTACCGACAACCGCAGCAGGTTCTCGCTGATCCCCATCAcc  <  4:62575/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
CACCAACGACGAAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGA‑GCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTACCGACAACCGCAGCAGGTTCTCGCTGATCCCCATCACC  >  NC_002947/5212346‑5212638

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