Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,212,488 C→A 7.9% intergenic (+52/+42) 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,4880CA 7.9% 150.7 / 4.2 51intergenic (+52/+42)PP_4593/PP_4594hypothetical protein/cystathionine gamma‑synthase
Reads supporting (aligned to +/- strand):  ref base C (22/25);  new base A (2/2);  total (24/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 = 3.00e-01

CACCAACGACGAAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGA‑GCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTACCGACAACCGCAGCAGGTTCTCGCTGA  >  NC_002947/5212346‑5212627
                                                                                                                                              |                                                                                                                                            
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                                                                <  6:312737/114‑1 (MQ=255)
                                                                                                          tCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGA‑GCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAg                                                                >  5:312737/1‑114 (MQ=255)
                                                                                                             aCCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGA‑GCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGGTGTCAACGCCCCAGCCAGGTCCTCGACCAGGTCCGGGGACTCGTCGATCGCTACCg                          >  8:118115/1‑149 (MQ=37)
                                                                                                              ccGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGA‑GCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTACCGa                         <  2:207150/149‑1 (MQ=255)
                                                                                                                      ttCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGA‑GCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTACCGa                         >  3:219635/1‑141 (MQ=255)
                                                                                                                      ttCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGA‑GCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTACCGa                         <  4:219635/141‑1 (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)
                                                                                                                                              |                                                                                                                                            
CACCAACGACGAAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGA‑GCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTACCGACAACCGCAGCAGGTTCTCGCTGA  >  NC_002947/5212346‑5212627

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