Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,212,487 A→G 7.6% intergenic (+51/+43) 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,4870AG 7.6% 145.5 / 2.9 53intergenic (+51/+43)PP_4593/PP_4594hypothetical protein/cystathionine gamma‑synthase
Reads supporting (aligned to +/- strand):  ref base A (24/24);  new base G (2/2);  total (27/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 = 1.84e-01

AACGGCACCAACGACGAAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTACCGACAACCGCAGCAGGTTCTCGCTGATC  >  NC_002947/5212341‑5212629
                                                                                                                                                  |                                                                                                                                              
aaCGGCACCAACGACGAAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCGGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACAcc                                                                                                                                              <  8:100188/149‑1 (MQ=255)
        cAACGACGAAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAg                                                                                                                                       >  6:238917/1‑148 (MQ=255)
         aaCGACGAAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAg                                                                                                                                       >  8:61088/1‑147 (MQ=255)
         aaCGACGAAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAg                                                                                                                                       >  8:175664/1‑147 (MQ=255)
          acgacgAAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAg                                                                                                                                    >  3:353219/1‑149 (MQ=255)
            gacgaAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCg                                                                                                                                  >  4:3154/1‑149 (MQ=255)
                        aaCCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTAcc                                                                                                                          >  4:222234/1‑145 (MQ=255)
                        aaCCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTAcc                                                                                                                          <  3:222234/145‑1 (MQ=255)
                         aCCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCcgcg                                                                                                                      >  3:156292/1‑148 (MQ=255)
                               tGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAg                                                                                                                                       <  2:357956/125‑1 (MQ=255)
                               tGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAg                                                                                                                                       >  1:357956/1‑125 (MQ=255)
                                 cgcAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGTAAACCCGCTCCTACGACGGTGTAGGAGCGGGCTTGCCCGCGaagaggc                                                                                                               <  1:147139/147‑6 (MQ=25)
                                 cgcAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCcgc                                                                                                                       >  7:57745/1‑139 (MQ=255)
                                 cgcAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCcgc                                                                                                                       <  8:57745/139‑1 (MQ=255)
                                  gcAACGAAACGTTGCTGATCAACGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTg                                                                                                            <  2:320263/149‑1 (MQ=255)
                                   cAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGAACATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGc                                                                                                           >  6:21442/1‑149 (MQ=255)
                                      ggAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGc                                                                                                        >  7:177320/1‑149 (MQ=255)
                                       gAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGTGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGcc                                                                                                       >  6:65074/1‑149 (MQ=255)
                                         aaGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATACACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGTGCCTGCGGCCCt                                                                                                     <  2:46585/149‑1 (MQ=255)
                                               gCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGaa                                                                                                                    <  8:90622/128‑1 (MQ=255)
                                               gCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGaa                                                                                                                    >  7:90622/1‑128 (MQ=255)
                                                   aTCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGc                                                                                             <  7:180350/147‑1 (MQ=255)
                                                   aTCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGa                                                                                           <  8:315286/149‑1 (MQ=255)
                                                      aaGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCa                                                                                        <  7:61088/149‑1 (MQ=255)
                                                      aaGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATCGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCCCACAGCGAc                                                                                          >  4:132509/1‑147 (MQ=255)
                                                         gACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAgcg                                                                                     >  2:11337/1‑149 (MQ=255)
                                                         gACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAACCCGCTCCTACGACGGTGTAGGAGCGGGCTTGcc                                                                                                                          <  1:35252/112‑1 (MQ=25)
                                                         gACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAACCCGCTCCTACGACGGTGTAGGAGCGGGCTTGcc                                                                                                                          >  2:35252/1‑112 (MQ=25)
                                                                tCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGAc                                                                                          >  8:36554/1‑137 (MQ=255)
                                                                tCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGAc                                                                                          <  7:36554/137‑1 (MQ=255)
                                                                 cGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCGCAGGCCCtttcgc                                                                                                >  5:252252/1‑125 (MQ=38)
                                                                 cGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCGCAGGCCCtttcgc                                                                                                <  6:252252/130‑6 (MQ=38)
                                                                   aaCGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGa                                                                                                                                     >  4:261028/1‑91 (MQ=255)
                                                                   aaCGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGa                                                                                                                                     <  3:261028/91‑1 (MQ=255)
                                                                   aaCGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCGTCACACAGCGACCAGCGCTTTCAACGc                                                                           <  1:11337/149‑1 (MQ=255)
                                                                   aaCGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCACACACAGCGACCAGCGCTTTCAACGc                                                                           <  4:156292/149‑1 (MQ=255)
                                                                            ttcatCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAg                                                                  <  5:238917/149‑1 (MQ=255)
                                                                                      cAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGAcc                                                        <  5:21442/149‑1 (MQ=255)
                                                                                               aCAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCAc                                                 <  4:353219/147‑1 (MQ=255)
                                                                                                   gCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAAt                                           >  4:215939/1‑149 (MQ=255)
                                                                                                     aTCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACTCCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATcc                                         >  1:256011/1‑149 (MQ=255)
                                                                                                        tATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCACCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTc                                       >  1:64716/1‑148 (MQ=255)
                                                                                                         aTTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCt                                        >  6:368810/1‑146 (MQ=255)
                                                                                                         aTTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCt                                        <  5:368810/146‑1 (MQ=255)
                                                                                                             aaTCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACAACGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGCCCAGGTCACGGTAATCATCGATcc                                  >  8:314989/1‑148 (MQ=255)
                                                                                                              aTCCACCGGCCTCTTCGCGGGGAAGCCCGCTCCTACACCGTCGTAGTAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCt                                <  3:202258/149‑1 (MQ=255)
                                                                                                              aTCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATccc                                 >  5:278016/1‑148 (MQ=255)
                                                                                                                 cACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTAcc                             <  7:175664/149‑1 (MQ=255)
                                                                                                                       ccTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCATCGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTACCGACAAc                       <  3:3154/149‑1 (MQ=255)
                                                                                                                                          gCTCCTACGACGGTGTAGGAGCGGGCTTGCCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTACCGACAACCGCAGCAGGTTCTCGCTGa    >  4:220718/1‑149 (MQ=39)
                                                                                                                                          gCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATc                                          <  5:332125/111‑1 (MQ=255)
                                                                                                                                          gCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAAAGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATc                                          >  6:332125/1‑111 (MQ=255)
                                                                                                                                              cTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTACCGACAACCGCAGCAGGTTCTCGCTGATc  <  2:256011/147‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                              
AACGGCACCAACGACGAAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTACCGACAACCGCAGCAGGTTCTCGCTGATC  >  NC_002947/5212341‑5212629

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